  Cooperation and defection may be considered as two extreme responses to a
social dilemma. Yet the reality is much less clear-cut. Between the two
extremes lies an interval of ambivalent choices, which may be captured
theoretically by means of continuous strategies defining the extent of the
contributions of each individual player to the common pool. If strategies are
chosen from the unit interval, where 0 corresponds to pure defection and 1
corresponds to the maximal contribution, the question is what is the
characteristic level of individual investments to the common pool that emerges
if the evolution is guided by different benefit functions. Here we consider the
steepness and the threshold as two parameters defining an array of generalized
benefit functions, and we show that in a structured population there exist
intermediate values of both at which the collective contributions are maximal.
However, as the cost-to-benefit ratio of cooperation increases the
characteristic threshold decreases, while the corresponding steepness
increases. Our observations remain valid if more complex sigmoid functions are
used, thus reenforcing the importance of carefully adjusted benefits for high
levels of public cooperation.


  Spontaneous symmetry breaking plays a fundamental role in many areas of
condensed matter and particle physics. A fundamental problem in ecology is the
elucidation of the mechanisms responsible for biodiversity and stability.
Neutral theory, which makes the simplifying assumption that all individuals
(such as trees in a tropical forest) --regardless of the species they belong
to-- have the same prospect of reproduction, death, etc., yields gross patterns
that are in accord with empirical data. We explore the possibility of birth and
death rates that depend on the population density of species while treating the
dynamics in a species-symmetric manner. We demonstrate that the dynamical
evolution can lead to a stationary state characterized simultaneously by both
biodiversity and spontaneously broken neutral symmetry.


  We suggest a new method to compute the asymptotic fitness distribution in the
Bak-Sneppen model of biological evolution. As applications we derive the full
asymptotic distribution in the four-species model, and give an explicit linear
recurrence relation for a set of coefficients determining the asymptotic
distribution in the five-species model.


  The rock-paper-scissor game -- which is characterized by three strategies
R,P,S, satisfying the non-transitive relations S excludes P, P excludes R, and
R excludes S -- serves as a simple prototype for studying more complex
non-transitive systems. For well-mixed systems where interactions result in
fitness reductions of the losers exceeding fitness gains of the winners,
classical theory predicts that two strategies go extinct. The effects of
spatial heterogeneity and dispersal rates on this outcome are analyzed using a
general framework for evolutionary games in patchy landscapes. The analysis
reveals that coexistence is determined by the rates at which dominant
strategies invade a landscape occupied by the subordinate strategy (e.g. rock
invades a landscape occupied by scissors) and the rates at which subordinate
strategies get excluded in a landscape occupied by the dominant strategy (e.g.
scissor gets excluded in a landscape occupied by rock). These invasion and
exclusion rates correspond to eigenvalues of the linearized dynamics near
single strategy equilibria. Coexistence occurs when the product of the invasion
rates exceeds the product of the exclusion rates. Provided there is sufficient
spatial variation in payoffs, the analysis identifies a critical dispersal rate
$d^*$ required for regional persistence. For dispersal rates below $d^*$, the
product of the invasion rates exceed the product of the exclusion rates and the
rock-paper-scissor metacommunities persist regionally despite being extinction
prone locally. For dispersal rates above $d^*$, the product of the exclusion
rates exceed the product of the invasion rates and the strategies are
extinction prone. These results highlight the delicate interplay between
spatial heterogeneity and dispersal in mediating long-term outcomes for
evolutionary games.


  Motivated by the observation that anomalous diffusion is a realistic feature
in the dynamics of biological populations, we investigate its implications in a
paradigmatic model for the evolution of a single species density $u(x,t)$. The
standard model includes growth and competition in a logistic expression, and
spreading is modeled through normal diffusion. Moreover, the competition term
is nonlocal, which has been shown to give rise to spatial patterns. We
generalize the diffusion term through the nonlinear form $\partial_t u(x,t) = D
\partial_{xx} u(x,t)^\nu$ (with $D, \nu>0$), encompassing the cases where the
state-dependent diffusion coefficient either increases ($\nu>1$) or decreases
($\nu<1$) with the density, yielding subdiffusion or superdiffusion,
respectively. By means of numerical simulations and analytical considerations,
we display how that nonlinearity alters the phase diagram. The type of
diffusion imposes critical values of the model parameters for the onset of
patterns and strongly influences their shape, inducing fragmentation in the
subdiffusive case. The detection of the main persistent mode allows analytical
prediction of the critical thresholds.


  We pursue the task of developing a finite population counterpart to Eigen's
