  Is there an overriding principle of nature, hitherto overlooked, that governs
all population behavior? A single principle that drives all the regimes
observed in nature - exponential-like growth, saturated growth, population
decline, population extinction, oscillatory behavior? In current orthodox
population theory, this diverse range of population behaviors is described by
many different equations - each with its own specific justification. The
signature of an overriding principle would be a differential equation which, in
a single statement, embraces all the panoply of regimes. A candidate such
governing equation is proposed. The principle from which the equation is
derived is this: The effect on the environment of a population's success is to
alter that environment in a way that opposes the success.


  Many approaches to three-dimensional constrained macromolecular chains at
thermal equilibrium, at about room temperatures, are based upon constrained
Classical Hamiltonian Dynamics (cCHDa). Quantum-mechanical approaches (QMa)
have also been treated by different researchers for decades. QMa address a
fundamental issue (constraints versus the uncertainty principle) and are
versatile: they also yield classical descriptions (which may not coincide with
those from cCHDa, although they may agree for certain relevant quantities).
Open issues include whether QMa have enough practical consequences which differ
from and/or improve those from cCHDa. We shall treat cCHDa briefly and deal
with QMa, by outlining old approaches and focusing on recent ones.


  The stochastic mutual repressor model is analysed using perturbation methods.
This simple model of a gene circuit consists of two genes and three promotor
states. Either of the two protein products can dimerize, forming a repressor
molecule that binds to the promotor of the other gene. When the repressor is
bound to a promotor, the corresponding gene is not transcribed and no protein
is produced. Either one of the promotors can be repressed at any given time or
both can be unrepressed, leaving three possible promotor states. This model is
analysed in its bistable regime in which the deterministic limit exhibits two
stable fixed points and an unstable saddle, and the case of small noise is
considered. On small time scales, the stochastic process fluctuates near one of
the stable fixed points, and on large time scales, a metastable transition can
occur, where fluctuations drive the system past the unstable saddle to the
other stable fixed point. To explore how different intrinsic noise sources
affect these transitions, fluctuations in protein production and degradation
are eliminated, leaving fluctuations in the promotor state as the only source
of noise in the system. Perturbation methods are then used to compute the
stability landscape and the distribution of transition times, or first exit
time density. To understand how protein noise affects the system, small
magnitude fluctuations are added back into the process, and the stability
landscape is compared to that of the process without protein noise. It is found
that significant differences in the random process emerge in the presence of
protein noise.


  Scaffolding is an important subproblem in "de novo" genome assembly in which
mate pair data are used to construct a linear sequence of contigs separated by
gaps. Here we present SLIQ, a set of simple linear inequalities derived from
the geometry of contigs on the line that can be used to predict the relative
positions and orientations of contigs from individual mate pair reads and thus
produce a contig digraph. The SLIQ inequalities can also filter out unreliable
mate pairs and can be used as a preprocessing step for any scaffolding
algorithm. We tested the SLIQ inequalities on five real data sets ranging in
complexity from simple bacterial genomes to complex mammalian genomes and
compared the results to the majority voting procedure used by many other
scaffolding algorithms. SLIQ predicted the relative positions and orientations
of the contigs with high accuracy in all cases and gave more accurate position
predictions than majority voting for complex genomes, in particular the human
genome. Finally, we present a simple scaffolding algorithm that produces linear
scaffolds given a contig digraph. We show that our algorithm is very efficient
compared to other scaffolding algorithms while maintaining high accuracy in
predicting both contig positions and orientations for real data sets.


  Many biological electron transfer (ET) reactions are mediated by metal
centres in proteins. NADH:ubiquinone oxidoreductase (complex I) contains an
intramolecular chain of seven iron-sulphur (FeS) clusters, one of the longest
chains of metal centres in biology and a test case for physical models of
intramolecular ET. In biology, intramolecular ET is commonly described as a
diffusive hopping process, according to the semi-classical theories of Marcus
and Hopfield. However, recent studies have raised the possibility that
non-trivial quantum mechanical effects play a functioning role in certain
biomolecular processes. Here, we extend the semi-classical model for biological
ET to incorporate both semi-classical and coherent quantum phenomena using a
quantum master equation based on the Holstein Hamiltonian. We test our model on
the structurally-defined chain of FeS clusters in complex I. By exploring a
wide range of realistic parameters we find that, when the energy profile for ET
along the chain is relatively flat, just a small coherent contribution can
provide a robust and significant increase in ET rate (above the semi-classical
diffusive-hopping rate), even at physiologically-relevant temperatures.
Conversely, when the on-site energies vary significantly along the chain the
coherent contribution is negligible. For complex I, a crucial respiratory
enzyme that is linked to many neuromuscular and degenerative diseases, our
results suggest a new contribution towards ensuring that intramolecular ET does
not limit the rate of catalysis. For the emerging field of quantum biology, our
model is intended as a basis for elucidating the general role of coherent ET in
biological ET reactions.


  We aim to clarify the relationship between interacting three-species models
and the two-species Lotka-Volterra (LV) model. We utilize mean-field theory and
Monte Carlo simulations on two-dimensional square lattices to explore the
temporal evolution characteristics of two different interacting three-species
predator-prey systems: (1) a cyclic rock-paper-scissors (RPS) model with
conserved total particle number but strongly asymmetric reaction rates that
lets the system evolve towards one corner of configuration space; (2) a
hierarchical food chain where an additional intermediate species is inserted
between the predator and prey in the LV model. For model variant (1), we
demonstrate that the evolutionary properties of both minority species in the
steady state of this stochastic spatial three-species corner RPS model are well
approximated by the LV system, with its emerging characteristic features of
localized population clustering, persistent oscillatory dynamics, correlated
spatio-temporal patterns, and fitness enhancement through quenched spatial
disorder in the predation rates. In contrast, we could not identify any regime
where the hierarchical model (2) would reduce to the two-species LV system. In
the presence of pair exchange processes, the system remains essentially
well-mixed, and we generally find the Monte Carlo simulation results for the
spatially extended model (2) to be consistent with the predictions from the
corresponding mean-field rate equations. If spreading occurs only through
nearest-neighbor hopping, small population clusters emerge; yet the requirement
of an intermediate species cluster obviously disrupts spatio-temporal
correlations between predator and prey, and correspondingly eliminates many of
the intriguing fluctuation phenomena that characterize the stochastic spatial
LV system.


  Protein collapse can be viewed as a dynamical phase transition, during which
new scales and collective variables become excited while the old ones recede
and fade away. This causes formidable computational bottle-necks in approaches
that are based on atomic scale scrutiny. Here we consider an effective
dynamical Landau theory to model the folding process at biologically relevant
time and distance scales. We reach both a substantial decrease in the execution
time and improvement in the accuracy of the final configuration, in comparison
to more conventional approaches. As an example we inspect the collapse of HP35
chicken villin headpiece subdomain, where there are detailed molecular dynamics
simulations to compare with. We start from a structureless, unbend and
untwisted initial configuration. In less than one second of wall-clock time on
a single processor personal computer we consistently reach the native state
with 0.5 Angstrom root mean square distance (RMSD) precision. We confirm that
our folding pathways are indeed akin those obtained in recent atomic level
molecular dynamics simulations. We conclude that our approach appears to have
the potential for a computationally economical method to accurately understand
theoretical aspects of protein collapse.


  The terminase motors of bacteriophages have been shown to be among the
strongest active machines in the biomolecular world, being able to package
several tens of kilobase pairs of viral genome into a capsid within minutes.
Yet these motors are hindered at the end of the packaging process by the
progressive build-up of a force resisting packaging associated with already
packaged DNA. In this experimental work, we raise the issue of what sets the
upper limit on the length of the genome that can be packaged by the terminase
motor of phage {\lambda} and still yield infectious virions, and the conditions
under which this can be efficiently performed. Using a packaging strategy
developed in our laboratory of building phage {\lambda} from scratch, together
with plaque assay monitoring, we have been able to show that the terminase
motor of phage {\lambda} is able to produce infectious particles with up to
110% of the wild-type (WT) {\lambda}-DNA length. However, the phage production
rate, and thus the infectivity, decreased exponentially with increasing DNA
length, and was a factor of 103 lower for the 110% {\lambda}-DNA phage.
Interestingly, our in vitro strategy was still efficient in fully packaging
phages with DNA lengths as high as 114% of the WT length, but these viruses
were unable to infect bacterial cells efficiently. Further, we demonstrated
that the phage production rate is modulated by the presence of multivalent
ionic species. The biological consequences of these finding are discussed.


  This paper investigates questions related to the modularity in discrete
models of biological interaction networks. We develop a theoretical framework
based on the analysis of their asymptotic dynamics. More precisely, we exhibit
formal conditions under which agents of interaction networks can be grouped
into modules. As a main result, we show that the usual decomposition in
strongly connected components fulfils the conditions of being a modular
organisation. Furthermore, we point out that our framework enables a finer
analysis providing a decomposition in elementary modules.


  We develop a theory of aggregation using statistical mechanical methods. An
example of a complicated aggregation system with several levels of structures
is peptide/protein self-assembly. The problem of protein aggregation is
important for the understanding and treatment of neurodegenerative diseases and
also for the development of bio-macromolecules as new materials. We write the
effective Hamiltonian in terms of interaction energies between protein
monomers, protein and solvent, as well as between protein filaments. The grand
partition function can be expressed in terms of a Zimm-Bragg-like transfer
matrix, which is calculated exactly and all thermodynamic properties can be
obtained. We start with two-state and three-state descriptions of protein
monomers using Potts models that can be generalized to include q-states, for
which the exactly solvable feature of the model remains. We focus on n X N
lattice systems, corresponding to the ordered structures observed in some real
fibrils. We have obtained results on nucleation processes and phase diagrams,
in which a protein property such as the sheet content of aggregates is
expressed as a function of the number of proteins on the lattice and
inter-protein or interfacial interaction energies. We have applied our methods
to A{\beta}(1-40) and Curli fibrils and obtained results in good agreement with
experiments.


  An individual-based model of the infectious disease spread among the urban
