  We investigate a game-theoretic model of a social system where both the rules
of the game and the interaction structure are shaped by the behavior of the
agents. We call this type of model, with several types of feedback couplings
from the behavior of the agents to their environment, a multiadaptive game. Our
model has a complex behavior with several regimes of different dynamic behavior
accompanied by different network topological properties. Some of these regimes
are characterized by heterogeneous, hierarchical interaction networks, where
cooperation and network topology coemerge from the dynamics.


  Eukaryotic DNA is packaged into chromatin: one-dimensional arrays of
nucleosomes separated by stretches of linker DNA are folded into 30-nm
chromatin fibers which in turn form higher-order structures. Each nucleosome,
the fundamental unit of chromatin, has 147 base pairs of DNA wrapped around a
histone octamer. In order to describe how chromatin fiber formation affects
nucleosome positioning and energetics, we have developed a thermodynamic model
of finite-size particles with effective nearest-neighbor interactions and
arbitrary DNA-binding energies. We show that both one- and two-body
interactions can be accurately extracted from one-particle density profiles
based on high-throughput maps of in vitro or in vivo nucleosome positions.
Although a simpler approach that neglects two-body interactions (even if they
are in fact present in the system) can be used to predict sequence determinants
of nucleosome positions, the full theory is required to disentangle one- and
two-body effects. Finally, we construct a minimal model in which nucleosomes
are positioned by steric exclusion and two-body interactions rather than
intrinsic histone-DNA sequence preferences. The model accurately reproduces
nucleosome occupancy patterns observed over transcribed regions in living
cells.


  In a previous paper we introduced a method called augmented sparse
reconstruction (ASR) that identifies links among nodes of ordinary differential
equation networks, given a small set of observed trajectories with various
initial conditions. The main purpose of that technique was to reconstruct
intracellular protein signaling networks. In this paper we show that a
recursive augmented sparse reconstruction generates artificial networks that
are homologous to a large, reference network, in the sense that kinase
inhibition of several reactions in the network alters the trajectories of a
sizable number of proteins in comparable ways for reference and reconstructed
networks. We show this result using a large in-silico model of the epidermal
growth factor receptor (EGF-R) driven signaling cascade to generate the data
used in the reconstruction algorithm. The most significant consequence of this
observed homology is that a nearly optimal combinatorial dosage of kinase
inhibitors can be inferred, for many nodes, from the reconstructed network, a
result potentially useful for a variety of applications in personalized
medicine.


  The transmission through contacts among MSM (men who have sex with men) is
one of the dominating contributors to HIV prevalence in industrialized
countries. In Amsterdam, the capital of the Netherlands, the MSM risk group has
been traced for decades. This has motivated studies which provide detailed
information about MSM's risk behavior statistically, psychologically and
sociologically. Despite the era of potent antiretroviral therapy, the incidence
of HIV among MSM increases. In the long term the contradictory effects of risk
behavior and effective therapy are still poorly understood. Using a previously
presented Complex Agent Network model, we describe steady and casual
partnerships to predict the HIV spreading among MSM. Behavior-related
parameters and values, inferred from studies on Amsterdam MSM, are fed into the
model; we validate the model using historical yearly incidence data.
Subsequently, we study scenarios to assess the contradictory effects of risk
behavior and effective therapy, by varying corresponding values of parameters.
Finally, we conduct quantitative analysis based on the resulting incidence
data. The simulated incidence reproduces the ACS historical incidence well and
helps to predict the HIV epidemic among MSM in Amsterdam. Our results show that
in the long run the positive influence of effective therapy can be outweighed
by an increase in risk behavior of at least 30% for MSM. Conclusion: We
recommend, based on the model predictions, that lowering risk behavior is the
prominent control mechanism of HIV incidence even in the presence of effective
therapy.


  The possibility of exploiting multiple resources is usually regarded as
positive from both the economic and the environmental point of view. However,
resource switching may also lead to unsustainable growth and, ultimately, to an
equilibrium condition which is worse than the one that could have been achieved
with a single resource. We applied a system dynamics model where users exploit
multiple resources and have different levels of preference among them. In this
setting, exploiting multiple resources leads to worse outcomes than the
single-resource case under a wide range of parameter configurations. Our
arguments are illustrated using two empirical situations, namely oil drilling
in the North Sea and whale hunting in the Antarctic.


  The method of flexi-Weighted Least Squares on evolutionary trees uses simple
polynomial or exponential functions of the evolutionary distance in place of
model-based variances. This has the advantage that unexpected deviations from
additivity can be modeled in a more flexible way. At present, only polynomial
weights have been used. However, a general family of exponential weights is
desirable to compare with polynomial weights and to potentially exploit recent
insights into fast least squares edge length estimation on trees. Here describe
families of weights that are multiplicative on trees, along with measures of
fit of data to tree. It is shown that polynomial, but also multiplicative
weights can approximate model-based variance of evolutionary distances well.
Both models are fitted to evolutionary data from yeast genomes and while the
polynomial weights model fits better, the exponential weights model can fit a
lot better than ordinary least squares. Iterated least squares is evaluated and
is seen to converge quickly and with minimal change in the fit statistics when
the data are in the range expected for the useful evolutionary distances and
simple Markov models of character change. In summary, both polynomial and
exponential weighted least squares work well and justify further investment
into developing the fastest possible algorithms for evaluating evolutionary
trees.


  Least squares trees, multi-dimensional scaling and Neighbor Nets are all
different and popular ways of visualizing multi-dimensional data. The method of
flexi-Weighted Least Squares (fWLS) is a powerful method of fitting
phylogenetic trees, when the exact form of errors is unknown. Here, both
polynomial and exponential weights are used to model errors. The exact same
models are implemented for multi-dimensional scaling to yield flexi-Weighted
MDS, including as special cases methods such as the Sammon Stress function.
Here we apply all these methods to population genetic data looking at the
relationships of "Abrahams Children" encompassing Arabs and now widely
dispersed populations of Jews, in relation to an African outgroup and a variety
of European populations. Trees, MDS and Neighbor Nets of this data are compared
within a common likelihood framework and the strengths and weaknesses of each
method are explored. Because the errors in this type of data can be complex,
for example, due to unexpected genetic transfer, we use a residual resampling
method to assess the robustness of trees and the Neighbor Net. Despite the
Neighbor Net fitting best by all criteria except BIC, its structure is ill
defined following residual resampling. In contrast, fWLS trees are favored by
BIC and retain considerable strong internal structure following residual
resampling. This structure clearly separates various European and Middle
Eastern populations, yet it is clear all of the models have errors much larger
than expected by sampling variance alone.


  The term unfoldome has been recently used to indicate the universe of
intrinsically disordered proteins. These proteins are characterized by an
ensemble of high-flexible interchangeable conformations and therefore they can
interact with many targets without requiring pre-existing stereo-chemical
complementarity. It has been suggested that intrinsically disordered proteins
are frequent in proteomes and disorder is widespread also in structured
proteins. However, several studies raise some doubt about these views. It this
paper we estimate the frequency of intrinsically disordered proteins in several
living organisms by using the ratio S between the likelihood, for a protein
sequence, of being composed mainly by order-promoting or disorder-promoting
residues. We scan several proteomes from Archaea, Bacteria and Eukarya. We find
the following figures: 1.63% for Archaea, 3.91% for Bacteria, 16.35% for
Eukarya. The frequencies we found can be considered an upper bound to the real
frequency of intrinsically disordered proteins in proteomes. Our estimates are
lower than those previously reported in several studies. A scanning of proteins
in the Protein Data Bank (PDB) searching for segments of non-observed residues
reveals that segments of non-observed residues longer than 30 amino acids, are
rare. Our observations support the idea that the spread of the unfoldome has
been often overestimated. If we exclude some exceptions, the structure-function
paradigm is generally valid and pre-existing stereo-chemical complementarity
among structures remains an important requisite for interactions between
biological macromolecules.


  We study the combined influence of amino acid composition and chain length on
the thermal stability of protein structures. A new parameterization of the
internal free energy is considered, as the sum of hydrophobic effect,
hydrogen-bond and de-hydration energy terms. We divided a non-redundant
selection of protein structures from the Protein Data Bank into three groups:
i) rich in order-promoting residues (OPR proteins); ii) rich in
disorder-promoting residues (DPR proteins); iii) belonging to a twilight zone
(TZ proteins). We observe a partition of PDB in several groups with different
internal free energies, amino acid compositions and protein lengths. Internal
free energy of 96% of the proteins analyzed ranges from -2 to -6.5 kJ/mol/res.
We found many DPR and OPR proteins with the same relative thermal stability.
Only OPR proteins with internal energy between -4 and -6.5 kJ/mol/res are
observed to have chains longer than 200 residues, with a high de-hydration
energy compensated by the hydrophobic effect. DPR and TZ proteins are shorter
than 200 residues and they have an internal energy above -4 kJ/mol/res, with a
few exceptions among TZ proteins. Hydrogen-bonds play an important role in the
stabilization of these DPR folds, often higher than contact energy. The new
parameterization of internal free energy let emerge a geography of thermal
stabilities of PDB structures. Amino acid composition per se is not sufficient
to determine the stability of protein folds, since. DPR and TZ proteins
generally have a relatively high internal free energy, and they are stabilized
by hydrogen-bonds. Long DPR proteins are not observed in the PDB, because their
low hydrophobicity cannot compensate the high de-hydration energy necessary to
accommodate residues within a highly packed globular fold.


  I examine the effect of exogenous spatial heterogeneity on the coexistence of