model. We consider the classical Wright-Fisher model describing the evolution
of a population of size $m$ of chromosomes of length $\ell$ over an alphabet of
cardinality $\kappa$. The mutation probability per locus is $q$. The
replication rate is $\sigma>1$ for the master sequence and $1$ for the other
sequences. We study the equilibrium distribution of the process in the regime
where $\ell\to+\infty$, $m\to+\infty$, $q\to0$, $\ell q\to a\in\,]0,+\infty[$,
$\frac{m}{\ell}\to\alpha\in [0,+\infty]$. We obtain an equation
$\alpha\psi(a)=\ln\kappa$ in the parameter space $(a,\alpha)$ separating the
regime where the equilibrium population is totally random from the regime where
a quasispecies is formed. We observe the existence of a critical population
size necessary for a quasispecies to emerge, and we recover the finite
population counterpart of the error threshold. The result is the twin brother
of the corresponding result for the Moran model. The proof is more complex, and
it relies on the Freidlin-Wentzell theory of random perturbations of dynamical
systems.


  The importance of statistical patterns of language has been debated over
decades. Although Zipf's law is perhaps the most popular case, recently,
Menzerath's law has begun to be involved. Menzerath's law manifests in
language, music and genomes as a tendency of the mean size of the parts to
decrease as the number of parts increases in many situations. This statistical
regularity emerges also in the context of genomes, for instance, as a tendency
of species with more chromosomes to have a smaller mean chromosome size. It has
been argued that the instantiation of this law in genomes is not indicative of
any parallel between language and genomes because (a) the law is inevitable and
(b) non-coding DNA dominates genomes. Here mathematical, statistical and
conceptual challenges of these criticisms are discussed. Two major conclusions
are drawn: the law is not inevitable and languages also have a correlate of
non-coding DNA. However, the wide range of manifestations of the law in and
outside genomes suggests that the striking similarities between non-coding DNA
and certain linguistics units could be anecdotal for understanding the
recurrence of that statistical law.


  Resistance to chemotherapies, particularly to anticancer treatments, is an
increasing medical concern. Among the many mechanisms at work in cancers, one
of the most important is the selection of tumor cells expressing resistance
genes or phenotypes. Motivated by the theory of mutation-selection in adaptive
evolution, we propose a model based on a continuous variable that represents
the expression level of a resistance gene (or genes, yielding a phenotype)
influencing in healthy and tumor cells birth/death rates, effects of
chemotherapies (both cytotoxic and cytostatic) and mutations. We extend
previous work by demonstrating how qualitatively different actions of
chemotherapeutic and cytostatic treatments may induce different levels of
resistance. The mathematical interest of our study is in the formalism of
constrained Hamilton-Jacobi equations in the framework of viscosity solutions.
We derive the long-term temporal dynamics of the fittest traits in the regime
of small mutations. In the context of adaptive cancer management, we also
analyse whether an optimal drug level is better than the maximal tolerated
dose.


  Group-based models arise in algebraic statistics while studying evolution
processes. They are represented by embedded toric algebraic varieties. Both
from the theoretical and applied point of view one is interested in determining
the ideals defining the varieties. Conjectural bounds on the degree in which
these ideals are generated were given by Sturmfels and Sullivant. We prove that
for the 3-Kimura model, corresponding to the group G=Z2xZ2, the projective
scheme can be defined by an ideal generated in degree 4. In particular, it is
enough to consider degree 4 phylogenetic invariants to test if a given point
belongs to the variety. We also investigate G-models, a generalization of
abelian group-based models. For any G-model, we prove that there exists a
constant $d$, such that for any tree, the associated projective scheme can be
defined by an ideal generated in degree at most d.


  A combined molecular docking and molecular structure in silico analysis on
the substrate and product of leukotriene A4 hydrolase (LTA4H) was performed.
The molecular structures of the substrate leukotriene A4 (LTA4) and product
leukotirene B4 (LTB4) were studied through Density Functional Theory (DFT)
calculations at the B3LYP/6-31+G(d) level of theory in both, gas and condensed
phases. The whole LTB4 molecule was divided into three fragments (hydrophobic
tail, triene motif, and a polar acidic group) which were subjected to a full
conformational study employing the most stable conformations of them to build
conformers of the complete molecule and geometry optimize further. LTA4
conformers structures were modeled from the LTB4 minimum energy conformers.
Both, protonated and deprotonated species of LTA4 and LTB4, were analyzed
according to pKa values founded in the literature. Finally, a binding model of
LTA4 with LTA4 hydrolase is proposed according to docking results which show
intermolecular interactions that position the protonated and deprotonated
ligand in the active site, in excellent agreement with experimental data.


  We introduce a stochastic model describing aggregation of misfolded proteins
and degradation by the protein quality control system in a single cell. In
analogy with existing literature, aggregates can grow, nucleate and fragment
stochastically. We assume that the quality control system acts as an enzyme
that can degrade aggregates at different stages of the growth process, with an
efficiency that decreases with the size of the aggregate. We show how this
stochastic dynamics, depending on the parameter choice, leads to two
qualitatively different behaviors: a homeostatic state, where the quality
control system is stable and aggregates of large sizes are not formed, and an
oscillatory state, where the quality control system periodically breaks down,
allowing for the formation of large aggregates. We discuss how these periodic
breakdowns may constitute a mechanism for the sporadic development of
neurodegenerative diseases.


  The goal of branch length estimation in phylogenetic inference is to estimate
the divergence time between a set of sequences based on compositional
differences between them. A number of software is currently available
facilitating branch lengths estimation for homogeneous and stationary
evolutionary models. Homogeneity of the evolutionary process imposes fixed
rates of evolution throughout the tree. In complex data problems this
assumption is likely to put the results of the analyses in question.
  In this work we propose an algorithm for parameter and branch lengths
inference in the discrete-time Markov processes on trees. This broad class of
nonhomogeneous models comprises the general Markov model and all its submodels,
including both stationary and nonstationary models.
  Here, we adapted the well-known Expectation-Maximization algorithm and
present a detailed performance study of this approach for a selection of
nonhomogeneous evolutionary models. We conducted an extensive performance
assessment on multiple sequence alignments simulated under a variety of
settings. We demonstrated high accuracy of the tool in parameter estimation and
branch lengths recovery, proving the method to be a valuable tool for
phylogenetic inference in real life problems. $\empar$ is an open-source C++
implementation of the methods introduced in this paper and is the first tool
designed to handle nonhomogeneous data.


  The atomistic characterization of the transition state is a fundamental step
to improve the understanding of the folding mechanism and the function of
proteins. From a computational point of view, the identification of the
conformations that build out the transition state is particularly cumbersome,
mainly because of the large computational cost of generating a
statistically-sound set of folding trajectories. Here we show that a biasing
algorithm, based on the physics of the ratchet-and-pawl, can be used to
identify efficiently the transition state. The basic idea is that the
algorithmic ratchet exerts a force on the protein when it is climbing the
free-energy barrier, while it is inactive when it is descending. The transition
state can be identified as the point of the trajectory where the ratchet
changes regime. Besides discussing this strategy in general terms, we test it
within a protein model whose transition state can be studied independently by
plain molecular dynamics simulations. Finally, we show its power in
explicit-solvent simulations, obtaining and characterizing a set of
transition--state conformations for ACBP and CI2.


  For decades, dimethyl sulfate (DMS) mapping has informed manual modeling of
RNA structure in vitro and in vivo. Here, we incorporate DMS data into
automated secondary structure inference using a pseudo-energy framework
developed for 2'-OH acylation (SHAPE) mapping. On six non-coding RNAs with
crystallographic models, DMS- guided modeling achieves overall false negative
and false discovery rates of 9.5% and 11.6%, comparable or better than
SHAPE-guided modeling; and non-parametric bootstrapping provides
straightforward confidence estimates. Integrating DMS/SHAPE data and including
CMCT reactivities give small additional improvements. These results establish
DMS mapping - an already routine technique - as a quantitative tool for
unbiased RNA structure modeling.


  The evolution and distribution of species body sizes for terrestrial mammals
is well-explained by a macroevolutionary tradeoff between short-term selective
advantages and long-term extinction risks from increased species body size,
unfolding above the 2g minimum size induced by thermoregulation in air. Here,
we consider whether this same tradeoff, formalized as a constrained
convection-reaction-diffusion system, can also explain the sizes of fully
aquatic mammals, which have not previously been considered. By replacing the
terrestrial minimum with a pelagic one, at roughly 7000g, the terrestrial
mammal tradeoff model accurately predicts, with no tunable parameters, the
observed body masses of all extant cetacean species, including the 175,000,000g
Blue Whale. This strong agreement between theory and data suggests that a
universal macroevolutionary tradeoff governs body size evolution for all
mammals, regardless of their habitat. The dramatic sizes of cetaceans can thus
be attributed mainly to the increased convective heat loss is water, which
shifts the species size distribution upward and pushes its right tail into
ranges inaccessible to terrestrial mammals. Under this macroevolutionary
tradeoff, the largest expected species occurs where the rate at which
smaller-bodied species move up into large-bodied niches approximately equals
the rate at which extinction removes them.


  We introduce and analyse an individual-based evolutionary model, in which a
population of genetically diverse organisms compete with each other for limited
resources. Through theoretical analysis and stochastic simulations, we show
that the model exhibits a pattern-forming instability which is highly amplified
by the effects of demographic noise, leading to the spontaneous formation of
genotypic clusters. This mechanism supports the thesis that stochasticity has a
central role in the formation and coherence of species.


  The linear noise approximation is commonly used to obtain intrinsic noise
statistics for biochemical networks. These estimates are accurate for networks
with large numbers of molecules. However it is well known that many biochemical
networks are characterized by at least one species with a small number of
molecules. We here describe version 0.3 of the software intrinsic Noise
Analyzer (iNA) which allows for accurate computation of noise statistics over
wide ranges of molecule numbers. This is achieved by calculating the next order
corrections to the linear noise approximation's estimates of variance and
covariance of concentration fluctuations. The efficiency of the methods is
significantly improved by automated just-in-time compilation using the LLVM
framework leading to a fluctuation analysis which typically outperforms that
obtained by means of exact stochastic simulations. iNA is hence particularly
well suited for the needs of the computational biology community.


  Many transport processes in nature take place on substrates, often considered
as unidimensional lanes. These unidimensional substrates are typically
non-static: affected by a fluctuating environment, they can undergo
conformational changes. This is particularly true in biological cells, where
the state of the substrate is often coupled to the active motion of
macromolecular complexes, such as motor proteins on microtubules or ribosomes
on mRNAs, causing new interesting phenomena. Inspired by biological processes
such as protein synthesis by ribosomes and motor protein transport, we
introduce the concept of localized dynamical sites coupled to a driven lattice
gas dynamics. We investigate the phenomenology of transport in the presence of
dynamical defects and find a novel regime characterized by an intermittent
current and subject to severe finite-size effects. Our results demonstrate the
impact of the regulatory role of the dynamical defects in transport, not only
in biology but also in more general contexts.


  We present an alternative domain concerning mathematics to investigate
universal evolution mechanisms by focusing on large cycles theory (LCT) - a
simplified version of well-known hamiltonian graph theory. LCT joins together a
number of $NP$-complete cycle problems in graph theory. $NP$-completeness is
the kay factor insuring (by conjecture of Cook) the generation of endless
developments and great diversity around large cycles problems. Originated about
60 years ago, the individuals (claims, propositions, lemmas, conjectures,
theorems, and so on) in LCT continually evolve and adapt to their environment
by an iterative process from primitive beginnings to best possible theorems
based on inductive reasoning. LCT evolves much more rapidly than biosphere and
has a few thousand pronounced species (theorems). Recall that life on earth
with more than 2 million species was originated about 3.7 billion years ago and
evolves extremely slowly. We show that all theorems in LCT have descended from
some common primitive propositions such as "every complete graph is
hamiltonian" or "every graph contains a cycle of length at least one" via
improvements, modifications and three kinds of generalizations - closing,
associating and extending. It is reasonable to review Darwinian mechanisms in
light of LCT evolution mechanisms (especially inductive reasoning) including
the origin and macroevolution disputable phenomena in the biosphere.


  A model with six mutually-exclusive compartments related to Dengue disease is
presented. In this model there are three vector control tools: insecticides
(larvicide and adulticide) and mechanical control. The problem is studied using
an Optimal Control (OC) approach. The human data for the model is based on the
Cape Verde Dengue outbreak. Some control measures are simulated and their
consequences analyzed.


  The mechanism by which silver staining of proteins in polyacrylamide gels
interferes with mass spectrometry of peptides produced by proteolysis has been
investigated. It was demonstrated that this interference increases with time
between silver staining and gel processing, although the silver image is
constant. This suggested an important role of the formaldehyde used in silver
staining development in this interference process. Consequently, a
formaldehyde-free staining protocol has been devised, using carbohydrazide as
the developing agent. This protocol showed much increased peptide coverage and
retained the sensitivity of silver staining. These results were however
obtained at the expense of an increased background in the stained gels and of a
reduced staining homogeneity.


  In recent years, significant progress in understanding the properties of
supercoiled DNA has been obtained due to nanotechniques that made stretching
and twisting of single molecules possible. Quantitative interpretation of such
experiments requires accurate knowledge of torques inside manipulated DNA. This
paper argues that it is not possible to transfer the entire magnitudes of
external torques to the twisting stress of the double helix, and that a
reducing torque transfer coefficient (TTC<1) should always be assumed. This
assertion agrees with simple physical intuition and is supported by the results
of all-atom molecular dynamics (MD) simulations. According to MD, the TTCs
around 0.8 are observed in nearly optimal conditions. Reaching higher values
requires special efforts and it should be difficult in practice. The TTC can be
partially responsible for the persistent discrepancies between the twisting
rigidity of DNA measured by different methods.


  In the framework of the paradigmatic prisoner's dilemma, we investigate the
evolutionary dynamics of social dilemmas in the presence of "cooperation
facilitators". In our model, cooperators and defectors interact as in the
classical prisoner's dilemma game, where selection favors defection. However,
here the presence of a small number of cooperation facilitators enhances the
fitness (reproductive potential) of cooperators, while it does not alter that
of defectors. In a finite population of size N, the dynamics of the prisoner's
dilemma with facilitators is characterized by the probability that cooperation
takes over (fixation probability) and by the mean times to reach the absorbing
states. These quantities are computed exactly and using Fokker-Planck
equations. Our findings, corroborated by stochastic simulations, demonstrate
that the influence of facilitators crucially depends on the difference between
their density z and the game's cost-to-benefit ratio r. When z>r, the fixation
of cooperators is likely in a large population and, under weak selection
pressure, invasion and replacement of defection by cooperation is favored by
selection if b(z-r)(1-z)>1/N, where 0<b<= 1 is the cooperation payoff benefit.
When z<r, the fixation probability of cooperators is exponentially enhanced by
the presence of facilitators but defection is the dominating strategy.


  In this paper, we inspect well-known population genetics and social dynamics
models. In these models, interacting individuals, while participating in a
self-organizing process, give rise to the emergence of complex behaviors and
patterns. While one main focus in population genetics is on the adaptive
behavior of a population, social dynamics is more often concerned with the
splitting of a connected array of individuals into a state of global
polarization, that is, the emergence of speciation. Applying computational and
mathematical tools we show that the way the mechanisms of selection,
interaction and replacement are constrained and combined in the modeling have
an important bearing on both adaptation and the emergence of speciation.
Differently (un)constraining the mechanism of individual replacement provides
the conditions required for either speciation or adaptation, since these
features appear as two opposing phenomena, not achieved by one and the same
model. Even though natural selection, operating as an external, environmental
mechanism, is neither necessary nor sufficient for the creation of speciation,
our modeling exercises highlight the important role played by natural selection
in the interplay of the evolutionary and the self-organization modeling
methodologies.


  Invasion fronts in ecology are well studied but very few mathematical results
concern the case with variable motility (possibly due to mutations). Based on
an apparently simple reaction-diffusion equation, we explain the observed
phenomena of front acceleration (when the motility is unbounded) as well as
other quantitative results, such as the selection of the most motile
individuals (when the motility is bounded). The key argument for the
construction and analysis of traveling fronts is the derivation of the
dispersion relation linking the speed of the wave and the spatial decay. When
the motility is unbounded we show that the position of the front scales as
$t^{3/2}$. When the mutation rate is low we show that the canonical equation
for the dynamics of the fittest trait should be stated as a PDE in our context.
It turns out to be a type of Burgers equation with source term.


  We present Quip, a lossless compression algorithm for next-generation
sequencing data in the FASTQ and SAM/BAM formats. In addition to implementing
reference-based compression, we have developed, to our knowledge, the first
assembly-based compressor, using a novel de novo assembly algorithm. A
probabilistic data structure is used to dramatically reduce the memory required
by traditional de Bruijn graph assemblers, allowing millions of reads to be
assembled very efficiently. Read sequences are then stored as positions within
the assembled contigs. This is combined with statistical compression of read
identifiers, quality scores, alignment information, and sequences, effectively
collapsing very large datasets to less than 15% of their original size with no
loss of information.
  Availability: Quip is freely available under the BSD license from
http://cs.washington.edu/homes/dcjones/quip.


  We develop an analytical method to calculate encounter times of two random
walkers in one dimension when each individual is segregated in its own spatial
domain and shares with its neighbor only a fraction of the available space,
finding very good agreement with numerically-exact calculations. We model a
population of susceptible and infected territorial individuals with this
spatial arrangement, and which may transmit an epidemic when they meet. We
apply the results on encounter times to determine analytically the macroscopic
propagation speed of the epidemic as a function of the microscopic
characteristics: the confining geometry, the animal diffusion constant, and the
infection transmission probability.


  We present a sequence-based probabilistic formalism that directly addresses
co-operative effects in networks of interacting positions in proteins,
providing significantly improved contact prediction, as well as accurate
quantitative prediction of free energy changes due to non-additive effects of
multiple mutations. In addition to these practical considerations, the
agreement of our sequence-based calculations with experimental data for both
structure and stability demonstrates a strong relation between the statistical
distribution of protein sequences produced by natural evolutionary processes,
and the thermodynamic stability of the structures to which these sequences
fold.


  Recent experimental results with humans involved in social dilemma games
suggest that cooperation may be a contagious phenomenon and that the selection
pressure operating on evolutionary dynamics (i.e., mimicry) is relatively weak.
I propose an evolutionary dynamics model that links these experimental findings
and evolution of cooperation. By assuming a small fraction of (imperfect)
zealous cooperators, I show that a large fraction of cooperation emerges in
evolutionary dynamics of social dilemma games. Even if defection is more
lucrative than cooperation for most individuals, they often mimic cooperation
of fellows unless the selection pressure is very strong. Then, zealous
cooperators can transform the population to be even fully cooperative under
standard evolutionary dynamics.


  This paper is concerned with front-like entire solutions for monostable
reactiondiffusion systems with cooperative and non-cooperative nonlinearities.
In the cooperative case, the existence and asymptotic behavior of spatially
independent solutions (SIS) are first proved. Combining a SIS and traveling
fronts with different wave speeds and directions, the existence and various
qualitative properties of entire solutions are then established using
comparison principle. In the non-cooperative case, we introduce two auxiliary
cooperative systems and establish some comparison arguments for the three
systems. The existence of entire solutions is then proved via the traveling
fronts and SIS of the auxiliary systems. Our results are applied to some
biological and epidemiological models. To the best of our knowledge, it is the
first work to study the entire solutions of non-cooperative reaction-diffusion
systems.


  In the present paper the macroscopic limits of the kinetic model for
inter-acting entities (individuals, organisms, cells) are studied. The kinetic
model is one-dimensional and entities are characterized by their position and
orientation (+/-) with swarming interaction controlled by the sensitivity
parameter. The macroscopic limits of the model are considered for solutions
close either to the diffusive (isotropic) or to the aligned (swarming)
equilibrium states for various sensitivity parameters. In the former case the
classical linear difusion equation results whereas in the latter a traveling
wave solution does both in the zeroth (`Euler') and frst (`Navier-Stokes')
order of approximation.


  Many transcription factors bind to DNA with a remarkable lack of specificity,
so that regulatory binding sites compete with an enormous number of
non-regulatory 'decoy' sites. For an auto-regulated gene, we show decoy sites
decrease noise in the number of unbound proteins to a Poisson limit that
results from binding and unbinding. This noise buffering is optimized for a
given protein concentration when decoys have a 1/2 probability of being
occupied. Decoys linearly increase the time to approach steady state and
exponentially increase the time to switch epigenetically between bistable
states.


  A central biological question is how natural organisms are so evolvable
(capable of quickly adapting to new environments). A key driver of evolvability
is the widespread modularity of biological networks--their organization as
functional, sparsely connected subunits--but there is no consensus regarding
why modularity itself evolved. While most hypotheses assume indirect selection
for evolvability, here we demonstrate that the ubiquitous, direct selection
pressure to reduce the cost of connections between network nodes causes the
emergence of modular networks. Experiments with selection pressures to maximize
network performance and minimize connection costs yield networks that are
significantly more modular and more evolvable than control experiments that
only select for performance. These results will catalyze research in numerous
disciplines, including neuroscience, genetics and harnessing evolution for
engineering purposes.


  We adopt an operatorial method based on the so-called creation, annihilation
and number operators in the description of different systems in which two
populations interact and move in a two-dimensional region. In particular, we
discuss diffusion processes modeled by a quadratic hamiltonian. This general
procedure will be adopted, in particular, in the description of migration
phenomena. With respect to our previous analogous results, we use here
fermionic operators since they automatically implement an upper bound for the
population densities.


  Biological structure and function depend on complex regulatory interactions
between many genes. A wealth of gene expression data is available from
high-throughput genome-wide measurement technologies, but effective gene
regulatory network inference methods are still needed. Model-based methods
founded on quantitative descriptions of gene regulation are among the most
promising, but many such methods rely on simple, local models or on ad hoc
inference approaches lacking experimental interpretability. We propose an
experimental design and develop an associated statistical method for inferring
a gene network by learning a standard quantitative, interpretable, predictive,
biophysics-based ordinary differential equation model of gene regulation. We
fit the model parameters using gene expression measurements from perturbed
steady-states of the system, like those following overexpression or knockdown
experiments. Although the original model is nonlinear, our design allows us to
transform it into a convex optimization problem by restricting attention to
steady-states and using the lasso for parameter selection. Here, we describe
the model and inference algorithm and apply them to a synthetic six-gene
system, demonstrating that the model is detailed and flexible enough to account
for activation and repression as well as synergistic and self-regulation, and
the algorithm can efficiently and accurately recover the parameters used to
generate the data.


  We investigate and quantify salient features of the charge distributions on
viral capsids. Our analysis combines the experimentally determined capsid
geometry with simple models for ionization of amino acids, thus yielding the
detailed description of spatial distribution for positive and negative charge
across the capsid wall. The obtained data is processed in order to extract the
mean radii of distributions, surface charge densities and dipole moment
densities. The results are evaluated and examined in light of previously
proposed models of capsid charge distributions, which are shown to have to some
extent limited value when applied to real viruses.


  A random network model which allows for tunable, quite general forms of
clustering, degree correlation and degree distribution is defined. The model is
an extension of the configuration model, in which stubs (half-edges) are paired
to form a network. Clustering is obtained by forming small completely connected
subgroups, and positive (negative) degree correlation is obtained by connecting
a fraction of the stubs with stubs of similar (dissimilar) degree. An SIR
(Susceptible -> Infective -> Recovered) epidemic model is defined on this
network. Asymptotic properties of both the network and the epidemic, as the
population size tends to infinity, are derived: the degree distribution, degree
correlation and clustering coefficient, as well as a reproduction number $R_*$,
the probability of a major outbreak and the relative size of such an outbreak.
The theory is illustrated by Monte Carlo simulations and numerical examples.
The main findings are that clustering tends to decrease the spread of disease,
the effect of degree correlation is appreciably greater when the disease is
close to threshold than when it is well above threshold and disease spread
broadly increases with degree correlation $\rho$ when $R_*$ is just above its
threshold value of one and decreases with $\rho$ when $R_*$ is well above one.


  Helically coiled filaments are a frequent motif in nature. In situations
commonly encountered in experiments coiled helices are squeezed flat onto two
dimensional surfaces. Under such 2-D confinement helices form "squeelices" -
peculiar squeezed conformations often resembling looped waves, spirals or
circles. Using theory and Monte-Carlo simulations we illuminate here the
mechanics and the unusual statistical mechanics of confined helices and show
that their fluctuations can be understood in terms of moving and interacting
discrete particle-like entities - the "twist-kinks". We show that confined
filaments can thermally switch between discrete topological twist quantized
states, with some of the states exhibiting dramatically enhanced
circularization probability while others displaying surprising
hyperflexibility.


  We introduce a sequence-dependent parametrization for a coarse-grained DNA
model [T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, J. Chem. Phys. 134,
085101 (2011)] originally designed to reproduce the properties of DNA molecules
with average sequences. The new parametrization introduces sequence-dependent
stacking and base-pairing interaction strengths chosen to reproduce the melting
temperatures of short duplexes. By developing a histogram reweighting
technique, we are able to fit our parameters to the melting temperatures of
thousands of sequences. To demonstrate the flexibility of the model, we study
the effects of sequence on: (a) the heterogeneous stacking transition of single
strands, (b) the tendency of a duplex to fray at its melting point, (c) the
effects of stacking strength in the loop on the melting temperature of
hairpins, (d) the force-extension properties of single strands and (e) the
structure of a kissing-loop complex. Where possible we compare our results with
experimental data and find a good agreement. A simulation code called oxDNA,
implementing our model, is available as free software.


  A critical step in the bacteriophage life cycle is genome ejection into host
bacteria. The ejection process for double-stranded DNA phages has been studied
thoroughly \textit{in vitro}, where after triggering with the cellular receptor
the genome ejects into a buffer. The experimental data have been interpreted in
terms of the decrease in free energy of the densely packed DNA associated with
genome ejection. Here we detail a simple model of genome ejection in terms of
the hydrostatic and osmotic pressures inside the phage, a bacterium, and a
buffer solution/culture medium. We argue that the hydrodynamic flow associated
with the water movement from the buffer solution into the phage capsid and
further drainage into the bacterial cytoplasm, driven by the osmotic gradient
between the bacterial cytoplasm and culture medium, provides an alternative
mechanism for phage genome ejection \textit{in vivo}; the mechanism is
perfectly consistent with phage genome ejection \textit{in vitro}.


  Reaction-diffusion systems may lead to the formation of steady state
heterogeneous spatial patterns, known as Turing patterns. Their mathematical
formulation is important for the study of pattern formation in general and play
central roles in many fields of biology, such as ecology and morphogenesis. In
the present study we focus on the role of Turing patterns in describing the
abundance distribution of predator and prey species distributed in patches in a
scale free network structure. We extend the original model proposed by Nakao
and Mikhailov by considering food chains with several interacting pairs of
preys and predators. We identify patterns of species distribution displaying
high degrees of apparent competition driven by Turing instabilities. Our
results provide further indication that differences in abundance distribution
among patches may be, at least in part, due to self organized Turing patterns,
and not necessarily to intrinsic environmental heterogeneity.


  The nonlinearities found in molecular networks usually prevent mathematical
analysis of network behaviour, which has largely been studied by numerical
simulation. This can lead to difficult problems of parameter determination.
However, molecular networks give rise, through mass-action kinetics, to
polynomial dynamical systems, whose steady states are zeros of a set of
polynomial equations. These equations may be analysed by algebraic methods, in
which parameters are treated as symbolic expressions whose numerical values do
not have to be known in advance. For instance, an "invariant" of a network is a
polynomial expression on selected state variables that vanishes in any steady
state. Invariants have been found that encode key network properties and that
discriminate between different network structures. Although invariants may be
calculated by computational algebraic methods, such as Gr\"obner bases, these
become computationally infeasible for biologically realistic networks. Here, we
exploit Chemical Reaction Network Theory (CRNT) to develop an efficient
procedure for calculating invariants that are linear combinations of
"complexes", or the monomials coming from mass action. We show how this
procedure can be used in proving earlier results of Horn and Jackson and of
Shinar and Feinberg for networks of deficiency at most one. We then apply our
method to enzyme bifunctionality, including the bacterial EnvZ/OmpR osmolarity
regulator and the mammalian
6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase glycolytic regulator,
whose networks have deficiencies up to four. We show that bifunctionality leads
to different forms of concentration control that are robust to changes in
initial conditions or total amounts. Finally, we outline a systematic procedure
for using complex-linear invariants to analyse molecular networks of any
deficiency.


  An intriguing unanswered question about the evolution of bilateral animals
with internal skeletons is how an internal skeleton evolved in the first place.
Computational modeling of the development of bilateral symmetric organisms
suggests an answer to this question. Our hypothesis is that an internal
skeleton may have evolved from a bilaterally symmetric ancestor with an
external skeleton. By growing the organism inside-out an external skeleton
becomes an internal skeleton. Our hypothesis is supported by a computational
theory of bilateral symmetry that allows us to model and simulate this process.
Inside-out development is achieved by an orientation switch. Given the
development of two bilateral founder cells that generate a bilateral organism,
a mutation that reverses the internal mirror orientation of those bilateral
founder cells leads to inside-out development. The new orientation is
epigenetically inherited by all progeny. A key insight is that each cell
contained in the newly evolved organism with the internal skeleton develops
according to the very same downstream developmental control network that
directs the development of its exoskeletal ancestor. The networks and their
genomes are are identical, but the interpretation is different because of the
cell's inverted orientation. The result is inside-out bilateral symmetric
development generating an inside-out organism with an internal skeleton.


  Complex systems, in many different scientific sectors, show coarse-grain
properties with simple growth laws with respect to fundamental microscopic
algorithms. We propose a classification scheme of growth laws which includes
human aging, tumor (and/or tissue) growth, logistic and generalized logistic
growth and the aging of technical devices. The proposed classification permits
to evaluate the aging/failure of combined new bio-technical "manufactured
products", where part of the system evolves in time according to
biological-mortality laws and part according to technical device behaviors.
Moreover it suggests a direct relation between the mortality leveling-off for
humans and technical devices and the observed small cure probability for large
tumors.


  Under certain circumstances such as lack of information or bounded
rationality, human players can take decisions on which strategy to choose in a
game on the basis of simple opinions. These opinions can be modified after each
round by observing own or others payoff results but can be also modified after
interchanging impressions with other players. In this way, the update of the
strategies can become a question that goes beyond simple evolutionary rules
based on fitness and become a social issue. In this work, we explore this
scenario by coupling a game with an opinion dynamics model. The opinion is
represented by a continuous variable that corresponds to the certainty of the
agents respect to which strategy is best. The opinions transform into actions
by making the selection of an strategy a stochastic event with a probability
regulated by the opinion. A certain regard for the previous round payoff is
included but the main update rules of the opinion are given by a model inspired
in social interchanges. We find that the dynamics fixed points of the coupled
model is different from those of the evolutionary game or the opinion models
alone. Furthermore, new features emerge such as the resilience of the fraction
of cooperators to the topology of the social interaction network or to the
presence of a small fraction of extremist players.


  Recent studies on epistatic networks of model organisms have unveiled a
certain type of modular property called monochromaticity in which the networks
are clusterable into functional modules that interact with each other through
the same type of epistasis. Here we propose and study three epistatic network
models that are inspired by the Duplication-Divergence mechanism to gain
insight into the evolutionary basis of monochromaticity and to test if it can
be explained as the outcome of a neutral evolutionary hypothesis. We show that
the epistatic networks formed by these stochastic evolutionary models have
monochromaticity conflict distributions that are centered close to zero and are
statistically significantly different from their randomized counterparts. In
particular, the last model we propose yields a strictly monochromatic solution.
Our results agree with the monochromaticity findings in real organisms and
point toward the possible role of a neutral mechanism in the evolution of this
phenomenon.


  The recent genealogical history of human populations is a complex mosaic
formed by individual migration, large-scale population movements, and other
demographic events. Population genomics datasets can provide a window into this
recent history, as rare traces of recent shared genetic ancestry are detectable
due to long segments of shared genomic material. We make use of genomic data
for 2,257 Europeans (the POPRES dataset) to conduct one of the first surveys of
recent genealogical ancestry over the past three thousand years at a
continental scale. We detected 1.9 million shared genomic segments, and used
the lengths of these to infer the distribution of shared ancestors across time
and geography. We find that a pair of modern Europeans living in neighboring
populations share around 10-50 genetic common ancestors from the last 1500
years, and upwards of 500 genetic ancestors from the previous 1000 years. These
numbers drop off exponentially with geographic distance, but since genetic
ancestry is rare, individuals from opposite ends of Europe are still expected
to share millions of common genealogical ancestors over the last 1000 years.
There is substantial regional variation in the number of shared genetic
ancestors: especially high numbers of common ancestors between many eastern
populations likely date to the Slavic and/or Hunnic expansions, while much
lower levels of common ancestry in the Italian and Iberian peninsulas may
indicate weaker demographic effects of Germanic expansions into these areas
and/or more stably structured populations. Recent shared ancestry in modern
Europeans is ubiquitous, and clearly shows the impact of both small-scale
migration and large historical events. Population genomic datasets have
considerable power to uncover recent demographic history, and will allow a much
fuller picture of the close genealogical kinship of individuals across the
world.


  The direct detection of haplotypes from short-read DNA sequencing data
requires changes to existing small-variant detection methods. Here, we develop
a Bayesian statistical framework which is capable of modeling multiallelic loci
in sets of individuals with non-uniform copy number. We then describe our
implementation of this framework in a haplotype-based variant detector,
FreeBayes.


  Coarse-grained (lattice-) models have a long tradition in aiding efforts to
decipher the physical or biological complexity of proteins. Despite the
simplicity of these models, however, numerical simulations are often
computationally very demanding and the quest for efficient algorithms is as old
as the models themselves. Expanding on our previous work [T. W\"ust and D. P.
Landau, Phys. Rev. Lett. 102, 178101 (2009)], we present a complete picture of
a Monte Carlo method based on Wang-Landau sampling in combination with
efficient trial moves (pull, bond-rebridging and pivot moves) which is
particularly suited to the study of models such as the hydrophobic-polar (HP)
lattice model of protein folding. With this generic and fully blind Monte Carlo
procedure, all currently known putative ground states for the most difficult
benchmark HP sequences could be found. For most sequences we could also
determine the entire energy density of states and, together with suitably
designed structural observables, explore the thermodynamics and intricate
folding behavior in the virtually inaccessible low-temperature regime. We
analyze the differences between random and protein-like heteropolymers for
sequence lengths up to 500 residues. Our approach is powerful both in terms of
robustness and speed, yet flexible and simple enough for the study of many
related problems in protein folding.


  We argue for a convergence of crystallography, materials science and biology,
that will come about through asking materials questions about biology and
biological questions about materials, illuminated by considerations of
information. The complex structures now being studied in biology and produced
in nanotechnology have outstripped the framework of classical crystallography,
and a variety of organizing concepts are now taking shape into a more modern
and dynamic science of structure, form and function. Absolute stability and
equilibrium are replaced by metastable structures existing in a flux of
energy-carrying information and moving within an energy landscape of complex
topology. Structures give place to processes and processes to systems. The
fundamental level is that of atoms. As smaller and smaller groups of atoms are
used for their physical properties, quantum effects become important; already
we see quantum computation taking shape. Concepts move towards those in life
with the emergence of specifically informational structures. We now see the
possibility of the artificial construction of a synthetic living system,
different from biological life, but having many or all of the same properties.
Interactions are essentially nonlinear and collective. Structures begin to have
an evolutionary history with episodes of symbiosis. Underlying all the
structures are constraints of time and space. Through hierarchization, a more
general principle than the periodicity of crystals, structures may be found
within structures on different scales. We must integrate unifying concepts from
dynamical systems and information theory to form a coherent language and
science of shape and structure beyond crystals. To this end, we discuss the
idea of categorizing structures based on information according to the
algorithmic complexity of their assembly.


  Incomplete lineage sorting (ILS) is a common source of gene tree incongruence
in multilocus analyses. A large number of methods have been developed to infer
species trees in the presence of ILS. Here we provide a mathematical analysis
of several coalescent-based methods. Our analysis is performed on a three-taxon
species tree and assumes that the gene trees are correctly reconstructed along
with their branch lengths.


  Haplotypes, the global patterns of DNA sequence variation, have important
implications for identifying complex traits. Recently, blocks of limited
haplotype diversity have been discovered in human chromosomes, intensifying the
research on modelling the block structure as well as the transitions or
co-occurrence of the alleles in these blocks as a way to compress the
variability and infer the associations more robustly. The haplotype block
structure analysis is typically complicated by the fact that the phase
information for each SNP is missing, i.e., the observed allele pairs are not
given in a consistent order across the sequence. The techniques for
circumventing this require additional information, such as family data, or a
more complex sequencing procedure. In this paper we present a hierarchical
statistical model and the associated learning and inference algorithms that
simultaneously deal with the allele ambiguity per locus, missing data, block
estimation, and the complex trait association. While the blo structure may
differ from the structures inferred by other methods, which use the pedigree
information or previously known alleles, the parameters we estimate, including
the learned block structure and the estimated block transitions per locus,
define a good model of variability in the set. The method is completely
datadriven and can detect Chron's disease from the SNP data taken from the
human chromosome 5q31 with the detection rate of 80% and a small error
variance.


  The classical Luria-Delbr\"uck model for fluctuation analysis is extended to
the case where cells can either divide or die at the end of their generation
time. This leads to a family of probability distributions generalizing the
Luria-Delbr\"uck family, and depending on three parameters: the expected number
of mutations, the relative fitness of normal cells compared to mutants, and the
death probability of mutants. The probabilistic treatment is similar to that of
the classical case; simulation and computing algorithms are provided. The
estimation problem is discussed: if the death probability is known, the two
other parameters can be reliably estimated. If the death probability is
unknown, the model can be identified only for large samples.


  To enhance the accuracy of protein-protein interaction function prediction, a
2-order graphic neighbor information feature extraction method based on
undirected simple graph is proposed in this paper, which extends the 1-order
graphic neighbor featureextraction method. And the chi-square test statistical
method is also involved in feature combination. To demonstrate the
effectiveness of our 2-order graphic neighbor feature, four logistic regression
models (logistic regression (abbrev. LR), diffusion kernel logistic regression
(abbrev. DKLR), polynomial kernel logistic regression (abbrev. PKLR), and
radial basis function (RBF) based kernel logistic regression (abbrev. RBF KLR))
are investigated on the two feature sets. The experimental results of protein
function prediction of Yeast Proteome Database (YPD) using the the
protein-protein interaction data of Munich Information Center for Protein
Sequences (MIPS) show that 2-order graphic neighbor information of proteins can
significantly improve the average overall percentage of protein function
prediction especially with RBF KLR. And, with a new 5-top chi-square feature
combination method, RBF KLR can achieve 99.05% average overall percentage on
2-order neighbor feature combination set.


  Gene transcriptional regulatory is an inherently noisy process. In this
paper, the study of fluctuations in a gene transcriptional regulatory system is
extended to the case of L\'evy noise, a kind of non-Gaussian noises which can
describe unpredictable jump changes of the random environment. The stationary
probability density is given to explore the key roles of L\'evy noise in the
gene regulatory networks. The results demonstrate that the parameters of L\'evy
noise, including noise intensity, stability index and skewness parameter can
induce switches between distinct gene-expression states. A further concern is
the switching time (from the high concentration state to the low concentration
one or from the low concentration state to the high concentration one), which
is a random variable and often referred to as the mean first passage time. The
effects of L\'evy noise in expression time and degradation time is studied by
computing the mean first passage time in two directions and a number of
different peculiarities of non-Gaussian L\'evy noise compared with Gaussian
noise are observed.


  Using all atom molecular dynamics simulations, we report spontaneous
unzipping and strong binding of small interfering RNA (siRNA) on graphene. Our
dispersion corrected density functional theory based calculations suggest that
nucleosides of RNA have stronger attractive interactions with graphene as
compared to DNA residues. These stronger interactions force the double stranded
siRNA to spontaneously unzip and bind to the graphene surface. Unzipping always
nucleates at one end of the siRNA and propagates to the other end after few
base-pairs get unzipped. While both the ends get unzipped, the middle part
remains in double stranded form because of torsional constraint. Unzipping
probability distributions fitted to single exponential function give unzipping
time (t) of the order of few nanoseconds which decrease exponentially with
temperature. From the temperature variation of unzipping time we estimate the
energy barrier to unzipping.


  The directed polymerization of actin networks is an essential element of many
biological processes, including cell migration. Different theoretical models
considering the interplay between the underlying processes of polymerization,
capping and branching have resulted in conflicting predictions. One of the main
reasons for this discrepancy is the assumption of a branching reaction that is
either first order (autocatalytic) or zeroth order in the number of existing
filaments. Here we introduce a unifying framework from which the two
established scenarios emerge as limiting cases for low and high filament
number. A smooth transition between the two cases is found at intermediate
conditions. We also derive a threshold for the capping rate, above which
autocatalytic growth is predicted at sufficiently low filament number. Below
the threshold, zeroth order characteristics are predicted to dominate the
dynamics of the network for all accessible filament numbers. Together, this
allows cells to grow stable actin networks over a large range of different
conditions.


  We study a stochastic linear discrete metapolulation model to understand the
effect of risk spreading by dispersion. We calculate analytically the stable
distribution of populations that live in different habitats. The result shows
that the simultaneous distribution of the populations has a complicated
self-similar structure, but a population at each habitat follows a log-normal
distribution.


  We consider a fitness-structured population model with competition and
migration between nearest neighbors. Under a combination of large population
and rare migration limits we are particularly interested in the asymptotic
behavior of the total population partition on supporting trait sites. For the
population without mutation on a finite trait space we obtain the equilibrium
configuration and characterize the right time scale for fixation. For the model
with mutation on an infinite trait space a jump process-trait substitution tree
model is established on a rarer mutation time scale against the rare migration
constrained in terms of a large population limit. Due to a change of the
fitness landscape provoked by a new mutant, some temporarily unfit types can be
recovered from time to time. In the end we shed light to illustrate sexual
reproduction in a diploid population on the genetic level.


  Protein mean and variance levels in a simple stochastic gene expression
circuit are controlled using proportional integral feedback. It is shown that
the protein mean level can be globally and robustly tracked to any desired
value using a simple PI controller that satisfies explicit sufficient
conditions. Controlling both the mean and variance on the other hand requires
the use of an additional control input, chosen here as the mRNA degradation
rate. Local robust tracking of mean and variance is proved to be achievable
using multivariable PI control, provided that the reference point satisfies
necessary conditions imposed by the system. Even more importantly, it is shown
that there exist PI controllers that locally, robustly and simultaneously
stabilize all the equilibrium points inside the admissible region. Simulation
examples illustrate the results.


  Experiments on the ultimatum game have revealed that humans are remarkably
fond of fair play. When asked to share an amount of money, unfair offers are
rare and their acceptance rate small. While empathy and spatiality may lead to
the evolution of fairness, thus far considered continuous strategies have
precluded the observation of solutions that would be driven by pattern
formation. Here we introduce a spatial ultimatum game with discrete strategies,
and we show that this simple alteration opens the gate to fascinatingly rich
dynamical behavior. Besides mixed stationary states, we report the occurrence
of traveling waves and cyclic dominance, where one strategy in the cycle can be
an alliance of two strategies. The highly webbed phase diagram, entailing
continuous and discontinuous phase transitions, reveals hidden complexity in
the pursuit of human fair play.


  A key goal of systems biology is the predictive mathematical description of
gene regulatory circuits. Different approaches are used such as deterministic
and stochastic models, models that describe cell growth and division explicitly
or implicitly etc. Here we consider simple systems of unregulated
(constitutive) gene expression and compare different mathematical descriptions
systematically to obtain insight into the errors that are introduced by various
common approximations such as describing cell growth and division by an
effective protein degradation term. In particular, we show that the population
average of protein content of a cell exhibits a subtle dependence on the
dynamics of growth and division, the specific model for volume growth and the
age structure of the population. Nevertheless, the error made by models with
implicit cell growth and division is quite small. Furthermore, we compare
various models that are partially stochastic to investigate the impact of
different sources of (intrinsic) noise. This comparison indicates that
different sources of noise (protein synthesis, partitioning in cell division)
contribute comparable amounts of noise if protein synthesis is not or only
weakly bursty. If protein synthesis is very bursty, the burstiness is the
dominant noise source, independent of other details of the model. Finally, we
discuss two sources of extrinsic noise: cell-to-cell variations in protein
content due to cells being at different stages in the division cycles, which we
show to be small (for the protein concentration and, surprisingly, also for the
protein copy number per cell) and fluctuations in the growth rate, which can
have a significant impact.


  Gene expression is inherently noisy as many steps in the read-out of the
genetic information are stochastic. To disentangle the effect of different
sources of stochasticity in such systems, we consider various models that
describe some processes as stochastic and others as deterministic. We review
earlier results for unregulated (constitutive) gene expression and present new
results for a gene controlled by negative autoregulation with cell growth
modeled by linear volume growth.


  The Conservation of Energy plays a pivotal part in the development of the
physical sciences. With the growth of computation and the study of other
discrete token based systems such as the genome, it is useful to ask if there
are conservation principles which apply to such systems and what kind of
functional behaviour they imply for such systems.
  Here I propose that the Conservation of Hartley-Shannon Information plays the
same over-arching role in discrete token based systems as the Conservation of
Energy does in physical systems. I will go on to prove that this implies
power-law behaviour in component sizes in software systems no matter what they
do or how they were built, and also implies the constancy of average gene
length in biological systems as reported for example by Lin Xu et al
(10.1093/molbev/msk019).
  These propositions are supported by very large amounts of experimental data
extending the first presentation of these ideas in Hatton (2011, IFIP / SIAM /
NIST Working Conference on Uncertainty Quantification in Scientific Computing,
Boulder, August 2011).


  The branching structure of biological evolution confers statistical
dependencies on phenotypic trait values in related organisms. For this reason,
comparative macroevolutionary studies usually begin with an inferred phylogeny
that describes the evolutionary relationships of the organisms of interest. The
probability of the observed trait data can be computed by assuming a model for
trait evolution, such as Brownian motion, over the branches of this fixed tree.
However, the phylogenetic tree itself contributes statistical uncertainty to
estimates of other evolutionary quantities, and many comparative evolutionary
biologists regard the tree as a nuisance parameter. In this paper, we present a
framework for analytically integrating over unknown phylogenetic trees in
comparative evolutionary studies by assuming that the tree arises from a
continuous-time Markov branching model called the Yule process. To do this, we
derive a closed-form expression for the distribution of phylogenetic diversity,
which is the sum of branch lengths connecting a set of taxa. We then present a
generalization of phylogenetic diversity which is equivalent to the expected
trait disparity in a set of taxa whose evolutionary relationships are generated
by a Yule process and whose traits evolve by Brownian motion. We derive
expressions for the distribution of expected trait disparity under a Yule tree.
Given one or more observations of trait disparity in a clade, we perform fast
likelihood-based estimation of the Brownian variance for unresolved clades. Our
method does not require simulation or a fixed phylogenetic tree. We conclude
with a brief example illustrating Brownian rate estimation for thirteen
families in the Mammalian order Carnivora, in which the phylogenetic tree for
each family is unresolved.


  The objective of the present study was to develop an understanding of short
single-stranded DNA (ssDNA) to assist the development of new DNA-based
biosensors. A ssDNA model containing twelve bases was constructed from the
130-145 codon sequence of the p53 gene. Various thermodynamic macroscopic
observables such as temperature, energy distributions, as well as root mean
square deviation (RMSD) of the nucleic acid backbone of the ssDNA were studied
using molecular dynamics (MD) simulations. The AMBER program was used for
building the structural model of the ssDNA, and atomistic MD simulations in
three different ensembles were carried out using the NAMD program. The
microcanonical (NVE), conical (NVT) and isobaric-isothermal (NPT) ensembles
were employed to compare the equilibrium characteristics of ssDNA in aqueous
solutions. Our results indicate that the conformational stability of the ssDNA
is dependent on the thermodynamic conditions.


  Motivation: Next Generation Sequencing technologies revolutionized many
fields in biology by enabling the fast and cheap sequencing of large amounts of
genomic data. The ever increasing sequencing capacities enabled by current
sequencing machines hold a lot of promise as for the future applications of
these technologies, but also create increasing computational challenges related
to the analysis and storage of these data. A typical sequencing data file may
occupy tens or even hundreds of gigabytes of disk space, prohibitively large
for many users. Raw sequencing data consists of both the DNA sequences (reads)
and per-base quality values that indicate the level of confidence in the
readout of these sequences. Quality values account for about half of the
required disk space in the commonly used FASTQ format and therefore their
compression can significantly reduce storage requirements and speed up analysis
and transmission of these data.
  Results: In this paper we present a framework for the lossy compression of
the quality value sequences of genomic read files. Numerical experiments with
reference based alignment using these quality values suggest that we can
achieve significant compression with little compromise in performance for
several downstream applications of interest, as is consistent with our
theoretical analysis. Our framework also allows compression in a regime - below
one bit per quality value - for which there are no existing compressors.


  The mechanisms of solid-state electron transport (ETp) via a monolayer of
immobilized Azurin (Az) was examined by conducting probe atomic force
microscopy (CP-AFM), both as function of temperature (248 - 373K) and of
applied tip force (6-12 nN). By varying both temperature and force in CP-AFM,
we find that the ETp mechanism can alter with a change in the force applied via
the tip to the proteins. As the applied force increases, ETp via Az changes
from temperature-independent to thermally activated at high temperatures. This
is in contrast to the Cu-depleted form of Az (apo-Az), where increasing the
applied force causes only small quantitative effects, that fit with a decrease
in electrode spacing. At low force ETp via holo-Az is temperature-independent
and thermally activated via apo-Az. This observation agrees with
macroscopic-scale measurements, thus confirming that the difference in ETp
dependence on temperature between holo- and apo-Az is an inherent one that may
reflect a difference in rigidity between the two forms. An important
implication of these results, which depend on CP-AFM measurements over a
significant temperature range, is that for ETp measurements on floppy systems,
such as proteins, the stress applied to the sample should be kept constant or,
at least controlled during measurement.


  Measuring electron transport (ETp) across proteins in the solid-state offers
a way to study electron transfer (ET) mechanism(s) that minimizes solvation
effects on the process. Solid state ETp is sensitive to any static
(conformational) or dynamic (vibrational) changes in the protein. Our
macroscopic measurement technique extends the use of ETp meas-urements down to
low temperatures and the concomitant lower current densities, because the
larger area still yields measurable currents. Thus, we reported previously a
surprising lack of temperature-dependence for ETp via the blue copper protein
azurin (Az), from 80K till denaturation, while ETp via apo-(Cu-free) Az was
found to be temperature de-pendent \geq 200K. H/D substitution (deuteration)
can provide a potentially powerful means to unravel factors that affect the ETp
mechanism at a molecular level. Therefore, we measured and report here the
kinetic deuterium isotope effect (KIE) on ETp through holo-Az as a function of
temperature (30-340K). We find that deuteration has a striking effect in that
it changes ETp from temperature independent to temperature dependent above
180K. This change is expressed in KIE values between 1.8 at 340K and 9.1 at
\leq 180K. These values are particularly remarkable in light of the previously
reported inverse KIE on the ET in Az in solution. The high values that we
obtain for the KIE on the ETp process across the protein monolayer are
consistent with a transport mechanism that involves
through-(H-containing)-bonds of the {\beta}-sheet structure of Az, likely those
of am-ide groups.


  Electrons can migrate via proteins over distances that are considered long
for non-conjugated systems. Proteins' nano-scale dimensions and the enormous
flexibility of their structures and chemistry makes them fascinating subjects
for investigating the mechanism of their electron transport (ETp) capacity. One
particular attractive research direction is that of tuning their ETp efficiency
by doping them with external small molecules. Here we report that solid-state
ETp across human serum albumin (HSA) increases by more than two orders of
magnitude upon retinoate (RA) binding to HSA. RA was chosen because optical
spectroscopy has provided evidence for the non-covalent binding of at least
three RA molecules to HSA and indications for their relative structural
positions. The temperature dependence of ETp shows that both the activation
energy and the distance-decay constant decrease with increasing RA binding to
HSA. Furthermore, the observed transition from temperature-activated ETp above
190K to temperature-independent ETp below this temperature suggests a change in
the ETp mechanism with temperature.


  In the last few decades, the human allosomes are engrossed in an intensive
attention among researchers. The allosomes are now already been sequenced and
found there are about 2000 and 78 genes in human X and Y chromosomes
respectively. The hemizygosity of the human X chromosome in males exposes
recessive disease alleles, and this phenomenon has prompted decades of
intensive study of X-linked disorders. By contrast, the small size of the human
Y chromosome, and its prominent long-arm heterochromatic region suggested
absence of function beyond sex determination. But the present problem is to
accomplish whether a given sequence of nucleotides i.e. a DNA is a Human X or Y
chromosomal genes or not, without any biological experimental support. In our
perspective, a proper quantitative understanding of these genes is required to
justify or nullify whether a given sequence is a Human X or Y chromosomal gene.
In this paper, some of the X and Y chromosomal genes have been quantified in
genomic and proteomic level through Fractal Geometric and Mathematical
Morphometric analysis. Using the proposed quantitative model, one can easily
make probable justification or deterministic nullification whether a given
sequence of nucleotides is a probable Human X or Y chromosomal gene or not,
without seeking any biological experiment. Of course, a further biological
experiment is essential to validate it as the probable Human X or Y chromosomal
gene homologue. This study would enable Biologists to understand these genes in
more quantitative manner instead of their qualitative features.


  We investigate the influences of the excluded volume of molecules on
biochemical reaction processes on 2-dimensional surfaces using a model of
signal transduction processes on biomembranes. We perform simulations of the
2-dimensional cell-based model, which describes the reactions and diffusion of
the receptors, signaling proteins, target proteins, and crowders on the cell
membrane. The signaling proteins are activated by receptors, and these
activated signaling proteins activate target proteins that bind autonomously
from the cytoplasm to the membrane, and unbind from the membrane if activated.
If the target proteins bind frequently, the volume fraction of molecules on the
membrane becomes so large that the excluded volume of the molecules for the
reaction and diffusion dynamics cannot be negligible. We find that such
excluded volume effects of the molecules induce non-trivial variations of the
signal flow, defined as the activation frequency of target proteins, as
follows. With an increase in the binding rate of target proteins, the signal
flow varies by i) monotonically increasing; ii) increasing then decreasing in a
bell-shaped curve; or iii) increasing, decreasing, then increasing in an
S-shaped curve. We further demonstrate that the excluded volume of molecules
influences the hierarchical molecular distributions throughout the reaction
processes. In particular, when the system exhibits a large signal flow, the
signaling proteins tend to surround the receptors to form receptor-signaling
protein clusters, and the target proteins tend to become distributed around
such clusters. To explain these phenomena, we analyze the stochastic model of
the local motions of molecules around the receptor.


  DNA has a well-defined structural transition -- the denaturation of its