population is considered. A system of stochastic equations, which describes
changes in quantities of four population groups, susceptible, exposed, infected
individuals and individuals in the state of remission, is built. The system of
equations of the model is supplemented with correlations, which consider
disease heaviness and duration for every infected individual. An algorithm and
a modelling program based on Monte-Carlo methods which allows to investigate
the group number dynamics is developed.


  Recent advances in anti-malware technologies have steered the security
industry away from maintaining vast signature databases and into newer defence
technologies such as behaviour blocking, application whitelisting and others.
Most would agree that the reasoning behind this is to keep up with the arms
race established between malware writers and the security community almost
three decades ago. Still, malware writers have not as yet created new
paradigms. Indeed, malicious code development is still largely limited to code
pattern changes utilizing polymorphic and metamorphic engines, as well as
executable packer and wrapper technologies. Each new malware instance retains
the exact same core functionality as its ancestor and only alters the way it
looks. What if, instead, malware were able to change its function or behaviour
autonomously? What if, in the absence of human intervention, computer viruses
resembled biological viruses in their ability to adapt to new defence
technologies as soon as they came into effect? In this paper, we will provide
the theoretical proof behind malware implementation that closely models
Darwinian evolution.


  Source-sink systems are metapopulations of patches that can be of variable
habitat quality. They can be seen as graphs, where vertices represent the
patches, and the weighted oriented edges give the probability of dispersal from
one patch to another. We consider either finite or source-transitive graphs,
i.e., graphs that are identical when viewed from a(ny) source. We assume
stochastic, individual-based, density-independent reproduction and dispersal.
By studying the path of a single random disperser, we are able to display
simple criteria for persistence, either necessary and sufficient, or just
sufficient. In case of persistence, we characterize the growth rate of the
population as well as the asymptotic occupancy frequencies of the line of
ascent of a random survivor. Our method allows to decouple the roles of
reproduction and dispersal. Finally, we extend our results to the case of
periodic or random environments, where some habitats can have variable growth
rates, autocorrelated in space and possibly in time. In the whole manuscript,
special attention is given to the example of regular graphs where each pair of
adjacent sources is separated by the same number of identical sinks. In the
case of a periodic and random environment, we also display examples where all
patches are sinks when forbidding dispersal but the metapopulation survives
with positive probability in the presence of dispersal, as previously known for
a two-patch mean-field model with parent-independent dispersal.


  We introduce a model for the evolution of species triggered by generation of
novel features and exhaustive combination with other available traits. Under
the assumption that innovations are rare, we obtain a bursty branching process
of speciations. Analysis of the trees representing the branching history
reveals structures qualitatively different from those of random processes. For
a tree with n leaves, the average distance of leaves from root scales as (log
n)^2 to be compared to log n for random branching. The mean values and standard
deviations for the tree shape indices depth (Sackin index) and imbalance
(Colless index) of the model are compatible with those of real phylogenetic
trees from databases. Earlier models, such as the Aldous' branching (AB) model,
show a larger deviation from data with respect to the shape indices.


  We introduce the nested canalyzing depth of a function, which measures the
extent to which it retains a nested canalyzing structure. We characterize the
structure of functions with a given depth and compute the expected activities
and sensitivities of the variables. This analysis quantifies how canalyzation
leads to higher stability in Boolean networks. It generalizes the notion of
nested canalyzing functions (NCFs), which are precisely the functions with
maximum depth. NCFs have been proposed as gene regulatory network models, but
their structure is frequently too restrictive and they are extremely sparse. We
find that functions become decreasingly sensitive to input perturbations as the
canalyzing depth increases, but exhibit rapidly diminishing returns in
stability. Additionally, we show that as depth increases, the dynamics of
networks using these functions quickly approach the critical regime, suggesting
that real networks exhibit some degree of canalyzing depth, and that NCFs are
not significantly better than functions of sufficient depth for many
applications of the modeling and reverse engineering of biological networks.


  By studying the literature about Tetracyclines (TCs), it becomes clearly
evident that TCs are very dynamic molecules. In some cases, their
structure-activity-relationship (SAR) are known, especially against bacteria,
while against other targets, they are virtually unknown. In other diverse
yields of research, such as neurology, oncology and virology the utility and
activity of the tetracyclines are being discovered and are also emerging as new
technological fronts. The first aim of this paper is classify the compounds
already used in therapy and prepare the schematic structure in which include
the next generation of TCs. The aim of this work is introduce a new framework
for the classification of old and new TCs, using a medicinal chemistry approach
to the structure of that drugs. A fully documented
Structure-Activity-Relationship (SAR) is presented with the analysis data of
antibacterial and nonantibacterial (antifungal, antiviral and anticancer)
tetracyclines. Lipophilicity of functional groups and conformations
interchangeably are determining rules in biological activities of TCs.


  We consider a generalized two-species population dynamic model and
analytically solve it for the amensalism and commensalism ecological
interactions. These two-species models can be simplified to a one-species model
with a time dependent extrinsic growth factor. With a one-species model with an
effective carrying capacity one is able to retrieve the steady state solutions
of the previous one-species model. The equivalence obtained between the
effective carrying capacity and the extrinsic growth factor is complete only
for a particular case, the Gompertz model. Here we unveil important aspects of
sigmoid growth curves, which are relevant to growth processes and population
dynamics.


  Pattern formation in microbial colonies of competing strains under purely
space-limited population growth has recently attracted considerable research
interest. We show that the reproduction time statistics of individuals has a
significant impact on the sectoring patterns. Generalizing the standard Eden
growth model, we introduce a simple one-parameter family of reproduction time
distributions indexed by the variation coefficient {\delta} \in [0, 1], which
includes deterministic ({\delta} = 0) and memoryless exponential distribution
({\delta} = 1) as extreme cases. We present convincing numerical evidence and
heuristic arguments that the generalized model is still in the KPZ universality
class, and the changes in patterns are due to changing prefactors in the
scaling relations, which we are able to predict quantitatively. At the example
of Saccharomyces cerevisiae, we show that our approach using the variation
coefficient also works for more realistic reproduction time distributions.


  Transcription factors (TFs) are regulatory proteins that bind DNA in promoter
regions of the genome and either promote or repress gene expression. Here we
predict analytically that enhanced homo-oligonucleotide sequence correlations,
such as poly(dA:dT) and poly(dC:dG) tracts, statistically enhance non-specific
TF-DNA binding affinity. This prediction is generic and qualitatively
independent of microscopic parameters of the model. We show that non-specific
TF binding affinity is universally controlled by the strength and symmetry of
DNA sequence correlations. We perform correlation analysis of the yeast genome
and show that DNA regions highly occupied by TFs exhibit stronger
homo-oligonucleotide sequence correlations, and thus higher propensity for
non-specific binding, as compared with poorly occupied regions. We suggest that
this effect plays the role of an effective localization potential enhancing the
quasi-one-dimensional diffusion of TFs in the vicinity of DNA, speeding up the
stochastic search process for specific TF binding sites. The predicted effect
also imposes an upper bound on the size of TF-DNA binding motifs.


  Summary: The ModuLand plug-in provides Cytoscape users an algorithm for
determining extensively overlapping network modules. Moreover, it identifies
several hierarchical layers of modules, where meta-nodes of the higher
hierarchical layer represent modules of the lower layer. The tool assigns
module cores, which predict the function of the whole module, and determines
key nodes bridging two or multiple modules. The plug-in has a detailed
JAVA-based graphical interface with various colouring options. The ModuLand
tool can run on Windows, Linux, or Mac OS. We demonstrate its use on protein
structure and metabolic networks. Availability: The plug-in and its user guide
can be downloaded freely from: http://www.linkgroup.hu/modules.php. Contact:
csermely.peter@med.semmelweis-univ.hu Supplementary information: Supplementary
information is available at Bioinformatics online.


  We investigate the folding mechanism of the WW domain Fip35 using a realistic
atomistic force field by applying the Dominant Reaction Pathways (DRP)
approach. We find evidence for the existence of two folding pathways, which
differ by the order of formation of the two hairpins. This result is consistent
with the analysis of the experimental data on the folding kinetics of WW
domains and with the results obtained from large-scale molecular dynamics (MD)
simulations of this system. Free-energy calculations performed in two
coarse-grained models support the robustness of our results and suggest that
the qualitative structure of the dominant paths are mostly shaped by the native
interactions. Computing a folding trajectory in atomistic detail only required
about one hour on 48 CPU's. The gain in computational efficiency opens the door
to a systematic investigation of the folding pathways of a large number of
globular proteins.


  Many complex networks from the World-Wide-Web to biological networks are
