  Biofilms are matrix-enclosed microbial communities that grow at interfaces.
They exhibit significant phenotypic changes that render them resistant to many
antibacterial agents that can kill their free-swimming counterparts. Biofilms
in the environment can consist of either a single or multiple microbial
species. Previous investigations in this group demonstrated that Bdellovibrio
bacteriovorus, a gram-negative predatory bacterium that consumes other
gram-negative bacteria, could eradicate a single-species E.coli biofilm in some
conditions1. These results suggest the potential value of B. bacteriovorus in
biofilm eradication in industrial, medical and environmental contexts.
  Here we carried out investigation of the potential of Bdellovibrio
bacteriovorus to interact with and remove multi-species biofilms, specifically
dual-species biofilms of E. coli (gram-negative prey) and M. luteus
(gram-positive decoy) at interfaces between liquid and air. We examined whether
this kind of dual-species biofilms is more, less, or equally susceptible to B.
bacteriovorus attack compared to their single-species counterparts. Different
research methods including bacterial culture, cell counting, crystal violet
staining, gram staining, optical microscopy, scanning electron microscopy, and
atomic force microscopy were exploited to gain an insight into Bdellovibrio's
interaction with biofilms. Our experiments showed that B. bacteriovorus not
only consumes E. coli but also weakens the attachment of M. luteus to the solid
surface, rendering the biofilms susceptible to removal. B. bacteriovorus keeps
in check not only prey but also decoy bacterial populations in the surrounding
media, the latter probably via competition for nutrients. These experiments
encourage us to consider how B. bacteriovorus might be used to control biofilms
in the environment.


  We present a Monte Carlo algorithm that provides efficient and unbiased
sampling of polymer melts consisting of two chains of equal length that jointly
visit all the sites of a cubic lattice with rod geometry L x L x rL and
non-periodic (hard wall) boundary conditions. Using this algorithm for chains
of length up to 40 000 monomers and aspect ratios 1 <= r <= 10, we show that in
the limit of a large lattice the two chains phase separate. This demixing
phenomenon is present already for r=1, and becomes more pronounced, albeit not
perfect, as r is increased.


  Evolutionary branching points are a paradigmatic feature of adaptive
dynamics, because they are potential starting points for adaptive
diversification. The antithesis to evolutionary branching points are
Continuously stable strategies (CSS's), which are convergent stable and
evolutionarily stable equilibrium points of the adaptive dynamics and hence are
thought to represent endpoints of adaptive processes. However, this assessment
is based on situations in which the invasion fitness function determining the
adaptive dynamics have non-zero second derivatives at a CSS. Here we show that
the scope of evolutionary branching can increase if the invasion fitness
function vanishes to higher than first order at a CSS. Using a class of
classical models for frequency-dependent competition, we show that if the
invasion fitness vanishes to higher orders, a CSS may be the starting point for
evolutionary branching, with the only additional requirement that mutant types
need to reach a certain threshold frequency, which can happen e.g. due to
demographic stochasticity. Thus, when invasion fitness functions vanish to
higher than first order at equilibrium points of the adaptive dynamics,
evolutionary diversification can occur even after convergence to an
evolutionarily stable strategy.


  The electronic structure of charybdotoxin (ChTX), a scorpion venom peptide
that is known to act as a potassium channel blocker, is investigated with the
aid of quantum mechanical calculations. The dipole moment vector (145 D) of
ChTX can be stirred by the full length KcsA potassium channel's macrodipole
(403 D) thereby assuming the proper orientation before binding the ion channel
on the cell surface. The localization of the frontier orbitals of ChTX has been
revealed for the first time. HOMO is localized on Trp14 while the three
lowest-energy MOs (LUMO, LUMO+1, and LUMO+2) are localized on the three
disulfide bonds that characterize this pepetide. An effective way to engineer
the HOMO-LUMO (H-L) gap of ChTX is that of replacing its Trp14 residue with
Ala14 whereas deletion of the LUMO-associated disulfide bond with the insertion
of a pair of L-alpha-aminobutyric acid residues does not affect the H-L energy
gap.


  Anthelmintic resistance is a major threat to current measures for helminth
control in humans and animals. The introduction of anthelmintic vaccines, as a
complement to or replacement for drug treatments, has been advocated as a
preventive measure. Here, a computer-based simulation, tracking the dynamics of
hosts, parasites and parasite-genes, shows that, depending on the degree of
host-population mixing, the frequency of totally recessive autosomes associated
with anthelmintic resistance can follow either a fast dynamical regime with a
low equilibrium point or a slow dynamical regime with a high equilibrium point.
For fully dominant autosomes, only one regime is predicted. The effectiveness
of anthelminthic vaccines against resistance is shown to be strongly influenced
by the underlying dynamics of resistant autosomes. Vaccines targeting adult
parasites, by decreasing helminth fecundity or lifespan, are predicted to be
more effective than vaccines targeting parasite larvae, by decreasing host
susceptibility to infection, in reducing the spread of resistance. These
results may inform new strategies to prevent, monitor and control the spread of
anthelmintic resistance, including the development of viable anthelmintic
vaccines.


  We present the latest advances of the multiscale approach to radiation damage
caused by irradiation of a tissue with energetic ions and report the most
recent advances in the calculations of complex DNA damage and the effects of
thermal spikes on biomolecules. The multiscale approach aims to quantify the
most important physical, chemical, and biological phenomena taking place during
and following irradiation with ions and provide a better means for
clinically-necessary calculations with adequate accuracy. We suggest a way of
quantifying the complex clustered damage, one of the most important features of
the radiation damage caused by ions. This method can be used for the
calculation of irreparable DNA damage. We include thermal spikes, predicted to
occur in tissue for a short time after ion's passage in the vicinity of the
ions' tracks in our previous work, into modeling of the thermal environment for
molecular dynamics analysis of ubiquitin and discuss the first results of these
simulations.


  We propose a mechanism that explains in a simple and natural form the
l-homochiralization of prebiotic aminoacids in a volume of water where a
geothermal gradient exists.


  The dynamics of interacting structured populations can be modeled by
$\frac{dx_i}{dt}= A_i (x)x_i$ where $x_i\in \R^{n_i}$, $x=(x_1,\dots,x_k)$, and
$A_i(x)$ are matrices with non-negative off-diagonal entries. These models are
permanent if there exists a positive global attractor and are robustly
permanent if they remain permanent following perturbations of $A_i(x)$.
Necessary and sufficient conditions for robust permanence are derived using
dominant Lyapunov exponents $\lambda_i(\mu)$ of the $A_i(x)$ with respect to
invariant measures $\mu$. The necessary condition requires $\max_i
\lambda_i(\mu)>0$ for all ergodic measures with support in the boundary of the
non-negative cone. The sufficient condition requires that the boundary admits a
Morse decomposition such that $\max_i \lambda_i(\mu)>0$ for all invariant
measures $\mu$ supported by a component of the Morse decomposition. When the
Morse components are Axiom A, uniquely ergodic, or support all but one
population, the necessary and sufficient conditions are equivalent.
Applications to spatial ecology, epidemiology, and gene networks are given.


  Lateral inhibition patterns mediated by the Notch-Delta signaling system
occur in diverse developmental contexts. These systems are based on an
intercellular feedback loop in which Notch activation leads to down-regulation
of Delta. However, even in relatively well-characterized systems, the pathway
leading from Notch activation to Delta repression often remains elusive. Recent
work has shown that cis-interactions between Notch and Delta lead to mutual
inactivation of both proteins. Here we show that this type of cis-interaction
enables a simpler and more direct mechanism for lateral inhibition feedback
than those proposed previously. In this mechanism, Notch signaling directly
up-regulates Notch expression, thereby inactivating Delta through the mutual
inactivation of Notch and Delta proteins. This mechanism, which we term
Simplest Lateral Inhibition by Mutual Inactivation (SLIMI), can implement
patterning without requiring any additional genes or regulatory interactions.
Moreover, the key interaction of Notch expression in response to Notch
signaling has been observed in some systems. Stability analysis and simulation
of SLIMI mathematical models show that this lateral inhibition circuit is
capable of pattern formation across a broad range of parameter values. These
results provide a simple and plausible explanation for lateral inhibition
pattern formation during development.


  Biodiversity widely observed in ecological systems is attributed to the
dynamical balance among the competing species. The time-varying populations of
the interacting species are often captured rather well by a set of
deterministic replicator equations in the evolutionary game theory. However,
intrinsic fluctuations arisen from the discreteness of populations lead to
stochastic derivations from the smooth evolution trajectories. The role of
these fluctuations is shown to be critical at causing extinction and
deteriorating the biodiversity of ecosystem. We use children's
rock-paper-scissors game to demonstrate how the intrinsic fluctuations arise
from the discrete populations and why the biodiversity of the ecosystem decays
exponentially, disregarding the detail parameters for competing mechanism and
initial distributions. The dissipative trend in biodiversity can be analogized
to the gradual erosion of kinetic energy of a moving particle due to air drag
or fluid viscosity. The dissipation-fluctuation theorem in statistical physics
seals the fate of these originally conserved quantities. This concept in
physics can be generalized to scrutinize the errors that might be incurred in
the ecological, biological, and quantitative economic modeling for which the
ingredients are all discrete in number.


  We investigate how scale-free (SF) and Erdos-Renyi (ER) topologies affect the
interplay between evolvability and robustness of model gene regulatory networks
with Boolean threshold dynamics. In agreement with Oikonomou and Cluzel (2006)
we find that networks with SFin topologies, that is SF topology for incoming
nodes and ER topology for outgoing nodes, are significantly more evolvable
towards specific oscillatory targets than networks with ER topology for both
incoming and outgoing nodes. Similar results are found for networks with SFboth
and SFout topologies. The functionality of the SFout topology, which most
closely resembles the structure of biological gene networks (Babu et al.,
2004), is compared to the ER topology in further detail through an extension to
multiple target outputs, with either an oscillatory or a non-oscillatory
nature. For multiple oscillatory targets of the same length, the differences
between SFout and ER networks are enhanced, but for non-oscillatory targets
both types of networks show fairly similar evolvability. We find that SF
networks generate oscillations much more easily than ER networks do, and this
may explain why SF networks are more evolvable than ER networks are for
oscillatory phenotypes. In spite of their greater evolvability, we find that
networks with SFout topologies are also more robust to mutations than ER
networks. Furthermore, the SFout topologies are more robust to changes in
initial conditions (environmental robustness). For both topologies, we find
that once a population of networks has reached the target state, further
neutral evolution can lead to an increase in both the mutational robustness and
the environmental robustness to changes in initial conditions.


  We use Fourier analysis and related techniques to investigate the question of
periodicities in fossil biodiversity. These techniques are able to identify
cycles superimposed on the long-term trends of the Phanerozoic. We review prior
results and analyze data previously reduced and published. Joint time-series
analysis of various reductions of the Sepkoski Data, Paleobiology Database, and
Fossil Record 2 indicate the same periodicity in biodiversity of marine animals
at 62 Myr. We have not found this periodicity in the terrestrial fossil record.
We have found that the signal strength decreases with time because of the
accumulation of apparently "resistant" long-lived genera. The existence of a
62-Myr periodicity despite very different treatment of systematic error,
particularly sampling-strength biases, in all three major databases strongly
argues for its reality in the fossil record.


  The set of regulatory interactions between genes, mediated by transcription
factors, forms a species' transcriptional regulatory network (TRN). By
comparing this network with measured gene expression data one can identify
functional properties of the TRN and gain general insight into transcriptional
control. We define the subnet of a node as the subgraph consisting of all nodes
topologically downstream of the node, including itself. Using a large set of
microarray expression data of the bacterium Escherichia coli, we find that the
gene expression in different subnets exhibits a structured pattern in response
to environmental changes and genotypic mutation. Subnets with less changes in
their expression pattern have a higher fraction of feed-forward loop motifs and
a lower fraction of small RNA targets within them. Our study implies that the
TRN consists of several scales of regulatory organization: 1) subnets with more
varying gene expression controlled by both transcription factors and
post-transcriptional RNA regulation, and 2) subnets with less varying gene
expression having more feed-forward loops and less post-transcriptional RNA
regulation.


  In this work we demonstrate how the first order phase transition in giant