double-stranded form into two single strands -- that strongly affects its
thermal transport properties. We show that, according to a widely implemented
model for DNA denaturation, one can engineer DNA "heattronic" devices that have
a rapidly increasing thermal conductance over a narrow temperature range across
the denaturation transition (~350 K). The origin of this rapid increase of
conductance, or "switching", is the softening of the lattice and suppression of
nonlinear effects as the temperature crosses the transition temperature and DNA
denatures. Most importantly, we demonstrate that DNA nanojunctions have a broad
range of thermal tunability due to varying the sequence and length, and
exploiting the underlying nonlinear behavior. We discuss the role of disorder
in the base sequence, as well as the relation to genomic DNA. These results set
the basis for developing thermal devices out of materials with nonlinear
structural dynamics, as well as understanding the underlying mechanisms of DNA
denaturation.


  Metabolism displays striking and robust regularities in the forms of
modularity and hierarchy, whose composition may be compactly described. This
renders metabolic architecture comprehensible as a system, and suggests the
order in which layers of that system emerged. Metabolism also serves as the
foundation in other hierarchies, at least up to cellular integration including
bioenergetics and molecular replication, and trophic ecology. The
recapitulation of patterns first seen in metabolism, in these higher levels,
suggests metabolism as a source of causation or constraint on many forms of
organization in the biosphere.
  We identify as modules widely reused subsets of chemicals, reactions, or