growing taking into account the heterogeneous features of the nodes. The
feature of a node might be a discrete quantity such as a classification of a
URL document as personal page, thematic website, news, blog, search engine,
social network, ect. or the classification of a gene in a functional module.
Moreover the feature of a node can be a continuous variable such as the
position of a node in the embedding space. In order to account for these
properties, in this paper we provide a generalization of growing network models
with preferential attachment that includes the effect of heterogeneous features
of the nodes. The main effect of heterogeneity is the emergence of an
"effective fitness" for each class of nodes, determining the rate at which
nodes acquire new links. The degree distribution exhibits a multiscaling
behaviour analogous to the the fitness model. This property is robust with
respect to variations in the model, as long as links are assigned through
effective preferential attachment. Beyond the degree distribution, in this
paper we give a full characterization of the other relevant properties of the
model. We evaluate the clustering coefficient and show that it disappears for
large network size, a property shared with the Barab\'asi-Albert model.
Negative degree correlations are also present in the studied class of models,
along with non-trivial mixing patterns among features. We therefore conclude
that both small clustering coefficients and disassortative mixing are outcomes
of the preferential attachment mechanism in general growing networks.


  Different extensions of the classical single-strain SIR model for the host
population, motivated by modeling dengue fever epidemiology, have reported a
rich dynamic structure including deterministic chaos which was able to explain
the large fluctuations of disease incidences. A comparison between the basic
two-strain dengue model, which already captures differences between primary and
secondary infections, with the four-strain dengue model, that introduces the
idea of competition of multiple strains in dengue epidemics shows that the
difference between first and secondary infections drives the rich dynamics more
than the detailed number of strains to be considered in the model structure.
Chaotic dynamics were found to happen at the same parameter region of interest,
for both the two and the four-strain models, able to explain the fluctuations
observed in empirical data and showing a qualitatively good agreement between
empirical data and model simulation. Since the law of parsimony favors the
simplest of two competing models, the two-strain model would be the better
candidate to be analyzed, giving the expected complex behavior to explain the
fluctuations observed in empirical data, and indeed the better option for
estimating all initial conditions as well as the few model parameters based on
the available incidence data.


  Many cellular networks rely on the regulated transport of their components to
transduce extracellular information into precise intracellular signals. The
dynamics of these networks is typically described in terms of compartmentalized
chemical reactions. There are many important situations, however, in which the
properties of the compartments change continuously in a way that cannot
naturally be described by chemical reactions. Here, we develop an approach
based on transport along a trafficking coordinate to precisely describe these
processes and we apply it explicitly to the TGF-{\beta} signal transduction
network, which plays a fundamental role in many diseases and cellular
processes. The results of this newly introduced approach accurately capture for
the first time the distinct TGF-{\beta} signaling dynamics of cells with and
without cancerous backgrounds and provide an avenue to predict the effects of
chemical perturbations in a way that closely recapitulates the observed
cellular behavior.


  Since its emergence in 1968, influenza A (H3N2) has evolved extensively in
genotype and antigenic phenotype. Antigenic evolution occurs in the context of
a two-dimensional 'antigenic map', while genetic evolution shows a
characteristic ladder-like genealogical tree. Here, we use a large-scale
individual-based model to show that evolution in a Euclidean antigenic space
provides a remarkable correspondence between model behavior and the
epidemiological, antigenic, genealogical and geographic patterns observed in
influenza virus. We find that evolution away from existing human immunity
results in rapid population turnover in the influenza virus and that this
population turnover occurs primarily along a single antigenic axis. Thus,
selective dynamics induce a canalized evolutionary trajectory, in which the
evolutionary fate of the influenza population is surprisingly repeatable and
hence, in theory, predictable.


  Understanding the emergence and evolution of multicellularity and cellular
differentiation is a core problem in biology. We develop a quantitative model
that shows that a multicellular form emerges from genetically identical
unicellular ancestors when the compartmentalization of poorly compatible
physiological processes into component cells of an aggregate produces a fitness
advantage. This division of labour between the cells in the aggregate occurs
spontaneously at the regulatory level due to mechanisms present in unicellular
ancestors and does not require any genetic pre-disposition for a particular
role in the aggregate or any orchestrated cooperative behaviour of aggregate
cells. Mathematically, aggregation implies an increase in the dimensionality of
phenotype space that generates a fitness landscape with new fitness maxima, and
in which the unicellular states of optimized metabolism become fitness saddle
points. Evolution of multicellularity is modeled as evolution of a hereditary
parameter, the propensity of cells to stick together, which determines the
fraction of time a cell spends in the aggregate form. Stickiness can increase
evolutionarily due to the fitness advantage generated by the division of labour
between cells in an aggregate.


  We examine the classic problem of homogeneous nucleation and growth by
deriving and analyzing a fully discrete stochastic master equation. Upon
comparison with results obtained from the corresponding mean-field
Becker-D\"{o}ring equations we find striking differences between the two
corresponding equilibrium mean cluster concentrations. These discrepancies
depend primarily on the divisibility of the total available mass by the maximum
allowed cluster size, and the remainder. When such mass incommensurability
arises, a single remainder particle can "emulsify" or "disperse" the system by
significantly broadening the mean cluster size distribution. This finite-sized
broadening effect is periodic in the total mass of the system and can arise
even when the system size is asymptotically large, provided the ratio of the
total mass to the maximum cluster size is finite. For such finite ratios we
show that homogeneous nucleation in the limit of large, closed systems is not
accurately described by classical mean-field mass-action approaches.


  A variety of genome-wide profiling techniques are available to probe
complementary aspects of genome structure and function. Integrative analysis of
heterogeneous data sources can reveal higher-level interactions that cannot be
detected based on individual observations. A standard integration task in
cancer studies is to identify altered genomic regions that induce changes in
the expression of the associated genes based on joint analysis of genome-wide
gene expression and copy number profiling measurements. In this review, we
provide a comparison among various modeling procedures for integrating
genome-wide profiling data of gene copy number and transcriptional alterations
and highlight common approaches to genomic data integration. A transparent
benchmarking procedure is introduced to quantitatively compare the cancer gene
prioritization performance of the alternative methods. The benchmarking
algorithms and data sets are available at http://intcomp.r-forge.r-project.org


  Neutral speciation mechanisms based on isolation by distance and sexual
selection, termed topopatric, have recently been shown to describe the observed
patterns of abundance distributions and species-area relationships. Previous
works have considered this type of process only in the context of hermaphrodic
populations. In this work we extend a hermaphroditic model of topopatric
speciation to populations where individuals are explicitly separated into males
and females. We show that for a particular carrying capacity speciation occurs
under similar conditions, but the number of species generated decreases as
compared to the hermaphroditic case. Evolution results in fewer species having
more abundant populations.


  Quantitative understanding of the principles regulating nucleosome occupancy
on a genome-wide level is a central issue in eukaryotic genomics. Here, we
address this question using budding yeast, Saccharomyces cerevisiae, as a model
organism. We perform a genome-wide computational analysis of nonspecific
transcription factor (TF)-DNA binding free energy landscape, and compare this
landscape with experimentally determined nucleosome binding preferences. We
show that DNA regions with enhanced nonspecific TF-DNA binding are
statistically significantly depleted of nucleosomes. We suggest therefore that
the competition between TFs with histones for nonspecific binding to genomic
sequences might be an important mechanism influencing nucleosome-binding
preferences in vivo. We also predict that poly(dA:dT) and poly(dC:dG) tracts
represent genomic elements with the strongest propensity for nonspecific TF-DNA
binding, thus allowing TFs to outcompete nucleosomes at these elements. Our
results suggest that nonspecific TF-DNA binding might provide a barrier for
statistical positioning of nucleosomes throughout the yeast genome. We predict
that the strength of this barrier increases with the concentration of DNA
binding proteins in a cell. We discuss the connection of the proposed mechanism
with the recently discovered pathway of active nucleosome reconstitution.


  We consider the problem of inferring the unknown parameters of a stochastic
biochemical network model from a single measured time-course of the
concentration of some of the involved species. Such measurements are available,
e.g., from live-cell fluorescence microscopy in image-based systems biology. In
addition, fluctuation time-courses from, e.g., fluorescence correlation
spectroscopy provide additional information about the system dynamics that can
be used to more robustly infer parameters than when considering only mean
concentrations. Estimating model parameters from a single experimental
trajectory enables single-cell measurements and quantification of cell--cell
variability. We propose a novel combination of an adaptive Monte Carlo sampler,
called Gaussian Adaptation, and efficient exact stochastic simulation
algorithms that allows parameter identification from single stochastic
trajectories. We benchmark the proposed method on a linear and a non-linear
reaction network at steady state and during transient phases. In addition, we
demonstrate that the present method also provides an ellipsoidal volume
estimate of the viable part of parameter space and is able to estimate the
physical volume of the compartment in which the observed reactions take place.


  Infectious disease remains, despite centuries of work to control and mitigate
its effects, a major problem facing humanity. This paper reviews the
mathematical modelling of infectious disease epidemics on networks, starting
from the simplest Erdos-Renyi random graphs, and building up structure in the
form of correlations, heterogeneity and preference, paying particular attention
to the links between random graph theory, percolation and dynamical systems
representing transmission. Finally, the problems posed by networks with a large
number of short closed looks are discussed.


  Most of the networks observed in real life obey power-law degree
distribution. It is hypothesized that the emergence of such a degree
distribution is due to preferential attachment of the nodes. Barabasi-Albert
model is a generative procedure that uses preferential attachment based on
degree and one can use this model to generate networks with power-law degree
distribution. In this model, the network is assumed to grow one node every time
step. After the evolution of such a network, it is impossible for one to
predict the exact order of node arrivals. We present in this article, a novel
strategy to partially predict the order of node arrivals in such an evolved
network. We show that our proposed method outperforms other centrality measure
based approaches. We bin the nodes and predict the order of node arrivals
between the bins with an accuracy of above 80%.


  Birth-death processes (BDPs) are continuous-time Markov chains that track the