unilamellar vesicles (GUVs) can function as a trigger for membrane fission.
When driven through their gel-fluid phase transition GUVs exhibit budding or
pearl formation. These buds remain connected to the mother vesicle presumably
by a small neck. Cooling these vesicles from the fluid phase (T>Tm) through the
phase transition into the gel state (T<Tm), leads to complete rupture and
fission of the neck, while the mother vesicle remains intact. Pearling tubes
which formed upon heating break-up and decay into multiple individual vesicles
which then diffuse freely. Finally we demonstrate that mimicking the
intracellular bulk viscosity by increasing the bulk viscosity to 40cP does not
affect the overall fission process, but leads to a significant decrease in size
of the released vesicles.


  This paper is concerned with SIR (susceptible $\to$ infected $\to$ removed)
household epidemic models in which the infection response may be either mild or
severe, with the type of response also affecting the infectiousness of an
individual. Two different models are analysed. In the first model, the
infection status of an individual is predetermined, perhaps due to partial
immunity, and in the second, the infection status of an individual depends on
the infection status of its infector and on whether the individual was infected
by a within- or between-household contact. The first scenario may be modelled
using a multitype household epidemic model, and the second scenario by a model
we denote by the infector-dependent-severity household epidemic model. Large
population results of the two models are derived, with the focus being on the
distribution of the total numbers of mild and severe cases in a typical
household, of any given size, in the event that the epidemic becomes
established. The aim of the paper is to investigate whether it is possible to
determine which of the two underlying explanations is causing the varying
response when given final size household outbreak data containing mild and
severe cases. We conduct numerical studies which show that, given data on
sufficiently many households, it is generally possible to discriminate between
the two models by comparing the Kullback-Leibler divergence for the two fitted
models to these data.


  The inference of gene predictors in the gene regulatory network has become an
important research area in the genomics and medical disciplines. Accurate
predicators are necessary for constructing the GRN model and to enable targeted
biological experiments that attempt to confirm or control the regulation
process. In this paper, we implement a SAT-based algorithm to determine the
gene predictor set from steady state gene expression data (attractor states).
Using the attractor states as input, the states are ordered into attractor
cycles. For each attractor cycle ordering, all possible predictors are
enumerated and a CNF expression is formulated which encodes these predictors
and their biological constraints. Each CNF is explored using a SAT solver to
find candidate predictor sets. Statistical analysis of the results selects the
most likely predictor set of the GRN corresponding to the attractor data. We
demonstrate our algorithm on attractor state data from a melanoma study, and
present our predictor set results.


  Suppose N is a phylogenetic network indicating a complicated relationship
among individuals and taxa. Often of interest is a much simpler network, for
example, a species tree T, that summarizes the most fundamental relationships.
The meaning of a species tree is made more complicated by the recent discovery
of the importance of hybridizations and lateral gene transfers. Hence it is
desirable to describe uniform well-defined procedures that yield a tree given a
network N. A useful tool toward this end is a connected surjective digraph
(CSD) map f from N to N' where N' is generally a much simpler network than N. A
set W of vertices in N is "restricted" if there is at most one vertex from
which there is an arc into W, thus yielding a bottleneck in N. A CSD map f from
N to N' is "restricted" if the inverse image of each vertex in N' is restricted
in N. This paper describes a uniform procedure that, given a network N, yields
a well-defined tree called the "restricted tree" of N. There is a restricted
CSD map from N to the restricted tree. Many relationships in the tree can be
proved to appear also in N.


  The polymerization of actin via branching at a cell membrane containing
nucleation-promoting factors is simulated using a stochastic-growth
methodology. The polymerized-actin distribution displays three types of
behavior: a) traveling waves, b) moving patches, and c) random fluctuations.
Increasing actin concentration causes a transition from patches to waves. The
waves and patches move by a treadmilling mechanism which does not require
myosin II. The effects of downregulation of key proteins on actin wave behavior
are evaluated.


  The accuracy of the recombination estimation method of Padhukasahasram et al.
2006 can be improved by including additional informative summary statistics in
the rejection scheme and by simulating datasets under a fixed segregating sites
model. A C++ program that outputs these summary statistics is freely available
from my website.


  A selective sweep describes the reduction of linked genetic variation due to
strong positive selection. If s is the fitness advantage of a homozygote for
the beneficial allele and h its dominance coefficient, it is usually assumed
that h=1/2, i.e. the beneficial allele is co-dominant. We complement existing
theory for selective sweeps by assuming that h is any value in [0,1]. We show
that genetic diversity patters under selective sweeps with strength s and
dominance 0<h<1 are similar to co-dominant sweeps with selection strength 2hs.
Moreover, we focus on the case h=0 of a completely recessive beneficial allele.
We find that the length of the sweep, i.e. the time from occurrence until
fixation of the beneficial allele, is of the order of sqrt(N/s) generations, if
N is the population size. Simulations as well as our results show that genetic
diversity patterns in the recessive case h=0 greatly differ from all other
cases.


  Species extinction occurs regularly and unavoidably in ecological systems.
The time scales for extinction can broadly vary and inform on the ecosystem's
stability. We study the spatio-temporal extinction dynamics of a paradigmatic
population model where three species exhibit cyclic competition. The cyclic
dynamics reflects the non-equilibrium nature of the species interactions. While
previous work focusses on the coarsening process as a mechanism that drives the
system to extinction, we found that unexpectedly the dynamics to extinction is
much richer. We observed three different types of dynamics. In addition to
coarsening, in the evolutionary relevant limit of large times, oscillating
traveling waves and heteroclinic orbits play a dominant role. The weight of the
different processes depends on the degree of mixing and the system size. By
analytical arguments and extensive numerical simulations we provide the full
characteristics of scenarios leading to extinction in one of the most
surprising models of ecology.


  Pandemics have the potential to cause immense disruption and damage to
communities and societies. In this paper, we propose a hybrid model to
determine how the pandemic spread through the world. The model combines the
SEIR-based model for local areas and the network model for global connection
between countries. We simulate the potential pandemic with different origins
and find how the difference of the origin of a pandemic influences the impact
in the world. We investigate the travelers network which is derived from real
data, and simulate 65 countries, and see how the pandemic spread through the
world from different 14 countries as origins of pandemic. We compare the
difference in terms of the impact in countries and the impact in the world. As
a result, the impact in the world increases when pandemic originates from the
United States, India, and China.


  Pandemics have the potential to cause immense disruption and damage to
communities and societies. In this paper, we model the Influenza Pandemic of
2009. We propose a hybrid model to determine how the pandemic spreads through
the world. The model considers both the SEIR-based model for local areas and
the network model for global connection between countries referring to data on
international travelers. Our interest is to reproduce the situation using the
data of early stage of pandemic and to predict the future transition by
extending the simulation cycle. Without considering the tendency of seasonal
flu, the simulation does not predict the second peak of the pandemic in the
real world. However, considering the seasonal tendency, the simulation result
predicts the next peak in winter. Thus we consider the seasonal tendency is an
important factor for the spreading of the pandemic.


  The Influenza Pandemic of 1918-1919, also called Spanish Flu Pandemic, was