functions, each with a conserved internal structure. At the small molecule
substrate level, module boundaries are generally associated with the most
complex reaction mechanisms and the most conserved enzymes. Cofactors form a
structurally and functionally distinctive control layer over the small-molecule
substrate. Complex cofactors are often used at module boundaries of the
substrate level, while simpler ones participate in widely used reactions.
Cofactor functions thus act as "keys" that incorporate classes of organic
reactions within biochemistry.
  The same modules that organize the compositional diversity of metabolism are
argued to have governed long-term evolution. Early evolution of core
metabolism, especially carbon-fixation, appears to have required few
innovations among a small number of conserved modules, to produce adaptations
to simple biogeochemical changes of environment. We demonstrate these features
of metabolism at several levels of hierarchy, beginning with the small-molecule
substrate and network architecture, continuing with cofactors and key conserved
reactions, and culminating in the aggregation of multiple diverse physical and
biochemical processes in cells.


  Southern and eastern African populations that speak non-Bantu languages with
click consonants are known to harbour some of the most ancient genetic lineages
in humans, but their relationships are poorly understood. Here, we report data
from 23 populations analyzed at over half a million single nucleotide
polymorphisms, using a genome-wide array designed for studying human history.
The southern African Khoisan fall into two genetic groups, loosely
corresponding to the northwestern and southeastern Kalahari, which we show
separated within the last 30,000 years. We find that all individuals derive at
least a few percent of their genomes from admixture with non-Khoisan
populations that began approximately 1,200 years ago. In addition, the east
African Hadza and Sandawe derive a fraction of their ancestry from admixture
with a population related to the Khoisan, supporting the hypothesis of an
ancient link between southern and eastern Africa


  Indirect reciprocity is a mechanism for cooperation in social dilemma