number of "particles" in a system over time. While widely used in population
biology, genetics and ecology, statistical inference of the instantaneous
particle birth and death rates remains largely limited to restrictive linear
BDPs in which per-particle birth and death rates are constant. Researchers
often observe the number of particles at discrete times, necessitating data
augmentation procedures such as expectation-maximization (EM) to find maximum
likelihood estimates. The E-step in the EM algorithm is available in
closed-form for some linear BDPs, but otherwise previous work has resorted to
approximation or simulation. Remarkably, the E-step conditional expectations
can also be expressed as convolutions of computable transition probabilities
for any general BDP with arbitrary rates. This important observation, along
with a convenient continued fraction representation of the Laplace transforms
of the transition probabilities, allows novel and efficient computation of the
conditional expectations for all BDPs, eliminating the need for approximation
or costly simulation. We use this insight to derive EM algorithms that yield
maximum likelihood estimation for general BDPs characterized by various rate
models, including generalized linear models. We show that our Laplace
convolution technique outperforms competing methods when available and
demonstrate a technique to accelerate EM algorithm convergence. Finally, we
validate our approach using synthetic data and then apply our methods to
estimation of mutation parameters in microsatellite evolution.


  Energetically responsive residues of the 173 amino acid N-terminal domain of
the cardiac Ryanodine receptor RyR2 are identified by a simple elastic net
model. Residues that respond in a correlated way to fluctuations of spatially
neighboring residues specify a hydrogen bonded path through the protein. The
evolutionarily conserved residues of the protein are all located on this path
or in its close proximity. All of the residues of the path are either located
on the two Mir domains of the protein or are hydrogen bonded them. Two calcium
binding residues, E171 and E173, are proposed as potential binding region,
based on insights gained from the elastic net analysis of another calcium
channel receptor, the inositol 1,4,5-triphosphate receptor, IP3R. Analysis of
the disease causing A77V mutated RyR2 showed that the path is disrupted by the
loss of energy responsiveness of certain residues.


  Complex hierarchical structures composed of simple nanoscale building blocks
form the basis of most biological materials. Here we demonstrate how analogies
between seemingly different fields enable the understanding of general
principles by which functional properties in hierarchical systems emerge,
similar to an analogy learning process. Specifically, natural hierarchical
materials like spider silk exhibit properties comparable to classical music in
terms of their hierarchical structure and function. As a comparative tool here
we apply hierarchical ontology logs (olog) that follow a rigorous mathematical
formulation based on category theory to provide an insightful system
representation by expressing knowledge in a conceptual map. We explain the
process of analogy creation, draw connections at several levels of hierarchy
and identify similar patterns that govern the structure of the hierarchical
systems silk and music and discuss the impact of the derived analogy for
nanotechnology.


  Boolean networks serve as discrete models of regulation and signaling in
biological cells. Identifying the key controllers of such processes is
important for understanding the dynamical systems and planning further
analysis. Here we quantify the dynamical impact of a node as the probability of
damage spreading after switching the node's state. We find that the leading
eigenvector of the adjacency matrix is a good predictor of dynamical impact in
the case of long-term spreading. This so-called eigenvector centrality is also
a good proxy measure of the influence a node's initial state has on the
attractor the system eventually arrives at. Quality of prediction is further
improved when eigenvector centrality is based on the weighted matrix of
activities rather than the unweighted adjacency matrix. Simulations are
performed with ensembles of random Boolean networks and a Boolean model of
signaling in fibroblasts. The findings are supported by analytic arguments from
a linear approximation of damage spreading.


  The consequences of recent experimental finding that hydrogen bonds of the
anti-parallel $\beta $-sheet in nonspecific binding site of serine proteases
become significantly shorter and stronger synchronously with the catalytic act
are examined. We investigate the effect of the transformation of an ordinary
hydrogen bond into a low-barrier one on the crankshaft motion a peptide group
in the anti-parallel $\beta $-sheet. For this purpose we make use of a
realistic model of the peptide chain with stringent microscopically derived
coupling interaction potential and effective on-site potential. The coupling
interaction characterizing the peptide chain rigidity is found to be
surprisingly weak and repulsive in character. The effective on-site potential
is found to be a hard one, i.e., goes more steep than a harmonic one. At
transformation of the ordinary hydrogen bond into the low-barrier one the
frequency of crankshaft motion of the corresponding peptide group in the
anti-parallel $\beta $-sheet is roughly doubled.


  Many natural populations are well modelled through time-inhomogeneous
stochastic processes. Such processes have been analysed in the physical
sciences using a method based on Lie algebras, but this methodology is not
widely used for models with ecological, medical and social applications. This
paper presents the Lie algebraic method, and applies it to three biologically
well motivated examples. The result of this is a solution form that is often
highly computationally advantageous.


  We present the Scalable Nucleotide Alignment Program (SNAP), a new short and
long read aligner that is both more accurate (i.e., aligns more reads with
fewer errors) and 10-100x faster than state-of-the-art tools such as BWA.
Unlike recent aligners based on the Burrows-Wheeler transform, SNAP uses a
simple hash index of short seed sequences from the genome, similar to BLAST's.
However, SNAP greatly reduces the number and cost of local alignment checks
performed through several measures: it uses longer seeds to reduce the false
positive locations considered, leverages larger memory capacities to speed
index lookup, and excludes most candidate locations without fully computing
their edit distance to the read. The result is an algorithm that scales well
for reads from one hundred to thousands of bases long and provides a rich error
model that can match classes of mutations (e.g., longer indels) that today's
fast aligners ignore. We calculate that SNAP can align a dataset with 30x
coverage of a human genome in less than an hour for a cost of $2 on Amazon EC2,
with higher accuracy than BWA. Finally, we describe ongoing work to further
improve SNAP.


  The tertiary structure of protein, as well as the local secondary structure
organization are fully determined by the angles of the peptidic bound. The
backbone dihedral angles not only determine the global fold of the protein, but
also the details of the local chain organization. Although a wealth of
structural information is available in different databases and numerous
structural biology softwares have been developed, rapid conformational
characterization remains challenging.
  We present here RamaDA, a program able to give a synthetic description of the
conformation of a protein. The RamaDA program is based on a model where the
Ramachadran plot is decomposed into seven conformational domains. Within the
framework of this model, each amino-acid of a given protein is assigned to one
of these domains.
  From this assignment secondary structure elements can be detected with an
accuracy equivalent to that of the DSSP program for helices and extended
strands, and with the added capability of detecting PolyProline II secondary
structures. Additionally, the determination of a z-score for each amino-acid of
the protein emphasizes any irregularities in the element.
  It is also possible to use this analysis to detect characteristic
conformational patterns. In the case of EF-hands, calcium-binding
helix-loop-helix domains, it is possible to design a strict consensus for the 9
amino-acids of the loop. 523 calcium binding protein files can be found into
the entire PDB with this pattern and only 2.7% false positive hits are
detected.
  The program RamaDA gathers several tools in one and is then able to give a
complete information on a protein structure, including loops and random coil
regions. Through the example of EF-hands, a promising approach of structural
biology is developed. RamaDA is freely available for download as well as online
usage at http://ramada.u-strasbg.fr


  Random heteropolymers do not display the typical equilibrium properties of
globular proteins, but are the starting point to understand the physics of
proteins and, in particular, to describe their non-native states. So far, they
have been studied only with mean-field models in the thermodynamic limit, or
with computer simulations of very small chains on lattice. After describing a
self-adjusting parallel-tempering technique to sample efficiently the
low-energy states of frustrated systems without the need of tuning the
system-dependent parameters of the algorithm, we apply it to random
heteropolymers moving in continuous space. We show that if the mean interaction
between monomers is negative, the usual description through the random energy
model is nearly correct, provided that it is extended to account for
non-compact conformations. If the mean interaction is positive, such a simple
description breaks out and the system behaves in a way more similar to Ising
spin glasses. The former case is a model for the denatured state of glob- ular
proteins, the latter of naturally-unfolded proteins, whose equilibrium
properties thus result qualitatively different.


  Drug resistance remains a major problem for the treatment of HIV. Resistance
can occur due to mutations that were present before treatment starts or due to
mutations that occur during treatment. The relative importance of these two
sources is unknown. We study three different situations in which HIV drug
resistance may evolve: starting triple-drug therapy, treatment with a single
dose of nevirapine and interruption of treatment. For each of these three cases
good data are available from literature, which allows us to estimate the
probability that resistance evolves from standing genetic variation. Depending
on the treatment we find probabilities of the evolution of drug resistance due
to standing genetic variation between 0 and 39%. For patients who start
triple-drug combination therapy, we find that drug resistance evolves from
standing genetic variation in approximately 6% of the patients. We use a
population-dynamic and population-genetic model to understand the observations
and to estimate important evolutionary parameters. We find that both, the
effective population size of the virus before treatment, and the fitness of the
resistant mutant during treatment, are key-parameters that determine the
probability that resistance evolves from standing genetic variation.
Importantly, clinical data indicate that both of these parameters can be
manipulated by the kind of treatment that is used.


  Compositional spectra (CS) analysis based on k-mer scoring of DNA sequences