competing species using a simple model of non-hierarchical competition for site
occupancy on a lattice. The sites on the lattice are divided into two types
representing two different habitats or spatial resources. The model features no
temporal variability, hierarchical competition, type-dependent interactions or
other features traditionally known to support more competing species than there
are resources. Nonetheless, stable coexistence of two habitat specialists and a
generalist is observed in this model for a range of parameter values. In the
spatially implicit mean field approximation of the model, such coexistence is
shown to be impossible, demonstrating that it indeed arises from the explicit
spatial structure.


  We consider a system of $N$ individuals consisting of $S$ species that
interact pairwise: $x_m+x_\ell \rightarrow 2x_m\,\,$ with arbitrary
probabilities $p_m^\ell $. With no spatial structure, the master equation
yields a simple set of rate equations in a mean field approximation, the focus
of this note. Generalizing recent findings of cyclically competing three- and
four-species models, we cast these equations in an appealingly simple form. As
a result, many general properties of such systems are readily discovered, e.g.,
the major difference between even and odd $S$ cases. Further, we find the
criteria for the existence of (subspaces of) fixed points and collective
variables which evolve trivially (exponentially or invariant). These apparently
distinct aspects can be traced to the null space associated with the
interaction matrix, $p_m^\ell $. Related to the left- and right-
zero-eigenvectors, these appear to be "dual" facets of the dynamics. We also
remark on how the standard Lotka-Volterra equations (which include birth/death
terms) can be regarded as a special limit of a pairwise interacting system.


  In this article the results of Waddell and Azad (2009) are extended. In
particular, the geometric percentage mean standard deviation measure of the fit
of distances to a phylogenetic tree is adjusted for the number of parameters
fitted to the model. The formulae are also presented in their general form for
any weight that is a function of the distance. The cell line gene expression
data set of Ross et al. (2000) is reanalyzed. It is shown that ordinary least
squares (OLS) is a much better fit to the data than a Neighbor Joining or BME
tree. Residual resampling shows that cancer cell lines do indeed fit a tree
fairly well and that the tree does have strong internal structure. Simulations
show that least squares tree building methods, including OLS, are strong
competitors with BME type methods for fitting model data, while real world
examples often suggest the same conclusion.


  Molecular biology explains function of molecules by their geometrical and
electronical structures that are mainly determined by utilization of quantum
effects in chemistry. However, further quantum effects are not thought to play
any significant role in the essential processes of life. On the contrary,
consideration of quantum circuits/protocols and organic molecules as software
and hardware of living systems that are co-optimized during evolution, may be
useful to overcome the difficulties raised by biochemical complexity and to
understand the physics of life. In this sense, we review quantum
information-theoretic approaches to the process of DNA replication and propose
a new model in which 1) molecular recognition of a nucleobase is assumed to
trigger an intrabase entanglement corresponding to a superposition of different
tautomer forms and 2) pairing of complementary nucleobases is described by
swapping intrabase entanglements with interbase entanglements. We examine
possible biochemical realizations of quantum circuits/protocols to be used to
obtain intrabase and interbase entanglements. We deal with the problem of
cellular decoherence by using the theory of decoherence-free subspaces and
subsystems. Lastly, we discuss feasibility of the computational or experimental
verification of the model and future research directions.


  The Dissertation is focused on the studies of associations between functional
elements in human genome and their nucleotide structure. The asymmetry in
nucleotide content (skew, bias) was chosen as the main feature for nucleotide
structure. A significant difference in nucleotide content asymmetry was found
for human exons vs. introns. Specifically, exon sequences display bias for
purines (i.e., excess of A and G over C and T), while introns exhibit
keto-amino skew (i.e. excess of G and T over A and C). The extents of these
biases depend upon gene expression patterns. The highest intronic keto-amino
skew is found in the introns of housekeeping genes. In the case of introns,
whose sequences are under weak repair system, the AT->GC and CG->TA
substitutions are preferentially accumulated. A comparative analysis of gene
sequences encoding cytochrome P450 2E1 of Homo sapiens and representative
mammals was done. The cladistic tree on the basis of coding sequences
similarity of the gene Cyp2e1 was constructed. A new programming tools of NCBI
database sequence mining and analysis was developed, resulting in construction
of a own database.


  We study spectra of directed networks with inhibitory and excitatory
couplings. We investigate in particular eigenvector localization properties of
various model networks for different value of correlation among their entries.
Spectra of random networks, with completely uncorrelated entries show a
circular distribution with delocalized eigenvectors, where as networks with
correlated entries have localized eigenvectors. In order to understand the
origin of localization we track the spectra as a function of connection
probability and directionality. As connections are made directed, eigenstates
start occurring in complex conjugate pairs and the eigenvalue distribution
combined with the localization measure shows a rich pattern. Moreover, for a
very well distinguished community structure, the whole spectrum is localized
except few eigenstates at boundary of the circular distribution. As the network
deviates from the community structure there is a sudden change in the
localization property for a very small value of deformation from the perfect
community structure. We search for this effect for the whole range of
correlation strengths and for different community configurations. Furthermore,
we investigate spectral properties of a metabolic network of zebrafish, and
compare them with those of the model networks.


  Malaria control strategies aiming at reducing disease transmission intensity
may impact both oocyst intensity and infection prevalence in the mosquito
vector. Thus far, mathematical models failed to identify a clear relationship
between Plasmodium falciparum gametocytes and their infectiousness to
mosquitoes. Natural isolates of gametocytes are genetically diverse and
biologically complex. Infectiousness to mosquitoes relies on multiple
parameters such as density, sex-ratio, maturity, parasite genotypes and host
immune factors. In this article, we investigated how density and genetic
diversity of gametocytes impact on the success of transmission in the mosquito
vector. We analyzed data for which the number of covariates plus attendant
interactions is at least of order of the sample size, precluding usage of
classical models such as general linear models. We then considered the variable
importance from random forests to address the problem of selecting the most
influent variables. The selected covariates were assessed in the zero inflated
negative binomial model which accommodates both over-dispersion and the sources
of non infected mosquitoes. We found that the most important covariates related
to infection prevalence and parasite intensity are gametocyte density and
multiplicity of infection.


  Mechanosensitive channels are ion channels activated by membrane tension. We
investigate the influence of bacterial mechanosensitive channels spatial
distribution on activation (gating). Based on elastic short-range interactions
we map this physical process onto an Ising-like model, which enables us to
predict the clustering of channels and the effects of clustering on their
gating. We conclude that the aggregation of channels and the consequent
interactions among them leads to a global cooperative gating behaviour with
potentially dramatic consequences for the cell.


  Two proteins, one belonging to the mainly alpha class and the other belonging
to the alpha/beta class, are selected to test a kinetic mechanism for protein
folding. Targeted molecular dynamics is applied to generate folding pathways
for those two proteins, starting from two well defined initial conformations: a
fully extended and a alpha-helical conformation. The results show that for both
proteins the alpha-helical initial conformation provides overall lower energy
pathways to the native state. For the alpha/beta protein, 30 % (40%) of the
pathways from an initial alpha-helix (fully extended) structure lead to
unentangled native folds, a success rate that can be increased to 85 % by the
introduction of a well-defined intermediate structure. These results open up a
new direction in which to look for a solution to the protein folding problem,
as detailed at the end.


  Empirical evidence is presented which contradicts the established
interpretation of the intraspecific genetic variability as the adaptive
potential of the species: the uniform evolutionary tempo across the life
kingdom, species discreteness, and absence of correlation between genetic
variability and prosperity of extant species testify against the ruling
paradigm. Consistent interpretation of the nature of intraspecific genetic
variability is based on recognizing the limited sensitivity of stabilising
selection, which allows for accumulation and persistence in the population of a
considerable amount of mutational substitutions, which, to some degree, erase
the meaningful genetic information of the species. The proposed interpretation
also provides solution to the inbreeding paradox in the invasive species.


  The construction of synthetic biochemical circuits is an essential step for
developing quantitative understanding of information processing in natural
organisms. Here, we report construction and analysis of an in vitro circuit
with positive autoregulation that consists of just four synthetic DNA strands
and three enzymes, bacteriophage T7 RNA polymerase, Escherichia coli
ribonuclease (RNase) H, and RNase R. The modularity of the DNA switch template
allowed a rational design of a synthetic DNA switch regulated by its RNA output
acting as a transcription activator. We verified that the thermodynamic and
kinetic constraints dictated by the sequence design criteria were enough to
experimentally achieve the intended dynamics: a transcription activator
configured to regulate its own production. Although only RNase H is necessary
to achieve bistability of switch states, RNase R is necessary to maintain
stable RNA signal levels and to control incomplete degradation products. A
simple mathematical model was used to fit ensemble parameters for the training
set of experimental results and was then directly applied to predict
time-courses of switch dynamics and sensitivity to parameter variations with
reasonable agreement. The positive autoregulation switches can be used to
provide constant input signals and store outputs of biochemical networks and
are potentially useful for chemical control applications.


  This paper is concerned with the dynamical properties of deterministically
modeled chemical reaction systems. Specifically, this paper provides a proof of
the Global Attractor Conjecture in the setting where the underlying reaction
diagram consists of a single linkage class, or connected component. The
conjecture dates back to the early 1970s and is the most well known and
important open problem in the field of chemical reaction network theory. The
resolution of the conjecture has important biological and mathematical
implications in both the deterministic and stochastic settings. One of our main
analytical tools, which is introduced here, will be a method for partitioning
the relevant monomials of the dynamical system along sequences of trajectory
points into classes with comparable growths. We will use this method to
conclude that if a trajectory converges to the boundary, then a whole family of
Lyapunov functions decrease along the trajectory. This will allow us to
overcome the fact that the usual Lyapunov functions of chemical reaction
network theory are bounded on the boundary of the positive orthant, which has
been the technical sticking point to a proof of the Global Attractor Conjecture
in the past.


  Competition is one of the most fundamental phenomena in physics, biology and