situations, in which an individual is motivated to help another to acquire a
good reputation and receive help from others afterwards. Ingroup favoritism is
another aspect of human cooperation, whereby individuals help members in their
own group more often than those in other groups. Ingroup favoritism is a puzzle
for the theory of cooperation because it is not easily evolutionarily stable.
In the context of indirect reciprocity, ingroup favoritism has been shown to be
a consequence of employing a double standard when assigning reputations to
ingroup and outgroup members; e.g., helping an ingroup member is regarded as
good, whereas the same action toward an outgroup member is regarded as bad. We
analyze a model of indirect reciprocity in which information sharing is
conducted groupwise. In our model, individuals play social dilemma games within
and across groups, and the information about their reputations is shared within
each group. We show that evolutionarily stable ingroup favoritism emerges even
if all the players use the same reputation assignment rule regardless of group
(i.e., a single standard). Two reputation assignment rules called simple
standing and stern judging yield ingroup favoritism. Stern judging induces much
stronger ingroup favoritism than does simple standing. Simple standing and
stern judging are evolutionarily stable against each other when groups
employing different assignment rules compete and the number of groups is
sufficiently large. In addition, we analytically show as a limiting case that
homogeneous populations of reciprocators that use reputations are unstable when
individuals independently infer reputations of individuals, which is consistent
with previously reported numerical results.


  Human immunodeficiency virus (HIV-1 or simply HIV) induces a persistent