was employed in this study for dot-plot comparison of human and primate
genomes. The detection of extended conserved synteny regions was based on
continuous fuzzy similarity rather than on chains of discrete anchors (genes or
highly conserved noncoding elements). In addition to the high correspondence
found in the comparisons of whole-genome sequences, a good similarity was also
found after masking gene sequences, indicating that CS analysis manages to
reveal phylogenetic signal in the organization of noncoding part of the genome
sequences, including repetitive DNA and the genome "dark matter". Obviously,
the possibility to reveal parallel ordering depends on the signal of common
ancestor sequence organization varying locally along the corresponding segments
of the compared genomes. We explored two sources contributing to this signal:
sequence composition (GC content) and sequence organization (abundances of
k-mers in the usual A,T,G,C or purine-pyrimidine alphabets). Whole-genome
comparisons based on GC distribution along the analyzed sequences indeed gives
reasonable results, but combining it with k-mer abundances dramatically
improves the ordering quality, indicating that compositional and organizational
heterogeneity comprise complementary sources of information on evolutionary
conserved similarity of genome sequences.


  We consider a stochastic version of the basic predator-prey differential
equation model. The model, which contains a parameter \omega which represents
the number of individuals for one unit of prey -- If x denotes the quantity of
prey in the differential equation model x = 1 means that there are \omega
individuals in the discontinuous one -- is derived from the classical birth and
death process. It is shown by the mean of simulations and explained by a
mathematical analysis based on results in singular perturbation theory (the so
called theory of Canards) that qualitative properties of the model like
persistence or extinction are dramatically sensitive to \omega. For instance,
in our example, if \omega = 107 we have extinction and if \omega = 108 we have
persistence. This means that we must be very cautious when we use continuous
variables in place of jump processes in dynamic population modeling even when
we use stochastic differential equations in place of deterministic ones.


  Actin networks, acting as an engine pushing against an external load, are
fundamentally important to cell motility. A measure of the effectiveness of an
engine is the velocity the engine is able to produce at a given force, the
force-velocity curve. One type of force-velocity curve, consisting of a concave
region where velocity is insensitive to increasing force followed by a decrease
in velocity, is indicative of an adaptive response. In contrast, an engine
whose velocity rapidly decays as a convex curve in response to increasing force
would indicate a lack of adaptive response. Even taken outside of a cellular
context, branching actin networks have been observed to exhibit both concave
and convex force-velocity curves. The exact mechanism that can explain both
force-velocity curves is not yet known. We carried out an agent-based
stochastic simulation to explore such a mechanism. Our results suggest that
upon loading, branching actin networks are capable of remodeling by increasing
the number filaments growing against the load. Our model provides a mechanism
that can account for both convex and concave force-velocity relationships
observed in branching actin networks. Finally, our model gives a potential
explanation to the experimentally observed force history dependence for actin
network velocity.


  A birth-death process is a continuous-time Markov chain that counts the
number of particles in a system over time. In the general process with $n$
current particles, a new particle is born with instantaneous rate $\lambda_n$
and a particle dies with instantaneous rate $\mu_n$. Currently no robust and
efficient method exists to evaluate the finite-time transition probabilities in
a general birth-death process with arbitrary birth and death rates. In this
paper, we first revisit the theory of continued fractions to obtain expressions
for the Laplace transforms of these transition probabilities and make explicit
an important derivation connecting transition probabilities and continued
fractions. We then develop an efficient algorithm for computing these
probabilities that analyzes the error associated with approximations in the
method. We demonstrate that this error-controlled method agrees with known
solutions and outperforms previous approaches to computing these probabilities.
Finally, we apply our novel method to several important problems in ecology,
evolution, and genetics.


  his paper describes and analyzes the spatial spread of tuberculosis (TB) on
complex metapopulation, that is, networks of populations connected by migratory
flows whose configurations are described in terms of connectivity distribution
of nodes (patches) and the conditional probabilities of connections among
classes of nodes sharing the same degree. The migration and transmission
processes occur simultaneously. For uncorrelated networks under the assumption
of standard incidence transmission, we compute the disease-free equilibrium and
the basic reproduction number, and show that the disease-free equilibrium is
locally asymptotically stable. Moreover, for uncorrelated networks and under
assumption of simple mass action transmission, we give a necessary and
sufficient conditions for the instability of the disease-free equilibrium.
  The existence of endemic equilibria is also discussed. Finally, the
prevalence of the TB infection across the metapopulation as a function of the
path connectivity is studied using numerical simulations.


  We analyze theoretically the problem of cargo transport along microtubules by
motors of two species with opposite polarities. We consider two different
one-dimensional models previously developed in the literature. On the one hand,
a quite widespread model which assumes equal force sharing, here referred to as
mean field model (MFM). On the other hand, a stochastic model (SM) which
considers individual motor-cargo links. We find that in generic situations the
MFM predicts larger cargo mean velocity, smaller mean run time and less
frequent reversions than the SM. These phenomena are found to be consequences
of the load sharing assumptions and can be interpreted in terms the
probabilities of the different motility states. We also explore the influence
of the viscosity in both models and the role of the stiffness of the
motor-cargo links within the SM. Our results show that the mean cargo velocity
is independent of the stiffness while the mean run time decreases with such a
parameter. We explore the case of symmetric forward and backward motors
considering kinesin- 1 parameters, and the problem of transport by kinesin-1
and cytoplasmic dyneins considering two different sets of parameters previously
proposed for dyneins.


  DNA based nanostructures built on a long single stranded DNA scaffold, known
as DNA origamis, offer the possibility to organize various molecules at the
nanometer scale in one pot experiments.
  The folding of the scaffold is guaranteed by the presence of short, single
stranded DNA sequences (staples), that hold together separate regions of the
scaffold. In this paper, we modelize the annealing-melting properties of these
DNA constructions. The model captures important features such as the hysteresis
between melting and annealing, as well as the dependence upon the topology of
the scaffold. We show that cooperativity between staples is critical to
quantitatively explain the folding process of DNA origamis.


  During this era of new drug designing, medicinal plants had become a very
interesting object of further research. Pharmacology screening of active
compound of medicinal plants would be time consuming and costly. Molecular
docking is one of the in silico method which is more efficient compare to in
vitro or in vivo method for its capability of finding the active compound in
medicinal plants. In this method, three-dimensional structure becomes very
important in the molecular docking methods, so we need a database that provides
information on three-dimensional structures of chemical compounds from
medicinal plants in Indonesia. Therefore, this study will prepare a database
which provides information of the three dimensional structures of chemical
compounds of medicinal plants. The database will be prepared by using MySQL
format and is designed to be placed in http://herbaldb.farmasi.ui.ac.id website
so that eventually this database can be accessed quickly and easily by users
via the Internet.


  Complex network dynamics have been analyzed with models of systems of coupled
switches or systems of coupled oscillators. However, many complex systems are
composed of components with diverse dynamics whose interactions drive the
system's evolution. We, therefore, introduce a new modeling framework that
describes the dynamics of networks composed of both oscillators and switches.
Both oscillator synchronization and switch stability are preserved in these
heterogeneous, coupled networks. Furthermore, this model recapitulates the
qualitative dynamics for the yeast cell cycle consistent with the hypothesized
dynamics resulting from decomposition of the regulatory network into dynamic
motifs. Introducing feedback into the cell-cycle network induces qualitative
dynamics analogous to limitless replicative potential that is a hallmark of
cancer. As a result, the proposed model of switch and oscillator coupling
provides the ability to incorporate mechanisms that underlie the synchronized
stimulus response ubiquitous in biochemical systems.


  Disease awareness in infection dynamics can be modeled with adaptive contact
networks whose rewiring rules reflect the attempt by susceptibles to avoid
infectious contacts. Simulations of this type of models show an active phase
with constant infected node density in which the interplay of disease dynamics
and link rewiring prompts the convergence towards a well defined degree
distribution, irrespective of the initial network topology. We develop a method
to study this dynamic equilibrium and give an analytic description of the
structure of the characteristic degree distributions and other network
measures. The method applies to a broad class of systems and can be used to
determine the steady-state topology of many other adaptive networks.


  To provide the Cytoscape users the possibility of integrating ITM Probe into
their workflows, we developed CytoITMprobe, a new Cytoscape plugin.
CytoITMprobe maintains all the desirable features of ITM Probe and adds
additional flexibility not achievable through its web service version. It
provides access to ITM Probe either through a web server or locally. The input,
consisting of a Cytoscape network, together with the desired origins and/or
destinations of information and a dissipation coefficient, is specified through
a query form. The results are shown as a subnetwork of significant nodes and
several summary tables. Users can control the composition and appearance of the
subnetwork and interchange their ITM Probe results with other software tools
through tab-delimited files.
  The main strength of CytoITMprobe is its flexibility. It allows the user to
specify as input any Cytoscape network, rather than being restricted to the
pre-compiled protein-protein interaction networks available through the ITM
Probe web service. Users may supply their own edge weights and
directionalities. Consequently, as opposed to ITM Probe web service,
CytoITMprobe can be applied to many other domains of network-based research
beyond protein-networks. It also enables seamless integration of ITM Probe
results with other Cytoscape plugins having complementary functionality for
data analysis.


  A toggle switch consists of two genes that mutually repress each other. This
regulatory motif is active during cell differentiation and is thought to act as
a memory device, being able to choose and maintain cell fate decisions. In this
contribution, we study the stability and dynamics of a two-stage gene
expression switch within a probabilistic framework inspired by the properties
of the Pu/Gata toggle switch in myeloid progenitor cells. We focus on low mRNA
numbers, high protein abundance and monomeric transcription factor binding.
Contrary to the expectation from a deterministic description, this switch shows
complex multi-attractor dy- namics without autoactivation and cooperativity.
Most importantly, the four attractors of the system, which only emerge in a
probabilistic two-stage description, can be identified with committed and
primed states in cell differentiation. We first study the dynamics of the
system and infer the mechanisms that move the system between attractors using
both the quasi-potential and the probability flux of the system. Second, we
show that the residence times of the system in one of the committed attractors
are geometrically distributed and provide an analytical expression of the
distribution. Most importantly we find that the mean residence time increases
linearly with the mean protein level. Finally, we study the implications of
this distribution for the stability of a switch and discuss the influence of
the stability on a specific cell differentiation mechanism. Our model explains
lineage priming and proposes the need of either high protein numbers or long
term modifications such as chromatin remodeling to achieve stable cell fate
decisions. Notably we present a system with high protein abundance that
nevertheless requires a probabilistic description to exhibit multistability,
complex switching dynamics and lineage priming.


  In a recent article, Behrens and Vingron (JCB 17, 12, 2010) compute waiting