economics. Recent studies of the competition between innovations have
highlighted the influence of switching costs and interaction networks, but the
problem is still puzzling. We introduce a model that reveals a novel
multi-percolation process, which governs the struggle of innovations trying to
penetrate a market. We find that innovations thrive as long as they percolate
in a population, and one becomes dominant when it is the only one that
percolates. Besides offering a theoretical framework to understand the
diffusion of competing innovations in social networks, our results are also
relevant to model other problems such as opinion formation, political
polarization, survival of languages and the spread of health behavior.


  In this paper we propose a malaria within-host model with k classes of age
for the parasitized red blood cells and n strains for the parasite. We provide
a global analysis for this model. A competitive exclusion principle holds. If
R0, the basic reproduction number, satisfies R0 </- 1, then the disease-free
equilibrium is globally asymptotically stable. On the contrary if R0 > 1, then
generically there is a unique endemic equilibrium which corresponds to the
endemic stabilization of the most virulent parasite strain and to the
extinction of all the other parasites strains. We prove that this equilibrium
is globally asymptotically stable on the positive orthant if a mild sufficient
condition is satisfied.


  The Red Queen's notion "It takes all the running you can do, to keep in the
same place" has been applied within evolutionary biology, politics and
economics. We find that a generalized version in which an adaptive Red Queen
(e.g. insurgency) sporadically edges ahead of a Blue King (e.g. military),
explains the progress curves for fatal insurgent attacks against the coalition
military within individual provinces in Afghanistan and Iraq. Remarkably
regular mathematical relations emerge which suggest a prediction formula for
the timing of the n'th future fatal day, and provide a common framework for
understanding how insurgents fight in different regions. Our findings are
consistent with a Darwinian selection hypothesis which favors a weak species
which can adapt rapidly, and establish an unexpected conceptual connection to
the physics of correlated walks.


  Organisms leave a distinctive chemical signature in their environment because
they synthesize those molecules that maximize their fitness. As a result, the
relative concentrations of related chemical monomers in life-bearing
environmental samples reflect, in part, those compounds' adaptive utility. In
contrast, rates of molecular synthesis in a lifeless environment are dictated
by reaction kinetics and thermodynamics, so concentrations of related monomers
in abiotic samples tend to exhibit specific patterns dominated by small, easily
formed, low-formation-energy molecules. We contend that this distinction can
serve as a universal biosignature: the measurement of chemical concentration
ratios that belie formation kinetics or equilibrium thermodynamics indicates
the likely presence of life. We explore the features of this biosignature as
observed in amino acids and carboxylic acids, using published data from
numerous studies of terrestrial sediments, abiotic (spark, UV, and high-energy
proton) synthesis experments, and meteorite bodies. We then compare these data
to the results of experimental studies of an evolving digital life system. We
observe the robust and repeatable evolution of an analogous biosignature in a
digital lifeform, suggesting that evolutionary selection necessarily constrains
organism composition and that the monomer abundance biosignature phenomenon is
universal to evolved biosystems.


  Actomyosin contractility is essential for biological force generation, and is
well understood in highly organized structures such as striated muscle.
Additionally, actomyosin bundles devoid of this organization are known to
contract both in vivo and in vitro, which cannot be described by standard
muscle models. To narrow down the search for possible contraction mechanisms in
these systems, we investigate their microscopic symmetries. We show that
contractile behavior requires non-identical motors that generate large enough
forces to probe the nonlinear elastic behavior of F-actin. This suggests a role
for filament buckling in the contraction of these bundles, consistent with
recent experimental results on reconstituted actomyosin bundles.


  The Michaelis-Menten equation has played a central role in our understanding
of biochemical processes. It has long been understood how this equation
approximates the dynamics of irreversible enzymatic reactions. However, a
similar approximation in the case of networks, where the product of one
reaction can act as an enzyme in another, has not been fully developed. Here we
rigorously derive such an approximation in a class of coupled enzymatic
networks where the individual interactions are of Michaelis-Menten type. We
show that the sufficient conditions for the validity of the total quasi steady
state assumption (tQSSA), obtained in a single protein case by Borghans, de
Boer and Segel can be extended to sufficient conditions for the validity of the
tQSSA in a large class of enzymatic networks. Secondly, we derive reduced
equations that approximate the network's dynamics and involve only protein
concentrations. This significantly reduces the number of equations necessary to
model such systems. We prove the validity of this approximation using geometric
singular perturbation theory and results about matrix differentiation. The
ideas used in deriving the approximating equations are quite general, and can
be used to systematize other model reductions.


  Mass spectrometry-based proteomics has become the tool of choice for
identifying and quantifying the proteome of an organism. Though recent years
have seen a tremendous improvement in instrument performance and the
computational tools used, significant challenges remain, and there are many
opportunities for statisticians to make important contributions. In the most
widely used "bottom-up" approach to proteomics, complex mixtures of proteins
are first subjected to enzymatic cleavage, the resulting peptide products are
separated based on chemical or physical properties and analyzed using a mass
spectrometer. The two fundamental challenges in the analysis of bottom-up
MS-based proteomics are as follows: (1) Identifying the proteins that are
present in a sample, and (2) Quantifying the abundance levels of the identified
proteins. Both of these challenges require knowledge of the biological and
technological context that gives rise to observed data, as well as the
application of sound statistical principles for estimation and inference. We
present an overview of bottom-up proteomics and outline the key statistical
issues that arise in protein identification and quantification.


  Inspired by dengue and yellow fever epidemics, we investigated the contact
process (CP) in a multiscale network constituted by one-dimensional chains
connected through a Barab\'asi-Albert scale-free network. In addition to the CP
dynamics inside the chains, the exchange of individuals between connected
chains (travels) occurs at a constant rate. A finite epidemic threshold and an
epidemic mean lifetime diverging exponentially in the subcritical phase,
concomitantly with a power law divergence of the outbreak's duration, were
found. A generalized scaling function involving both regular and SF components
was proposed for the quasistationary analysis and the associated critical
exponents determined, demonstrating that the CP on this hybrid network and
nonvanishing travel rates establishes a new universality class.


  We develop a model to describe the force generated by the polymerization of
an array of parallel biofilaments. The filaments are assumed to be coupled only
through mechanical contact with a movable barrier. We calculate the filament
density distribution and the force-velocity relation with a mean-field approach
combined with simulations. We identify two regimes: a non-condensed regime at
low force in which filaments are spread out spatially, and a condensed regime
at high force in which filaments accumulate near the barrier. We confirm a
result previously known from other related studies, namely that the stall force
is equal to N times the stall force of a single filament. In the model studied
here, the approach to stalling is very slow, and the velocity is practically
zero at forces significantly lower than the stall force.


  Genetic interactions have been widely used to define functional relationships
between proteins and pathways. In this study, we demonstrated that yeast
synthetic lethal genetic interactions can be explained by the genetic
interactions between domains of those proteins. The domain genetic interactions
rarely overlap with the domain physical interactions from iPfam database and
provide a complementary view about domain relationships. Moreover, we found
that domains in multidomain yeast proteins contribute to their genetic
interactions differently. The domain genetic interactions help more precisely
define the function related to the synthetic lethal genetic interactions, and
then help understand how domains contribute to different functionalities of
multidomain proteins. Using the probabilities of domain genetic interactions,
we were able to predict novel yeast synthetic lethal genetic interactions.
Furthermore, we had also identified novel compensatory pathways from the
predicted synthetic lethal genetic interactions. Our study significantly
improved the understanding of yeast mulitdomain proteins, the synthetic lethal
genetic interactions and the functional relationships between proteins and
pathways.


  The Yule model and the coalescent model are two neutral stochastic models for