infection, which in the absence of treatment leads to AIDS and death in almost
all infected individuals. HIV infection elicits a vigorous immune response
starting about 2-3 weeks post infection that can lower the amount of virus in
the body, but which cannot eradicate the virus. How HIV establishes a chronic
infection in the face of a strong immune response remains poorly understood. It
has been shown that HIV is able to rapidly change its proteins via mutation to
evade recognition by virus-specific cytotoxic T lymphocytes (CTLs). Typically,
an HIV-infected patient will generate 4-12 CTL responses specific for parts of
viral proteins called epitopes. Such CTL responses lead to strong selective
pressure to change the viral sequences encoding these epitopes so as to avoid
CTL recognition. Here we review experimental data on HIV evolution in response
to CTL pressure, mathematical models developed to explain this evolution, and
highlight problems associated with the data and previous modeling efforts. We
show that estimates of the strength of the epitope-specific CTL response depend
on the method used to fit models to experimental data and on the assumptions
made regarding how mutants are generated during infection. We illustrate that
allowing CTL responses to decay over time may improve the fit to experimental
data and provides higher estimates of the killing efficacy of HIV-specific
CTLs. We also propose a novel method for simultaneously estimating the killing
efficacy of multiple CTL populations specific for different epitopes of HIV
using stochastic simulations. Lastly, we show that current estimates of the
efficacy at which HIV-specific CTLs clear virus-infected cells can be improved
by more frequent sampling of viral sequences and by combining data on sequence
evolution with experimentally measured CTL dynamics.


  Rationale: In a shotgun proteomics experiment with data-dependent