times for k-mers to appear during DNA evolution under the assumption that the
considered k-mers do not occur in the initial DNA sequence, an issue arising
when studying the evolution of regulatory DNA sequences with regard to
transcription factor (TF) binding site emergence. The mathematical analysis
underlying their computation assumes that occurrences of words under interest
do not overlap. We relax here this assumption by use of an automata approach.
In an alphabet of size 4 like the DNA alphabet, most words have no or a low
autocorrelation; therefore, globally, our results confirm those of Behrens and
Vingron. The outcome is quite different when considering highly autocorrelated
k-mers; in this case, the autocorrelation pushes down the probability of
occurrence of these k-mers at generation 1 and, consequently, increases the
waiting time for apparition of these k-mers up to 40%. An analysis of existing
TF binding sites unveils a significant proportion of k-mers exhibiting
autocorrelation. Thus, our computations based on automata greatly improve the
accuracy of predicting waiting times for the emergence of TF binding sites to
appear during DNA evolution. We do the computation in the Bernoulli or M0
model; computations in the M1 model, a Markov model of order 1, are more costly
in terms of time and memory but should produce similar results. While Behrens
and Vingron considered specifically promoters of length 1000, we extend the
results to promoters of any size; we exhibit the property that the probability
that a k-mer occurs at generation time 1 while being absent at time 0 behaves
linearly with respect to the length of the promoter, which induces a hyperbolic
behaviour of the waiting time of any k-mer with respect to the length of the
promoter.


  Cellular signal transduction usually involves activation cascades, the
sequential activation of a series of proteins following the reception of an
input signal. Here we study the classic model of weakly activated cascades and
obtain analytical solutions for a variety of inputs. We show that in the
special but important case of optimal-gain cascades (i.e., when the
deactivation rates are identical) the downstream output of the cascade can be
represented exactly as a lumped nonlinear module containing an incomplete gamma
function with real parameters that depend on the rates and length of the
cascade, as well as parameters of the input signal. The expressions obtained
can be applied to the non-identical case when the deactivation rates are random
to capture the variability in the cascade outputs. We also show that cascades
can be rearranged so that blocks with similar rates can be lumped and
represented through our nonlinear modules. Our results can be used both to
represent cascades in computational models of differential equations and to fit
data efficiently, by reducing the number of equations and parameters involved.
In particular, the length of the cascade appears as a real-valued parameter and
can thus be fitted in the same manner as Hill coefficients. Finally, we show
how the obtained nonlinear modules can be used instead of delay differential
equations to model delays in signal transduction.


  Biopolymer looping is a dynamic process that occurs ubiquitously in cells for
gene regulation, protein folding, etc. In cellular environments, biopolymers
are often subject to tensions which are either static, or temporally
fluctuating far away from equilibrium. We study the dynamics of semiflexible
polymer looping in the presence of such tensions by using Brownian dynamics
simulation combined with an analytical theory. We show a minute tension
dramatically changes the looping time, especially for long chains. Considering
a dichotomically flipping noise as a simple example of the nonequilibrium
tension, we find the phenomenon of resonant activation, where the looping time
can be the minimum at an optimal flipping time. We discuss our results in
connection with recent experiments.


  The simple and partial Mantel tests are routinely used in many areas of
evolutionary biology to assess the significance of the association between two
or more matrices of distances relative to the same pairs of individuals or
demes. Partial Mantel tests rather than simple Mantel tests are widely used to
assess the relationship between two variables displaying some form of
structure.
  We show that contrarily to a widely shared belief, partial Mantel tests are
not valid in this case, and their bias remains close to that of the simple
Mantel test.
  We confirm that strong biases are expected under a sampling design and
spatial correlation parameter drawn from an actual study.
  The Mantel tests should not be used in case auto-correlation is suspected in
both variables compared under the null hypothesis. We outline alternative
strategies. The R code used for our computer simulations is distributed as
supporting material.


  In this paper we explore the life expectancy at birth in the Netherlands by
based on a recent theory and a new methodology but also a classical theory of
fitting and forecasting. We use the data from 1850 to 2006 provided by the
Human Mortality Database (HMD) for the annual deaths per year of age and the
structure of the population per year of age. We apply the IM first exit time
model which includes also the infant mortality by using the appropriate
non-linear regression analysis program that we have developed in Excel. We
provide all the related material from the website: http://www.cmsim.net


  Using a simple example of biological macromolecules which are partitioned
between bulk solution and membrane, we investigate T.L. Hill's phenomenological
nanothermodynamics for small systems. By introducing a {\em systems size
dependent} equilibrium constant for the bulk-membrane partition, we obtain
Hill's results on differential and integral chemical potentials $\mu$ and
$\hat{\mu}$ from computations based on standard Gibbsian equilibrium
statistical mechanics. It is shown that their difference can be understood from
an equilibrium re-partitioning between bulk and membrane fractions upon a
change in system's size; it is closely related to the system's fluctuations and
inhomogeneity. These results provide a better understanding of the
nanothermodynamics and clarify its logical relation with the theory of
statistical mechanics.


  In cellular environment, confinement and macromulecular crowding play an
important role on thermal stability and folding kinetics of a protein. We have
resorted to a generalized version of the Wako-Saito-Munoz-Eaton model for
protein folding to study the behavior of six different protein structures
confined between two walls. Changing the distance 2R between the walls, we
found, in accordance with previous studies, two confinement regimes: starting
from large R and decreasing R, confinement first enhances the stability of the
folded state as long as this is compact and until a given value of R; then a
further decrease of R leads to a decrease of folding temperature and folding
rate. We found that in the low confinement regime both unfolding temperatures
and logarithm of folding rates scale as R-{\gamma} where {\gamma} values lie in
between 1.42 and 2.35.


  We consider exact enumerations and probabilistic properties of ranked trees
when generated under the random coalescent process. Using a new approach, based
on generating functions, we derive several statistics such as the exact
probability of finding k cherries in a ranked tree of fixed size n. We then
extend our method to consider also the number of pitchforks. We find a
recursive formula to calculate the joint and conditional probabilities of
cherries and pitch- forks when the size of the tree is fixed.


  Most biological systems are formed by component parts that to some degree are
inter-related. Groups of parts that are more associated among themselves and
are relatively autonomous from others are called modules. One of the
consequences of modularity is that biological systems usually present an
unequal distribution of the genetic variation among variables. Estimating the
covariance matrix that describes these systems is a difficult problem due to a
number of factors such as poor sample sizes and measurement errors. We show
that this problem will be exacerbated whenever matrix inversion is required, as
in directional selection reconstruction analysis. We explore the consequences
of varying degrees of modularity and signal-to-noise ratio on selection
reconstruction. We then present and test the efficiency of available methods
for controlling noise in matrix estimates. In our simulations, controlling
matrices for noise vastly improves the reconstruction of selection gradients.
We also perform an analysis of selection gradients reconstruction over a New
World Monkeys skull database in order to illustrate the impact of noise on such
analyses. Noise- controlled estimates render far more plausible interpretations
that are in full agreement with previous results.


  We investigate the possibility that prebiotic homochirality can be achieved
exclusively through chiral-selective reaction rate parameters without any other
explicit mechanism for chiral bias. Specifically, we examine an open network of
polymerization reactions, where the reaction rates can have chiral-selective
values. The reactions are neither autocatalytic nor do they contain explicit
enantiomeric cross-inhibition terms. We are thus investigating how rare a set
of chiral-selective reaction rates needs to be in order to generate a
reasonable amount of chiral bias. We quantify our results adopting a
statistical approach: varying both the mean value and the rms dispersion of the
relevant reaction rates, we show that moderate to high levels of chiral excess
can be achieved with fairly small chiral bias, below 10%. Considering the
various unknowns related to prebiotic chemical networks in early Earth and the
dependence of reaction rates to environmental properties such as temperature
and pressure variations, we argue that homochirality could have been achieved
from moderate amounts of chiral selectivity in the reaction rates.


  Metastasis is one of the most enigmatic aspects of cancer pathogenesis and is
a major cause of cancer-associated mortality. Secondary bone cancer (SBC) is a
complex disease caused by metastasis of tumor cells from their primary site and
is characterized by intricate interplay of molecular interactions.
Identification of targets for multifactorial diseases such as SBC, the most
frequent complication of breast and prostate cancers, is a challenge. Towards
achieving our aim of identification of targets specific to SBC, we constructed
a 'Cancer Genes Network', a representative protein interactome of cancer genes.
Using graph theoretical methods, we obtained a set of key genes that are
relevant for generic mechanisms of cancers and have a role in biological
essentiality. We also compiled a curated dataset of 391 SBC genes from
published literature which serves as a basis of ontological correlates of
secondary bone cancer. Building on these results, we implement a strategy based
on generic cancer genes, SBC genes and gene ontology enrichment method, to
obtain a set of targets that are specific to bone metastasis. Through this
study, we present an approach for probing one of the major complications in
cancers, namely, metastasis. The results on genes that play generic roles in
cancer phenotype, obtained by network analysis of 'Cancer Genes Network', have
broader implications in understanding the role of molecular regulators in
mechanisms of cancers. Specifically, our study provides a set of potential
targets that are of ontological and regulatory relevance to secondary bone
cancer.


  Chargaff once said that "I saw before me in dark contours the beginning of a