generating trees in phylogenetics and population genetics, respectively.
Although these models are quite different, they lead to identical distributions
concerning the probability that pre-specified groups of taxa form monophyletic
groups (clades) in the tree. We extend earlier work to derive exact formulae
for the probability of finding one or more groups of taxa as clades in a rooted
tree, or as `clans' in an unrooted tree. Our findings are relevant for
calculating the statistical significance of observed monophyly and reciprocal
monophyly in phylogenetics.


  We consider the problem of detecting multiple changepoints in large data
sets. Our focus is on applications where the number of changepoints will
increase as we collect more data: for example in genetics as we analyse larger
regions of the genome, or in finance as we observe time-series over longer
periods. We consider the common approach of detecting changepoints through
minimising a cost function over possible numbers and locations of changepoints.
This includes several established procedures for detecting changing points,
such as penalised likelihood and minimum description length. We introduce a new
method for finding the minimum of such cost functions and hence the optimal
number and location of changepoints that has a computational cost which, under
mild conditions, is linear in the number of observations. This compares
favourably with existing methods for the same problem whose computational cost
can be quadratic or even cubic. In simulation studies we show that our new
method can be orders of magnitude faster than these alternative exact methods.
We also compare with the Binary Segmentation algorithm for identifying
changepoints, showing that the exactness of our approach can lead to
substantial improvements in the accuracy of the inferred segmentation of the
data.


  In recent experiments by Richardson et al. ((2010), PLoS ONE 5(3): e9621.
doi:10.1371/ journal.pone.0009621) ant motion out of the nest is shown to be a
non-stationary process intriguingly similar to the dynamics encountered in
\emph{physical aging} of glassy systems. Specifically, exit events can be
described as a Poisson process in logarithmic time, or, for short, a
log-Poisson process. Nouvellet et al.(J. Theor. Biol. 266, 573. (2010))
criticized these conclusions and performed new experiments where the exit
process could more simply be described by standard Poisson statistics. In their
reply, (J Theor. Biol. 269, 356-358 (2011)) Richardson et al. stressed that the
two sets of experiments were performed under very different conditions and
claimed that this was the likely source of the discrepancy. The focal point of
this work is whether log-Poisson and Poisson statistics both are possible under
different external conditions. To this end, a model is introduced where
interacting ants move in a stochastic fashion from one site to a neighboring
site on a finite 2D lattice. The probability of each move is determined by the
ensuing changes of a utility function which is a sum of pairwise interactions
between ants, weighted by distance. Depending on how the interactions are
defined and on a control parameter dubbed 'degree of stochasticity' (DS), the
dynamics either quickly converges to a stationary state, where movements are a
standard Poisson process, or may enter a non-stationary regime, where exits can
be described as suggested by Richardson et al.


  Integrin signaling regulates cell migration and plays a pivotal role in
developmental processes and cancer metastasis. Integrin signaling has been
studied extensively and much data is available on pathway components and
interactions. Yet the data is fragmented and an integrated model is missing. We
use a rule-based modeling approach to integrate available data and test
biological hypotheses regarding the role of talin, Dok1 and PIPKI in integrin
activation. The detailed biochemical characterization of integrin signaling
provides us with measured values for most of the kinetics parameters. However,
measurements are not fully accurate and the cellular concentrations of
signaling proteins are largely unknown and expected to vary substantially
across different cellular conditions. By sampling model behaviors over the
physiologically realistic parameter range we find that the model exhibits only
two different qualitative behaviours and these depend mainly on the relative
protein concentrations, which offers a powerful point of control to the cell.
Our study highlights the necessity to characterize model behavior not for a
single parameter optimum, but to identify parameter sets that characterize
different signaling modes.


  We analyze the timescales of conflict decision-making in a primate society.
We present evidence for multiple, periodic timescales associated with social
decision-making and behavioral patterns. We demonstrate the existence of
periodicities that are not directly coupled to environmental cycles or known
ultraridian mechanisms. Among specific biological and socially-defined
demographic classes, periodicities span timescales between hours and days, and
many are not driven by exogenous or internal regularities. Our results indicate
that they are instead driven by strategic responses to social interaction
patterns. Analyses also reveal that a class of individuals, playing a critical
functional role, policing, have a signature timescale on the order of one hour.
We propose a classification of behavioral timescales analogous to those of the
nervous system, with high-frequency, or $\alpha$-scale, behavior occurring on
hour-long scales, through to multi-hour, or $\beta$-scale, behavior, and,
finally $\gamma$ periodicities observed on a timescale of days.


  This article aims to clarify the case and the mechanism where sanction or
punishment by institution can deliver the evolution of cooperation. Compared to
peer sanctioning, institutional sanctioning may be sensitive to players'
attitude toward players who do not pre-commit punishment. Departed from former
studies based on the punisher who always acts cooperatively, we assume that the
punishing player is skeptical in that she cooperates in proportion to how many
same types join in her team. Relying on stochastic adaptive dynamics, we show
that institutional sanctioning coupled with skeptical punisher can make
cooperation evolve for the case where peer sanctioning may not.


  The white-rot fungus, Physisporinus vitreus, degrades the membranes of
bordered pits in tracheids and consequently increases the permeability of wood,
which is a process that can be used by the wood industry to improve the uptake
of wood preservatives and environmentally benign wood modification substances
to enhance the use and sustainability of native conifer wood species. To
understand and apply this process requires an understanding of how a complex
system (fungus-wood) interacts under defined conditions. We present a
three-dimensional fungal growth model (FGM) of the hyphal growth of P. vitreus
in the heartwood of Norway spruce. The model considers hyphae and nutrients as
discrete structures and links the microscopic interactions between fungus and
wood (e.g. degradation rate and degree of opening of pits) with macroscopic
system properties, such penetration depth of the fungus, biomass and
distribution of destroyed pits in early- and latewood. Simulations were
compared with experimental data. The growth of P. vitreus is characterized by a
stepwise capture of the substrate and the effect of this on wood according to
different model parameters is discussed.


  In this paper we explore the life expectancy limits by based on the
stochastic modeling of mortality and applying the first exit or hitting time
theory of a stochastic process. The main assumption is that the health state or
the "vitality", according to Strehler and Mildvan, of an individual is a
stochastic variable and thus it was introduced and applied a first exit time
density function to mortality data. The model is used to estimate the
development of mortality rates in the late stages of the human life span, to
make better fitting to population mortality data including the infant
mortality, to compare it with the classical Gompertz curve, and to make
comparisons between the Carey med-fly data and the population mortality data
estimating the health state or "vitality" functions. Furthermore, we apply the
model to the life table data of Italy, France, USA, Canada, Sweden, Norway and
Japan, and we analyze the characteristic parameters of the model and make
forecasts.


  We describe a new general connection between $\Lambda$-coalescents and
genealogies of continuous-state branching processes. This connection is based
on the construction of an explicit coupling using a particle representation
inspired by the lookdown process of Donnelly and Kurtz. This coupling has the
property that the coalescent comes down from infinity if and only if the
branching process becomes extinct, thereby answering a question of Bertoin and
Le Gall. The coupling also offers new perspective on the speed of coming down
from infinity and allows us to relate power-law behavior for $N^{\Lambda}(t)$
to the classical upper and lower indices arising in the study of pathwise
properties of L\'{e}vy processes.


  Random Boolean networks (RBNs) have been a popular model of genetic
regulatory networks for more than four decades. However, most RBN studies have
been made with random topologies, while real regulatory networks have been
found to be modular. In this work, we extend classical RBNs to define modular
RBNs. Statistical experiments and analytical results show that modularity has a
strong effect on the properties of RBNs. In particular, modular RBNs have more
attractors and are closer to criticality when chaotic dynamics would be
expected, compared to classical RBNs.


  This contribution is concerned with mathematical models for the dynamics of
the genetic composition of populations evolving under recombination.
Recombination is the genetic mechanism by which two parent individuals create
the mixed type of their offspring during sexual reproduction. The corresponding
models are large, nonlinear dynamical systems (for the deterministic treatment
that applies in the infinite-population limit), or interacting particle systems
(for the stochastic treatment required for finite populations). We review
recent progress on these difficult problems. In particular, we present a closed
solution of the deterministic continuous-time system, for the important special
case of single crossovers; we extract an underlying linearity; we analyse how
this carries over to the corresponding stochastic setting; and we provide a
solution of the analogous deterministic discrete-time dynamics, in terms of its
generalised eigenvalues and a simple recursion for the corresponding
coefficients.


  We show that two important problems that have applications in computational
biology are ASP-complete, which implies that, given a solution to a problem, it
is NP-complete to decide if another solution exists. We show first that a
variation of Betweenness, which is the underlying problem of questions related
to radiation hybrid mapping, is ASP-complete. Subsequently, we use that result
to show that Quartet Compatibility, a fundamental problem in phylogenetics that
asks whether a set of quartets can be represented by a parent tree, is also
ASP-complete. The latter result shows that Steel's \sc Quartet Challenge, which
asks whether a solution to Quartet Compatibility is unique, is coNP-complete.


  Regulation of intrinsic noise in gene expression is essential for many
cellular functions. Correspondingly, there is considerable interest in
understanding how different molecular mechanisms of gene expression impact
variations in protein levels across a population of cells. In this work, we
analyze a stochastic model of bursty gene expression which considers general
waiting-time distributions governing arrival and decay of proteins. By mapping
the system to models analyzed in queueing theory, we derive analytical
expressions for the noise in steady-state protein distributions. The derived
results extend previous work by including the effects of arbitrary probability
distributions representing the effects of molecular memory and bursting. The
analytical expressions obtained provide insight into the role of
transcriptional, post-transcriptional and post-translational mechanisms in
controlling the noise in gene expression.


  Large scale motion of cytoplasm called cytoplasmic streaming occurs in some
large eukaryotic cells to stir the cell's constituents. In Drosophila oocytes,
microtubules have been observed to undergo undulating motion, curving to form
travelling waves during cytoplasmic streaming. Here we show that this wave-like
motion can be understood physically as due to the hydrodynamic drag of
streaming impellers attached to kinesin motors moving toward the plus-ends of
microtubules whose minus ends are anchored to the cell cortex. The tangential
forces applied to such microtubules by kinesin give rise to bending and leads
to chiral symmetry breaking causing the microtubules to propagate long
travelling waves. The waves are reminiscent of those seen in flagellar motion
but of a much longer time scale and by a different physical mechanism. We show
how kinesin movement can produce a bulk flow of cytoplasm surrounding a
microtubule with the range of flow greatly enhanced by the effect of
hydrodynamic coupling between impellers. That is, a relatively small number of
motors can move a large amount of fluid. The chaotic nature of the fluid motion
of cytoplasm caused by kinesin movement along constantly changing microtubule
trajectories is important as it greatly enhances the efficiency of mixing.
Existing data on in vitro microtubule gliding assays also show this chiral
instability in two dimensions and an analysis of this gives quantitative
estimates for the forces exerted by motors and the drag coefficient.


  Computer experiments are performed to investigate why protein contact
networks (networks induced by spatial contacts between amino acid residues of a
protein) do not have shorter average shortest path lengths in spite of their
importance to protein folding. We find that shorter average inter-nodal
distances is no guarantee of finding a global optimum more easily. Results from
the experiments also led to observations which parallel an existing view that
neither short-range nor long-range interactions dominate the protein folding
process. Nonetheless, runs where there was a slight delay in the use of
long-range interactions yielded the best search performance. We incorporate
this finding into the optimization function by giving more weight to
short-range links. This produced results showing that randomizing long-range
links does not yield better search performance than protein contact networks au
natural even though randomizing long-range links significantly reduces average
path lengths and retains much of the clustering and positive degree-degree
correlation inherent in protein contact networks. Hence there can be
explanations, other than the excluded volume argument, beneath the topological
limits of protein contact networks.


  Biological cells sense external chemical stimuli in their environment using
cell-surface receptors. To increase the sensitivity of sensing, receptors often
cluster, most noticeably in bacterial chemotaxis, a paradigm for signaling and
sensing in general. While amplification of weak stimuli is useful in absence of
noise, its usefulness is less clear in presence of extrinsic input noise and
intrinsic signaling noise. Here, exemplified on bacterial chemotaxis, we
combine the allosteric Monod-Wyman- Changeux model for signal amplification by
receptor complexes with calculations of noise to study their
interconnectedness. Importantly, we calculate the signal-to-noise ratio,
describing the balance of beneficial and detrimental effects of clustering for
the cell. Interestingly, we find that there is no advantage for the cell to
build receptor complexes for noisy input stimuli in absence of intrinsic
signaling noise. However, with intrinsic noise, an optimal complex size arises
in line with estimates of the sizes of chemoreceptor complexes in bacteria and
protein aggregates in lipid rafts of eukaryotic cells.


  The problem of the helix-coil transition of biopolymers in explicit solvents,
like water, with the ability for hydrogen bonding with solvent is addressed
analytically using a suitably modified version of the Generalized Model of
Polypeptide Chains. Besides the regular helix-coil transition, an additional
coil-helix or reentrant transition is also found at lower temperatures. The
reentrant transition arises due to competition between polymer-polymer and
polymer-water hydrogen bonds. The balance between the two types of hydrogen
bonding can be shifted to either direction through changes not only in
temperature, but also by pressure, mechanical force, osmotic stress or other
external influences. Both polypeptides and polynucleotides are considered
within a unified formalism. Our approach provides an explanation of the
experimental difficulty of observing the reentrant transition with pressure;
and underscores the advantage of pulling experiments for studies of DNA.
Results are discussed and compared with those reported in a number of recent
publications with which a significant level of agreement is obtained.


  Agent-based stochastic models for finite populations have recently received
much attention in the game theory of evolutionary dynamics. Both the ultimate
fixation and the pre-fixation transient behavior are important to a full
understanding of the dynamics. In this paper, we study the transient dynamics
of the well-mixed Moran process through constructing a landscape function. It
is shown that the landscape playing a central theoretical "device" that
integrates several lines of inquiries: the stable behavior of the replicator
dynamics, the long-time fixation, and continuous diffusion approximation
associated with asymptotically large population. Several issues relating to the
transient dynamics are discussed: (i) multiple time scales phenomenon
associated with intra- and inter-attractoral dynamics; (ii) discontinuous
transition in stochastically stationary process akin to Maxwell construction in
equilibrium statistical physics; and (iii) the dilemma diffusion approximation
facing as a continuous approximation of the discrete evolutionary dynamics. It
is found that rare events with exponentially small probabilities, corresponding
to the uphill movements and barrier crossing in the landscape with multiple
wells that are made possible by strong nonlinear dynamics, plays an important
role in understanding the origin of the complexity in evolutionary, nonlinear
biological systems.


  RNA viruses exist as genetically diverse populations displaying different
phenotypes, including diverse degrees of virulence. The evolution of virulence
in viral populations is, however, poorly understood. Based on the experimental
observation of an RNA virus clone in cell culture diversifying into two
subpopulations of different virulence, we study the dynamics of heterogeneous
virus populations and the evolution of virulence. We introduce a
competition-colonization trade-off into standard mathematical models of
intra-host viral infection. Colonizers are fast spreading, virulent strains,
whereas competitors are less virulent variants that are more successful within
coinfected cells. We observe biphasic dynamics of the population: Early in the
infection the population is dominated by colonizers, which later will be
outcompeted by competitors. The simulations suggest the existence of a steady
state where few low-virulence variants coexist. This equilibrium implies
collective virulence attenuation in the population, in contrast to previous
models predicting development of the population towards increased virulence.
Nevertheless, the attenuation effect disappears if we include a highly
simplified immune response in our models. Thus, the competition-colonization
trade-off indicates a role for virulence in the modulation of the viral
population diversity. The evolution of virulence is a dynamic feature of the
population shaped by interactions between individuals and by the structure of
the patchy habitat.


  In this note we review recent progress about fluctuation relations and their
applicability to free energy recovery in single molecule experiments. We
underline the importance of the operational definition for the mechanical work
and the non-invariance of fluctuation relations under Galilean transformations,
both aspects currently amenable to experimental test. Finally we describe a
generalization of the Crooks fluctuation relation useful to recover free
energies of partially equilibrated states and thermodynamic branches.


  This paper has been withdrawn by the author due to a crucial error in
mathematical derivation and copyrights.


  Biophysical models are often populated by a large number of input parameters
that are difficult to predict or measure experimentally. The validity and
robustness of a given model can be evaluated by a sensitivity test to its input
parameters. In this study, we performed local (based on a Taylor-like method)
and global sensitivity (based on Monte Carlo filtering techniques) analyses of
a previously derived PAR1-mediated activation model of endothelial cells. This
activation model previously demonstrated that peptide-activated PAR1 has a
different receptor/G-protein binding affinity that favors Galpha_q activation
over Galpha_12/13 by approximately 800-fold. Interestingly, the present study
shows that the parameter regulating the binding rate of activated PAR1 to
Galpha_12/13 is indeed important to obtain the expected RhoGTP response.
Moreover, we show that the parameters representing the rate of PAR1
deactivation and the rate of PAR1 binding to G_q, are the most important
parameters in the system. Finally, we illustrate that the kinetic model
considered in this study is robust and we provide complementary insights into
the biological meaning and importance of its kinetic parameters.


  The diversity of cell populations is regulated by extracellular and
intracellular variability. The latter includes genetic, epigenetic and
stochastic variability, all contributing to the experimentally observed
heterogeneity in response to external death-inducing stimuli. Studies of
sources and regulation of variability in commitment to apoptotic cancer cell
death are likely to identify the fundamental features of apoptotic protein
networks that are responsible for determining the ultimate cell fate. Systems
biology approaches, involving computer simulations of the biochemical reactions
accompanied, if possible, by experimental verification of selected components
of the model, are proving useful in determining the origins of cell-to-cell
variability in response to external stress stimuli. Here we summarize our
current understanding of the origins of stochastic variability in cells'
commitment to apoptosis, and its implications in the field on cancer therapy.


  Genetic interaction measures how different genes collectively contribute to a
phenotype, and can reveal functional compensation and buffering between
pathways under genetic perturbations. Recently, genome-wide screening for
genetic interactions has revealed genetic interaction networks that provide
novel insights either when analyzed by themselves or when integrated with other
functional genomic datasets. For higher eukaryotes such as human, the above
reverse-genetics approaches are not straightforward since the phenotypes of
interest for higher eukaryotes are difficult to study in a cell based assay. We
propose a general framework for constructing and analyzing human genetic
interaction networks from genome-wide single nucleotide polymorphism (SNP) data
used for case-control studies on complex diseases. Specifically, the approach
contains three major steps: (1) estimating SNP-SNP genetic interactions, (2)
identifying linkage disequilibrium (LD) blocks and mapping SNP-SNP interactions
to block-block interactions, and (3) functional mapping for LD blocks. We
performed two sets of functional analyses for each of the six datasets used in
the paper, and demonstrated that (i) the constructed genetic interaction
networks are supported by functional evidence from independent biological
databases, and (ii) the network can be used to discover pairs of compensatory
gene modules (between-pathway models) in their joint association with a disease
phenotype. The proposed framework should provide novel insights beyond existing
approaches that either ignore interactions between SNPs or model different
SNP-SNP pairs with genetic interactions separately. Furthermore, our study
provides evidence that some of the core properties of genetic interaction
networks based on reverse genetics in model organisms like yeast are also
present in genetic interactions revealed by natural variation in human
populations.


  A number of methods have been proposed over the last decade for encoding
information using deoxyribonucleic acid (DNA), giving rise to the emerging area
of DNA data embedding. Since a DNA sequence is conceptually equivalent to a
sequence of quaternary symbols (bases), DNA data embedding (diversely called
DNA watermarking or DNA steganography) can be seen as a digital communications
problem where channel errors are tantamount to mutations of DNA bases.
Depending on the use of coding or noncoding DNA hosts, which, respectively,
denote DNA segments that can or cannot be translated into proteins, DNA data
embedding is essentially a problem of communications with or without side
information at the encoder. In this paper the Shannon capacity of DNA data
embedding is obtained for the case in which DNA sequences are subject to
substitution mutations modelled using the Kimura model from molecular evolution
studies. Inferences are also drawn with respect to the biological implications
of some of the results presented.


  Gene Set Enrichment Analysis (GSEA) and its variations aim to discover
collections of genes that show moderate but coordinated differences in
expression. However, such techniques may be ineffective if many individual
genes in a phenotype-related gene set have weak discriminative power. A
potential solution is to search for combinations of genes that are highly
differentiating even when individual genes are not. Although such techniques
have been developed, these approaches have not been used with GSEA to any
significant degree because of the large number of potential gene combinations
and the heterogeneity of measures that assess the differentiation provided by
gene groups of different sizes.
  To integrate the search for differentiating gene combinations and GSEA, we
propose a general framework with two key components: (A) a procedure that
reduces the number of scores to be handled by GSEA to the number of genes by
summarizing the scores of the gene combinations involving a particular gene in
a single score, and (B) a procedure to integrate the heterogeneous scores from
combinations of different sizes and from different gene combination measures by
mapping the scores to p-values. Experiments on four gene expression data sets
demonstrate that the integration of GSEA and gene combination search can
enhance the power of traditional GSEA by discovering gene sets that include
genes with weak individual differentiation but strong joint discriminative
power. Also, gene sets discovered by the integrative framework share several
common biological processes and improve the consistency of the results among
three lung cancer data sets.


  Inferring genetic networks from gene expression data is one of the most
challenging work in the post-genomic era, partly due to the vast space of
possible networks and the relatively small amount of data available. In this
field, Gaussian Graphical Model (GGM) provides a convenient framework for the
discovery of biological networks. In this paper, we propose an original
approach for inferring gene regulation networks using a robust biological prior
on their structure in order to limit the set of candidate networks.
  Pathways, that represent biological knowledge on the regulatory networks,
will be used as an informative prior knowledge to drive Network Inference. This
approach is based on the selection of a relevant set of genes, called the
"molecular signature", associated with a condition of interest (for instance,
the genes involved in disease development). In this context, differential
expression analysis is a well established strategy. However outcome signatures
are often not consistent and show little overlap between studies. Thus, we will
dedicate the first part of our work to the improvement of the standard process
of biomarker identification to guarantee the robustness and reproducibility of
the molecular signature.
  Our approach enables to compare the networks inferred between two conditions
of interest (for instance case and control networks) and help along the
biological interpretation of results. Thus it allows to identify differential
regulations that occur in these conditions. We illustrate the proposed approach
by applying our method to a study of breast cancer's response to treatment.


  Individuals interact with conspecifics in a number of behavioural contexts or
dimensions. Here, we formalise this by considering a social network between n
individuals interacting in b behavioural dimensions as a nxnxb multidimensional
object. In addition, we propose that the topology of this object is driven by
individual needs to reduce uncertainty about the outcomes of interactions in
one or more dimension. The proposal grounds social network dynamics and
evolution in individual selection processes and allows us to define the
uncertainty of the social network as the joint entropy of its constituent
interaction networks. In support of these propositions we use simulations and
natural 'knock-outs' in a free-ranging baboon troop to show (i) that such an
object can display a small-world state and (ii) that, as predicted, changes in
interactions after social perturbations lead to a more certain social network,
in which the outcomes of interactions are easier for members to predict. This
new formalisation of social networks provides a framework within which to
predict network dynamics and evolution under the assumption that it is driven
by individuals seeking to reduce the uncertainty of their social environment.


  On the basis of first-principles GW calculations, we study the quasiparticle
properties of the guanine, adenine, cytosine, thymine, and uracil DNA and RNA
nucleobases. Beyond standard G0W0 calculations, starting from Kohn-Sham
eigenstates obtained with (semi)local functionals, a simple self-consistency on
the eigenvalues allows to obtain vertical ionization energies and electron
affinities within an average 0.11 eV and 0.18 eV error respectively as compared
to state-of-the-art coupled-cluster and multi-configurational perturbative
quantum chemistry approaches. Further, GW calculations predict the correct \pi
-character of the highest occupied state, thanks to several level crossings
between density functional and GW calculations. Our study is based on a recent
gaussian-basis implementation of GW with explicit treatment of dynamical
screening through contour deformation techniques.


  We investigate the mutation-selection dynamics for an evolutionary
computation model based on Turing Machines that we introduced in a previous
article.
  The use of Turing Machines allows for very simple mechanisms of code growth
and code activation/inactivation through point mutations. To any value of the
point mutation probability corresponds a maximum amount of active code that can
be maintained by selection and the Turing machines that reach it are said to be
at the error threshold. Simulations with our model show that the Turing
machines population evolve towards the error threshold.
  Mathematical descriptions of the model point out that this behaviour is due
more to the mutation-selection dynamics than to the intrinsic nature of the
Turing machines. This indicates that this result is much more general than the
model considered here and could play a role also in biological evolution.


  Genetic regulatory networks enable cells to respond to the changes in
internal and external conditions by dynamically coordinating their gene
expression profiles. Our ability to make quantitative measurements in these
biochemical circuits has deepened our understanding of what kinds of
computations genetic regulatory networks can perform and with what reliability.
These advances have motivated researchers to look for connections between the
architecture and function of genetic regulatory networks. Transmitting
information between network's inputs and its outputs has been proposed as one
such possible measure of function, relevant in certain biological contexts.
Here we summarize recent developments in the application of information theory
to gene regulatory networks. We first review basic concepts in information
theory necessary to understand recent work. We then discuss the functional
complexity of gene regulation which arrises from the molecular nature of the
regulatory interactions. We end by reviewing some experiments supporting the
view that genetic networks responsible for early development of multicellular
organisms might be maximizing transmitted 'positional' information.


  All biological processes are controlled by complex systems of enzymatic
chemical reactions. Although the majority of enzymatic networks have very
elaborate structures, there are many experimental observations indicating that
some turnover rates still follow a simple Michaelis-Menten relation with a
hyperbolic dependence on a substrate concentration. The original
Michaelis-Menten mechanism has been derived as a steady-state approximation for
a single-pathway enzymatic chain. The validity of this mechanism for many
complex enzymatic systems is surprising. To determine general conditions when
this relation might be observed in experiments, enzymatic networks consisting
of coupled parallel pathways are investigated theoretically. It is found that
the Michaelis-Menten equation is satisfied for specific relations between
chemical rates, and it also corresponds to the situation with no fluxes between
parallel pathways. Our results are illustrated for simple models. The
importance of the Michaelis-Menten relationship and derived criteria for
single-molecule experimental studies of enzymatic processes are discussed.


  Understanding design principles of biomolecular recognition is a key question
of molecular biology. Yet the enormous complexity and diversity of biological
molecules hamper the efforts to gain a predictive ability for the free energy
of protein-protein, protein-DNA, and protein-RNA binding. Here, using a variant
of the Derrida model, we predict that for a large class of biomolecular
interactions, it is possible to accurately estimate the relative free energy of
binding based on the fluctuation properties of their energy spectra, even if a
finite number of the energy levels is known. We show that the free energy of
the system possessing a wider binding energy spectrum is almost surely lower
compared with the system possessing a narrower energy spectrum. Our predictions
imply that low-affinity binding scores, usually wasted in protein-protein and
protein-DNA docking algorithms, can be efficiently utilized to compute the free
energy. Using the results of Rosetta docking simulations of protein-protein
interactions from Andre et al., Proc. Natl. Acad. Sci. U.S.A. 105, 16148
(2008), we demonstrate the power of our predictions.


  We present a rigidity analysis on a large number of X-ray crystal structures
of the enzyme HIV-1 protease using the 'pebble game' algorithm of the software
FIRST. We find that although the rigidity profile remains similar across a
comprehensive set of high resolution structures, the profile changes
significantly in the presence of an inhibitor. Our study shows that the action
of the inhibitors is to restrict the flexibility of the beta-hairpin flaps
which allow access to the active site. The results are discussed in the context
of full molecular dynamics simulations as well as data from NMR experiments.


  External information propagates in the cell mainly through signaling cascades
and transcriptional activation, allowing it to react to a wide spectrum of
environmental changes. High throughput experiments identify numerous molecular
components of such cascades that may, however, interact through unknown
partners. Some of them may be detected using data coming from the integration
of a protein-protein interaction network and mRNA expression profiles. This
inference problem can be mapped onto the problem of finding appropriate optimal
connected subgraphs of a network defined by these datasets. The optimization
procedure turns out to be computationally intractable in general. Here we
present a new distributed algorithm for this task, inspired from statistical
physics, and apply this scheme to alpha factor and drug perturbations data in
yeast. We identify the role of the COS8 protein, a member of a gene family of
previously unknown function, and validate the results by genetic experiments.
The algorithm we present is specially suited for very large datasets, can run
in parallel, and can be adapted to other problems in systems biology. On
renowned benchmarks it outperforms other algorithms in the field.


  Conventional kinesin is a two-headed homodimeric motor protein, which is able
to walk along microtubules processively by hydrolyzing ATP. Its neck linkers,
which connect the two motor domains and can undergo a docking/undocking
transition, are widely believed to play the key role in the coordination of the
chemical cycles of the two motor domains and, consequently, in force production
and directional stepping. Although many experiments, often complemented with
partial kinetic modeling of specific pathways, support this idea, the ultimate
test of the viability of this hypothesis requires the construction of a
complete kinetic model. Considering the two neck linkers as entropic springs
that are allowed to dock to their head domains and incorporating only the few
most relevant kinetic and structural properties of the individual heads, here
we develop the first detailed, thermodynamically consistent model of kinesin
that can (i) explain the cooperation of the heads (including their gating
mechanisms) during walking and (ii) reproduce much of the available
experimental data (speed, dwell time distribution, randomness, processivity,
hydrolysis rate, etc.) under a wide range of conditions (nucleotide
concentrations, loading force, neck linker length and composition, etc.).
Besides revealing the mechanism by which kinesin operates, our model also makes
it possible to look into the experimentally inaccessible details of the
mechanochemical cycle and predict how certain changes in the protein affect its
motion.


  We analyze the dynamics of evolutionary games in which fitness is defined as
an affine function of the expected payoff and a constant contribution. The
resulting inhomogeneous replicator equation has an homogeneous equivalent with
modified payoffs. The affine terms also influence the stochastic dynamics of a
two-strategy Moran model of a finite population. We then apply the affine
fitness function in a model for tumor-normal cell interactions to determine
which are the most successful tumor strategies. In order to analyze the
dynamics of concurrent strategies within a tumor population, we extend the
model to a three-strategy game involving distinct tumor cell types as well as
normal cells. In this model, interaction with normal cells, in combination with
an increased constant fitness, is the most effective way of establishing a
population of tumor cells in normal tissue.


  The immense increase in the generation of genomic scale data poses an unmet
analytical challenge, due to a lack of established methodology with the
required flexibility and power. We propose a first principled approach to
statistical analysis of sequence-level genomic information. We provide a
growing collection of generic biological investigations that query pairwise
relations between tracks, represented as mathematical objects, along the
genome. The Genomic HyperBrowser implements the approach and is available at
http://hyperbrowser.uio.no.


  We employ Monte Carlo simulations to numerically study the temporal evolution
and transient oscillations of the population densities, the associated
frequency power spectra, and the spatial correlation functions in the
(quasi-)steady state in two-dimensional stochastic May--Leonard models of
mobile individuals, allowing for particle exchanges with nearest-neighbors and
hopping onto empty sites. We therefore consider a class of four-state
three-species cyclic predator-prey models whose total particle number is not
conserved. We demonstrate that quenched disorder in either the reaction or in
the mobility rates hardly impacts the dynamical evolution, the emergence and
structure of spiral patterns, or the mean extinction time in this system. We
also show that direct particle pair exchange processes promote the formation of
regular spiral structures. Moreover, upon increasing the rates of mobility, we
observe a remarkable change in the extinction properties in the May--Leonard
system (for small system sizes): (1) As the mobility rate exceeds a threshold
that separates a species coexistence (quasi-)steady state from an absorbing
state, the mean extinction time as function of system size N crosses over from
a functional form ~ e^{cN} / N (where c is a constant) to a linear dependence;
(2) the measured histogram of extinction times displays a corresponding
crossover from an (approximately) exponential to a Gaussian distribution. The
latter results are found to hold true also when the mobility rates are randomly
distributed.


  Working together in groups may be beneficial if compared to isolated efforts.
Yet this is true only if all group members contribute to the success. If not,
group efforts may act detrimentally on the fitness of their members. Here we
study the evolution of cooperation in public goods games on scale-free networks
that are subject to deletion of links that are connected to the highest-degree
individuals, i.e., on networks that are under attack. We focus on the case
where all groups a player belongs to are considered for the determination of
payoffs; the so-called multi-group public goods games. We find that the effect
of link deletions on the evolution of cooperation is predominantly detrimental,
although there exist regions of the multiplication factor where the existence
of an "optimal" number of removed links for deterioration of cooperation can
also be demonstrated. The findings are explained by means of wealth
distributions and analytical approximations, confirming that socially diverse
states are crucial for the successful evolution of cooperation.


  A theory for thermomechanical behavior of homogeneous DNA at thermal
equilibrium predicts critical temperatures for denaturation under torque and
stretch, phase diagrams for stable B--DNA, supercoiling, optimally stable
torque, and the overstretching transition as force-induced DNA melting.
Agreement with available single molecule manipulation experiments is excellent.


  \textit{Vibrio harveyi} and \textit{Vibrio cholerae} have quorum sensing
pathways with similar design and highly homologous components including
multiple small RNAs (sRNAs). However, the associated luminescence phenotypes of
strains with sRNA deletions differ dramatically: in \textit{V. harveyi}, the
sRNAs act additively; however, in \textit{V. cholerae}, the sRNAs act
redundantly. Furthermore, there are striking differences in the luminescence
phenotypes for different pathway mutants in \textit{V. harveyi} and \textit{V.
cholerae}. However these differences have not been connected with the observed
differences for the sRNA deletion strains in these bacteria. In this work, we
present a model for quorum sensing induced luminescence phenotypes focusing on
the interactions of multiple sRNAs with target mRNA. Within our model, we find
that one key parameter -- the fold-change in protein concentration necessary
for luminescence activation -- can control whether the sRNAs appear to act
additively or redundantly. For specific parameter choices, we find that
differences in this key parameter can also explain hitherto unconnected
luminescence phenotypes differences for various pathway mutants in \textit{V.
harveyi} and \textit{V. cholerae}. The model can thus provide a unifying
explanation for observed differences in luminescence phenotypes and can also be
used to make testable predictions for future experiments.


  We demonstrate a technique based on noise measurements which can be utilized
to study dynamical processes in protein assembly. Direct visualization of
dynamics in membrane protein system such as bacteriorhodopsin (bR) upon
photostimulation are quite challenging. bR represents a model system where the
stimulus-triggered structural dynamics and biological functions are directly
correlated. Our method utilizes a pump-probe near field microscopy method in
the transmission mode and involves analyzing the transmittance fluctuations
from a finite size of molecular assembly. Probability density distributions
indicating the effects of finite size and statistical correlations appear as a
characteristic frequency distribution in the noise spectra of bR whose origin
can be traced to photocycle kinetics. Valuable insight into the molecular
processes were obtained from the noise studies of bR and its mutant D96N as a
function of external parameters such as temperature, humidity or presence of an
additional pump source.


  The sequence diversity of viral populations within individual hosts is the
starting material for selection and subsequent evolution of RNA viruses such as
foot-and-mouth disease virus (FMDV). Using next-generation sequencing (NGS)
performed on a Genome Analyzer platform (Illumina), this study compared the
viral populations within two bovine epithelial samples (foot lesions) from a
single animal with the Inoculum used to initiate experimental infection.
Genomic sequences were determined in duplicate sequencing runs, and the
consensus sequence determined by NGS, for the Inoculum, was identical to that
previously determined using the Sanger method. However, NGS reveals the fine
polymorphic sub-structure of the viral population, from nucleotide variants
present at just below 50% frequency to those present at fractions of 1%. Some
of the higher frequency polymorphisms identified encoded changes within codons
associated with heparan sulphate binding and were present in both feet lesions
revealing intermediate stages in the evolution of a tissue-culture adapted
virus replicating within a mammalian host. We identified 2,622, 1,434 and 1,703
polymorphisms in the Inoculum, and in the two foot lesions respectively: most
of the substitutions occurred only in a small fraction of the population and
represent the progeny from recent cellular replication prior to onset of any
selective pressures. We estimated an upper limit for the genome-wide mutation
rate of the virus within a cell to be 7.8 x 10-4 per nt. The greater depth of
detection, achieved by NGS, demonstrates that this method is a powerful and
valuable tool for the dissection of FMDV populations within-hosts.


  To understand the effect of assortative mating on the genetic evolution of a
population, we consider a finite population in which each individual has a
type, determined by a sequence of n diallelic loci. We assume that the
population evolves according to a Moran model with weak assortative mating,
strong recombination and low mutation rates. With an appropriate rescaling of
time, we obtain that the evolution of the genotypic frequencies in a large
population can be approximated by the evolution of the product of the allelic
frequencies at each locus, and the vector of the allelic frequencies is
approximately governed by a diffusion. We present some features of the limiting
diffusions (in particular their boundary behaviour and conditions under which
the allelic frequencies at different loci evolve independently). If mutation
rates are strictly positive then the limiting diffusion is reversible and,
under some assumptions, the critical points of the stationary density can be
characterised.


  First shells of hydration and bulk solvent plays a crucial role in the
folding of proteins. Here, the role of water in the dynamics of proteins has
been investigated using a theoretical protein-solvent model and a statistical
physics approach. We formulate a hydration model where the hydrogen bonds
between water molecules pertaining to the first shell of the protein
conformation may be either mainly formed or broken. At thermal equilibrium,
hydrogen bonds are formed at low temperature and are broken at high
temperature. To explore the solvent effect, we follow the folding of a large
sampling of protein chains, using a master-equation evolution. The dynamics
shows a clear mechanism. Above the glass-transition temperature, a large ratio
of chains fold very rapidly into the native structure irrespective of the
temperature, following pathways of high transition rates through structures
surrounded by the solvent with broken hydrogen bonds. Although these states
have an infinitesimal probability, they act as strong dynamical attractors and
fast folding proceeds along these routes rather than pathways with small
transition rates between configurations of much higher equilibrium
probabilities. At a given low temperature, a broad jump in the folding times is
observed. Below this glass temperature, the pathways where hydrogen bonds are
mainly formed become those of highest rates although with conformational
changes of huge relaxation times. The present results reveal that folding obeys
a double-funnel mechanism.


  The usefulness of a `total-evidence' approach to human population genetics