acquisition, real-time analysis of a precursor scan results in selection of a
handful of peaks for subsequent isolation, fragmentation and secondary
scanning. This peak selection protocol typically focuses on the most abundant
peaks in the precursor scan, while attempting to avoid re-sampling the same m/z
values in rapid succession. The protocol does not, however, incorporate
analysis of previous fragmentation scans into the peak selection procedure.
  Methods: In this work, we investigate the feasibility and utility of
incorporating analysis of previous fragmentation scans into the peak selection
protocol. We demonstrate that real-time identification of fragmentation spectra
is feasible in principle, and we investigate, via simulations, several
strategies to make use of the resulting peptide identifications during peak
selection.
  Results: Our simulations fail to provide evidence that peptide
identifications can provide a large improvement in the total number of peptides
identified by a shotgun proteomics experiment.
  Conclusions: These results are significant because they point out the
feasibility of using peptide identifications during peak selection, and because
our experiments may provide a starting point for others working in this
direction.


  Recent findings raise the possibility of PARP inhibitor therapy in colorectal
cancers(CRCs). However, the extent of PARP-1 protein expression in clinical
specimens of CRC is not known. Using immunohistochemistry we assessed PARP-1
protein expression in tissue microarrays of 151 CRCs and its association with
the patient's age, sex, Astler-Coller stage, grade and site of the tumor. High
PARP nuclear immunoreactivity was found in 68.2% (103/151) of all cases. In
turn, 31.8% (48/151)of tumors showed low PARP expression, including 9 (6%)
PARP-1 negative CRCs. There was a significant association of PARP-1 expression
with the site of CRC and Astler-Coller stage. A high PARP expression was noted
in 79.1% of colon vs. 53.9% of rectal tumors (p = 0.001). The mean PARP-1 score
was 1.27 times higher in colon vs. rectal cancers (p = 0.009) and it was higher
in stage B2 vs. stage C of CRCs (p = 0.018). In conclusion, the level of PARP-1
protein nuclear expression is associated with the tumor site and heterogeneous
across clinical specimens of CRC, with the majority of CRCs expressing a high
level but minority - low or no PARP-1 expression. These findings may have a
clinical significance because the assessment of PARP-1 expression in tumor
samples may improve selection of patients with CRC for PARP inhibitor therapy.


  Existing state-of-the-art methods that take a single RNA sequence and predict