grammar of Biology". In linguistics, "grammar" is the set of natural language
rules, but we do not know for sure what Chargaff meant by "grammar" of Biology.
Nevertheless, assuming the metaphor, Chargaff himself started a "grammar of
Biology" discovering the so called Chargaff's rules. In this work, we further
develop his grammar. Using new concepts, we were able to discovery new genomic
rules that seem to be invariant across a large set of organisms, and show a
fractal-like property, since no matter the scale, the same pattern is observed
(self-similarity). We hope that these new invariant genomic rules may be used
in different contexts since short read data bias detection to genome assembly
quality assessment.


  We analyse metapopulation dynamics in terms of an individual-based,
stochastic model of a finite metapopulation. We suggest a new approach, using
the number of patches in the population as a large parameter. This approach
does not require that the number of individuals per patch is large, neither is
it necessary to assume a time-scale separation between local population
dynamics and migration. Our approach makes it possible to accurately describe
the dynamics of metapopulations consisting of many small patches. We focus on
metapopulations on the brink of extinction. We estimate the time to extinction
and describe the most likely path to extinction. We find that the logarithm of
the time to extinction is proportional to the product of two vectors, a vector
characterising the distribution of patch population sizes in the quasi-steady
state, and a vector -- related to Fisher's reproduction vector -- that
quantifies the sensitivity of the quasi-steady state distribution to
demographic fluctuations. We compare our analytical results to stochastic
simulations of the model, and discuss the range of validity of the analytical
expressions. By identifying fast and slow degrees of freedom in the
metapopulation dynamics, we show that the dynamics of large metapopulations
close to extinction is approximately described by a deterministic equation
originally proposed by Levins (1969). We were able to compute the rates in
Levins' equation in terms of the parameters of our stochastic, individual-based
model. It turns out, however, that the interpretation of the dynamical variable
depends strongly on the intrinsic growth rate and carrying capacity of the
patches. Only when the growth rate and the carrying capacity are large does the
slow variable correspond to the number of patches, as envisaged by Levins. Last
but not least, we discuss how our findings relate to other, widely used
metapopulation models.


  We introduce a nonlinear structured population model with diffusion in the
state space. Individuals are structured with respect to a continuous variable
which represents a pathogen load. The class of uninfected individuals
constitutes a special compartment that carries mass, hence the model is
equipped with generalized Wentzell (or dynamic) boundary conditions. Our model
is intended to describe the spread of infection of a vertically transmitted
disease, for example Wolbachia in a mosquito population. Therefore the
(infinite dimensional) nonlinearity arises in the recruitment term. First we
establish global existence of solutions and the Principle of Linearised
Stability for our model. Then, in our main result, we formulate simple
conditions, which guarantee the existence of non-trivial steady states of the
model. Our method utilizes an operator theoretic framework combined with a
fixed point approach. Finally, in the last section we establish a sufficient
condition for the local asymptotic stability of the positive steady state.


  Game-theoretical models where the rules of the game and the interaction
structure both coevolves with the game dynamics -- multiadaptive games --
capture very flexible situations where cooperation among selfish agents can
emerge. In this work, we will discuss a multiadaptive model presented in a
recent Letter [Phys. Rev. Lett. 106, 028702 (2011)], and generalizations of it.
The model captures a non-equilibrium situation where social unrest increases
the incentive to cooperate and, simultaneously, agents are partly free to
influence with whom they interact. First, we investigate the details of how the
feedback from the behavior of agents determines the emergence of cooperation
and hierarchical contact structures. We also study the stability of the system
to different types of noise, and find that different regions of parameter space
show very different response. Some types of noise can destroy an all-cooperator
state. If, on the other hand, hubs are stable, then so is the all-C state.
Finally, we investigate the dependence of the ratio between the timescales of
strategy updates and the evolution of the interaction structure. We find that a
comparatively fast strategy dynamics is a prerequisite for the emergence of
cooperation.


  The purpose of the present essay is to suggest a possible model to describe
the worldwide healthcare crisis, where diseases that have been considered to be
eradicated or under our control are re-emerging today.


  Mitochondrial hypervariable region I genetic data from ancient populations at