was assessed through a clustering analysis of combined genome-wide SNP
datasets. The combination contained only 3146 SNPs. Detailed examination of the
results nonetheless enables the extraction of relevant clues about the history
of human populations, some pertaining to events as ancient as the first
migration out of Africa. The results are mostly coherent with what is known
from history, linguistics, and previous genetic analyses. These promising
results suggest that cross-studies data confrontation have the potential to
yield interesting new hypotheses about human population history.


  We consider control of reaction and population systems by deterministically
imposed transitions between the states with different numbers of particles or
individuals. Even where the imposed transitions are significantly less frequent
than spontaneous transitions, they can exponentially strongly modify the rates
of rare events, including switching between metastable states or population
extinction. We also study optimal control of rare events, and specifically,
optimal control of disease extinction for a limited vaccine supply. A
comparison is made with control of rare events by modulating the rates of
elementary transitions rather than imposing transitions. It is found that,
unexpectedly, for the same mean control parameters, controlling the transitions
rates can be more efficient.


  The structure and stability of strongly charged complex fibers formed by
complexation of a single long semi-flexible polyelectrolyte (PE) chain and many
oppositely charged spherical macroions are investigated numerically at the
ground-state level using a chain-sphere cell model. The model takes into
account chain elasticity as well as electrostatic interactions between charged
spheres and chain segments. Using a numerical optimization method based on a
periodically repeated unit cell, we obtain fiber configurations that minimize
the total energy. The optimal configurations exhibit a variety of helical
structures for the arrangement of macroions including zig-zag, solenoidal and
beads-on-a-string patterns. These structures are determined by a competition
between attraction between spheres and the PE chain (which favors chain
wrapping around the spheres), chain bending and electrostatic repulsion between
chain segments (which favor unwrapping of the chain), and the interactions
between neighboring sphere-chain complexes which can be attractive or repulsive
depending on the system parameters such as medium salt concentration, macroion
charge and chain length per macroion (linker size). At about physiological salt
concentration, dense zig-zag patterns are found to be energetically most stable
when parameters appropriate for the DNA-histone system in a chromatin fiber are
adopted. In fact, the predicted fiber diameter in this regime is found to be
around 30nm, which appears to agree with the thickness observed in in vitro
experiments on chromatin. We also find a macroion density of 5-6 per 11nm which
agrees with the zig-zag or cross-linker models of chromatin. Since our study
deals primarily with a generic model, these findings suggest that
chromatin-like structures should also be observable for PE-macroion complexes
formed in solutions of DNA and synthetic nano-colloids of opposite charge.


  Multistationarity in biological systems is a mechanism of cellular decision
