  Using Wang-Landau sampling with suitable Monte Carlo trial moves (pull moves
and bond-rebridging moves combined) we have determined the density of states
and thermodynamic properties for a short sequence of the HP protein model. For
free chains these proteins are known to first undergo a collapse "transition"
to a globule state followed by a second "transition" into a native state. When
placed in the proximity of an attractive surface, there is a competition
between surface adsorption and folding that leads to an intriguing sequence of
"transitions". These transitions depend upon the relative interaction strengths
and are largely inaccessible to "standard" Monte Carlo methods.


  A classical result, fundamental to evolutionary biology, states that an
edge-weighted tree $T$ with leaf set $X$, positive edge weights, and no
vertices of degree 2 can be uniquely reconstructed from the set of leaf-to-leaf
distances between any two elements of $X$. In biology, $X$ corresponds to a set
of taxa (e.g. extant species), the tree $T$ describes their phylogenetic
relationships, the edges correspond to earlier species evolving for a time
until splitting in two or more species by some speciation/bifurcation event,
and their length corresponds to the genetic change accumulating over that time
in such a species. In this paper, we investigate which subsets of
$\binom{X}{2}$ suffice to determine (`lasso') a tree from the leaf-to-leaf
distances induced by that tree. The question is particularly topical since
reliable estimates of genetic distance - even (if not in particular) by modern
mass-sequencing methods - are, in general, available only for certain
combinations of taxa.


  The draft sequence of the human genome became available almost a decade ago
but the encoded proteome is not being explored to its fullest. Our bibliometric
analysis of several large protein families, including those known to be
"druggable", reveals that, even today, most papers focus on proteins that were
known prior to 2000. It is evident that one or more aspects of the biomedical
research system severely limits the exploration of the proteins in the 'dark
matter' of the proteome, despite unbiased genetic approaches that have pointed
to their functional relevance. It is perhaps not surprising that relatively few
genome-derived targets have led to approved drugs.


  We present experimental results regarding the electrotransfer of plasmid DNA
into phosphatidylcholine giant unilamellar vesicles (GUVs). Our observations
indicate that a direct entry is the predominant mechanism of electrotransfer. A
quantitative analysis of the DNA concentration increments inside the GUVs is
also performed, and we find that our experimental data are very well described
by a simple theoretical model in which DNA entry is mostly driven by
electrophoresis. Our theoretical framework allows for the prediction of the
amount of transfered DNA as a function of the electric field parameters, and
thus paves the way towards a novel method for encapsulating with high
efficiency not only DNA, but any negatively charged macromolecule into GUVs.


  Electroporation is a physical method to induce the uptake of therapeutic
drugs and DNA, by eukaryotic cells and tissues. The phenomena behind
electro-mediated membrane permeabilization to plasmid DNA have been shown to be
significantly more complex than those for small molecules. Small molecules
cross the permeabilized membrane by diffusion whereas plasmid DNA first
interacts with the electropermeabilized part of the cell surface, forming
localized aggregates. The dynamics of this process is still poorly understood
because direct observations have been limited to scales of the order of
seconds. Here, cells are electropermeabilized in the presence of plasmid DNA
and monitored with a temporal resolution of 2 ms. This allows us to show that
during the first pulse application, plasmid complexes, or aggregates, start to
form at distinct sites on the cell membrane. FRAP measurements show that the
positions of these sites are remarkably immobile during the application of
further pluses. A theoretical model is proposed to explain the appearance of
distinct interaction sites, the quantitative increase in DNA and also their
immobility leading to a tentative explanation for the success of
electro-mediated gene delivery.


  The efficiency of institutionalized punishment is studied by evaluating the
stationary states in the spatial public goods game comprising unconditional
defectors, cooperators, and cooperating pool-punishers as the three competing
strategies. Fine and cost of pool-punishment are considered as the two main
parameters determining the stationary distributions of strategies on the square
lattice. Each player collects its payoff from five five-person public goods
games, and the evolution of strategies is subsequently governed by imitation
based on pairwise comparisons at a low level of noise. The impact of
pool-punishment on the evolution of cooperation in structured populations is
significantly different from that reported previously for peer-punishment.
Representative phase diagrams reveal remarkably rich behavior, depending also
on the value of the synergy factor that characterizes the efficiency of
investments payed into the common pool. Besides traditional single and
two-strategy stationary states, a rock-paper-scissors type cyclic dominance can
emerge in strikingly different ways.


  The structural origin of the size of the 11 nm nucleosomal disc is addressed.
On the nanometer length-scale the organization of DNA as chromatin in the
chromosomes involves a coiling of DNA around the histone core of the
nucleosome. We suggest that the size of the nucleosome core particle is
dictated by the fulfillment of two criteria: One is optimizing the volume
fraction of the DNA double helix; this requirement for close-packing has its
root in optimizing atomic and molecular interactions. The other criterion being
that of having a zero strain-twist coupling; being a zero-twist structure is a
necessity when allowing for transient tensile stresses during the
reorganization of DNA, e.g., during the reposition, or sliding, of a nucleosome
along the DNA double helix. The mathematical model we apply is based on a
tubular description of double helices assuming hard walls. When the base-pairs
of the linker-DNA is included the estimate of the size of an ideal nucleosome
is in close agreement with the experimental numbers. Interestingly, the size of
the nucleosome is shown to be a consequence of intrinsic properties of the DNA
double helix.


  Evolutionary dynamics in finite populations is known to fixate eventually in
the absence of mutation. We here show that a similar phenomenon can be found in
stochastic game dynamical batch learning, and investigate fixation in learning
processes in a simple 2x2 game, for two-player games with cyclic interaction,
and in the context of the best-shot network game. The analogues of finite
populations in evolution are here finite batches of observations between
strategy updates. We study when and how such fixation can occur, and present
results on the average time-to-fixation from numerical simulations. Simple
cases are also amenable to analytical approaches and we provide estimates of
the behaviour of so-called escape times as a function of the batch size. The
differences and similarities with escape and fixation in evolutionary dynamics
are discussed.


  In this paper, we consider approximability issues of the following four
problems: triangle packing, full sibling reconstruction, maximum profit
coverage and 2-coverage. All of them are generalized or specialized versions of
set-cover and have applications in biology ranging from full-sibling
reconstructions in wild populations to biomolecular clusterings; however, as
this paper shows, their approximability properties differ considerably. Our
inapproximability constant for the triangle packing problem improves upon the
previous results; this is done by directly transforming the inapproximability
gap of Haastad for the problem of maximizing the number of satisfied equations
for a set of equations over GF(2) and is interesting in its own right. Our
approximability results on the full siblings reconstruction problems answers
questions originally posed by Berger-Wolf et al. and our results on the maximum
profit coverage problem provides almost matching upper and lower bounds on the
approximation ratio, answering a question posed by Hassin and Or.


  Mass-action chemical reaction systems are frequently used in Computational
Biology. The corresponding polynomial dynamical systems are often large
(consisting of tens or even hundreds of ordinary differential equations) and
poorly parametrized (due to noisy measurement data and a small number of data
points and repetitions). Therefore, it is often difficult to establish the
existence of (positive) steady states or to determine whether more complicated
phenomena such as multistationarity exist. If, however, the steady state ideal
of the system is a binomial ideal, then we show that these questions can be
answered easily. The focus of this work is on systems with this property, and
we say that such systems have toric steady states. Our main result gives
sufficient conditions for a chemical reaction system to have toric steady
states. Furthermore, we analyze the capacity of such a system to exhibit
positive steady states and multistationarity. Examples of systems with toric
steady states include weakly-reversible zero-deficiency chemical reaction
systems. An important application of our work concerns the networks that
describe the multisite phosphorylation of a protein by a kinase/phosphatase
pair in a sequential and distributive mechanism.


  We present a comprehensive new framework for handling biologically accurate
models of molecular evolution. This model provides a systematic framework for
studying models of molecular evolution that implement heterogeneous rates,
conservation of reading frame, differing rates of insertion and deletion,
customizable parametrization of the probabilities and types of substitutions,
insertions, and deletions, as well as neighboring dependencies. We have stated
the model in terms of an infinite state Markov chain in order to maximize the
number of applicable theorems useful in the analysis of the model. We use such
theorems to develop an alignment-free parameter estimation technique. This
alignment-free technique circumvents many of the nuanced issues related to
alignment-dependent estimation. We then apply an implementation of our model to
reproduce (in a completely alignment-free fashion) some observed results of
Zhang and Gerstein (2003) regarding indel length distribution in human
ribosomal protein pseudogenes.


  Evolutionary game theory has been an important tool for describing economic
and social behaviour for decades. Approximate mean value equations describing
the time evolution of strategy concentrations can be derived from the players'
microscopic update rules. We show that they can be generalized to a learning
process. As an example, we compare a restricted imitation process, in which
unused parts of the role model's meta-strategy are hidden from the imitator,
with the widely used imitation rule that allows the imitator to adopt the
entire meta-strategy of the role model. This change in imitation behaviour
greatly affects dynamics and stationary states in the iterated prisoner
dilemma. Particularly we find Grim Trigger to be a more successful strategy
than Tit-For-Tat especially in the weak selection regime.


  We investigate stall force and polymerization kinetics of rigid
protofilaments in a microtubule or interacting filaments in bundles under an
external load force in the framework of a discrete growth model. We introduce
the concecpt of polymerization cycles to describe the stochastic growth
kinetics, which allows us to derive an exact expression for the stall force. We
find that the stall force is independent of ensemble geometry and load
distribution. Furthermore, the stall force is proportional to the number of
filaments and increases linearly with the strength of lateral filament
interactions. These results are corroborated by simulations, which also show a
strong influence of ensemble geometry on growth kinetics below the stall force.


  The polymerization of actin filaments is coupled to the hydrolysis of
adenosine triphosphate (ATP), which involves both the cleavage of ATP and the
release of inorganic phosphate. We describe hydrolysis by a reduced two-state
model with a cooperative cleavage mechanism, where the cleavage rate depends on
the state of the neighboring actin protomer in a filament. We obtain
theoretical predictions of experimentally accessible steady state quantities
such as the size of the ATP-actin cap, the size distribution of ATP-actin
islands, and the cleavage flux for cooperative cleavage mechanisms.


  Traditionally, resource limitation in evolutionary game theory is assumed
just to impose a constant population size. Here we show that resource
limitations may generate dynamical payoffs able to alter an original prisoner's
dilemma, and to allow for the stable coexistence between unconditional
cooperators and defectors in well-mixed populations. This is a consequence of a
self-organizing process that turns the interaction payoff matrix into
evolutionary neutral, and represents a resource-based control mechanism
preventing the spread of defectors. To our knowledge, this is the first example
of coexistence in well-mixed populations with a game structure different from a
snowdrift game.


  FoF1-ATP synthase is the ubiquitous membrane-bound enzyme in mitochondria,
chloroplasts and bacteria which provides the 'chemical energy currency'
adenosine triphosphate (ATP) for cellular processes. In Escherichia coli ATP
synthesis is driven by a proton motive force (PMF) comprising a proton
concentration difference {\Delta}pH plus an electric potential {\Delta}{\Psi}
across the lipid membrane. Single-molecule in vitro experiments have confirmed
that proton-driven subunit rotation within FoF1-ATP synthase is associated with
ATP synthesis. Based on intramolecular distance measurements by single-molecule
fluorescence resonance energy transfer (FRET) the kinetics of subunit rotation
and the step sizes of the different rotor parts have been unraveled. However,
these experiments were accomplished in the presence of a PMF consisting of a
maximum {\Delta}pH ~ 4 and an unknown {\Delta}{\Psi}. In contrast, in living
bacteria the maximum {\Delta}pH across the plasma membrane is likely 0.75, and
{\Delta}{\Psi} has been measured between -80 and -140 mV. Thus the problem of
in vivo catalytic turnover rates, or the in vivo rotational speed in single
FoF1-ATP synthases, respectively, has to be solved. In addition, the absolute
number of functional enzymes in a single bacterium required to maintain the
high ATP levels has to be determined. We report our progress of measuring
subunit rotation in single FoF1-ATP synthases in vitro and in vivo, which was
enabled by a new labeling approach for single-molecule FRET measurements.


  We formulate general plant-herbivore interaction models with monotone plant
growth functions (rates). We study the impact of monotone plant growth
functions in general plant-herbivore models on their dynamics. Our study shows
that all monotone plant growth models generate a unique interior equilibrium
and they are uniform persistent under certain range of {parameters} values.
However, if the attacking rate of herbivore is too small or the quantity of
plant is not enough, then herbivore goes extinct. Moreover, these models lead
to noise sensitive bursting which can be identified as a dynamical mechanism
for almost periodic outbreaks of the herbivore infestation. Monotone and
non-monotone plant growth models are contrasted with respect to bistability and
crises of chaotic attractors.


  In order to understand the role of space in ecological communities where each
species produces a certain type of resource and has varying abilities to
exploit the resources produced by its own species and by the other species, we
carry out a comparative study of an interacting particle system and its
mean-field approximation. For a wide range of parameter values, we show both
analytically and numerically that the spatial model results in predictions that
significantly differ from its nonspatial counterpart, indicating that the use
of the mean-field approach to describe the evolution of communities in which
individuals only interact locally is invalid. In two-species communities, the
disagreements between the models appear either when both species compete by
producing resources that are more beneficial for their own species or when both
species cooperate by producing resources that are more beneficial for the other
species. In particular, while both species coexist if and only if they
cooperate in the mean-field approximation, the inclusion of space in the form
of local interactions may prevent coexistence even in cooperative communities.
Introducing additional species, cooperation is no longer the only mechanism
that promotes coexistence. We prove that, in three-species communities,
coexistence results either from a global cooperative behavior, or from
rock-paper-scissors type interactions, or from a mixture of these dynamics,
which excludes in particular all cases in which two species compete.


  We propose a simple mathematical model by applying Michaelis-Menton equations
of enzyme kinetics to study the mutualistic interaction between the leaf cutter
ant and its fungus garden at the early stage of colony expansion. We derive the
sufficient conditions on the extinction and coexistence of these two species.
In addition, we give a region of initial condition that leads to the extinction
of two species when the model has an interior attractor. Our global analysis
indicates that the division of labor by workers ants and initial conditions are
two important {factors} that determine whether leaf cutter ants colonies and
their fungus garden survive and grow can exist or not. We validate the model by
doing the comparing between model simulations and data on fungal and ant colony
growth rates under laboratory conditions. We perform sensitive analysis and
parameter estimation of the model based on the experimental data to gain more
biological insights on the ecological interactions between leaf cutter ants and
their fungus garden. Finally, we give conclusions and {discuss} potential
future {work}.


  Large scale databases are available that contain homologous gene families
constructed from hundreds of complete genome sequences from across the three
domains of Life. Here we discuss approches of increasing complexity aimed at
extracting information on the pattern and process of gene family evolution from
such datasets. In particular, we consider models that invoke processes of gene
birth (duplication and transfer) and death (loss) to explain the evolution of
gene families. First, we review birth-and-death models of family size evolution
and their implications in light of the universal features of family size
distribution observed across different species and the three domains of life.
Subsequently, we proceed to recent developments on models capable of more
completely considering information in the sequences of homologous gene families
through the probabilistic reconciliation of the phylogenetic histories of
individual genes with the phylogenetic history of the genomes in which they
have resided. To illustrate the methods and results presented, we use data from
the HOGENOM database, demonstrating that the distribution of homologous gene
family sizes in the genomes of the Eukaryota, Archaea and Bacteria exhibit
remarkably similar shapes. We shown that these distributions are best described
by models of gene family size evolution where for individual genes the death
(loss) rate is larger than the birth (duplication and transfer) rate, but new
families are continually supplied to the genome by a process of origination.
Finally, we use probabilistic reconciliation methods to take into consideration
additional information from gene phylogenies, and find that, for prokaryotes,
the majority of birth events are the result of transfer.


  Fifteen years ago Monique Tirion showed that the low-frequency normal modes
of a protein are not significantly altered when non-bonded interactions are
replaced by Hookean springs, for all atom pairs whose distance is smaller than
a given cutoff value. Since then, it has been shown that coarse-grained
versions of Tirion's model are able to provide fair insights on many dynamical
properties of biological macromolecules. In this text, theoretical tools
required for studying these so-called Elastic Network Models are described,
focusing on practical issues and, in particular, on possible artifacts. Then,
an overview of some typical results that have been obtained by studying such
models is given.


  Fluorescent detection of proteins is a popular method of detection allying
sensitivity, linearity and compatibility with mass spectrometry. Among the
numerous methods described in the literature, staining with ruthenium II
tris(bathophenanthroline disulfonate) is particularly cost-effective, but
slightly cumbersome owing to difficulties in the preparation of the complex and
complexity of staining protocols. We describe here the modifications on both
aspects that allow to perform a higher contrast staining and offer a more
robust method of complex preparation, thereby maximizing the advantages of the
method.


  Patterns of biodiversity predicted by the neutral theory rely on a simple
phenomenological model of speciation. To further investigate the effect of
speciation on neutral biodiversity, we analyze a spatially-explicit neutral
model based on population genetics. We define the metacommunity as a system of
populations exchanging migrants and we use this framework to introduce
speciation with little or no gene flow (allopatric and parapatric speciation).
We find that with realistic mutation rates, our metacommunity model driven by
neutral processes cannot support more than a few species. Adding natural
selection in the population genetics of speciation increases the number of
species in the metacommunity but the level of diversity found in Barro Colorado
Island is difficult to reach.


  Native oyster populations in Chesapeake Bay have been the focus of three
decades of restoration attempts, which have generally failed to rebuild the
populations and oyster reef structure. Recent restoration successes and field
experiments suggest that high-relief reefs offset heavy sedimentation and
promote oyster survival, disease resistance and growth, in contrast to
low-relief reefs which degrade in just a few years. These findings suggest the
existence of alternative stable states in oyster reef populations. We developed
a mathematical model consisting of three differential equations that represent
volumes of live oysters, dead oyster shells (= accreting reef), and sediment.
Bifurcation analysis and numerical simulations demonstrated that multiple
nonnegative equilibria can exist for live oyster, accreting reef and sediment
volume at an ecologically reasonable range of parameter values; the initial
height of oyster reefs determined which equilibrium was reached. This
investigation thus provides a conceptual framework for alternative stable
states in native oyster populations, and can be used as a tool to improve the
likelihood of success in restoration efforts.


  We report the results of molecular dynamics simulations of an off-lattice
protein model featuring a physical force-field and amino-acid sequence. We show
that localized modes of nonlinear origin (discrete breathers) emerge naturally
as continuations of a subset of high-frequency normal modes residing at
specific sites dictated by the native fold. In the case of the small
$\beta$-barrel structure that we consider, localization occurs on the turns
connecting the strands. At high energies, discrete breathers stabilize the
structure by concentrating energy on few sites, while their collapse marks the
onset of large-amplitude fluctuations of the protein. Furthermore, we show how
breathers develop as energy-accumulating centres following perturbations even
at distant locations, thus mediating efficient and irreversible energy
transfers. Remarkably, due to the presence of angular potentials, the breather
induces a local static distortion of the native fold. Altogether, the
combination of this two nonlinear effects may provide a ready means for
remotely controlling local conformational changes in proteins.


  Cosmic radiation backgrounds are a constraint on life, and their distribution
will affect the Galactic Habitable Zone. Life on Earth has developed in the
context of these backgrounds, and characterizing event rates will elaborate the
important influences. This in turn can be a base for comparison with other
potential life-bearing planets. In this review we estimate the intensities and
rates of occurrence of many kinds of strong radiation bursts by astrophysical
entities ranging from gamma-ray bursts at cosmological distances to the Sun
itself. Many of these present potential hazards to the biosphere: on timescales
long compared with human history, the probability of an event intense enough to
disrupt life on the land surface or in the oceans becomes large. We enumerate
the known sources of radiation and characterize their intensities at the Earth
and rates or upper limits on these quantities. When possible, we estimate a
"lethal interval", our best estimate of how often a major extinction-level
event is probable given the current state of knowledge; we base these estimates
on computed or expected depletion of stratospheric ozone. In general, moderate
level events are dominated by the Sun, but the far more severe infrequent
events are probably dominated by gamma-ray bursts and supernovae. We note for
the first time that so-called "short-hard" gamma-ray bursts are a substantial
threat, comparable in magnitude to supernovae and greater than that of the
higher-luminosity long bursts considered in most past work. Given their
precursors, short bursts may come with little or no warning.


  The simplest, and most common, stochastic model for population processes,
including those from biochemistry and cell biology, are continuous time Markov
chains. Simulation of such models is often relatively straightforward as there
are easily implementable methods for the generation of exact sample paths.
However, when using ensemble averages to approximate expected values, the
computational complexity can become prohibitive as the number of computations
per path scales linearly with the number of jumps of the process. When such
methods become computationally intractable, approximate methods, which
introduce a bias, can become advantageous. In this paper, we provide a general
framework for understanding the weak error, or bias, induced by different
numerical approximation techniques in the current setting. The analysis takes
into account both the natural scalings within a given system and the step-size
of the numerical method. Examples are provided to demonstrate the main
analytical results as well as the reduction in computational complexity
achieved by the approximate methods.


  Cells use genetic switches to shift between alternate gene expression states,
e.g., to adapt to new environments or to follow a developmental pathway. Here,
we study the dynamics of switching in a generic-feedback on/off switch. Unlike
protein-only models, we explicitly account for stochastic fluctuations of mRNA,
which have a dramatic impact on switch dynamics. Employing the WKB theory to
treat the underlying chemical master equations, we obtain accurate results for
the quasi-stationary distributions of mRNA and protein copy numbers and for the
mean switching time, starting from either state. Our analytical results agree
well with Monte Carlo simulations. Importantly, one can use the approach to
study the effect of varying biological parameters on switch stability.


  A dynamical model is presented for chiral change in DNA molecules. The model
is an extension of the conventional elastic model which incorporates the
structure of base pairs and uses a spinor representation for the DNA
configuration together with a gauge principle. Motivated by a recent experiment
reporting chiral transitions between right-handed B-DNA and left-handed Z-DNA
[M. Lee, et. al., Proc. Natl. Acad. Sci. (USA) 107, 4985 (2010)], we analyze
the free energy for the particular case of linear DNA with an externally
applied torque. The model shows that there exists, at low temperature, a rapid
structural change depending on the torque exerted on the DNA, which causes
switching in B and Z domain sizes. This can explain the frequent switches of
DNA extension observed in experiments.


  We report on electrical properties of the two sensing proteins:
bacteriorhodopsin and rat olfactory receptor OR-I7. As relevant transport
parameters we consider the small-signal impedance spectrum and the static
current-voltage characteristics. Calculations are compared with available
experimental results and the model predictability is tested for future
perspectives.


  A recent survey of 17 134 proteins has identified a new class of proteins
which are expected to yield stretching induced force-peaks in the range of 1
nN. Such high force peaks should be due to forcing of a slip-loop through a
cystine ring, i.e. by generating a cystine slipknot. The survey has been
performed in a simple coarse grained model. Here, we perform all-atom steered
molecular dynamics simulations on 15 cystine knot proteins and determine their
resistance to stretching. In agreement with previous studies within a coarse
grained structure based model, the level of resistance is found to be
substantially higher than in proteins in which the mechanical clamp operates
through shear. The large stretching forces arise through formation of the
cystine slipknot mechanical clamp and the resulting steric jamming. We
elucidate the workings of such a clamp in an atomic detail. We also study the
behavior of five top strength proteins with the shear-based mechanostability in
which no jamming is involved. We show that in the atomic model, the jamming
state is relieved by moving one amino acid at a time and there is a choice in
the selection of the amino acid that advances the first. In contrast, the
coarse grained model also allows for a simultaneous passage of two amino acids.


  Evolutionary dynamics have been traditionally studied in the context of
homogeneous populations, mainly described my the Moran process. Recently, this
approach has been generalized in \cite{LHN} by arranging individuals on the
nodes of a network. Undirected networks seem to have a smoother behavior than
directed ones, and thus it is more challenging to find suppressors/amplifiers
of selection. In this paper we present the first class of undirected graphs
which act as suppressors of selection, by achieving a fixation probability that
is at most one half of that of the complete graph, as the number of vertices
increases. Moreover, we provide some generic upper and lower bounds for the
fixation probability of general undirected graphs. As our main contribution, we
introduce the natural alternative of the model proposed in \cite{LHN}, where
all individuals interact simultaneously and the result is a compromise between
aggressive and non-aggressive individuals. That is, the behavior of the
individuals in our new model and in the model of \cite{LHN} can be interpreted
as an "aggregation" vs. an "all-or-nothing" strategy, respectively. We prove
that our new model of mutual influences admits a \emph{potential function},
which guarantees the convergence of the system for any graph topology and any
initial fitness vector of the individuals. Furthermore, we prove fast
convergence to the stable state for the case of the complete graph, as well as
we provide almost tight bounds on the limit fitness of the individuals. Apart
from being important on its own, this new evolutionary model appears to be
useful also in the abstract modeling of control mechanisms over invading
populations in networks. We demonstrate this by introducing and analyzing two
alternative control approaches, for which we bound the time needed to stabilize
to the "healthy" state of the system.


  We propose a novel stochastic method to generate paths conditioned to start
in an initial state and end in a given final state during a certain time
$t_{f}$. These paths are weighted with a probability given by the overdamped
Langevin dynamics. We show that these paths can be exactly generated by a
non-local stochastic differential equation. In the limit of short times, we
show that this complicated non-solvable equation can be simplified into an
approximate stochastic differential equation. For longer times, the paths
generated by this approximate equation can be reweighted to generate the
correct statistics. In all cases, the paths generated by this equation are
statistically independent and provide a representative sample of transition
paths. In case the reaction takes place in a solvent (e.g. protein folding in
water), the explicit solvent can be treated. The method is illustrated on the
one-dimensional quartic oscillator.


  Two giants of evolutionary theory, Sewall Wright and R. A. Fisher, fought
bitterly for over thirty years. The Wright-Fisher controversy forms a
cornerstone of the history and philosophy of biology. I argue that the standard
interpretations of the Wright-Fisher controversy do not accurately represent
the ideas and arguments of these two key historical figures. The usual account
contrasts the major slogans attached to each name: Wright's adaptive landscape
and shifting balance theory of evolution versus Fisher's fundamental theorem of
natural selection. These alternative theories are in fact incommensurable.
Wright's theory is a detailed dynamical model of evolutionary change in actual
populations. Fisher's theory is an abstract invariance and conservation law
that, like all physical laws, captures essential features of a system but does
not account for all aspects of dynamics in real examples. This key contrast
between embodied theories of real cases and abstract laws is missing from prior
analyses of Wright versus Fisher. They never argued about this contrast.
Instead, the issue at stake in their arguments concerned the actual dynamics of
real populations. Both agreed that fluctuations of nonadditive (epistatic) gene
combinations play a central role in evolution. Wright emphasized stochastic
fluctuations of gene combinations in small, isolated populations. By contrast,
Fisher believed that fluctuating selection in large populations was the main
cause of fluctuation in nonadditive gene combinations. Close reading shows that
widely cited views attributed to Fisher mostly come from what Wright said about
Fisher, whereas Fisher's own writings clearly do not support such views.


  Whole-genome sequencing in an isolated population with few founders directly
ascertains variants from the population bottleneck that may be rare elsewhere.
In such populations, shared haplotypes allow imputation of variants in
unsequenced samples without resorting to statistical methods, as in studies of
outbred cohorts. We focus on an isolated population cohort from the Pacific
Island of Kosrae, Micronesia, where we previously collected SNP array and rich
phenotype data for the majority of the population. We report identification of
long regions with haplotypes co-inherited between pairs of individuals and
methodology to leverage such shared genetic content for imputation. Our
estimates show that sequencing as few as 40 personal genomes allows for
imputation in up to 60% of the 3,000-person cohort at the average locus. We
ascertained a pilot data-set of whole-genome sequences from seven Kosraean
individuals, with average 5X coverage. This dataset identified 5,735,306 unique
sites of which 1,212,831 were previously unknown. Additionally, these Kosraen
variants are unusually enriched for alleles that are rare in other populations
when compared to geographic neighbors. We were able to use the presence of
shared haplotypes between the seven individuals to estimate imputation accuracy
of known and novel variants and achieved levels of 99.6% and 97.3%,
respectively. This study presents the first whole-genome analysis of a
homogenous isolate population with emphasis on rare variant inference.


  Understanding design principles of molecular interaction networks is an
important goal of molecular systems biology. Some insights have been gained
into features of their network topology through the discovery of graph
theoretic patterns that constrain network dynamics. This paper contributes to
the identification of patterns in the mechanisms that govern network dynamics.
The control of nodes in gene regulatory, signaling, and metabolic networks is
governed by a variety of biochemical mechanisms, with inputs from other network
nodes that act additively or synergistically. This paper focuses on a certain
type of logical rule that appears frequently as a regulatory pattern. Within
the context of the multistate discrete model paradigm, a rule type is
introduced that reduces to the concept of nested canalyzing function in the
Boolean network case. It is shown that networks that employ this type of
multivalued logic exhibit more robust dynamics than random networks, with few
attractors and short limit cycles. It is also shown that the majority of
regulatory functions in many published models of gene regulatory and signaling
networks are nested canalyzing.


  A quantum theory on conformation-electron system is presented. Protein
folding is regarded as the quantum transition between torsion states on
polypeptide chain, and the folding rate is calculated by nonadiabatic operator
method. The theory is used to study the temperature dependences of folding rate
of 15 proteins and their non-Arrhenius behavior can all be deduced in a natural
way. A general formula on the rate-temperature dependence has been deduced
which is in good accordance with experimental data. These temperature
dependences are further analyzed in terms of torsion potential parameters. Our
results show it is necessary to move outside the realm of classical physics
when the temperature dependence of protein folding is studied quantitatively.


  In this work we explore the effects of human mobility on the dispersion of a
vector borne disease. We combine an already presented stochastic model for
dengue with a simple representation of the daily motion of humans on a
schematic city of 20x20 blocks with 100 inhabitants in each block. The pattern
of motion of the individuals is described in terms of complex networks in which
links connect different blocks and the link length distribution is in
accordance with recent findings on human mobility. It is shown that human
mobility can turn out to be the main driving force of the disease dispersal.


  A computational challenge to validate the candidate disease genes identified
in a high-throughput genomic study is to elucidate the associations between the
set of candidate genes and disease phenotypes. The conventional gene set
enrichment analysis often fails to reveal associations between disease
phenotypes and the gene sets with a short list of poorly annotated genes,
because the existing annotations of disease causative genes are incomplete. We
propose a network-based computational approach called rcNet to discover the
associations between gene sets and disease phenotypes. Assuming coherent
associations between the genes ranked by their relevance to the query gene set,
and the disease phenotypes ranked by their relevance to the hidden target
disease phenotypes of the query gene set, we formulate a learning framework
maximizing the rank coherence with respect to the known disease phenotype-gene
associations. An efficient algorithm coupling ridge regression with label
propagation, and two variants are introduced to find the optimal solution of
the framework. We evaluated the rcNet algorithms and existing baseline methods
with both leave-one-out cross-validation and a task of predicting recently
discovered disease-gene associations in OMIM. The experiments demonstrated that
the rcNet algorithms achieved the best overall rankings compared to the
baselines. To further validate the reproducibility of the performance, we
applied the algorithms to identify the target diseases of novel candidate
disease genes obtained from recent studies of GWAS, DNA copy number variation
analysis, and gene expression profiling. The algorithms ranked the target
disease of the candidate genes at the top of the rank list in many cases across
all the three case studies. The rcNet algorithms are available as a webtool for
disease and gene set association analysis at
http://compbio.cs.umn.edu/dgsa_rcNet.


  Regulatory genes called small RNAs (sRNAs) are known to play critical roles
in cellular responses to changing environments. For several sRNAs, regulation
is effected by coupled stoichiometric degradation with messenger RNAs (mRNAs).
The nonlinearity inherent in this regulatory scheme indicates that exact
analytical solutions for the corresponding stochastic models are intractable.
Here, we present a variational approach to analyze a well-studied stochastic
model for regulation by sRNAs via coupled degradation. The proposed approach is
efficient and provides accurate estimates of mean mRNA levels as well as higher
order terms. Results from the variational ansatz are in excellent agreement
with data from stochastic simulations for a wide range of parameters, including
regions of parameter space where mean-field approaches break down. The proposed
approach can be applied to quantitatively model stochastic gene expression in
complex regulatory networks.


  De novo discovery of over-represented DNA motifs is one of the major
challenges in computational biology. Although numerous tools have been
available for de novo motif discovery, many of these tools are subject to local
optima phenomena, which may hinder detection of multiple distinct motifs. A
greedy algorithm based tool named dMotifGreedy was developed. dMotifGreedy
begins by searching for candidate motifs from pair-wise local alignments of
input sequences and then computes an optimal global solution for each candidate
motif through a greedy algorithm. dMotifGreedy has competitive performance in
detecting a true motif and greatly enhanced performance in detecting multiple
distinct true motifs. dMotifGreedy is freely available via a stand-alone
program at http://lambchop.ads.uga.edu/dmotifgreedy/download.php.


  Certain cellular processes take the form of activity pulses that can be
interpreted in terms of noise-driven excitable dynamics. Here we present an
overview of different gene circuit architectures that exhibit excitable pulses
of protein expression, when subject to molecular noise. Different types of
excitable dynamics can occur depending on the bifurcation structure leading to
the specific excitable phase-space topology. The bifurcation structure is not,
however, linked to a particular circuit architecture. Thus a given gene circuit
design can sustain different classes of excitable dynamics depending on the
system parameters.


  Two species of Endangered Leaf frogs, Agalychnis moreletii and Agalychnis
annae, belonging to the tree frog Subfamily Phyllomedusinae, Genus Agalychnis,
were hybridized for the first time whilst being maintained in captivity.
Previous to this, these allopatric Central American species were considered as
being distinctly separate. Crossbreeding following genetic analysis reveals
that the two species are extremely closely related, and the hybrid of A.
moreletii and A. annae is presented for the first time. The importance of
identifying degrees of genetic variation between species and different
populations of the same species, for conservation purposes, is highlighted and
discussed.


  Many biological processes involve the action of molecular motors that
interact with the cell cytoskeleton. Some processes, such as the transport of
cargoes is achieved mainly by the action of individual motors. Other, such as
cell motility and division, require the cooperative work of many motors.
Collective motor dynamics can be quite complex and unexpected. One beautiful
example is the bidirectional ("back and forth") motion of filaments which is
induced when the motors within a group exert forces in opposite directions.
This review tackles the puzzle emerging from a recent experimental work in
which it has been shown that the characteristic reversal times of the
bidirectional motion are practically independent of the number of motors. This
result is in a striking contradiction with existing theoretical models that
predict an exponential growth of the reversal times with the size of the
system. We argue that the solution to this puzzle may be the crosstalk between
the motors which is mediated by the elastic tensile stress that develops in the
cytoskeleton track. The crosstalk does not directly correlate the attachment
and detachment of the motors, which work independently of each other. However,
it highly accelerates their detachments by making the detachment rates system
size dependent.


  Discriminative patterns are association patterns that occur with
disproportionate frequency in some classes versus others, and have been studied
under names such as emerging patterns and contrast sets. Such patterns have
demonstrated considerable value for classification and subgroup discovery, but
a detailed understanding of the types of interactions among items in a
discriminative pattern is lacking. To address this issue, we propose to
categorize discriminative patterns according to four types of item interaction:
(i) driver-passenger, (ii) coherent, (iii) independent additive and (iv)
synergistic beyond independent additive. Either of the last three is of
practical importance, with the latter two representing a gain in the
discriminative power of a pattern over its subsets. Synergistic patterns are
most restrictive, but perhaps the most interesting since they capture a
cooperative effect. For domains such as genetic research, differentiating among
these types of patterns is critical since each yields very different biological
interpretations. For general domains, the characterization provides a novel
view of the nature of the discriminative patterns in a dataset, which yields
insights beyond those provided by current approaches that focus mostly on
pattern-based classification and subgroup discovery. This paper presents a
comprehensive discussion that defines these four pattern types and investigates
their properties and their relationship to one another. In addition, these
ideas are explored for a variety of datasets (ten UCI datasets, one gene
expression dataset and two genetic-variation datasets). The results demonstrate
the existence, characteristics and statistical significance of the different
types of patterns. They also illustrate how pattern characterization can
provide novel insights into discriminative pattern mining and the
discriminative structure of different datasets.


  In this paper we explore relevant electrical properties of two olfactory
receptors (ORs), one from rat OR I7 and the other from human OR 17-40, which
are of interest for the realization of smell nanobiosensors. The investigation
compares existing experiments, coming from electrochemical impedance
spectroscopy, with the theoretical expectations obtained from an impedance
network protein analogue, recently developed. The changes in the response due
to the sensing action of the proteins are correlated with the conformational
change undergone by the single protein. The satisfactory agreement between
theory and experiments points to a promising development of a new class of
nanobiosensors based on the electrical properties of sensing proteins.


  The comparison of computer generated protein structural models is an
important element of protein structure prediction. It has many uses including
model quality evaluation, selection of the final models from a large set of
candidates or optimisation of parameters of energy functions used in
template-free modelling and refinement. Although many protein comparison
methods are available online on numerous web servers, they are not well suited
for large scale model comparison: (1) they operate with methods designed to
compare actual proteins, not the models of the same protein, (2) majority of
them offer only a single pairwise structural comparison and are unable to scale
up to a required order of thousands of comparisons. To bridge the gap between
the protein and model structure comparison we have developed the Protein Models
Comparator (pm-cmp). To be able to deliver the scalability on demand and handle
large comparison experiments the pm-cmp was implemented "in the cloud".
  Protein Models Comparator is a scalable web application for a fast
distributed comparison of protein models with RMSD, GDT TS, TM-score and
Q-score measures. It runs on the Google App Engine (GAE) cloud platform and is
a showcase of how the emerging PaaS (Platform as a Service) technology could be
used to simplify the development of scalable bioinformatics services. The
functionality of pm-cmp is accessible through API which allows a full
automation of the experiment submission and results retrieval. Protein Models
Comparator is free software released on the Affero GNU Public Licence and is
available with its source code at: http://www.infobiotics.org/pm-cmp
  This article presents a new web application addressing the need for a
large-scale model-specific protein structure comparison and provides an insight
into the GAE (Google App Engine) platform and its usefulness in scientific
computing.


  One of the classical questions in evolutionary biology is how evolutionary
processes are coupled at the gene and species level. With this motivation, we
compare the topological properties (mainly the depth scaling, as a
characterization of balance) of a large set of protein phylogenies with a set
of species phylogenies. The comparative analysis shows that both sets of
phylogenies share remarkably similar scaling behavior, suggesting the
universality of branching rules and of the evolutionary processes that drive
biological diversification from gene to species level. In order to explain such
generality, we propose a simple model which allows us to estimate the
proportion of evolvability/robustness needed to approximate the scaling
behavior observed in the phylogenies, highlighting the relevance of the
robustness of a biological system (species or protein) in the scaling
properties of the phylogenetic trees. Thus, the rules that govern the
incapability of a biological system to diversify are equally relevant both at
the gene and at the species level.


  In spatial evolutionary games the fitness of each individual is traditionally
determined by the payoffs it obtains upon playing the game with its neighbors.
Since defection yields the highest individual benefits, the outlook for
cooperators is gloomy. While network reciprocity promotes collaborative
efforts, chances of averting the impending social decline are slim if the
temptation to defect is strong. It is therefore of interest to identify viable
mechanisms that provide additional support for the evolution of cooperation.
Inspired by the fact that the environment may be just as important as
inheritance for individual development, we introduce a simple switch that
allows a player to either keep its original payoff or use the average payoff of
all its neighbors. Depending on which payoff is higher, the influence of either
option can be tuned by means of a single parameter. We show that, in general,
taking into account the environment promotes cooperation. Yet coveting the
fitness of one's neighbors too strongly is not optimal. In fact, cooperation
thrives best only if the influence of payoffs obtained in the traditional way
is equal to that of the average payoff of the neighborhood. We present results
for the prisoner's dilemma and the snowdrift game, for different levels of
uncertainty governing the strategy adoption process, and for different
neighborhood sizes. Our approach outlines a viable route to increased levels of
cooperative behavior in structured populations, but one that requires a
thoughtful implementation.


  In this paper we present geometric features of group based models. We focus
on the 3-Kimura model. We present a precise geometric description of the
variety associated to any tree on a Zariski open set. In particular this set
contains all biologically meaningful points. Our motivation is a conjecture of
Sturmfels and Sullivant on the degree in which the ideal associated to 3-Kimura
model is generated.


  Resource competition is a fundamental interaction in natural
communities.However little is known about competition in spatial environments
where organisms are able to regulate resource distributions. Here, we analyze
the competition of two consumers for two resources in a one-dimensional habitat
in which the resources are supplied from opposite sides. We show that the
success of an invading species crucially depends on the slope of the resource
gradients shaped by the resident. Our analysis reveals that parameter
combinations which lead to coexistence in a uniform environment may favor
alternative stable states in a spatial system, and vice versa. Furthermore,
differences in growth rate, mortality or dispersal abilities allow a consumer
to coexist stationarily with - or even outcompete - a competitor with lower
resource requirements. Applying our theory to a phytoplankton model, we explain
shifts in the community structure that are induced by environmental changes.


  The understanding of molecular cell biology requires insight into the
structure and dynamics of networks that are made up of thousands of interacting
molecules of DNA, RNA, proteins, metabolites, and other components. One of the
central goals of systems biology is the unraveling of the as yet poorly
characterized complex web of interactions among these components. This work is
made harder by the fact that new species and interactions are continuously
discovered in experimental work, necessitating the development of adaptive and
fast algorithms for network construction and updating. Thus, the
"reverse-engineering" of networks from data has emerged as one of the central
concern of systems biology research.
  A variety of reverse-engineering methods have been developed, based on tools
from statistics, machine learning, and other mathematical domains. In order to
effectively use these methods, it is essential to develop an understanding of
the fundamental characteristics of these algorithms. With that in mind, this
chapter is dedicated to the reverse-engineering of biological systems.
  Specifically, we focus our attention on a particular class of methods for
reverse-engineering, namely those that rely algorithmically upon the so-called
"hitting-set" problem, which is a classical combinatorial and computer science
problem, Each of these methods utilizes a different algorithm in order to
obtain an exact or an approximate solution of the hitting set problem. We will
explore the ultimate impact that the alternative algorithms have on the
inference of published in silico biological networks.


  Social, biological and economic networks grow and decline with occasional
fragmentation and re-formation, often explained in terms of external
perturbations. We show that these phenomena can be a direct consequence of
simple imitation and internal conflicts between 'cooperators' and 'defectors'.
We employ a game-theoretic model of dynamic network formation where successful
individuals are more likely to be imitated by newcomers who adopt their
strategies and copy their social network. We find that, despite using the same
mechanism, cooperators promote well-connected highly prosperous networks and
defectors cause the network to fragment and lose its prosperity; defectors are
unable to maintain the highly connected networks they invade. Once the network
is fragmented it can be reconstructed by a new invasion of cooperators, leading
to the cycle of formation and fragmentation seen, for example, in bacterial
communities and socio-economic networks. In this endless struggle between
cooperators and defectors we observe that cooperation leads to prosperity, but
prosperity is associated with instability. Cooperation is prosperous when the
network has frequent formation and fragmentation.


  We analyze the nonlinear waves that propagate on a microtubule that is
tethered at its minus end due to kinesin walking on it, as is seen during the
fluid mixing caused by cytoplasmic streaming in Drosophila oocytes.The model we
use assumes that the microtubule can be modeled as an elastic string in a
viscous medium. The effect of the kinesin is to apply a force tangential to the
microtubule and we also consider the addition of a uniform cytoplasmic velocity
field. We show that travelling wave solutions exist and analyze their
properties. There exist scale invariant families of solutions and solutions can
exist that are flat or helical. The relationship between the period and
wavelength is obtained by both analytic and numerical means. Numerical
implementation of the equation of motion verifies our analytical predictions.


  The actin cytoskeleton of adherent tissue cells often condenses into filament
bundles contracted by myosin motors, so-called stress fibers, which play a
crucial role in the mechanical interaction of cells with their environment.
Stress fibers are usually attached to their environment at the endpoints, but
possibly also along their whole length. We introduce a theoretical model for
such contractile filament bundles which combines passive viscoelasticity with
active contractility. The model equations are solved analytically for two
different types of boundary conditions. A free boundary corresponds to stress
fiber contraction dynamics after laser surgery and results in good agreement
with experimental data. Imposing cyclic varying boundary forces allows us to
calculate the complex modulus of a single stress fiber.


  Creatures as varied as mammals, fish, insects, reptiles, and migratory birds
have an intriguing `sixth' sense that allows them to distinguish north from
south by using the Earth's intrinsic magnetic field. Yet despite decades of
study, the physical basis of this magnetic sense remains elusive. A likely
mechanism is furnished by magnetically sensitive radical pair reactions
occurring in the retina, the light-sensitive part of the eyes. A photoreceptor,
cryptochrome, has been suggested to endow birds with magnetoreceptive abilities
as the protein has been shown to exhibit the biophysical properties required
for an animal magnetoreceptor to operate properly. Here, we propose a concrete
light-driven reaction cycle in cryptochrome that lets a magnetic field
influence the signaling state of the photoreceptor. The reaction cycle ties
together transient absorption and electron-spin-resonance observations with
known facts on avian magnetoreception. Our analysis establishes the feasibility
of cryptochrome to act as a geomagnetic compass in birds, gives insight into
structural and dynamic design features required for optimal detection of the
geomagnetic field direction, and suggests further theoretical and experimental
studies.


  In an implicit combinatorial optimization problem, the constraints are not
enumerated explicitly but rather stated implicitly through equations, other
constraints or auxiliary algorithms. An important subclass of such problems is
the implicit set cover (or, equivalently, hitting set) problem in which the
sets are not given explicitly but rather defined implicitly For example, the
well-known minimum feedback arc set problem is such a problem. In this paper,
we consider such a cover problem that arises in the study of wild populations
in biology in which the sets are defined implicitly via the Mendelian
constraints and prove approximability results for this problem.


  Understanding functional organization of genetic information is a major
challenge in modern biology. Following the initial publication of the human
genome sequence in 2001, advances in high-throughput measurement technologies
and efficient sharing of research material through community databases have
opened up new views to the study of living organisms and the structure of life.
In this thesis, novel computational strategies have been developed to
investigate a key functional layer of genetic information, the human
transcriptome, which regulates the function of living cells through protein
synthesis. The key contributions of the thesis are general exploratory tools
for high-throughput data analysis that have provided new insights to
cell-biological networks, cancer mechanisms and other aspects of genome
function.
  A central challenge in functional genomics is that high-dimensional genomic
observations are associated with high levels of complex and largely unknown
sources of variation. By combining statistical evidence across multiple
measurement sources and the wealth of background information in genomic data
repositories it has been possible to solve some the uncertainties associated
with individual observations and to identify functional mechanisms that could
not be detected based on individual measurement sources. Statistical learning
and probabilistic models provide a natural framework for such modeling tasks.
Open source implementations of the key methodological contributions have been
released to facilitate further adoption of the developed methods by the
research community.


  Combined with all density-dependent factors, the per capita growth rate of a
species may be non-monotonic. One important consequence is that species may
suffer from weak Allee effects or strong Allee effects. In this paper, we study
the permanence of a discrete-time two-species-interaction model with
non-monotonic per capita growth rates for the first time. By using the average
Lyapunov functions and extending the ecological concept of the relative
nonlinearity, we find a simple sufficient condition for guaranteeing the
permanence of systems that can model complicated two-species interactions. The
extended relative nonlinearity allows us to fully characterize the effects of
nonlinearities in the per capita growth functions with non-monotonicity. These
results are illustrated with specific two species competition and predator-prey
models of generic forms with non-monotone per capita growth rates.


  Background: Human populations during the last 10,000 years have undergone
rapid decreases in average brain size as measured by endocranial volume or as
estimated from linear measurements of the cranium. A null hypothesis to explain
the evolution of brain size is that reductions result from genetic correlation
of brain size with body mass or stature.
  Results: The absolute change of endocranial volume in the study samples was
significantly greater than would be predicted from observed changes in body
mass or stature.
  Conclusions: The evolution of smaller brains in many recent human populations
must have resulted from selection upon brain size itself or on other features
more highly correlated with brain size than are gross body dimensions. This
selection may have resulted from energetic or nutritional demands in Holocene
populations, or to life history constraints on brain development.


  We develop a transfer matrix formalism to visualize the framing of discrete
piecewise linear curves in three dimensional space. Our approach is based on
the concept of an intrinsically discrete curve, which enables us to more
effectively describe curves that in the limit where the length of line segments
vanishes approach fractal structures in lieu of continuous curves. We verify
that in the case of differentiable curves the continuum limit of our discrete
equation does reproduce the generalized Frenet equation. As an application we
consider folded proteins, their Hausdorff dimension is known to be fractal. We
explain how to employ the orientation of $C_\beta$ carbons of amino acids along
a protein backbone to introduce a preferred framing along the backbone. By
analyzing the experimentally resolved fold geometries in the Protein Data Bank
we observe that this $C_\beta$ framing relates intimately to the discrete
Frenet framing. We also explain how inflection points can be located in the
loops, and clarify their distinctive r\^ole in determining the loop structure
of foldel proteins.


  We investigate the evolutionary dynamics of an idealised model for the robust
self-assembly of two-dimensional structures called polyominoes. The model
includes rules that encode interactions between sets of square tiles that drive
the self-assembly process. The relationship between the model's rule set and
its resulting self-assembled structure can be viewed as a genotype-phenotype
map and incorporated into a genetic algorithm. The rule sets evolve under
selection for specified target structures. The corresponding, complex fitness
landscape generates rich evolutionary dynamics as a function of parameters such
as the population size, search space size, mutation rate, and method of
recombination. Furthermore, these systems are simple enough that in some cases
the associated model genome space can be completely characterised, shedding
light on how the evolutionary dynamics depends on the detailed structure of the
fitness landscape. Finally, we apply the model to study the emergence of the
preference for dihedral over cyclic symmetry observed for homomeric protein
tetramers.


  The green area of economy is the key of healthy living. It is necessary to
convene economic and ecologic framework to establish a market attentive to
drastic reduction of emissions damaging our climate and landscapes in rural
areas, to the protection of biological diversity of the planet, to stop
producing nuclear waste, etc. This paper tries to demonstrate human concern for
a waste recycling economy that will provide new jobs, will create economic and
social stability and will ensure a healthier and cleaner environment. Green
Economy and its support system (planetary ecosystem) won't be in conflict
anymore. Green Economy will be able to support economic progress for future.


  This article presents a theoretical model for a dynamic system based on
sustainable development. Due to the relatively absence of theoretical studies
and practical issues in the area of sustainable development, Romania aspires to
the principles of sustainable development. Based on the concept as a process in
which economic, social, political and natural environment are combined in order
to sustain planet management, our goal is to promote an economic tool for
Romanian decision-makers in order to evaluate scenarios and planning options.


  We review models of biological evolution in which the population frequency
changes deterministically with time. If the population is self-replicating,
although the equations for simple prototypes can be linearised, nonlinear
equations arise in many complex situations. For sexual populations, even in the
simplest setting, the equations are necessarily nonlinear due to the mixing of
the parental genetic material. The solutions of such nonlinear equations
display interesting features such as multiple equilibria and phase transitions.
We mainly discuss those models for which an analytical understanding of such
nonlinear equations is available.


  We apply the Wako-Saito-Munoz-Eaton model to the study of Myotrophin, a small
ankyrin repeat protein, whose folding equilibrium and kinetics have been
recently characterized experimentally. The model, which is a native-centric
with binary variables, provides a finer microscopic detail than the Ising
model, that has been recently applied to some different repeat proteins, while
being still amenable for an exact solution. In partial agreement with the
experiments, our results reveal a weakly three-state equilibrium and a
two-state-like kinetics of the wild type protein despite the presence of a
non-trivial free-energy profile. These features appear to be related to a
careful "design" of the free-energy landscape, so that mutations can alter this
picture, stabilizing some intermediates and changing the position of the
rate-limiting step. Also the experimental findings of two alternative pathways,
an N-terminal and a C-terminal one, are qualitatively confirmed, even if the
variations in the rates upon the experimental mutations cannot be
quantitatively reproduced. Interestingly, folding and unfolding pathway appear
to be different, even if closely related: a property that is not generally
considered in the phenomenological interpretation of the experimental data.


  In a community-structured population, public goods games (PGG) occur both
within and between communities. Such type of PGG is referred as multilevel
public goods games (MPGG). We propose a minimalist evolutionary model of the
MPGG and analytically study the evolution of cooperation. We demonstrate that
in the case of sufficiently large community size and community number, if the
imitation strength within community is weak, i.e., an individual imitates
another one in the same community almost randomly, cooperation as well as
punishment are more abundant than defection in the long run; if the imitation
strength between communities is strong, i.e., the more successful strategy in
two individuals from distinct communities is always imitated, cooperation and
punishment are also more abundant. However, when both of the two imitation
intensities are strong, defection becomes the most abundant strategy in the
population. Our model provides insight into the investigation of the
large-scale cooperation in public social dilemma among contemporary
communities.


  One approach to estimating a species tree from a collection of gene trees is
to first estimate probabilities of clades from the gene trees, and then to
construct the species tree from the estimated clade probabilities. While a
greedy consensus algorithm, which consecutively accepts the most probable
clades compatible with previously accepted clades, can be used for this second
stage, this method is known to be statistically inconsistent under the
multispecies coalescent model. This raises the question of whether it is
theoretically possible to reconstruct the species tree from known probabilities
of clades on gene trees. We investigate clade probabilities arising from the
multispecies coalescent model, with an eye toward identifying features of the
species tree. Clades on gene trees with probability greater than 1/3 are shown
to reflect clades on the species tree, while those with smaller probabilities
may not. Linear invariants of clade probabilities are studied both
computationally and theoretically, with certain linear invariants giving
insight into the clade structure of the species tree. For species trees with
generic edge lengths, these invariants can be used to identify the species tree
topology. These theoretical results both confirm that clade probabilities
contain full information on the species tree topology and suggest future
directions of study for developing statistically consistent inference methods
from clade frequencies on gene trees.


  We investigate the force generation by polymerizing bundles of filaments,
which form because of short-range attractive filament interactions. We show
that bundles can generate forces by a zipping mechanism, which is not limited
by buckling and operates in the fully buckled state. The critical zipping
force, i.e. the maximal force that a bundle can generate, is given by the
adhesive energy gained during bundle formation. For opposing forces larger than
the critical zipping force, bundles undergo a force-induced unbinding
transition. For larger bundles, the critical zipping force depends on the
initial configuration of the bundles. Our results are corroborated by Monte
Carlo simulations.


  We establish a framework for assessing whether the transition state location
of a biopolymer, which can be inferred from single molecule pulling
experiments, corresponds to the ensemble of structures that have equal
probability of reaching either the folded or unfolded states (Pfold = 0.5).
Using results for the forced-unfolding of a RNA hairpin, an exactly soluble
model and an analytic theory, we show that Pfold is solely determined by s, an
experimentally measurable molecular tensegrity parameter, which is a ratio of
the tensile force and a compaction force that stabilizes the folded state.
Applications to folding landscapes of DNA hairpins and leucine zipper with two
barriers provide a structural interpretation of single molecule experimental
data. Our theory can be used to assess whether molecular extension is a good
reaction coordinate using measured free energy profiles.


  A detailed description of the conformational plasticity of double stranded
DNA (ds) is a necessary framework for understanding protein-DNA interactions.
Until now, however structure and kinetics of the transition from the basic
conformation of ds-DNA (B state) to the 1.7 times longer and partially unwound
conformation (S state) have not been defined. The force-extension relation of
the ds-DNA of lambda-phage is measured here with unprecedented resolution using
a dual laser optical tweezers that can impose millisecond force steps of 0.5-2
pN (25 C). This approach reveals the kinetics of the transition between
intermediate states of ds-DNA and uncovers the load-dependence of the rate
constant of the unitary reaction step. DNA overstretching transition results
essentially a two-state reaction composed of 5.85 nm steps, indicating
cooperativity of ~25 base pairs. This mechanism increases the free energy for
the unitary reaction to ~94 kBT, accounting for the stability of the basic
conformation of DNA, and explains the absence of hysteresis in the
force-extension relation at equilibrium. The novel description of the kinetics
and energetics of the B-S transition of ds-DNA improves our understanding the
biological role of the S state in the interplay between mechanics and
enzymology of the DNA-protein machinery.


  A simple method, inspired by procedures used in physics of nuclear
multifragmentation, allows to establish order of precedence and age of pairs of
haplotypes separated by one mutation. For both haplotypes of the pair, searches
for existing haplotypes, differing by increasing number of mutations, are
carried out using a database. The resulting ratios of frequences of haplotypes,
found at given mutation distances, are compared to calculated probability
ratios. The order of precedence and age of the pair of haplotypes can be
deduced when the resulting ratios follow hyperbolic dependence. Method can be
used with relatively small and not necessarily complete samples, using publicly
accessible databases.


  Motivation: Modeling biological signaling networks using ordinary
differential equations (ODEs) has proven to be a powerful technique for
generating insight into cellular dynamics, but it typically requires estimating
rate parameters based on experimentally observed concentrations. New
measurement methods can measure concentrations for all molecular species in a
pathway, which creates a new opportunity to decompose the optimization of rate
parameters. Results: In contrast with conventional methods that minimize the
disagreement between simulated and observed concentrations, the BPPE method
fits a spline curve through the observed concentration points, and then matches
the derivatives of the spline-curve to the production and consumption of each
species. Whereas traditional methods follow the ODEs exactly and then attempt
to match the data, BPPE follows the data exactly and then attempts to match the
ODEs. The new objective function is an extreme decomposition of the problem
because each factor in the function is enforcing the equality of one ODE at one
timeslice. A "loopy belief propagation" algorithm solves this factorized
approximation of the parameter estimation problem providing systematic coverage
of the search space and unique asymptotic behavior; the run time is polynomial
in the number of molecules and timepoints, but exponential in the degree of the
biochemical network. The implementation is a global-local hybrid optimization,
and we compare with the performance of local, global, and hybrid methods. BPPE
is demonstrated for a novel model of Akt activation dynamics including
redox-mediated inactivation of PTEN. Availability: Software and supplementary
information are available at http://webbppe.nus.edu.sg:8080/opal2/WebBPPE .
Contact: LisaTK@nus.edu.sg . Keywords: probabilistic graphical models,
physico-chemical modeling, systems biology, signal transduction.


  We apply theoretically open-loop quantum optimal control techniques to
provide methods for the verification of various quantum coherent transport
mechanisms in natural and artificial light-harvesting complexes under realistic
experimental constraints. We demonstrate that optimally shaped laser pulses
allow to faithfully prepare the photosystem in specified initial states (such
as localized excitation or coherent superposition, i.e. propagating and
non-propagating states) and to probe efficiently the dynamics. These results
provide a path towards the discrimination of the different transport pathways
and to the characterization of environmental properties, enhancing our
understanding of the role that coherent processes may play in biological
complexes.


  Here we introduce a new design framework for synthetic biology that exploits
the advantages of Bayesian model selection. We will argue that the difference
between inference and design is that in the former we try to reconstruct the
system that has given rise to the data that we observe, while in the latter, we
seek to construct the system that produces the data that we would like to
observe, i.e. the desired behavior. Our approach allows us to exploit methods
from Bayesian statistics, including efficient exploration of models spaces and
high-dimensional parameter spaces, and the ability to rank models with respect
to their ability to generate certain types of data. Bayesian model selection
furthermore automatically strikes a balance between complexity and (predictive
or explanatory) performance of mathematical models. In order to deal with the
complexities of molecular systems we employ an approximate Bayesian computation
scheme which only requires us to simulate from different competing models in
order to arrive at rational criteria for choosing between them. We illustrate
the advantages resulting from combining the design and modeling (or in-silico
prototyping) stages currently seen as separate in synthetic biology by
reference to deterministic and stochastic model systems exhibiting adaptive and
switch-like behavior, as well as bacterial two-component signaling systems.


  Robust advances in interactome analysis demand comprehensive, non-redundant
and consistently annotated datasets. By non-redundant, we mean that the
accounting of evidence for every interaction should be faithful: each
independent experimental support is counted exactly once, no more, no less.
While many interactions are shared among public repositories, none of them
contains the complete known interactome for any model organism. In addition,
the annotations of the same experimental result by different repositories often
disagree. This brings up the issue of which annotation to keep while
consolidating evidences that are the same. The iRefIndex database, including
interactions from most popular repositories with a standardized protein
nomenclature, represents a significant advance in all aspects, especially in
comprehensiveness. However, iRefIndex aims to maintain all
information/annotation from original sources and requires users to perform
additional processing to fully achieve the aforementioned goals.
  To address issues with iRefIndex and to achieve our goals, we have developed
ppiTrim, a script that processes iRefIndex to produce non-redundant,
consistently annotated datasets of physical interactions. Our script proceeds
in three stages: mapping all interactants to gene identifiers and removing all
undesired raw interactions, deflating potentially expanded complexes, and
reconciling for each interaction the annotation labels among different source
databases. As an illustration, we have processed the three largest organismal
datasets: yeast, human and fruitfly. While ppiTrim can resolve most apparent
conflicts between different labelings, we also discovered some unresolvable
disagreements mostly resulting from different annotation policies among
repositories.
  URL: http://www.ncbi.nlm.nih.gov/CBBresearch/Yu/downloads/ppiTrim.html


  Exponential distribution is ubiquitous in the framework of multi-agent
systems. Usually, it appears as an equilibrium state in the asymptotic time
evolution of statistical systems. It has been explained from very different
perspectives. In statistical physics, it is obtained from the principle of
maximum entropy. In the same context, it can also be derived without any
consideration about information theory, only from geometrical arguments under
the hypothesis of equiprobability in phase space. Also, several multi-agent
economic models based on mappings, with random, deterministic or chaotic
interactions, can give rise to the asymptotic appearance of the exponential
wealth distribution. An alternative approach to this problem in the framework
of iterations in the space of distributions has been recently presented.
Concretely, the new iteration given by $ f_{n+1}(x) =
\int\int_{u+v>x}{f_n(u)f_n(v)\over u+v} dudv.$. It is found that the
exponential distribution is a stable fixed point of the former functional
iteration equation. From this point of view, it is easily understood why the
exponential wealth distribution (or by extension, other kind of distributions)
is asymptotically obtained in different multi-agent economic models.


  Population dynamics in systems composed of cyclically competing species has
been of increasing interest recently. Here, we investigate a system with four
or more species. Using mean field theory, we study in detail the trajectories
in configuration space of the population fractions. We discover a variety of
orbits, shaped like saddles, spirals, and straight lines. Many of their
properties are found explicitly. Most remarkably, we identify a collective
variable which evolves simply as an exponential: $\mathcal{Q}% \propto
e^{\lambda t}$, where $\lambda$ is a function of the reaction rates. It
provides information on the state of the system for late times (as well as for
$t\rightarrow -\infty $). We discuss implications of these results for the
evolution of a finite, stochastic system. A generalization to an arbitrary
number of cyclically competing species yields valuable insights into universal
properties of such systems.


  Ants are known to be able to find paths of minimal length between the nest
and food sources. The deposit of pheromones while they search for food and
their chemotactical response to them has been proposed as a crucial element in
the mechanism for finding minimal paths. We investigate both individual and
collective behavior of ants in some simple networks representing basic mazes.
The character of the graphs considered is such that it allows a fully rigorous
mathematical treatment via analysis of some markovian processes in terms of
which the evolution can be represented. Our analytical and computational
results show that in order for the ants to follow shortest paths between nest
and food, it is necessary to superimpose to the ants' random walk the
chemotactic reinforcement. It is also needed a certain degree of persistence so
that ants tend to move preferably without changing their direction much. It is
also important the number of ants, since we will show that the speed for
finding minimal paths increases very fast with it.


  Food-web perturbations stemming from climate change, overexploitation,
invasive species, and habitat degradation often cause an initial loss of
species that results in a cascade of secondary extinctions, posing considerable
challenges to ecosystem conservation efforts. Here we devise a systematic
network-based approach to reduce the number of secondary extinctions using a
predictive modeling framework. We show that the extinction of one species can
often be compensated by the concurrent removal or population suppression of
other specific species, which is a counterintuitive effect not previously
tested in complex food webs. These compensatory perturbations frequently
involve long-range interactions that are not evident from local predator-prey
relationships. In numerous cases, even the early removal of a species that
would eventually be extinct by the cascade is found to significantly reduce the
number of cascading extinctions. These compensatory perturbations only exploit
resources available in the system, and illustrate the potential of human
intervention combined with predictive modeling for ecosystem management.


  One of the hallmarks of biological organisms is their ability to integrate
disparate information sources to optimize their behavior in complex
environments. How this capability can be quantified and related to the
functional complexity of an organism remains a challenging problem, in
particular since organismal functional complexity is not well-defined. We
present here several candidate measures that quantify information and
integration, and study their dependence on fitness as an artificial agent
("animat") evolves over thousands of generations to solve a navigation task in
a simple, simulated environment. We compare the ability of these measures to
predict high fitness with more conventional information-theoretic processing
measures. As the animat adapts by increasing its "fit" to the world,
information integration and processing increase commensurately along the
evolutionary line of descent. We suggest that the correlation of fitness with
information integration and with processing measures implies that high fitness
requires both information processing as well as integration, but that
information integration may be a better measure when the task requires memory.
A correlation of measures of information integration (but also information
processing) and fitness strongly suggests that these measures reflect the
functional complexity of the animat, and that such measures can be used to
quantify functional complexity even in the absence of fitness data.


  Rooted phylogenetic networks are often used to represent conflicting
phylogenetic signals. Given a set of clusters, a network is said to represent
these clusters in the "softwired" sense if, for each cluster in the input set,
at least one tree embedded in the network contains that cluster. Motivated by
parsimony we might wish to construct such a network using as few reticulations
as possible, or minimizing the "level" of the network, i.e. the maximum number
of reticulations used in any "tangled" region of the network. Although these
are NP-hard problems, here we prove that, for every fixed k >= 0, it is
polynomial-time solvable to construct a phylogenetic network with level equal
to k representing a cluster set, or to determine that no such network exists.
However, this algorithm does not lend itself to a practical implementation. We
also prove that the comparatively efficient Cass algorithm correctly solves
this problem (and also minimizes the reticulation number) when input clusters
are obtained from two not necessarily binary gene trees on the same set of taxa
but does not always minimize level for general cluster sets. Finally, we
describe a new algorithm which generates in polynomial-time all binary
phylogenetic networks with exactly r reticulations representing a set of input
clusters (for every fixed r >= 0).


  Kinesin-5, also known as Eg5 in vertebrates is a processive motor with 4
heads, which moves on filamentous tracks called microtubules. The basic
function of Kinesin-5 is to slide apart two anti-parallel microtubules by
simultaneously walking on both the microtubules. We develop an analytical
expression for the steady-state relative velocity of this sliding in terms of
the rates of attachments and detachments of motor heads with the ATPase sites
on the microtubules. We first analyse the motion of one pair of motor heads on
one microtubule and then couple it to the motion of the other pair of motor
heads of the same motor on the second microtubule to get the relative velocity
of sliding.


  For a given set $\mathcal{L}$ of species and a set $\mathcal{T}$ of triplets
on $\mathcal{L}$, one wants to construct a phylogenetic network which is
consistent with $\mathcal{T}$, i.e which represents all triplets of
$\mathcal{T}$. The level of a network is defined as the maximum number of
hybrid vertices in its biconnected components. When $\mathcal{T}$ is dense,
there exist polynomial time algorithms to construct level-$0,1,2$ networks (Aho
et al. 81, Jansson et al. 04, Iersel et al. 08). For higher levels, partial
answers were obtained by Iersel et al. 2008 with a polynomial time algorithm
for simple networks. In this paper, we detail the first complete answer for the
general case, solving a problem proposed by Jansson et al. 2004: for any $k$
fixed, it is possible to construct a minimum level-$k$ network consistent with
$\mathcal{T}$, if there is any, in time
$O(|\mathcal{T}|^{k+1}n^{\lfloor\frac{4k}{3}\rfloor+1})$. This is an improved
result of our preliminary version presented at CPM'2009.


  Materials in biology span all the scales from Angstroms to meters and
typically consist of complex hierarchical assemblies of simple building blocks.
Here we describe an application of category theory to describe structural and
resulting functional properties of biological protein materials by developing
so-called ologs. An olog is like a "concept web" or "semantic network" except
that it follows a rigorous mathematical formulation based on category theory.
This key difference ensures that an olog is unambiguous, highly adaptable to
evolution and change, and suitable for sharing concepts with other olog. We
consider simple cases of alpha-helical and amyloid-like protein filaments
subjected to axial extension and develop an olog representation of their
structural and resulting mechanical properties. We also construct a
representation of a social network in which people send text-messages to their
nearest neighbors and act as a team to perform a task. We show that the olog
for the protein and the olog for the social network feature identical
category-theoretic representations, and we proceed to precisely explicate the
analogy or isomorphism between them. The examples presented here demonstrate
that the intrinsic nature of a complex system, which in particular includes a
precise relationship between structure and function at different hierarchical
levels, can be effectively represented by an olog. This, in turn, allows for
comparative studies between disparate materials or fields of application, and
results in novel approaches to derive functionality in the design of de novo
hierarchical systems. We discuss opportunities and challenges associated with
the description of complex biological materials by using ologs as a powerful
tool for analysis and design in the context of materiomics, and we present the
potential impact of this approach for engineering, life sciences, and medicine.


  Recent whole genome polymerase binding assays have shown that a large
proportion of unexpressed genes have pre-assembled RNA pol II transcription
initiation complex stably bound to their promoters. Some such promoter
proximally paused genes are regulated at transcription elongation rather than
at initiation; it has been proposed that this difference allows these genes to
both express faster and achieve more synchronous expression across populations
of cells, thus overcoming molecular "noise" arising from low copy number
factors. It has been established experimentally that genes which are regulated
at elongation tend to express faster and more synchronously; however, it has
not been shown directly whether or not it is the change in the regulated step
{\em per se} that causes this increase in speed and synchrony. We investigate
this question by proposing and analyzing a continuous-time Markov chain model
of polymerase complex assembly regulated at one of two steps: initial
polymerase association with DNA, or release from a paused, transcribing state.
Our analysis demonstrates that, over a wide range of physical parameters,
increased speed and synchrony are functional consequences of elongation
control. Further, we make new predictions about the effect of elongation
regulation on the consistent control of total transcript number between cells,
and identify which elements in the transcription induction pathway are most
sensitive to molecular noise and thus may be most evolutionarily constrained.
Our methods produce symbolic expressions for quantities of interest with
reasonable computational effort and can be used to explore the interplay
between interaction topology and molecular noise in a broader class of
biochemical networks. We provide general-purpose code implementing these
methods.


  Functional effects of different mutations are known to combine to the total
effect in highly nontrivial ways. For the trait under evolutionary selection
(`fitness'), measured values over all possible combinations of a set of
mutations yield a fitness landscape that determines which mutational states can
be reached from a given initial genotype. Understanding the accessibility
properties of fitness landscapes is conceptually important in answering
questions about the predictability and repeatability of evolutionary
adaptation. Here we theoretically investigate accessibility of the globally
optimal state on a wide variety of model landscapes, including landscapes with
tunable ruggedness as well as neutral `holey' landscapes. We define a
mutational pathway to be accessible if it contains the minimal number of
mutations required to reach the target genotype, and if fitness increases in
each mutational step. Under this definition accessibility is high, in the sense
that at least one accessible pathwayexists with a substantial probability that
approaches unity as the dimensionality of the fitness landscape (set by the
number of mutational loci) becomes large. At the same time the number of
alternative accessible pathways grows without bound. We test the model
predictions against an empirical 8-locus fitness landscape obtained for the
filamentous fungus \textit{Aspergillus niger}. By analyzing subgraphs of the
full landscape containing different subsets of mutations, we are able to probe
the mutational distance scale in the empirical data. The predicted effect of
high accessibility is supported by the empirical data and very robust, which we
argue to reflect the generic topology of sequence spaces.


  We undertake a systematic numerical exploration of self-organized states in a
deterministic model of interacting self-propelled particles in two dimensions.
In the process, we identify various types of collective motion, namely,
disordered swarms, rings, and droplets. We construct a "phase diagram," which
summarizes our results as it delineates phase transitions ?all discontinuous?
between disordered swarms and vortical flocks on one hand and bound vortical
flocks and expanding formations on the other. One of the transition lines is
found to have a close analogy with the temperature-driven gas-liquid
transition, in finite clusters with the same interparticle potential.
Furthermore, we report on a type of flocking which takes place in the presence
of a uniform external driver. Altogether, our results set a rather firm stage
for experimental refinement and/or falsification of this class of models.


  Many online marketplaces enjoy great success. Buyers and sellers in
successful markets carry out cooperative transactions even if they do not know
each other in advance and a moral hazard exists. An indispensable component
that enables cooperation in such social dilemma situations is the reputation
system. Under the reputation system, a buyer can avoid transacting with a
seller with a bad reputation. A transaction in online marketplaces is better
modeled by the trust game than other social dilemma games, including the
donation game and the prisoner's dilemma. In addition, most individuals
participate mostly as buyers or sellers; each individual does not play the two
roles with equal probability. Although the reputation mechanism is known to be
able to remove the moral hazard in games with asymmetric roles, competition
between different strategies and population dynamics of such a game are not
sufficiently understood. On the other hand, existing models of reputation-based
cooperation, also known as indirect reciprocity, are based on the symmetric
donation game. We analyze the trust game with two fixed roles, where trustees
(i.e., sellers) but not investors (i.e., buyers) possess reputation scores. We
study the equilibria and the replicator dynamics of the game. We show that the
reputation mechanism enables cooperation between unacquainted buyers and
sellers under fairly generous conditions, even when such a cooperative
equilibrium coexists with an asocial equilibrium in which buyers do not buy and
sellers cheat. In addition, we show that not many buyers may care about the
seller's reputation under cooperative equilibrium. Buyers' trusting behavior
and sellers' reputation-driven cooperative behavior coevolve to alleviate the
social dilemma.