two sites from Asia, Linzi in Shandong (northern China) and Egyin Gol in
Mongolia, were reanalyzed to detect population affinities. Data from a total of
51 modern populations were used to generate distance measures (Fst's) to the
two ancient populations. The tests first analyzed relationships at the regional
level, and then compiled the top regional matches for an overall comparison to
the two probe populations. The reanalysis showed that the Egyin Gol and Linzi
populations have clear distinctions in genetic affinity. The Egyin Gol
population as a whole appears to bear close affinities with modern populations
of northern East Asia. The Linzi population does seem to have some genetic
affinities with the West as suggested by the original analysis, though the
original attribution of "European-like" seems to be misleading. This study
suggests that the Linzi individuals are potentially related to early Iranians,
who are thought to have been widespread in parts of Central Eurasia and the
steppe regions in the first millennium BC, though some significant admixture
between a number of populations of varying origin cannot be ruled out. The
study also examines the effect of sequence length on this type of genetic data
analysis and provides analysis and explanation for the results of previous
studies on the Linzi sample as compared to this one.


  Mitochondrial DNA from fourteen archaeological samples at the Ural State
University in Yekaterinburg, Russia was extracted to test the feasibility of
ancient DNA work on their collection. These samples come from a number of sites
that fall into two groupings. Seven samples are from three sites that belong to
a northern group of what are thought to be Ugrians dating to the 8th-12th
century AD, who lived along the Ural Mountains in northwestern Siberia. The
remaining seven samples are from two sites that belong to a southern group
representing the Sargat culture, dating between roughly the 5th century BC and
the 5th century AD, from southwestern Siberia near the Ural Mountains and the
present-day Kazakhstan border. The samples derived from several burial types,
including kurgan burials. They also represented a number of different skeletal
elements, as well as a range of observed preservation. The northern sites
repeatedly failed to amplify after multiple extraction and amplification
attempts, but the samples from the southern sites were successfully extracted
and amplified. The sequences obtained from the southern sites support the
hypothesis that the Sargat culture was a potential zone of intermixture between
native Ugrian and/or Siberian populations and steppe peoples from the South,
possibly early Iranian or Indo-Iranian, which has been previously suggested by
archaeological analysis.


  We introduce a class of weighted graphs whose properties are meant to mimic
the topological features of idiotypic networks, namely the interaction networks
involving the B-core of the immune system. Each node is endowed with a
bit-string representing the idiotypic specificity of the corresponding B cell
and a proper distance between any couple of bit-strings provides the coupling
strength between the two nodes. We show that a biased distribution of the
entries in bit-strings can yield fringes in the (weighted) degree distribution,
small-worlds features, and scaling laws, in agreement with experimental
findings. We also investigate the role of ageing, thought of as a progressive
increase in the degree of bias in bit-strings, and we show that it can possibly
induce mild percolation phenomena, which are investigated too.


  We study complexity of rearrangement problems in the generalized breakpoint
model and settle several open questions. The model was introduced by Tannier et
al. (2009) who showed that the median problem is solvable in polynomial time in
the multichromosomal circular and mixed breakpoint models. This is intriguing,
since in most other rearrangement models (DCJ, reversal, unichromosomal or
multilinear breakpoint models), the problem is NP-hard. The complexity of the
small or even the large phylogeny problem under the breakpoint distance
remained an open problem. We improve the algorithm for the median problem and
show that it is equivalent to the problem of finding maximum cardinality
non-bipartite matching (under linear reduction). On the other hand, we prove
that the more general small phylogeny problem is NP-hard. Surprisingly, we show
that it is already NP-hard (or even APX-hard) for 4 species (a quartet
phylogeny). In other words, while finding an ancestor for 3 species is easy,
finding two ancestors for 4 species is already hard. We also show that, in the
unichromosomal and the multilinear breakpoint model, the halving problem is
NP-hard, thus refuting the conjecture of Tannier et al. Interestingly, this is
the first problem which is harder in the breakpoint model than in the DCJ or
reversal models.


  I consider the microscopic mechanisms by which a particular left-right (L/R)
asymmetry is generated at the organism level from the microscopic handedness of
cytoskeletal molecules. In light of a fundamental symmetry principle, the
typical pattern-formation mechanisms of diffusion plus regulation cannot
implement the "right-hand rule"; at the microscopic level, the cell's
cytoskeleton of chiral filaments seems always to be involved, usually in
collective states driven by polymerization forces or molecular motors. It seems
particularly easy for handedness to emerge in a shear or rotation in the
background of an effectively two-dimensional system, such as the cell membrane
or a layer of cells, as this requires no pre-existing axis apart from the layer
normal. I detail a scenario involving actin/myosin layers in snails and in C.
elegans, and also one about the microtubule layer in plant cells. I also survey
the other examples that I am aware of, such as the emergence of handedness such
as the emergence of handedness in neurons, in eukaryote cell motility, and in
non-flagellated bacteria.


  Channels are Catalysts for Diffusion and in that sense are Enzymes. This idea
is useful for the design and interpretation of experiments.


  Established already in the Biblical times, the Matthew effect stands for the
fact that in societies rich tend to get richer and the potent even more
powerful. Here we investigate a game theoretical model describing the evolution
of cooperation on structured populations where the distribution of public goods
is driven by the reproductive success of individuals. Phase diagrams reveal
that cooperation is promoted irrespective of the uncertainty by strategy
adoptions and the type of interaction graph, yet the complete dominance of
cooperators is elusive due to the spontaneous emergence of super-persistent
defectors that owe their survival to extremely rare microscopic patterns. This
indicates that success-driven mechanisms are crucial for effectively harvesting
benefits from collective actions, but that they may also account for the
observed persistence of maladaptive behavior.


  This work concerns the identification of the allelic sequences of the DQB1
gene of the human leukocyte antigen system conferring susceptibility to the
development of insulin-dependent diabetes mellitus (IDDM) in DNA samples with
no need of PCR amplification. Our method is based on the time-resolved analysis
of a F\"orster energy-transfer mechanism that occurs in a dual-labeled
fluorescent probe specific for the base sequence of the allelic variant of
interest. Such an oligonucleotide probe is labeled, at the two ends, by a pair
of chromophores that operate as donor and acceptor in a F\"orster resonant
energy-transfer. The donor fluorescence is quenched with an efficiency that is
strongly dependent on the donor-to-acceptor distance, hence on the
configuration of the probe after hybridization with the DNA containing or not
the selected allelic sequence. By time-correlated single-photon counting,
performed with an excitation/detection system endowed with 30-ps resolution, we
measure the time-resolved fluorescence decay of the donor and discriminate, by
means of the decay time value, the DNA bearing the allele conferring
susceptibility to IDDM from the DNAs bearing any other sequence in the same
region of the DQB1 gene.


  The evolution of the hemagglutinin amino acids sequences of Influenza A virus
is studied by a method based on an informational metrics, originally introduced
by Rohlin for partitions in abstract probability spaces. This metrics does not
require any previous functional or syntactic knowledge about the sequences and
it is sensitive to the correlated variations in the characters disposition. Its
efficiency is improved by algorithmic tools, designed to enhance the detection
of the novelty and to reduce the noise of useless mutations. We focus on the
USA data from 1993/94 to 2010/2011 for A/H3N2 and on USA data from 2006/07 to
2010/2011 for A/H1N1 . We show that the clusterization of the distance matrix
gives strong evidence to a structure of domains in the sequence space, acting
as weak attractors for the evolution, in very good agreement with the
epidemiological history of the virus. The structure proves very robust with
respect to the variations of the clusterization parameters, and extremely
coherent when restricting the observation window. The results suggest an
efficient strategy in the vaccine forecast, based on the presence of
"precursors" (or "buds") populating the most recent attractor.


  We consider SIR and SIRS models with differential mortality. Global stability
of equilibria is established by using Lyapunov's method


  This paper descibes a new method for deriving incidence rates of a chronic
disease from prevalence data. It is based on a new ordinary differential
equation, which relates the change in the age-specific prevalence to the
agespecific incidence and mortality rates. The method allows the extraction of
longtudinal information from cross-sectional studies. Applicability of the
method is tested in the prevalence of dementia in Germany. The derived
age-specific incidence is in good agreement with published values.


  George Williams defined an evolutionary unit as hereditary information for
which the selection bias between competing units dominates the informational
decay caused by imperfect transmission. In this article, I extend Williams'
approach to show that the ratio of selection bias to transmission bias provides
a unifying framework for diverse biological problems. Specific examples include
Haldane and Lande's mutation-selection balance, Eigen's error threshold and
quasispecies, Van Valen's clade selection, Price's multilevel formulation of
group selection, Szathmary and Demeter's evolutionary origin of primitive
cells, Levin and Bull's short-sighted evolution of HIV virulence, Frank's
timescale analysis of microbial metabolism, and Maynard Smith and Szathmary's
major transitions in evolution. The insights from these diverse applications
lead to a deeper understanding of kin selection, group selection, multilevel
evolutionary analysis, and the philosophical problems of evolutionary units and
individuality.


  Global transposable characteristics in the complete DNA sequence of the
Saccharomyces cevevisiae yeast is determined by using the metric representation
and recurrence plot methods. In the form of the correlation distance of
nucleotide strings, 16 chromosome sequences of the yeast, which are divided
into 5 groups, display 4 kinds of the fundamental transposable characteristics:
a short period increasing, a long quasi-period increasing, a long major value
and hardly relevant.


  The Luria-Delbr\"uck mutation model has a long history and has been
mathematically formulated in several different ways. Here we tackle the problem
in the case of a continuous distribution using some mathematical tools from
nonlinear statistical physics. Starting from the classical formulations we
derive the corresponding differential models and show that under a suitable
mean field scaling they correspond to generalized Fokker-Planck equations for
the mutants distribution whose solutions are given by the corresponding
Luria-Delbr\"uck distribution. Numerical results confirming the theoretical
analysis are also presented.


  We study the accuracy of Bayesian supervised method used to cluster
individuals into genetically homogeneous groups on the basis of dominant or
codominant molecular markers. We provide a formula relating an error criterion
the number of loci used and the number of clusters. This formula is exact and
holds for arbitrary number of clusters and markers. Our work suggests that
dominant markers studies can achieve an accuracy similar to that of codominant
markers studies if the number of markers used in the former is about 1.7 times
larger than in the latter.


  Background. Organisms use a variety of mechanisms to protect themselves
against perturbations. For example, repair mechanisms fix damage, feedback
loops keep homeostatic systems at their setpoints, and biochemical filters
distinguish signal from noise. Such buffering mechanisms are often discussed in
terms of robustness, which may be measured by reduced sensitivity of
performance to perturbations. Methodology/Principal Findings. I use a
mathematical model to analyze the evolutionary dynamics of robustness in order
to understand aspects of organismal design by natural selection. I focus on two
characters: one character performs an adaptive task; the other character
buffers the performance of the first character against perturbations. Increased
perturbations favor enhanced buffering and robustness, which in turn decreases
sensitivity and reduces the intensity of natural selection on the adaptive
character. Reduced selective pressure on the adaptive character often leads to
a less costly, lower performance trait. Conclusions/Significance. The paradox
of robustness arises from evolutionary dynamics: enhanced robustness causes an
evolutionary reduction in the adaptive performance of the target character,
leading to a degree of maladaptation compared to what could be achieved by
natural selection in the absence of robustness mechanisms. Over evolutionary
time, buffering traits may become layered on top of each other, while the
underlying adaptive traits become replaced by cheaper, lower performance
components. The paradox of robustness has widespread implications for
understanding organismal design.


  Electrical signaling via voltage-gated ion channels depends upon the function
of the voltage sensor (VS), identified with the S1-S4 domain of voltage-gated K
channels. Here we investigate some physical aspects of the sliding-helix model
of the VS using simulations based on VS charges, linear dielectrics and
whole-body motion. Model electrostatics in voltage-clamped boundary conditions
are solved using a boundary element method. The statistical mechanical
consequences of the electrostatic configurational energy are computed to gain
insight into the sliding-helix mechanism and to predict experimentally measured
ensemble properties such as gating charge displaced by an applied voltage.
Those consequences and ensemble properties are investigated for variations of:
S4 configuration ({\alpha}- and 3(10)-helical), intrinsic counter-charges,
protein polarizability, geometry of the gating canal, screening of S4 charges
by the baths, and protein charges located at the bath interfaces. We find that
the sliding helix VS has an inherent electrostatic stability and its function
as a VS is robust in the parameter space explored. Maximal charge displacement
is limited by geometry, specifically the range of movement where S4 charges and
counter-charges overlap in the region of weak dielectric. The steepness of
charge rearrangement in the physiological voltage range is sensitive to the
landscape of electrostatic energy: energy differences of <2 kT have substantial
consequences. Such variations of energy landscape are produced by all
variations of model features tested. The amount of free energy per unit voltage
that a sliding-helix VS can deliver to other parts of the channel (conductance
voltage sensitivity) is limited by both the maximal displaced charge and the
steepness of charge redistribution by voltage (sensor voltage sensitivity).


  We study the evolution of cooperation in public goods games on the square
lattice, focusing on the effects that are brought about by different sizes of
groups where individuals collect their payoffs and search for potential
strategy donors. We find that increasing the group size does not necessarily
lead to mean-field behavior, as is traditionally observed for games governed by
pairwise interactions, but rather that public cooperation may be additionally
promoted by means of enhanced spatial reciprocity that sets in for very large
groups. Our results highlight that the promotion of cooperation due to spatial
interactions is not rooted solely in having restricted connections amongst
players, but also in individuals having the opportunity to collect payoffs
separately from their direct opponents. Moreover, in large groups the presence
of a small number of defectors is bearable, which makes the mixed phase region
expand with increasing group size. Having a chance of exploiting distant
players, however, offers defectors a new way to break the phalanx of
cooperators, and even to resurrect from small numbers to eventually completely
invade the population.


  In biology, the evolution of increasingly cooperative groups has shaped the
history of life. Genes collaborate in the control of cells; cells efficiently
divide tasks to produce cohesive multicellular individuals; individual members
of insect colonies cooperate in integrated societies. Biological cooperation
provides a foundation on which to understand human behavior. Conceptually, the
economics of efficient allocation and the game-like processes of strategy are
well understood in biology; we find the same essential processes in many
successful theories of human sociality. Historically, the trace of biological
evolution informs in two ways. First, the evolutionary transformations in
biological cooperation provide insight into how economic and strategic
processes play out over time--a source of analogy that, when applied
thoughtfully, aids analysis of human sociality. Second, humans arose from
biological history--a factual account of the past that tells us much about the
material basis of human behavior.