making. In particular, signaling pathways regulated by protein phosphorylation
display features that facilitate a variety of responses to different biological
inputs. The features that lead to multistationarity are of particular interest
to determine as well as the stability properties of the steady states. In this
paper we determine conditions for the emergence of multistationarity in small
motifs without feedback that repeatedly occur in signaling pathways. We derive
an explicit mathematical relationship between the concentration of a chemical
species at steady state and a conserved quantity of the system such as the
total amount of substrate available. We show that the relation determines the
number of steady states and provides a necessary condition for a steady state
to be stable, that is, to be biologically attainable. Further, we identify
characteristics of the motifs that lead to multistationarity, and extend the
view that multistationarity in signaling pathways arises from multisite
phosphorylation. Our approach relies on mass-action kinetics and the
conclusions are drawn in full generality without resorting to simulations or
random generation of parameters. The approach is extensible to other systems.


  We propose and study a class-expansion/innovation/loss model of genome
evolution taking into account biological roles of genes and their constituent
domains. In our model numbers of genes in different functional categories are
coupled to each other. For example, an increase in the number of metabolic
enzymes in a genome is usually accompanied by addition of new transcription
factors regulating these enzymes. Such coupling can be thought of as a
proportional "recipe" for genome composition of the type "a spoonful of sugar
for each egg yolk". The model jointly reproduces two known empirical laws: the
distribution of family sizes and the nonlinear scaling of the number of genes
in certain functional categories (e.g. transcription factors) with genome size.
In addition, it allows us to derive a novel relation between the exponents
characterising these two scaling laws, establishing a direct quantitative
connection between evolutionary and functional categories. It predicts that
functional categories that grow faster-than-linearly with genome size to be
characterised by flatter-than-average family size distributions. This relation
is confirmed by our bioinformatics analysis of prokaryotic genomes. This proves
that the joint quantitative trends of functional and evolutionary classes can
be understood in terms of evolutionary growth with proportional recipes.


  A wealth of new research has highlighted the critical roles of small RNAs