one of the severest pandemics in history. It is thought that the First World
War much influenced the spread of the pandemic. In this paper, we model the
pandemic considering both civil and military traffic. We propose a hybrid model
to determine how the pandemic spread through the world. Our approach considers
both the SEIR-based model for local areas and the network model for global
connection between countries. First, we reproduce the situation in 12
countries. Then, we simulate another scenario: there was no military traffic
during the pandemic, to determine the influence of the influenced of the war on
the pandemic. By considering the simulation results, we find that the influence
of the war varies in countries; in countries which were deeply involved in the
war, the infections were much influenced by the war, while in countries which
were not much engaged in the war, the infections were not influenced by the
war.


  Influenza Pandemic of 1957-1958, also called Asian Flu Pandemic, was one of
the most widespread pandemics in history. In this paper, we model the pandemic,
considering the effect of the Cold War. There were some restrictions between
Western and Eastern nations due to the Cold War during the pandemic. We expect
that such restrictions influenced the spread of the pandemic. We propose a
hybrid model to determine how the pandemic spread through the world. The model
combines the SEIR-based model for local areas and the network model for global
connection between countries. First, we reproduce the situation in 19
countries. Then, we run another experiment to find the influence of the war in
the spread of the pandemic; simulation considering international relationships
in different years. The simulation results show that the impact of the pandemic
in each country was much influenced by international relationships. This study
indicates that if there was less effect of the Cold War, Western nations would
have larger number of death cases, Eastern nations would have smaller number of
death cases, and the world impact would be increased somewhat.


  Thanks to recent technological advances, it is now possible to track with an
unprecedented precision and for long periods of time the movement patterns of
many living organisms in their habitat. The increasing amount of data available
on single trajectories offers the possibility of understanding how animals move
and of testing basic movement models. Random walks have long represented the
main description for micro-organisms and have also been useful to understand
the foraging behaviour of large animals. Nevertheless, most vertebrates, in
particular humans and other primates, rely on sophisticated cognitive tools
such as spatial maps, episodic memory and travel cost discounting. These
properties call for other modeling approaches of mobility patterns. We propose
a foraging framework where a learning mobile agent uses a combination of
memory-based and random steps. We investigate how advantageous it is to use
memory for exploiting resources in heterogeneous and changing environments. An
adequate balance of determinism and random exploration is found to maximize the
foraging efficiency and to generate trajectories with an intricate
spatio-temporal order. Based on this approach, we propose some tools for
analysing the non-random nature of mobility patterns in general.


  The second law of thermodynamics implies that no macroscopic system may
oscillate indefinitely without consuming energy. The question of the number of
possible oscillations and the coherent quality of these oscillations remain
unanswered. This paper proves the upper-bounds on the number and quality of
such oscillations when the system in question is homogeneously driven and has a
discrete network of states. In a closed system, the maximum number of
oscillations is bounded by the number of states in the network. In open
systems, the size of the network bounds the quality factor of oscillation. This
work also explores how the quality factor of macrostate oscillations, such as
would be observed in chemical reactions, are bounded by the smallest equivalent
loop of the network, not the size of the entire system. The consequences of
this limit are explored in the context of chemical clocks and limit cycles.


  The coat proteins of many viruses spontaneously form icosahedral capsids
around nucleic acids or other polymers. Elucidating the role of the packaged
polymer in capsid formation could promote biomedical efforts to block viral
replication and enable use of capsids in nanomaterials applications. To this
end, we perform Brownian dynamics on a coarse-grained model that describes the
dynamics of icosahedral capsid assembly around a flexible polymer. We identify
several mechanisms by which the polymer plays an active role in its
encapsulation, including cooperative polymer-protein motions. These mechanisms
are related to experimentally controllable parameters such as polymer length,
protein concentration, and solution conditions. Furthermore, the simulations
demonstrate that assembly mechanisms are correlated to encapsulation
efficiency, and we present a phase diagram that predicts assembly outcomes as a
function of experimental parameters. We anticipate that our simulation results
will provide a framework for designing in vitro assembly experiments on
single-stranded RNA virus capsids.


  Darwin's book, Origin of the Species has been a source of public controversy
for more than hundred and fifty years. Court cases and mountains of words have
not dispelled this controversy. In this paper, a quantitative approach using
simple mathematics shows that the concept of evolution by natural selection
using only random choice of variables does work. The procedure applied to the
optical equations forming the phenotype of a spider eye produces an eye design
modeled after the measurements by Land and Barth.


  I describe an analytical approximation for calculating the short-term
probability of loss of a chromosome under the neutral Wright-Fisher model with
recombination. I also present an upper and lower bound for this probability.
Exact analytical calculation of this quantity is difficult and computationally
expensive because the number of different ways in which a chromosome can be
lost, grows very large in the presence of recombination. Simulations indicate
that the probabilities obtained using my approximate formula are always
comparable to the true expectations provided that the number of generations
remains small. These results are useful in the context of an algorithm that we
recently developed for simulating Wright-Fisher populations forward in time.
C++ programs that can efficiently calculate these formulas are available on
request.


  Ribosome is a molecular machine that moves on a mRNA track while,
simultaneously, polymerizing a protein using the mRNA also as the corresponding
template. We define, and analytically calculate, two different measures of the
efficiency of this machine. However, we arugue that its performance is
evaluated better in terms of the translational fidelity and the speed with
which it polymerizes a protein. We define both these quantities and calculate
these analytically. Fidelity is a measure of the quality of the products while
the total quantity of products synthesized in a given interval depends on the
speed of polymerization. We show that for synthesizing a large quantity of
proteins, it is not necessary to sacrifice the quality. We also explore the
effects of the quality control mechanism on the strength of mechano-chemical
coupling. We suggest experiments for testing some of the ideas presented here.


  We analyze the time-series data for the onset of A(H1N1)v influenza pandemic
in India during the period June 1- September 30, 2009. Using a variety of
statistical fitting procedures, we obtain a robust estimate of the exponential
growth rate $\langle \lambda \rangle \simeq 0.15$. This corresponds to a basic
reproductive number $R_0 \simeq 1.45$ for influenza A(H1N1)v in India, a value
which lies towards the lower end of the range of values reported for different
countries affected by the pandemic.


  Cells are regulated by networks of controllers having many targets, and
targets affected by many controllers, but these "many-to-many" combinatorial
control systems are poorly understood. Here we analyze distinct cellular
networks (transcription factors, microRNAs, and protein kinases) and a
drug-target network. Certain network properties seem universal across systems
and species, suggesting the existence of common control strategies in biology.
The number of controllers is ~8% of targets and the density of links is 2.5%
\pm 1.2%. Links per node are predominantly exponentially distributed, implying
conservation of the average, which we explain using a mathematical model of
robustness in control networks. These findings suggest that optimal
pharmacological strategies may benefit from a similar, many-to-many
combinatorial structure, and molecular tools are available to test this
approach.


  We use analytical techniques based on an expansion in the inverse system size
to study the stochastic evolutionary dynamics of finite populations of players
interacting in a repeated prisoner's dilemma game. We show that a mechanism of
amplification of demographic noise can give rise to coherent oscillations in
parameter regimes where deterministic descriptions converge to fixed points
with complex eigenvalues. These quasi-cycles between co-operation and defection
have previously been observed in computer simulations; here we provide a
systematic and comprehensive analytical characterization of their properties.
We are able to predict their power spectra as a function of the mutation rate
and other model parameters, and to compare the relative magnitude of the cycles
induced by different types of underlying microscopic dynamics. We also extend
our analysis to the iterated prisoner's dilemma game with a win-stay lose-shift
strategy, appropriate in situations where players are subject to errors of the
trembling-hand type.


  We investigate the effects of heterogeneous and clustered contact patterns on
the timescale and final size of infectious disease epidemics. The abundance of
transitive relationships (the number of 3 cliques) in a network and the
variance of the degree distribution are shown to have large effects on the
number ultimately infected and how quickly the epidemic propagates. The network
model is based on a simple generalization of the configuration model, and
epidemic dynamics are modeled with a low dimensional system of ordinary
differential equations. Because of the simplicity of this model, we are able to
explore a large parameter space and characterize dynamics over a wide range of
network topologies. We find that the interaction between clustering and the
degree distribution is complex, and that clustering always slows down an
epidemic, but that simultaneously increasing clustering and variance of the
degree distribution can potentially increase final epidemic size. In contrast
to solutions for unclustered configuration model networks, we find that bond
percolation solutions for the final epidemic size are potentially biased if
they do not take variable infectious periods into account.


  Inferential summaries of tree estimates are useful in the setting of