the corresponding RNA secondary-structure are thermodynamic methods. These
predict the most stable RNA structure, but do not consider the process of
structure formation. We have by now ample experimental and theoretical
evidence, however, that sequences in vivo fold while being transcribed and that
the process of structure formation matters. We here present a conceptually new
method for predicting RNA secondary-structure, called CoFold, that combines
thermodynamic with kinetic considerations. Our method significantly improves
the state-of-art in terms of prediction accuracy, especially for long sequences
of more than a thousand nucleotides length such as ribosomal RNAs.


  This article characterizes certain small multistationary chemical reaction
networks. We consider the set of fully open networks, those for which all
chemical species participate in inflow and outflow, containing one non-flow
(reversible or irreversible) reaction. We show that such a network admits
multiple positive mass-action steady states if and only if the stoichiometric
coefficients in the non-flow reaction satisfy a certain simple arithmetic
relation. The multistationary fully open one-reaction networks are identified
with the chemical process of autocatalysis. Using the notion of `embedded
network' defined recently by Joshi and Shiu, we provide new sufficient
conditions for establishing multistationarity of fully open networks,
applicable well beyond the one-reaction setting.


  Here we present a new fixed parameter tractable algorithm to compute the
hybridization number r of two rooted, not necessarily binary phylogenetic trees
on taxon set X in time (6^r.r!).poly(n)$, where n=|X|. The novelty of this
approach is its use of terminals, which are maximal elements of a natural
partial order on X, and several insights from the softwired clusters
literature. This yields a surprisingly simple and practical bounded-search
algorithm and offers an alternative perspective on the underlying combinatorial
structure of the hybridization number problem.


  An extension of the Truscott-Brindley model (Bull. Math. Biol. {\bf 56}, 981
(1994)) is derived to account for the effect of demographic fluctuations. In
the presence of seasonal forcing, and sufficiently shallow water conditions,
the fluctuations induced by the discreteness of the zooplankton component
appear sufficient to cause switching between the bloom and no-bloom cycle
predicted at the mean-field level by the model.The destabilization persists in
the thermodynamic limit of a water basin infinitely extended in the horizontal
direction.


  In the framework of the three-party constrained voter model, where voters of
two radical parties (A and B) interact with "centrists" (C and Cz), we study
the competition between a persuasive majority and a committed minority. In this
model, A's and B's are incompatible voters that can convince centrists or be
swayed by them. Here, radical voters are more persuasive than centrists, whose
sub-population consists of susceptible agents C and a fraction zeta of centrist
zealots Cz. Whereas C's may adopt the opinions A and B with respective rates
1+delta_A and 1+delta_B (with delta_A>=delta_B>0), Cz's are committed
individuals that always remain centrists. Furthermore, A and B voters can
become (susceptible) centrists C with a rate 1. The resulting competition
between commitment and persuasion is studied in the mean field limit and for a
finite population on a complete graph. At mean field level, there is a
continuous transition from a coexistence phase when
zeta<Delta_c=delta_A/(1+delta_A) to a phase where centrism prevails when zeta>=
Delta_c. In a finite population of size N, demographic fluctuations lead to
centrism consensus and the dynamics is characterized by the mean consensus time
tau. Because of the competition between commitment and persuasion, here
consensus is reached much slower (zeta<Delta_c) or faster (zeta>=Delta_c) than
in the absence of zealots (when tau\simN). In fact, when zeta<Delta_c and there
is an initial minority of centrists, the mean consensus time asymptotically
grows as tau\simN^{-1/2} e^{N gamma}, where gamma is determined. The dynamics
is thus characterized by a metastable state where the most persuasive voters
and centrists coexist when delta_A>delta_B, whereas all species coexist when
delta_A=delta_B. When zeta>=Delta_c and the initial density of centrists is
low, one finds tau\simln N (when N>>1). Our analytical findings are
corroborated by stochastic simulations.


  We address the problem of the joint statistical inference of phylogenetic
trees and multiple sequence alignments from unaligned molecular sequences. This
problem is generally formulated in terms of string-valued evolutionary
processes along the branches of a phylogenetic tree. The classical evolutionary
process, the TKF91 model, is a continuous-time Markov chain model comprised of
insertion, deletion and substitution events. Unfortunately this model gives
rise to an intractable computational problem---the computation of the marginal
likelihood under the TKF91 model is exponential in the number of taxa. In this
work, we present a new stochastic process, the Poisson Indel Process (PIP), in
which the complexity of this computation is reduced to linear. The new model is
closely related to the TKF91 model, differing only in its treatment of
insertions, but the new model has a global characterization as a Poisson
process on the phylogeny. Standard results for Poisson processes allow key
computations to be decoupled, which yields the favorable computational profile
of inference under the PIP model. We present illustrative experiments in which
Bayesian inference under the PIP model is compared to separate inference of
phylogenies and alignments.


  A wide range of infectious diseases are both vertically and horizontally
transmitted. Such diseases are spatially transmitted via multiple species in
heterogeneous environments, typically described by complex meta-population
models. The reproduction number is a critical metric predicting whether the
disease can invade the meta-population system. This paper presents the
reproduction number for a generic disease vertically and horizontally
transmitted among multiple species in heterogeneous networks, where nodes are
locations, and links reflect outgoing or incoming movement flows. The
metapopulation model for vertically and horizontally transmitted diseases is
gradually formulated from two species, two-node network models. We derived an
explicit expression of the reproduction number, which is the spectral radius of
a matrix reduced in size with respect to the original next generation matrix.
The reproduction number is shown to be a function of vertical and horizontal
transmission parameters, and the lower bound is the reproduction number for
horizontal transmission. As an application, the reproduction number and its
bounds for the Rift Valley fever zoonosis, where livestock, mosquitoes, and
humans are the involved species are derived. By computing the reproduction
number for different scenarios through numerical simulations, we found the
reproduction number is affected by livestock movement rates only when
parameters are heterogeneous across nodes. To summarize, our study contributes
the reproduction number for vertically and horizontally transmitted diseases in
heterogeneous networks. This explicit expression is easily adaptable to
specific infectious diseases, affording insights into disease evolution.


  In experiments and in simulations, the free energy of a state of a system can
be determined from the probability that the state is occupied. However, it is
often necessary to impose a biasing potential on the system so that high energy
states are sampled with sufficient frequency. The unbiased energy is typically
obtained from the data using the weighted histogram analysis method (WHAM).
Here we present differential energy surface analysis (DESA), in which the
gradient of the energy surface, dE/dx, is extracted from data taken with a
series of harmonic biasing potentials. It is shown that DESA produces a maximum
likelihood estimate of the folding landscape gradient. DESA is demonstrated by
analyzing data from a simulated system as well as data from a single-molecule
unfolding experiment in which the end-to-end distance of a DNA hairpin is
measured. It is shown that the energy surface obtained from DESA is
indistinguishable from the energy surface obtained when WHAM is applied to the
same data. Two criteria are defined which indicate whether the DESA results are
self-consistent. It is found that these criteria can detect a situation where
the energy is not a single-valued function of the measured reaction coordinate.
The criteria were found to be satisfied for the experimental data analyzed,
confirming that end-to-end distance is a good reaction coordinate for the
experimental system. The combination of DESA and the optical trap assay in
which a structure is disrupted under harmonic constraint facilitates an
extremely accurate measurement of the folding energy surface.


  Single molecule experiments on B-DNA stretching have revealed one or two
structural transitions, when increasing the external force. They are
characterized by a sudden increase of DNA contour length and a decrease of the
bending rigidity. It has been proposed that the first transition, at forces of
60--80 pN, is a transition from B to S-DNA, viewed as a stretched duplex DNA,
while the second one, at stronger forces, is a strand peeling resulting in
single stranded DNAs (ssDNA), similar to thermal denaturation. But due to
experimental conditions these two transitions can overlap, for instance for
poly(dA-dT). We derive analytical formula using a coupled discrete worm like
chain-Ising model. Our model takes into account bending rigidity, discreteness
of the chain, linear and non-linear (for ssDNA) bond stretching. In the limit
of zero force, this model simplifies into a coupled model already developed by
us for studying thermal DNA melting, establishing a connexion with previous
fitting parameter values for denaturation profiles. We find that: (i) ssDNA is
fitted, using an analytical formula, over a nanoNewton range with only three
free parameters, the contour length, the bending modulus and the monomer size;
(ii) a surprisingly good fit on this force range is possible only by choosing a
monomer size of 0.2 nm, almost 4 times smaller than the ssDNA nucleobase
length; (iii) mesoscopic models are not able to fit B to ssDNA (or S to ss)
transitions; (iv) an analytical formula for fitting B to S transitions is
derived in the strong force approximation and for long DNAs, which is in
excellent agreement with exact transfer matrix calculations; (v) this formula
fits perfectly well poly(dG-dC) and $\lambda$-DNA force-extension curves with
consistent parameter values; (vi) a coherent picture, where S to ssDNA
transitions are much more sensitive to base-pair sequence than the B to S one,
emerges.


  We consider a stochastic evolutionary model for a phenotype developing
amongst n related species with unknown phylogeny. The unknown tree is modelled
by a Yule process conditioned on n contemporary nodes. The trait value is
assumed to evolve along lineages as an Ornstein-Uhlenbeck process. As a result,
the trait values of the n species form a sample with dependent observations. We
establish three limit theorems for the sample mean corresponding to three
domains for the adaptation rate. In the case of fast adaptation, we show that
for large $n$ the normalized sample mean is approximately normally distributed.
Using these limit theorems, we develop novel confidence interval formulae for
the optimal trait value.