(sRNAs) in diverse processes such as quorum sensing and cellular responses to
stress. The pathways controlling these processes often have a central motif
comprising of a master regulator protein whose expression is controlled by
multiple sRNAs. However, the regulation of stochastic gene expression of a
single target gene by multiple sRNAs is currently not well understood. To
address this issue, we analyze a stochastic model of regulation of gene
expression by multiple sRNAs. For this model, we derive exact analytic results
for the regulated protein distribution including compact expressions for its
mean and variance. The derived results provide novel insights into the roles of
multiple sRNAs in fine-tuning the noise in gene expression. In particular, we
show that, in contrast to regulation by a single sRNA, multiple sRNAs provide a
mechanism for independently controlling the mean and variance of the regulated
protein distribution.


  During the last decade, network approaches became a powerful tool to describe
protein structure and dynamics. Here we review the links between disordered
proteins and the associated networks, and describe the consequences of local,
mesoscopic and global network disorder on changes in protein structure and
dynamics. We introduce a new classification of protein networks into
cumulus-type, i.e., those similar to puffy (white) clouds, and stratus-type,
i.e., those similar to flat, dense (dark) low-lying clouds, and relate these
network types to protein disorder dynamics and to differences in energy
transmission processes. In the first class, there is limited overlap between
the modules, which implies higher rigidity of the individual units; there the
conformational changes can be described by an energy transfer mechanism. In the
second class, the topology presents a compact structure with significant
overlap between the modules; there the conformational changes can be described
by multi-trajectories; that is, multiple highly populated pathways. We further
propose that disordered protein regions evolved to help other protein segments
reach rarely visited but functionally-related states. We also show the role of
disorder in spatial games of amino acids; highlight the effects of
intrinsically disordered proteins (IDPs) on cellular networks and list some
possible studies linking protein disorder and protein structure networks.


  Unsupervised two-view learning, or detection of dependencies between two