evolutionary biology, where phylogenetic trees have been built from DNA data
since the 1960's. In bioinformatics, psychometrics and data mining,
hierarchical clustering techniques output the same mathematical objects, and
practitioners have similar questions about the stability and `generalizability'
of these summaries. This paper provides an implementation of the geometric
distance between trees developed by Billera, Holmes and Vogtmann (2001) [BHV]
equally applicable to phylogenetic trees and hieirarchical clustering trees,
and shows some of the applications in statistical inference for which this
distance can be useful. In particular, since BHV have shown that the space of
trees is negatively curved (a CAT(0) space), a natural representation of a
collection of trees is a tree. We compare this representation to the Euclidean
approximations of treespace made available through Multidimensional Scaling of
the matrix of distances between trees. We also provide applications of the
distances between trees to hierarchical clustering trees constructed from
microarrays. Our method gives a new way of evaluating the influence both of
certain columns (positions, variables or genes) and of certain rows (whether
species, observations or arrays).


  DNA supercoiling plays an important role in a variety of cellular processes.
The torsional stress related with supercoiling may be also involved in gene
regulation through the local structure and dynamics of the double helix. To
check this possibility steady torsional stress was applied to DNA in the course
of all-atom molecular dynamics simulations. It is found that small static
untwisting significantly reduces the torsional persistence length ($l_t$) of
GC-alternating DNA. For the AT-alternating sequence a smaller effect of the
opposite sign is observed. As a result, the measured $l_t$ values are similar
under zero stress, but diverge with untwisting. The effect is traced to
sequence-specific asymmetry of local torsional fluctuations, and it should be
small in long random DNA due to compensation. In contrast, the stiffness of
special short sequences can vary significantly, which gives a simple
possibility of gene regulation via probabilities of strong fluctuations. These
results have important implications for the role of local DNA twisting in
complexes with transcription factors.


  Nanoparticles introduced in living cells are capable of strongly promoting
the aggregation of peptides and proteins. We use here molecular dynamics
simulations to characterise in detail the process by which nanoparticle
surfaces catalyse the self- assembly of peptides into fibrillar structures. The
simulation of a system of hundreds of peptides over the millisecond timescale
enables us to show that the mechanism of aggregation involves a first phase in
which small structurally disordered oligomers assemble onto the nanoparticle
and a second phase in which they evolve into highly ordered beta-sheets as
their size increases.


  The intrinsic property of proteins to form structural motifs such as
alpha-helices and beta-sheets leads to a complex phase behavior in which
proteins can assemble into various types of aggregates including crystals,
liquidlike phases of unfolded or natively folded proteins, and amyloid fibrils.
Here we use a coarse-grained protein model that enables us to perform Monte
Carlo simulations for determining the phase diagram of natively folded
alpha-helical and unfolded beta-sheet forming peptides. The simulations reveal
the existence of various metastable peptide phases. The liquidlike phases are
metastable with respect to the fibrillar phases, and there is a hierarchy of
metastability.


  We study the solutions of Von Neumann's expanding model with reversible
processes for an infinite reaction network. We show that, contrary to the
irreversible case, the solution space need not be convex in contracting phases
(i.e. phases where the concentrations of reagents necessarily decrease over
time). At optimality, this implies that, while multiple dynamical paths of
global contraction exist, optimal expansion is achieved by a unique time
evolution of reaction fluxes. This scenario is investigated in a statistical
mechanics framework by a replica symmetric theory. The transition from a
non-convex to a convex solution space, which turns out to be well described by
a phenomenological order parameter (the fraction of unused reversible
reactions) is analyzed numerically.


  We investigate the presence of structural collective motions on a picosecond
time scale for the heme protein, cytochrome c, as a function of oxidation and
hydration, using terahertz (THz) time-domain spectroscopy and molecular
dynamics simulations. The THz response dramatically increases with oxidation,
with the largest increase for lowest hydrations and highest frequencies. For
both oxidation states the THz response rapidly increases with hydration
saturating above ~25% (g H2O/g protein). Quasi-harmonic vibrational modes and
dipole-dipole correlation functions are calculated from molecular dynamics
trajectories. The collective mode density of states alone reproduces the
measured hydration dependence providing strong evidence of the existence of
these motions. The large oxidation dependence is reproduced only by the
dipole-dipole correlation function, indicating the contrast arises from
diffusive motions consistent with structural changes occurring in the vicinity
of a buried internal water molecule.


  The assumption of linear response of protein molecules to thermal noise or
structural perturbations, such as ligand binding or detachment, is broadly used
in the studies of protein dynamics. Conformational motions in proteins are
traditionally analyzed in terms of normal modes and experimental data on
thermal fluctuations in such macromolecules is also usually interpreted in
terms of the excitation of normal modes. We have chosen two important protein
motors - myosin V and kinesin KIF1A - and performed numerical investigations of
their conformational relaxation properties within the coarse-grained elastic
network approximation. We have found that the linearity assumption is deficient
for ligand-induced conformational motions and can even be violated for
characteristic thermal fluctuations. The deficiency is particularly pronounced
in KIF1A where the normal mode description fails completely in describing
functional mechanochemical motions. These results indicate that important
assumptions of the theory of protein dynamics may need to be reconsidered.
Neither a single normal mode, nor a superposition of such modes yield an
approximation of strongly nonlinear dynamics.


  We present investigations on volume regulation and beading shape transitions
in PC12 neurites conducted using a flow-chamber technique. By disrupting the
cell cytoskeleton with specific drugs we investigate the role of its individual
components in the volume regulation response. We find that microtubule
disruption increases both swelling rate and maximum volume attained, but does
not affect the ability of the neurite to recover its initial volume. In
addition, investigation of axonal beading --also known as pearling
instability-- provides additional clues on the mechanical state of the neurite.
We conclude that the initial swelling phase is mechanically slowed down by
microtubules, while the volume recovery is driven by passive diffusion of
osmolites. Our experiments provide a framework to investigate the role of
cytoskeletal mechanics in volume homeostasis.


  The coalescent revolutionised theoretical population genetics, simplifying,
or making possible for the first time, many analyses, proofs, and derivations,
and offering crucial insights about the way in which the structure of data in
samples from populations depends on the demographic history of the population.
However statistical inference under the coalescent model is extremely
challenging, effectively because no explicit expressions are available for key
sampling probabilities. This led initially to approximation of these
probabilities by ingenious application of modern computationally-intensive
statistical methods. A key breakthrough occurred when Li and Stephens
introduced a different model, similar in spirit to the coalescent, for which
efficient calculations are feasible. In turn, the Li and Stephens model has
changed statistical inference for the wealth of data now available which
documents molecular genetic variation within populations. We briefly review the
coalescent and associated measure-valued diffusions, describe the Li and
Stephens model, and introduce and apply a generalisation of it for inference of
population structure in the presence of linkage disequilibrium.


  The establishment and spreading of biological populations depends crucially
on population growth at low densities. The Allee effect is a problem in those
populations where the per-capita growth rate at low densities is reduced. We
examine stochastic spatial models in which the reproduction rate changes across
a gradient g so that the population undergoes a 2D-percolation transition.
Without the Allee effect, the transition is continuous and the width w of the
hull scales as in conventional (i.e., uncorrelated) gradient percolation,
proportional to g^(-0.57). However, with a strong Allee effect the transition
is first order and w is proportional to g^(-0.26).


  Computational Intelligence (CI) is a sub-branch of Artificial Intelligence
paradigm focusing on the study of adaptive mechanisms to enable or facilitate
intelligent behavior in complex and changing environments. There are several
paradigms of CI [like artificial neural networks, evolutionary computations,
swarm intelligence, artificial immune systems, fuzzy systems and many others],
each of these has its origins in biological systems [biological neural systems,
natural Darwinian evolution, social behavior, immune system, interactions of
organisms with their environment]. Most of those paradigms evolved into
separate machine learning (ML) techniques, where probabilistic methods are used
complementary with CI techniques in order to effectively combine elements of
learning, adaptation, evolution and Fuzzy logic to create heuristic algorithms
that are, in some sense, intelligent. The current trend is to develop consensus
techniques, since no single machine learning algorithms is superior to others
in all possible situations. In order to overcome this problem several
meta-approaches were proposed in ML focusing on the integration of results from
different methods into single prediction. We discuss here the Landau theory for
the nonlinear equation that can describe the adaptive integration of
information acquired from an ensemble of independent learning agents. The
influence of each individual agent on other learners is described similarly to
the social impact theory. The final decision outcome for the consensus system
is calculated using majority rule in the stationary limit, yet the minority
solutions can survive inside the majority population as the complex
intermittent clusters of opposite opinion.


  Under favourable conditions, many proteins can assemble into macroscopically
large aggregate's, Parkinson's and other neurological and systemic diseases.
The overall process of protein aggregation is characterized by initial lag time
during which no detectable aggregation occurs in the solution and by maximal
aggregation rate at which the dissolved protein converts into aggregates. In
this study, the correlation between the lag time and the maximal rate of
protein aggregation is analyzed. It is found that the product of these two
quantities depends on a single numerical parameter, the kinetic index of the
curve quantifying the time evolution of the fraction of protein aggregated. As
this index depends relatively little on the conditions and/or system studied,
our finding provides insight into why for many experiments the values of the
product of the lag time and the maximal aggregation rate are often equal or
quite close to each other. It is shown how the kinetic index is related to a
basic kinetic parameter of a recently proposed theory of protein aggregation.


  We consider nucleation of amyloid fibrils in the case when the process occurs
by the mechanism of direct polymerization of practically fully extended protein
segments, i.e. beta-strands, into beta-sheets. Applying the classical
nucleation theory, we derive a general expression for the work to form a
nanosized amyloid fibril (protofilament) constituted of successively layered
beta-sheets. Analysis of this expression reveals that with increasing its size,
the fibril transforms from one-dimensional into two-dimensional aggregate in
order to preserve the equilibrium shape corresponding to minimal formation
work. We determine the size of the fibril nucleus, the fibril nucleation work
and the fibril nucleation rate as explicit functions of the concentration and
temperature of the protein solution. The results obtained are applicable to
homogeneous nucleation which occurs when the solution is sufficiently pure
and/or strongly supersaturated.


  This study presents a method to discover an outbreak of an infectious disease
in a region for which data are missing, but which is at work as a disease
spreader. Node discovery for the spread of an infectious disease is defined as
discriminating between the nodes which are neighboring to a missing disease
spreader node, and the rest, given a dataset on the number of cases. The spread
is described by stochastic differential equations. A perturbation theory
quantifies the impact of the missing spreader on the moments of the number of
cases. Statistical discriminators examine the mid-body or tail-ends of the
probability density function, and search for the disturbance from the missing
spreader. They are tested with computationally synthesized datasets, and
applied to the SARS outbreak and flu pandemic.


  The quite recent technological rise in molecular biology allowed single
molecule manipulation experiments, where molecule stretching plays a primary
role. In order to understand the experimental data, it is felt the urge of some
physical and mathematical models to quantitatively express the mechanical
properties of the observed molecules. In this paper we reconsider a simple
phenomenological model which reproduces the behaviour of single and double
stranded DNA under tension. The problem is easily solved via the cavity method
both in the small forces range and in presence of overstretching transition, so
that some properties such as bending stiffness and elasticity of DNA emerge in
a very clear manner. Our theoretical findings are successfully fitted to real
measurements and compared to Monte Carlo simulations, confirming the quality of
the approach.


  A replicator equation with explicit space and global regulation is
considered. This model provides a natural framework to follow frequencies of
species that are distributed in the space. For this model, analogues to
classical notions of the Nash equilibrium and evolutionary stable strategy are
provided. A sufficient condition for a uniform stationary state to be a
spatially distributed evolutionary stable state is presented and illustrated
with examples.


  We live in a time where climate models predict future increases in
environmental variability and biological invasions are becoming increasingly
frequent. A key to developing effective responses to biological invasions in
increasingly variable environments will be estimates of their rates of spatial
spread and the associated uncertainty of these estimates. Using stochastic,
stage-structured, integro-difference equation models, we show analytically that
invasion speeds are asymptotically normally distributed with a variance that
decreases in time. We apply our methods to a simple juvenile-adult model with
stochastic variation in reproduction and an illustrative example with published
data for the perennial herb, \emph{Calathea ovandensis}. These examples
buttressed by additional analysis reveal that increased variability in vital
rates simultaneously slow down invasions yet generate greater uncertainty about
rates of spatial spread. Moreover, while temporal autocorrelations in vital
rates inflate variability in invasion speeds, the effect of these
autocorrelations on the average invasion speed can be positive or negative
depending on life history traits and how well vital rates ``remember'' the
past.


  In this paper we provide some sharp asymptotic results for a stochastic model
of species survival recently proposed by Guiol, Marchado, and Schinazi.


  In this work we study the process of mutant invasion on an example of a
consumer-resource system with annual character of the behavior. Namely,
individuals are active during seasons of fixed length separated by winter
periods. All individuals die at the end of the season and the size of the next
generation is determined by the number of offspring produced during the past
season. The rate at which the consumers produce immature offspring depends on
their internal energy which can be increased by feeding. The reproduction of
the resource simply occurs at a constant rate. At the beginning, we consider a
population of consumers maximizing their common fitness, all consumers being
individuals having the same goal function and acting for the common good. We
suppose that a small fraction of the consumer population may appear at the
beginning of one season and start to behave as mutants in the main population.
We study how such invasion occurs.


  We consider a population of identical individuals preying on an exhaustible
resource. The individuals in the population choose a strategy that defines how
they use their available time over the course of their life for feeding, for
reproducing (say laying eggs), or split their energy in between these two
activities. We here suppose that their life lasts a full season, so that the
chosen strategy results in the production of a certain number of eggs laid over
the season. This number then helps define the long-run evolution of the
population since the eggs constitute the basis for the population for the next
season. However, in this paper, we strictly concentrate on what occurs within a
season, by considering two possible strategies: the collective optimum and the
uninvadable strategy. The (collective) optimal strategy involves a singular
arc, which has some resource saving feature and may, on the long-term, lead to
an equilibrium. However, it is susceptible to invasion by a greedy mutant which
free-rides on this resource saving strategy. Therefore, the "optimal" strategy
is not evolutionarily stable. We thus look for an evolutionarily stable
strategy, using the first ESS condition ("Nash-Wardrop condition"). This yields
a different singular arc and a strategy that conserves less the resource.
Unfortunately, we show in a forthcoming paper that, though it is an ESS, it is
not long-term (interseasonal) sustainable. This is an instance of the classical
"tragedy of the commons". In this paper we also investigate whether the mixed
strategy on the singular arcs may be interpreted in terms of population shares
using pure strategies. The answer is negative.


  The advent of high-throughput sequencing technologies constituted a major
advance in genomic studies, offering new prospects in a wide range of
applications. We propose a rigorous and flexible algorithmic solution to
mapping SOLiD color-space reads to a reference genome. The solution relies on
an advanced method of seed design that uses a faithful probabilistic model of
read matches and, on the other hand, a novel seeding principle especially
adapted to read mapping. Our method can handle both lossy and lossless
frameworks and is able to distinguish, at the level of seed design, between
SNPs and reading errors. We illustrate our approach by several seed designs and
demonstrate their efficiency.


  Autocatalysis is a fundamental concept, used in a wide range of domains. From
the most general definition of autocatalysis, that is a process in which a
chemical compound is able to catalyze its own formation, several different
systems can be described. We detail the different categories of autocatalyses,
and compare them on the basis of their mechanistic, kinetic, and dynamic
properties. It is shown how autocatalytic patterns can be generated by
different systems of chemical reactions. The notion of autocatalysis covering a
large variety of mechanistic realisations with very similar behaviors, it is
proposed that the key signature of autocatalysis is its kinetic pattern
expressed in a mathematical form.


  We investigate the dynamical behavior of simple modules composed of two genes
with two or three regulating connections. Continuous dynamics for mRNA and
protein concentrations is compared to a Boolean model for gene activity. Using
a generalized method, we study within a single framework different continuous
models and different types of regulatory functions, and establish conditions
under which the system can display stable oscillations. These conditions
concern the time scales, the degree of cooperativity of the regulating
interactions, and the signs of the interactions. Not all models that show
oscillations under Boolean dynamics can have oscillations under continuous
dynamics, and vice versa.


  Cellular functions are based on the complex interplay of proteins, therefore
the structure and dynamics of these protein-protein interaction (PPI) networks
are the key to the functional understanding of cells. In the last years,
large-scale PPI networks of several model organisms were investigated.
Methodological improvements now allow the analysis of PPI networks of multiple
organisms simultaneously as well as the direct modeling of ancestral networks.
This provides the opportunity to challenge existing assumptions on network
evolution. We utilized present-day PPI networks from integrated datasets of
seven model organisms and developed a theoretical and bioinformatic framework
for studying the evolutionary dynamics of PPI networks. A novel filtering
approach using percolation analysis was developed to remove low confidence
interactions based on topological constraints. We then reconstructed the
ancient PPI networks of different ancestors, for which the ancestral proteomes,
as well as the ancestral interactions, were inferred. Ancestral proteins were
reconstructed using orthologous groups on different evolutionary levels. A
stochastic approach, using the duplication-divergence model, was developed for
estimating the probabilities of ancient interactions from today's PPI networks.
The growth rates for nodes, edges, sizes and modularities of the networks
indicate multiplicative growth and are consistent with the results from
independent static analysis. Our results support the duplication-divergence
model of evolution and indicate fractality and multiplicative growth as general
properties of the PPI network structure and dynamics.


  Predicting protein structure from amino acid sequence is one of the most
important unsolved problems of molecular biology and biophysics.Not only would
a successful prediction algorithm be a tremendous advance in the understanding
of the biochemical mechanisms of proteins, but, since such an algorithm could
conceivably be used to design proteins to carry out specific
functions.Prediction of the secondary structure of a protein (alpha-helix,
beta-sheet, coil) is an important step towards elucidating its three
dimensional structure as well as its function. In this research, we use
different Hidden Markov models for protein secondary structure prediction. In
this paper we have proposed an algorithm for predicting protein secondary
structure. We have used Hidden Markov model with sliding window for secondary
structure prediction.The secondary structure has three regular forms, for each
secondary structural element we are using one Hidden Markov Model.


  We study implications of the dynamical and spatial contact structure between
Brazilian escorts and sex-buyers for the spreading of sexually transmitted
infections (STI). Despite a highly skewed degree distribution diseases
spreading in this contact structure have rather well-defined epidemic
thresholds. Temporal effects create a broad distribution of outbreak sizes even
if the transmission probability is taken to the hypothetical value of 100%.
Temporal correlations speed up outbreaks, especially in the early phase,
compared to randomized contact structures. The time-ordering and the network
topology, on the other hand, slow down the epidemics. Studying compartmental
models we show that the contact structure can probably not support the spread
of HIV, not even if individuals were sexually active during the acute
infection. We investigate hypothetical means of containing an outbreak and find
that travel restrictions are about as efficient as removal of the vertices of
highest degree. In general, the type of commercial sex we study seems not like
a major factor in STI epidemics.


  The evolution of cooperation has been a perennial problem for evolutionary
biology because cooperation is undermined by selfish cheaters (or "free
riders") that profit from cooperators but do not invest any resources
themselves. In a purely "selfish" view of evolution, those cheaters should be
favored. Evolutionary game theory has been able to show that under certain
conditions, cooperation nonetheless evolves stably. One of these scenarios
utilizes the power of punishment to suppress free riders, but only if players
interact in a structured population where cooperators are likely to be
surrounded by other cooperators. Here we show that cooperation via punishment
can evolve even in well-mixed populations that play the "public goods" game, if
the synergy effect of cooperation is high enough. As the synergy is increased,
populations transition from defection to cooperation in a manner reminiscent of
a phase transition. If punishment is turned off, the critical synergy is
significantly higher, illustrating that (as shown before) punishment aids in
establishing cooperation. We also show that the critical point depends on the
mutation rate so that higher mutation rates discourage cooperation, as has been
observed before in the Prisoner's Dilemma.


  Evolutionary adaptation is the process that increases the fit of a population
to the fitness landscape it inhabits. As a consequence, evolutionary dynamics
is shaped, constrained, and channeled, by that fitness landscape. Much work has
been expended to understand the evolutionary dynamics of adapting populations,
but much less is known about the structure of the landscapes. Here, we study
the global and local structure of complex fitness landscapes of interacting
loci that describe protein folds or sets of interacting genes forming pathways
or modules. We find that in these landscapes, high peaks are more likely to be
found near other high peaks, corroborating Kauffman's "Massif Central"
hypothesis. We study the clusters of peaks as a function of the ruggedness of
the landscape and find that this clustering allows peaks to form interconnected
networks. These networks undergo a percolation phase transition as a function
of minimum peak height, which indicates that evolutionary trajectories that
take no more than two mutations to shift from peak to peak can span the entire
genetic space. These networks have implications for evolution in rugged
landscapes, allowing adaptation to proceed after a local fitness peak has been
ascended.


  The four-year oscillations of the number of spawning sockeye salmon
(Oncorhynchus nerka) that return to their native stream within the Fraser River
basin in Canada are a striking example of population oscillations. The period
of the oscillation corresponds to the dominant generation time of these fish.
Various - not fully convincing - explanations for these oscillations have been
proposed, including stochastic influences, depensatory fishing, or genetic
effects. Here, we show that the oscillations can be explained as a stable
dynamical attractor of the population dynamics, resulting from a strong
resonance near a Neimark Sacker bifurcation. This explains not only the
long-term persistence of these oscillations, but also reproduces correctly the
empirical sequence of salmon abundance within one period of the oscillations.
Furthermore, it explains the observation that these oscillations occur only in
sockeye stocks originating from large oligotrophic lakes, and that they are
usually not observed in salmon species that have a longer generation time.


  RNA-RNA binding is an important phenomenon observed for many classes of
non-coding RNAs and plays a crucial role in a number of regulatory processes.
Recently several MFE folding algorithms for predicting the joint structure of
two interacting RNA molecules have been proposed. Here joint structure means
that in a diagram representation the intramolecular bonds of each partner are
pseudoknot-free, that the intermolecular binding pairs are noncrossing, and
that there is no so-called ``zig-zag'' configuration. This paper presents the
combinatorics of RNA interaction structures including their generating
function, singularity analysis as well as explicit recurrence relations. In
particular, our results imply simple asymptotic formulas for the number of
joint structures.


  The evolution and maintenance of cooperation fascinated researchers for
several decades. Recently, theoretical models and experimental evidence show
that costly punishment may facilitate cooperation in human societies, but may
not be used by winners. The puzzle how the costly punishment behaviour evolves
can be solved under voluntary participation. Could the punishers emerge if
participation is compulsory? Is the punishment inevitably a selfish behaviour
or an altruistic behaviour? The motivations behind punishment are still an
enigma. Based on public goods interactions, we present a model in which just a
certain portion of the public good is divided equally among all members. The
other portion is distributed to contributors when paying a second cost.
Contributors who are willing to pay a second cost can be costly (and then
altruistic) punishers, but they can also flourish or dominate the population,
in this case we may call them "advantageous punishers". We argue that most of
successful cooperators and punishers in nature are advantageous punishers, and
costly punishment mostly happens in humans. This indicates a universal
surviving rule: contributing more and gaining more. Our models show
theoretically that the original motivation behind punishment is to retrieve
deserved payoff from their own contributions, a selfish incentive.


  Evolutionary game dynamics in finite populations is typically subject to
noise, inducing effects which are not present in deterministic systems,
including fixation and extinction. In the first part of this paper we
investigate the phenomenon of drift reversal in finite populations, taking into
account that drift is a local quantity in strategy space. Secondly, we study a
simple imitation dynamics, and show that it can lead to fixation at internal
mixed-strategy fixed points even in finite populations. Imitation in infinite
populations is adequately described by conventional replicator dynamics, and
these equations are known to have internal fixed points. Internal absorption in
finite populations on the other hand is a novel dynamic phenomenon. Due to an
outward drift in finite populations this type of dynamic arrest is not found in
other commonly studied microscopic dynamics, not even in those with the same
deterministic replicator limit as imitation.


  Phylogenetic trees and networks are leaf-labelled graphs that are used to
describe evolutionary histories of species. The Tree Containment problem asks
whether a given phylogenetic tree is embedded in a given phylogenetic network.
Given a phylogenetic network and a cluster of species, the Cluster Containment
problem asks whether the given cluster is a cluster of some phylogenetic tree
embedded in the network. Both problems are known to be NP-complete in general.
In this article, we consider the restriction of these problems to several
well-studied classes of phylogenetic networks. We show that Tree Containment is
polynomial-time solvable for normal networks, for binary tree-child networks,
and for level-$k$ networks. On the other hand, we show that, even for
tree-sibling, time-consistent, regular networks, both Tree Containment and
Cluster Containment remain NP-complete.


  Karlin's (1982) Theorem 5.2 shows that linear systems alternating between
growth and mixing phases have lower asymptotic growth with greater mixing. Here
this result is extended to linear differential equations that combine
site-specific growth or decay rates, and mixing between sites, showing that the
spectral abscissa of a matrix D + m A decreases with m, where D does-not-equal
c I is a real diagonal matrix, A is an irreducible matrix with non-negative
off-diagonal elements (an ML- or essentially non-negative matrix), and m >= 0.
The result is based on the inequality: u' A v < r(A), where u and v are the
left and right Perron vectors of the matrix D + A, and r(A) is the spectral
abscissa and Perron root of A. The result gives an analytic solution to prior
work that relied on two-site or numerical simulation of models of growth and
mixing, such as source and sink ecological models, or multiple tissue
compartment models of microbe growth. The result has applications to the
Lyapunov stability of perturbations in nonlinear systems.


  Background: Transposable elements (TEs) have played an important role in the
diversification and enrichment of mammalian transcriptomes through various
mechanisms such as exonization and intronization (the birth of new
exons/introns from previously intronic/exonic sequences, respectively), and
insertion into first and last exons. However, no extensive analysis has
compared the effects of TEs on the transcriptomes of mammalian, non-mammalian
vertebrates and invertebrates. Results: We analyzed the influence of TEs on the
transcriptomes of five species, three invertebrates and two non-mammalian
vertebrates. Compared to previously analyzed mammals, there were lower levels
of TE introduction into introns, significantly lower numbers of exonizations
originating from TEs and a lower percentage of TE insertion within the first
and last exons. Although the transcriptomes of vertebrates exhibit a
significant level of exonizations of TEs, only anecdotal cases were found in
invertebrates. In vertebrates, as in mammals, the exonized TEs are mostly
alternatively spliced, indicating selective pressure maintains the original
mRNA product generated from such genes. Conclusions: Exonization of TEs is
wide-spread in mammals, less so in non- mammalian vertebrates, and very low in
invertebrates. We assume that the exonization process depends on the length of
introns. Vertebrates, unlike invertebrates, are characterized by long introns
and short internal exons. Our results suggest that there is a direct link
between the length of introns and exonization of TEs and that this process
became more prevalent following the appearance of mammals.


  Insertion of transposed elements within mammalian genes is thought to be an
important contributor to mammalian evolution and speciation. Insertion of
transposed elements into introns can lead to their activation as alternatively
spliced cassette exons, an event called exonization. Elucidation of the
evolutionary constraints that have shaped fixation of transposed elements
within human and mouse protein coding genes and subsequent exonization is
important for understanding of how the exonization process has affected
transcriptome and proteome complexities. Here we show that exonization of
transposed elements is biased towards the beginning of the coding sequence in
both human and mouse genes. Analysis of single nucleotide polymorphisms (SNPs)
revealed that exonization of transposed elements can be population-specific,
implying that exonizations may enhance divergence and lead to speciation. SNP
density analysis revealed differences between Alu and other transposed
elements. Finally, we identified cases of primate-specific Alu elements that
depend on RNA editing for their exonization. These results shed light on TE
fixation and the exonization process within human and mouse genes.


  We study the mechanical unfolding pathways of the $FnIII_{10}$ domain of
fibronectin by means of an Ising--like model, using both constant force and
constant velocity protocols. At high forces and high velocities our results are
consistent with experiments and previous computational studies. Moreover, the
simplicity of the model allows us to probe the biologically relevant low force
regime, where we predict the existence of two intermediates with very close
elongations. The unfolding pathway is characterized by stochastic transitions
between these two intermediates.


  The Natural Immune System (NIS) is a distributed system that solves
challenging search and response problems while operating under constraints
imposed by physical space and resource availability. Remarkably, NIS search and
response times do not scale appreciably with the physical size of the animal in
which its search is conducted. Many distributed systems are engineered to solve
analogous problems, and the NIS demonstrates how such engineered systems can
achieve desirable scalability. We hypothesize that the architecture of the NIS,
composed of a hierarchical decentralized detection network of lymph nodes (LN)
facilitates efficient search and response. A sub-modular architecture in which
LN numbers and size both scale with organism size is shown to efficiently
balance tradeoffs between local antigen detection and global antibody
production, leading to nearly scale-invariant detection and response. We
characterize the tradeoffs as balancing local and global communication and show
that similar tradeoffs exist in distributed systems like LN inspired artificial
immune system (AIS) applications and peer-to-peer (P2P) systems. Taking
inspiration from the architecture of the NIS, we propose a modular RADAR
(Robust Adaptive Decentralized search with Automated Response) strategy for
distributed systems. We demonstrate how two existing distributed systems (a LN
inspired multi-robot control application and a P2P system) can be improved by a
modular RADAR strategy. Such a sub-modular architecture is shown to balance the
tradeoffs between local communication (within artificial LNs and P2P clusters)
and global communication (between artificial LNs and P2P clusters), leading to
efficient search and response.


  The last few years have seen remarkably fast progress in the understanding of
statistics and epidemic dynamics of various clustered networks. This paper
considers a class of networks based around a concept (the locale) that allows
asymptotically exact results to be derived for epidemic dynamics. While there
is no restriction on the motifs that can be found in such graphs, each node
must be uniquely assigned to a generally clustered subgraph to obtain analytic
traction.


  We define catalytic networks as chemical reaction networks with an
essentially catalytic reaction pathway: one which is on in the presence of
certain catalysts and off in their absence. We show that examples of catalytic
networks include synthetic DNA molecular circuits that have been shown to
perform signal amplification and molecular logic. Recall that a critical siphon
is a subset of the species in a chemical reaction network whose absence is
forward invariant and stoichiometrically compatible with a positive point. Our
main theorem is that all weakly-reversible networks with critical siphons are
catalytic. Consequently, we obtain new proofs for the persistence of atomic
event-systems of Adleman et al., and normal networks of Gnacadja. We define
autocatalytic networks, and conjecture that a weakly-reversible reaction
network has critical siphons if and only if it is autocatalytic.


  We study the effect of quenched randomness in the arc-length dependent
spontaneous curvature of a wormlike chain under tension. In the weakly bending
approximation in two dimensions, we obtain analytic results for the
force-elongation curve and the width of transverse fluctuations. We compare
quenched and annealed disorder and conclude that the former cannot always be
reduced to a simple change in the stiffness of the pure system. We also discuss
the effect of a random transverse force on the stretching response of a
wormlike chain without spontaneous curvature.


  Variation in high energy cosmic rays (HECRs) has been proposed to explain a
62 My periodicity in terrestrial fossil biodiversity. It has been suggested
that the infall of our galaxy toward the Virgo cluster could generate an
extragalactic shock, accelerating charged particles and exposing the earth to a
flux of high energy cosmic rays (HECRs). The oscillation of the Sun
perpendicular to the galactic plane could induce 62 My periodicity in the HECR
flux on the Earth, with a magnitude much higher than the Galactic cosmic ray
change we see in a solar cycle. This mechanism could potentially explain the
observed 62 My periodicity in terrestrial biodiversity over the past 500 My. In
addition to direct effects on life from secondaries, HECRs induced air showers
ionize the atmosphere leading to changes in atmospheric chemistry and
microphysical processes that can lead to cloud formation including low altitude
cloud cover. An increase in ionization changes the global electric circuit
which could enhance the formation of cloud condensation nuclei (CCN) through
microphysical processes such as electroscavenging and ion mediated nucleation,
leading to an increase in the cloud cover. This could increase the albedo and
reduce the solar flux reaching the ground, reducing the global temperature.
Using an existing model, we have calculated the enhancement in atmospheric
ionization at low altitudes resulting from exposure to HECRs. We use a
conservative model to estimate the change in low altitude cloud cover from this
increased ionization.


  In Templeton (2010), the Approximate Bayesian Computation (ABC) algorithm
(see, e.g., Pritchard et al., 1999, Beaumont et al., 2002, Marjoram et al.,
2003, Ratmann et al., 2009) is criticised on mathematical and logical grounds:
"the [Bayesian] inference is mathematically incorrect and formally illogical".
Since those criticisms turn out to be bearing on Bayesian foundations rather
than on the computational methodology they are primarily directed at, we
endeavour to point out in this note the statistical errors and inconsistencies
in Templeton (2010), refering to Beaumont et al. (2010) for a reply that is
broader in scope since it also covers the phylogenetic aspects of nested clade
versus a model-based approach.


  The network of 5823 cities of Mexico with a population more than 5000
inhabitants is studied. Our analysis is focused to the spectral properties of
the adjacency matrix, the small-world properties of the network, the
distribution of the clustering coefficients and the degree distribution of the
vertices. The connection of these features with the spread of epidemics on this
network is also discussed.


  We have developed a two dimensional stochastic molecular dynamics model for
the description of intra cellular collective motion of bio motors, in
particular Kinesins, on a microtubular track. The model is capable or
reproducing the hand-over-hand mechanism of the directed motion along the
microtubule. The model gives the average directed velocity and the current of
Kinesins along the microtubule. It is shown that beyond a certain density of
Kinesins, the average velocity and current undergo notable decrease which is
due to formation of traffic jams in the system.


  Recently several minimum free energy (MFE) folding algorithms for predicting
the joint structure of two interacting RNA molecules have been proposed. Their
folding targets are interaction structures, that can be represented as diagrams
with two backbones drawn horizontally on top of each other such that (1)
intramolecular and intermolecular bonds are noncrossing and (2) there is no
"zig-zag" configuration. This paper studies joint structures with arc-length at
least four in which both, interior and exterior stack-lengths are at least two
(no isolated arcs). The key idea in this paper is to consider a new type of
shape, based on which joint structures can be derived via symbolic enumeration.
Our results imply simple asymptotic formulas for the number of joint structures
with surprisingly small exponential growth rates. They are of interest in the
context of designing prediction algorithms for RNA-RNA interactions.


  Pre-calculated libraries of molecular fragment configurations have previously
been used as a basis for both equilibrium sampling (via "library-based Monte
Carlo") and for obtaining absolute free energies using a polymer-growth
formalism. Here, we combine the two approaches to extend the size of systems
for which free energies can be calculated. We study a series of all-atom
poly-alanine systems in a simple dielectric "solvent" and find that precise
free energies can be obtained rapidly. For instance, for 12 residues, less than
an hour of single-processor is required. The combined approach is formally
equivalent to the "annealed importance sampling" algorithm; instead of
annealing by decreasing temperature, however, interactions among fragments are
gradually added as the molecule is "grown." We discuss implications for future
binding affinity calculations in which a ligand is grown into a binding site.


  The intrinsic stochasticity of gene expression can lead to large variability
of protein levels across a population of cells. Variability (or noise) in
protein distributions can be modulated by cellular mechanisms of gene
regulation; in particular, there is considerable interest in understanding the
role of post-transcriptional regulation. To address this issue, we propose and
analyze a stochastic model for post-transcriptional regulation of gene
expression. The analytical solution of the model provides insight into the
effects of different mechanisms of post-transcriptional regulation on the noise
in protein distributions. The results obtained also demonstrate how different
sources of intrinsic noise in gene expression can be discriminated based on
observations of regulated protein distributions.


  This paper proposed a new method to build the large scale DNA sequences
search system based on web search engine technology. We give a very brief
introduction for the methods used in search engine first. Then how to build a
DNA search system like Google is illustrated in detail. Since there is no local
alignment process, this system is able to provide the ms level search services
for billions of DNA sequences in a typical server.


  A series of simulations aimed at elucidating the self-assembly dynamics of
spherical virus capsids is described. This little-understood phenomenon is a
fascinating example of the complex processes that occur in the simplest of
organisms. The fact that different viruses adopt similar structural forms is an
indication of a common underlying design, motivating the use of simplified,
low-resolution models in exploring the assembly process. Several versions of a
molecular dynamics approach are described. Polyhedral shells of different sizes
are involved, the assembly pathways are either irreversible or reversible, and
an explicit solvent is optionally included. Model design, simulation
methodology and analysis techniques are discussed. The analysis focuses on the
growth pathways and the nature of the intermediate states, properties that are
hard to access experimentally. Among the key observations are that efficient
growth proceeds by means of a cascade of highly reversible stages, and that
while there are a large variety of possible partial assemblies, only a
relatively small number of strongly bonded configurations are actually
encountered.


  MYH9 has been proposed as a major genetic risk locus for a spectrum of
non-diabetic end stage kidney disease (ESKD). We use recently released
sequences from the 1000 Genomes Project to identify two western African
specific missense mutations (S342G and I384M) in the neighbouring APOL1 gene,
and demonstrate that these are more strongly associated with ESKD than
previously reported MYH9 variants. We also show that the distribution of these
risk variants in African populations is consistent with the pattern of African
ancestry ESKD risk previously attributed to the MYH9 gene. Additional
associations were also found among other members of the APOL gene family, and
we propose that ESKD risk is caused by western African variants in members of
the APOL gene family, which evolved to confer protection against pathogens,
such as Trypanosoma.


  These are notes for a set of 7 two-hour lectures given at the 2010 Summer
School on Quantitative Evolutionary and Comparative Genomics at OIST, Okinawa,
Japan. The emphasis is on understanding how biological systems process
information. We take a physicist's approach of looking for simple
phenomenological descriptions that can address the questions of biological
function without necessarily modeling all (mostly unknown) microscopic details;
the example that is developed throughout the notes is transcriptional
regulation in genetic regulatory networks. We present tools from information
theory and statistical physics that can be used to analyze noisy nonlinear
biological networks, and build generative and predictive models of regulatory
processes.


  We present an application of optimal control theory to Dengue epidemics. This
epidemiologic disease is an important theme in tropical countries due to the
growing number of infected individuals. The dynamic model is described by a set
of nonlinear ordinary differential equations, that depend on the dynamic of the
Dengue mosquito, the number of infected individuals, and the people's
motivation to combat the mosquito. The cost functional depends not only on the
costs of medical treatment of the infected people but also on the costs related
to educational and sanitary campaigns. Two approaches to solve the problem are
considered: one using optimal control theory, another one by discretizing first
the problem and then solving it with nonlinear programming. The results
obtained with OC-ODE and IPOPT solvers are given and discussed. We observe that
with current computational tools it is easy to obtain, in an efficient way,
better solutions to Dengue problems, leading to a decrease of infected
mosquitoes and individuals in less time and with lower costs.


  Although an exchange of genetic information by recombination plays an
important role in the evolution of viruses, it is not clear how it generates
diversity. {\it Geminiviruses} are plant viruses which have ambisense
single-stranded circular DNA genomes and one of the most economically important
plant viruses in agricultural production. Small circular single-stranded DNA
satellites, termed DNA-$\beta$, have recently been found associated with some
geminivirus infections. In this paper we analyze a satellite molecule
DNA-$\beta$ of geminiviruses for recombination events using phylogenetic and
statistical analysis and we find that one strain from ToLCMaB has a
recombination pattern and is possibly recombinant molecule between two strains
from two species, PaLCuB-[IN:Chi:05] (major parent) and ToLCB-[IN:CP:04] (minor
parent).


  We numerically study the phase diagram of a three-state host-parasite model
on the square lattice motivated by population biology. The model is an
extension of the contact process, and the three states correspond to an empty
site, a host, and a parasite. We determine the phase diagram of the model by
scaling analysis. In agreement with previous results, three phases are
identified: the phase in which both hosts and parasites are extinct (S_{0}),
the phase in which hosts survive but parasites are extinct (S_{01}), and the
phase in which both hosts and parasites survive (S_{012}). We argue that both
the S_{0}-S_{01} and S_{01}-S_{012} boundaries belong to the directed
percolation class. In this model, it has been suggested that an excessively
large reproduction rate of parasites paradoxically extinguishes hosts and
parasites and results in S_{0}. We show that this paradoxical extinction is a
finite size effect; the corresponding parameter region is likely to disappear
in the limit of infinite system size.


  We adapt existing statistical modeling techniques for social networks to
study consumption data observed in trophic food webs. These data describe the
feeding volume (non-negative) among organisms grouped into nodes, called
trophic species, that form the food web. Model complexity arises due to the
extensive amount of zeros in the data, as each node in the web is predator/prey
to only a small number of other trophic species. Many of the zeros are regarded
as structural (non-random) in the context of feeding behavior. The presence of
basal prey and top predator nodes (those who never consume and those who are
never consumed, with probability 1) creates additional complexity to the
statistical modeling. We develop a special statistical social network model to
account for such network features. The model is applied to two empirical food
webs; focus is on the web for which the population size of seals is of concern
to various commercial fisheries.


  A parameterisation of generalised network clustering, in the form of
four-motif prevalences, is presented. This involves three real parameters that
are conditional on one- two- and three-motif prevalences. Interpretations of
these real parameters are presented that motivate a set of rewiring schemes to
create appropriately clustered networks. Finally, the dynamical implications of
higher order structure, as parameterised, for a contact process are considered.


  The prevalence of common chronic non-communicable diseases (CNCDs) far
overshadows the prevalence of both monogenic and infectious diseases combined.
All CNCDs, also called complex genetic diseases, have a heritable genetic
component that can be used for pre-symptomatic risk assessment. Common single
nucleotide polymorphisms (SNPs) that tag risk haplotypes across the genome
currently account for a non-trivial portion of the germ-line genetic risk and
we will likely continue to identify the remaining missing heritability in the
form of rare variants, copy number variants and epigenetic modifications. Here,
we describe a novel measure for calculating the lifetime risk of a disease,
called the genetic composite index (GCI), and demonstrate its predictive value
as a clinical classifier. The GCI only considers summary statistics of the
effects of genetic variation and hence does not require the results of
large-scale studies simultaneously assessing multiple risk factors. Combining
GCI scores with environmental risk information provides an additional tool for
clinical decision-making. The GCI can be populated with heritable risk
information of any type, and thus represents a framework for CNCD
pre-symptomatic risk assessment that can be populated as additional risk
information is identified through next-generation technologies.


  We present a theoretical model of facilitated diffusion of proteins in the
cell nucleus. This model, which takes into account the successive
binding/unbinding events of proteins to DNA, relies on a fractal description of
the chromatin which has been recently evidenced experimentally. Facilitated
diffusion is shown quantitatively to be favorable for a fast localization of a
target locus by a transcription factor, and even to enable the minimization of
the search time by tuning the affinity of the transcription factor with DNA.
This study shows the robustness of the facilitated diffusion mechanism, invoked
so far only for linear conformations of DNA.


  We perform a quantum mechanical study of the peptides that are part of the
LH2 complex from Rhodopseudomonas acidophila, a non-sulfur purple bacteria that
has the ability of producing chemical energy from photosynthesis. The
electronic structure calculations indicate that the transmembrane helices of
these peptides are characterized by dipole moments with a magnitude of ~150 D.
When the full nonamer assembly made of eighteen peptides is considered, then a
macrodipole of magnitude 704 D is built up from the vector sum of each monomer
dipole. The macrodipole is oriented normal to the membrane plane and with the
positive tip toward the cytoplasm thereby indicating that the electronic charge
of the protein scaffold is polarized toward the periplasm. The results obtained
here suggest that the asymmetric charge distribution of the protein scaffold
contributes an anisotropic electrostatic environment which differentiates the
absorption properties of the bacteriochlorophyll pigments, B800 and B850,
embedded in the LH2 complex.


  During outbreaks of serious infectious diseases many individuals closely
follow media reports of the outbreak and take steps, including self-isolation,
to protect themselves from infection and possibly death. Self-isolation can
take many forms including restricting local and long-distance travel, using
face masks, and choosing to receive a vaccine. We use mathematical modeling to
show that public health agencies working together with the media can
significantly decrease the severity of an outbreak by providing timely and
accurate accounts of the numbers of new infections and deaths. Our model also
shows that although providing such information beginning as early as possible
is best, even starting to provide it well into the course of an outbreak can
significantly reduce the severity of the outbreak. We illustrate our results
with a simulated outbreak of Ebola Hemorrhagic Fever in the town of Huntington,
West Virginia (population about 50, 000).


  The internal dynamics of macro-molecular systems is characterized by widely
separated time scales, ranging from fraction of ps to ns. In ordinary molecular
dynamics simulations, the elementary time step dt used to integrate the
equation of motion needs to be chosen much smaller of the shortest time scale,
in order not to cut-off important physical effects. We show that, in systems
obeying the over-damped Langevin Eq., the fast molecular dynamics which occurs
at time scales smaller than dt can be analytically integrated out and gives
raise to a time-dependent correction to the diffusion coefficient, which we
rigorously compute. The resulting effective Langevin equation describes by
construction the same long-time dynamics, but has a lower time resolution
power, hence it can be integrated using larger time steps dt. We illustrate and
validate this method by studying the diffusion of a point-particle in a
one-dimensional toy-model and the denaturation of a protein.


  We will introduce the biological motivation of the $\gamma$- food-limited
model with variable parameters. New criteria are established for the existence
and global stability of positive periodic solutions. To prove the existence of
steady-state solutions, we used the upper-lower solution method where the
existence of at least one positive periodic solution is obtained by
constructing a pair of upper and lower solutions and application of the
Friedreichs Theorem. Numerical simulations illustrate effects of periodic
variation in the values of the basic biological and environmental parameters
and how the adaptive harvesting strategies affect fishing stocks.


  Conflict destabilizes social interactions and impedes cooperation at multiple
scales of biological organization. Of fundamental interest are the causes of
turbulent periods of conflict. We analyze conflict dynamics in a monkey society
model system. We develop a technique, Inductive Game Theory, to extract
directly from time-series data the decision-making strategies used by
individuals and groups. This technique uses Monte Carlo simulation to test
alternative causal models of conflict dynamics. We find individuals base their
decision to fight on memory of social factors, not on short timescale
ecological resource competition. Furthermore, the social assessments on which
these decisions are based are triadic (self in relation to another pair of
individuals), not pairwise. We show that this triadic decision making causes
long conflict cascades and that there is a high population cost of the large
fights associated with these cascades. These results suggest that individual
agency has been over-emphasized in the social evolution of complex aggregates,
and that pair-wise formalisms are inadequate. An appreciation of the empirical
foundations of the collective dynamics of conflict is a crucial step towards
its effective management.


  Organisms are composed of multiple chemical elements such as carbon,
nitrogen, and phosphorus. The scarcity of any of these elements can severely
restrict organismal and population growth. However, many trophic interaction
models only consider carbon limitation via energy flow. In this paper, we
construct an algal growth model with the explicit incorporation of light and
nutrient availability to characterize both carbon and phosphorus limitations.
We provide a global analysis of this model to illustrate how light and nutrient
availability regulate algal dynamics.


  Recently we have introduced a simplified model of ecosystem assembly (Capitan
et al., 2009) for which we are able to map out all assembly pathways generated
by external invasions in an exact manner. In this paper we provide a deeper
analysis of the model, obtaining analytical results and introducing some
approximations which allow us to reconstruct the results of our previous work.
In particular, we show that the population dynamics equations of a very general
class of trophic-level structured food-web have an unique interior equilibrium
point which is globally stable. We show analytically that communities found as
end states of the assembly process are pyramidal and we find that the
equilibrium abundance of any species at any trophic level is approximately
inversely proportional to the number of species in that level. We also find
that the per capita growth rate of a top predator invading a resident community
is key to understand the appearance of complex end states reported in our
previous work. The sign of these rates allows us to separate regions in the
space of parameters where the end state is either a single community or a
complex set containing more than one community. We have also built up
analytical approximations to the time evolution of species abundances that
allow us to determine, with high accuracy, the sequence of extinctions that an
invasion may cause. Finally we apply this analysis to obtain the communities in
the end states. To test the accuracy of the transition probability matrix
generated by this analytical procedure for the end states, we have compared
averages over those sets with those obtained from the graph derived by
numerical integration of the Lotka-Volterra equations. The agreement is
excellent.