paired data sets, is typically done by some variant of canonical correlation
analysis (CCA). CCA searches for a linear projection for each view, such that
the correlations between the projections are maximized. The solution is
invariant to any linear transformation of either or both of the views; for
tasks with small sample size such flexibility implies overfitting, which is
even worse for more flexible nonparametric or kernel-based dependency discovery
methods. We develop variants which reduce the degrees of freedom by assuming
constraints on similarity of the projections in the two views. A particular
example is provided by a cancer gene discovery application where chromosomal
distance affects the dependencies between gene copy number and activity levels.
Similarity constraints are shown to improve detection performance of known
cancer genes.


  The evolution of organismal populations is not typically thought of in terms
of classical mechanics. However, many of the conceptual models used to
approximate evolutionary trajectories have implicit parallels to dynamic
physical systems. The parallels between currently-used evolutionary models and
a type of model related to Lagrangian Coherent Structures (LCS) will be
explored. The limits of evolvability in a population can be treated in a way
analogous to fronts, waves, and other aggregate formations observed in fluid
dynamics. Various measures and architectural features will be introduced.
Relevant scenarios include so-called evolvable boundaries and related scenarios
involving evolutionary neutrality, such as migrations, demographic bottlenecks,
and island biogeography. The LCS-like model introduced here could eventually be
applied to a wide range of problems that normally utilize forms of evolutionary
modeling.


  We demonstrate how niche theory and Hubbell's original formulation of neutral
theory can be blended together into a general framework modeling the combined
effects of selection, drift, speciation, and dispersal on community dynamics.
This framework connects many seemingly unrelated ecological population models,
and allows for quantitative predictions to be made about the impact of niche
stabilizing and destabilizing forces on population extinction times and
abundance distributions. In particular, the existence of niche stabilizing
forces in our blended framework can simultaneously resolve two major problems
with the dynamics of neutral theory, namely predictions of species lifetimes
that are too short and species ages that are too long.


  In this work we monitor the catalytic mechanism of P-glycoprotein (Pgp) using
single-molecule fluorescence resonance energy transfer (FRET). Pgp, a member of
the ATP binding cassette family of transport proteins, is found in the plasma
membrane of animal cells where it is involved in the ATP hydrolysis driven
export of hydrophobic molecules. When expressed in the plasma membrane of
cancer cells, the transport activity of Pgp can lead to the failure of
chemotherapy by excluding the mostly hydrophobic drugs from the interior of the
cell. Despite ongoing effort, the catalytic mechanism by which Pgp couples
MgATP binding and hydrolysis to translocation of drug molecules across the
lipid bilayer is poorly understood. Using site directed mutagenesis, we have
introduced cysteine residues for fluorescence labeling into different regions
of the nucleotide binding domains (NBDs) of Pgp. Double-labeled single Pgp
molecules showed fluctuating FRET efficiencies during drug stimulated ATP
hydrolysis suggesting that the NBDs undergo significant movements during
catalysis. Duty cycle-optimized alternating laser excitation (DCO-ALEX) is
applied to minimize FRET artifacts and to select the appropriate molecules. The
data show that Pgp is a highly dynamic enzyme that appears to fluctuate between
at least two major conformations during steady state turnover.


  We show that the Eigen model and the asexual Wright-Fisher model can be
obtained as different limit cases of a unique stochastic model. This derivation
makes clear which are the exact differences between these two models.
  The two key concepts introduced with the Eigen model, the error threshold and
the quasispecies, are not affected by these differences, so that they are
naturally present also in population genetics models. According to this fact,
in the last part of the paper, we use the classical diploid mutation-selection
equation and the single peak fitness approximation to obtain the error
threshold for sexual diploids. Finally, we compare the results with the asexual
case.


  The paper uses mesoscopic, non-linear lattice dynamics based
(Peyrard-Bishop-Dauxois, PBD) modeling to describe thermal properties of DNA
below and near the denaturation temperature. Computationally efficient notation
is introduced for the relevant statistical mechanics. Computed melting profiles
of long and short heterogeneous sequences are presented, using a recently
introduced reparametrization of the PBD model, and critically discussed. The
statistics of extended open bubbles and bound clusters is formulated and
results are presented for selected examples.


  We consider a simple mathematical model of gradual Darwinian evolution in
continuous time and continuous trait space, due to intraspecific competition
for common resource in an asexually reproducing population in constant
environment, while far from evolutionary stable equilibrium. The model admits
exact analytical solution. In particular, Gaussian distribution of the trait
emerges from generic initial conditions.