  Linear bio-molecular motors move unidirectionally along a track by
coordinating several different processes, such as fuel (ATP) capture,
hydrolysis, conformational changes, binding and unbinding from a track, and
center-of-mass diffusion. A better understanding of the interdependencies
between these processes, which take place over a wide range of different time
scales, would help elucidate the general operational principles of molecular
motors. Artificial molecular motors present a unique opportunity for such a
study because motor structure and function are a priori known. Here we describe
use of a Master equation approach, integrated with input from Langevin and
molecular dynamics modeling, to stochastically model a molecular motor across
many time scales. We apply this approach to a specific concept for an
artificial protein motor, the Tumbleweed.


  Statistical consistency in phylogenetics has traditionally referred to the
accuracy of estimating phylogenetic parameters for a fixed number of species as
we increase the number of characters. However, as sequences are often of fixed
length (e.g. for a gene) although we are often able to sample more taxa, it is
useful to consider a dual type of statistical consistency where we increase the
number of species, rather than characters. This raises some basic questions:
what can we learn about the evolutionary process as we increase the number of
species? In particular, does having more species allow us to infer the
ancestral state of characters accurately? This question is particularly
relevant when sequence site evolution varies in a complex way from character to
character, as well as for reconstructing ancestral sequences. In this paper, we
assemble a collection of results to analyse various approaches for inferring
ancestral information with increasing accuracy as the number of taxa increases.


  Fuzzy c-means clustering is widely used to identify cluster structures in
high-dimensional data sets, such as those obtained in DNA microarray and
quantitative proteomics experiments. One of its main limitations is the lack of
a computationally fast method to determine the two parameters fuzzifier and
cluster number. Wrong parameter values may either lead to the inclusion of
purely random fluctuations in the results or ignore potentially important data.
The optimal solution has parameter values for which the clustering does not
yield any results for a purely random data set but which detects cluster
formation with maximum resolution on the edge of randomness. Estimation of the
optimal parameter values is achieved by evaluation of the results of the
clustering procedure applied to randomized data sets. In this case, the optimal
value of the fuzzifier follows common rules that depend only on the main
properties of the data set. Taking the dimension of the set and the number of
objects as input values instead of evaluating the entire data set allows us to
propose a functional relationship determining its value directly. This result
speaks strongly against setting the fuzzifier equal to 2 as typically done in
many previous studies. Validation indices are generally used for the estimation
of the optimal number of clusters. A comparison shows that the minimum distance
between the centroids provides results that are at least equivalent or better
than those obtained by other computationally more expensive indices.


  The coupling between the depolymerization of microtubules (MTs) and the
motion of the Dam1 ring complex is now thought to play an important role in the
generation of forces during mitosis. Our current understanding of this motion
is based on a number of detailed computational models. Although these models
realize possible mechanisms for force transduction, they can be extended by
variation of any of a large number of poorly measured parameters and there is
no clear strategy for determining how they might be distinguished
experimentally. Here we seek to identify and analyze two distinct mechanisms
present in the computational models. In the first the splayed protofilaments at
the end of the depolymerizing MT physically prevent the Dam1 ring from falling
off the end, in the other an attractive binding secures the ring to the
microtubule. Based on this analysis, we discuss how to distinguish between
competing models that seek to explain how the Dam1 ring stays on the MT. We
propose novel experimental approaches that could resolve these models for the
first time, either by changing the diffusion constant of the Dam1 ring (e.g.,
by tethering a long polymer to it) or by using a time varying load.


  Network epidemiology often assumes that the relationships defining the social
network of a population are static. The dynamics of relationships is only taken
indirectly into account, by assuming that the relevant information to study
epidemic spread is encoded in the network obtained by considering numbers of
partners accumulated over periods of time roughly proportional to the
infectious period of the disease at hand. On the other hand, models explicitly
including social dynamics are often too schematic to provide a reasonable
representation of a real population, or so detailed that no general conclusions
can be drawn from them. Here we present a model of social dynamics that is
general enough that its parameters can be obtained by fitting data from surveys
about sexual behaviour, but that can still be studied analytically, using mean
field techniques. This allows us to obtain some general results about epidemic
spreading. We show that using accumulated network data to estimate the static
epidemic threshold leads to a significant underestimation of it. We also show
that, for a dynamic network, the relative epidemic threshold is an increasing
function of the infectious period of the disease, implying that the static
value is a lower bound to the real threshold.


  Metapopulation models provide the theoretical framework for describing
disease spread between different populations connected by a network. In
particular, these models are at the basis of most simulations of pandemic
spread. They are usually studied at the mean-field level by neglecting
fluctuations. Here we include fluctuations in the models by adopting fully
stochastic descriptions of the corresponding processes. This level of
description allows to address analytically, in the SIS and SIR cases, problems
such as the existence and the calculation of an effective threshold for the
spread of a disease at a global level. We show that the possibility of the
spread at the global level is described in terms of (bond) percolation on the
network. This mapping enables us to give an estimate (lower bound) for the
pandemic threshold in the SIR case for all values of the model parameters and
for all possible networks.


  In Aldous and Shields (1988), a model for a rooted, growing random binary
tree was presented. For some c>0, an external vertex splits at rate c^(-i) (and
becomes internal) if its distance from the root (depth) is i. For c>1, we
reanalyse the tree profile, i.e. the numbers of external vertices in depth
i=1,2,.... Our main result are concrete formulas for the expectation and
covariance-structure of the profile. In addition, we present the application of
the model to cellular ageing. Here, we assume that nodes in depth h+1 are
senescent, i.e. do not split. We obtain a limit result for the proportion of
non-senescent vertices for large h.


  Using a recently developed mesoscopic theory of protein dielectrics, we have
calculated the salt bridge energies, total residue electrostatic potential
energies, and transfer energies into a low dielectric amyloid-like phase for 12
species and mutants of the prion protein. Salt bridges and self energies play
key roles in stabilizing secondary and tertiary structural elements of the
prion protein. The total electrostatic potential energy of each residue was
found to be invariably stabilizing. Residues frequently found to be mutated in
familial prion disease were among those with the largest electrostatic
energies. The large barrier to charged group desolvation imposes regional
constraints on involvement of the prion protein in an amyloid aggregate,
resulting in an electrostatic amyloid recruitment profile that favours regions
of sequence between alpha helix 1 and beta strand 2, the middles of helices 2
and 3, and the region N-terminal to alpha helix 1. We found that the
stabilization due to salt bridges is minimal among the proteins studied for
disease-susceptible human mutants of prion protein.


  Background: A metabolic genotype comprises all chemical reactions an organism
can catalyze via enzymes encoded in its genome. A genotype is viable in a given
environment if it is capable of producing all biomass components the organism
needs to survive and reproduce. Previous work has focused on the properties of
individual genotypes while little is known about how genome-scale metabolic
networks with a given function can vary in their reaction content. Results: We
here characterize spaces of such genotypes. Specifically, we study metabolic
genotypes whose phenotype is viability in minimal chemical environments that
differ in their sole carbon sources. We show that regardless of the number of
reactions in a metabolic genotype, the genotypes of a given phenotype typically
form vast, connected, and unstructured sets -- genotype networks -- that nearly
span the whole of genotype space. The robustness of metabolic phenotypes to
random reaction removal in such spaces has a narrow distribution with a high
mean. Different carbon sources differ in the number of metabolic genotypes in
their genotype network; this number decreases as a genotype is required to be
viable on increasing numbers of carbon sources, but much less than if metabolic
reactions were used independently across different chemical environments.
Conclusions: Our work shows that phenotype-preserving genotype networks have
generic organizational properties and that these properties are insensitive to
the number of reactions in metabolic genotypes.


  The one-dimensional problem of selecting the triple helix with the highest
volume fraction is solved and hence the condition for a helix to be
close-packed is obtained. The close-packed triple helix is shown to have a
pitch angle of $v_{CP} =43.3 ^\circ$. Contrary to the conventional notion, we
suggest that close packing form the underlying principle behind the structure
of collagen, and the implications of this suggestion are considered. Further,
it is shown that the unique zero-twist structure with no strain-twist coupling
is practically identical to the close-packed triple helix. Some of the
difficulties for the current understanding of the structure of collagen are
reviewed: The ambiguity in assigning crystal structures for collagen-like
peptides, and the failure to satisfactorily calculate circular dichroism
spectra. Further, the proposed new geometrical structure for collagen is better
packed than both the 10/3 and the 7/2 structure. A feature of the suggested
collagen structure is the existence of a central channel with negatively
charged walls. We find support for this structural feature in some of the early
x-ray diffraction data of collagen. The central channel of the structure
suggests the possibility of a one-dimensional proton lattice. This geometry can
explain the observed magic angle effect seen in NMR studies of collagen. The
central channel also offers the possibility of ion transport and may cast new
light on various biological and physical phenomena, including
biomineralization.


  We develop a spatially explicit biophysical model of the hippocampal CA3-CA1
presynaptic bouton to study local calcium dynamics leading to vesicle fusion. A
kinetic model with two calcium sensors is formulated specifically for the
CA3-CA1 synapse. The model includes a sensor for fast synchronous release that
lasts a few tens of milliseconds and a sensor for slow asynchronous release
that lasts a few hundred milliseconds. We show that a variety of extant data on
CA3-CA1 synapse can be accounted for consistently only when a refractory period
of the order of few milliseconds between releases is introduced. Including a
second sensor for asynchronous release that has a slow unbinding site and
therefore an embedded long memory, is shown to play a role in short-term
plasticity by facilitating release. For synchronous release mediated by
Synaptotagmin II a third time scale is revealed in addition to the fast and
slow release. This third time scale corresponds to "stimulus-correlated
super-fast" neurotransmitter release. Our detailed spatial simulation indicates
that all three-time scales of neurotransmitter release are an emergent property
of the calcium sensor and independent of synaptic ultrastructure. Furthermore,
it allows us to identify features of synaptic transmission that are universal
and those that are modulated by structure.


  The observed cooperation on the level of genes, cells, tissues, and
individuals has been the object of intense study by evolutionary biologists,
mainly because cooperation often flourishes in biological systems in apparent
contradiction to the selfish goal of survival inherent in Darwinian evolution.
In order to resolve this paradox, evolutionary game theory has focused on the
Prisoner's Dilemma (PD), which incorporates the essence of this conflict. Here,
we encode strategies for the iterated Prisoner's Dilemma (IPD) in terms of
conditional probabilities that represent the response of decision pathways
given previous plays. We find that if these stochastic strategies are encoded
as genes that undergo Darwinian evolution, the environmental conditions that
the strategies are adapting to determine the fixed point of the evolutionary
trajectory, which could be either cooperation or defection. A transition
between cooperative and defective attractors occurs as a function of different
parameters such a mutation rate, replacement rate, and memory, all of which
affect a player's ability to predict an opponent's behavior.


  Realistic 3D-conformations of protein structures can be embedded in a cubic
lattice using exclusively integer numbers, additions, subtractions and boolean
operations.


  A series of coarse-grained models have been developed for the study of the
molecular dynamics of RNA nanostructures. The models in the series have one to
three beads per nucleotide and include different amounts of detailed structural
information. Such a treatment allows us to reach, for the systems of thousands
of nucleotides, a time scale of microseconds (i.e. by three orders of magnitude
longer than in the full atomistic modelling) and thus to enable simulations of
large RNA polymers in the context of bionanotechnology. We find that the
3-beads-per-nucleotide models, described by a set of just a few universal
parameters, are able to describe different RNA conformations and are comparable
in structural precision to the models where detailed values of the backbone
P-C4' dihedrals taken from a reference structure are included. These findings
are discussed in the context of the RNA conformation classes.


  An active microbiota, reaching up to 10 E+7 cells/g, was found to inhabit a
naturally occurring asphalt lake characterized by low water activity and
elevated temperature. Geochemical and molecular taxonomic approaches revealed
novel and deeply branching microbial assemblages mediating anaerobic
hydrocarbon degradation, metal respiration and C1 utilization pathways. These
results open a window into the origin and adaptive evolution of microbial life
within recalcitrant hydrocarbon matrices, and establish the site as a useful
analog for the liquid hydrocarbon environments on Saturn's moon Titan.


  We study the evolution of cooperation under the assumption that the
collective benefits of group membership can only be harvested if the fraction
of cooperators within the group, i.e. their critical mass, exceeds a threshold
value. Considering structured populations, we show that a moderate fraction of
cooperators can prevail even at very low multiplication factors if the critical
mass is minimal. For larger multiplication factors, however, the level of
cooperation is highest at an intermediate value of the critical mass. The
latter is robust to variations of the group size and the interaction network
topology. Applying the optimal critical mass threshold, we show that the
fraction of cooperators in public goods games is significantly larger than in
the traditional linear model, where the produced public good is proportional to
the fraction of cooperators within the group.


  Balanced minimum evolution (BME) is a statistically consistent distance-based
method to reconstruct a phylogenetic tree from an alignment of molecular data.
In 2000, Pauplin showed that the BME method is equivalent to optimizing a
linear functional over the BME polytope, the convex hull of the BME vectors
obtained from Pauplin's formula applied to all binary trees. The BME method is
related to the Neighbor Joining (NJ) algorithm, now known to be a greedy
optimization of the BME principle. Further, the NJ and BME algorithms have been
studied previously to understand when the NJ Algorithm returns a BME tree for
small numbers of taxa. In this paper we aim to elucidate the structure of the
BME polytope and strengthen knowledge of the connection between the BME method
and NJ Algorithm. We first prove that any subtree-prune-regraft move from a
binary tree to another binary tree corresponds to an edge of the BME polytope.
Moreover, we describe an entire family of faces parametrized by disjoint
clades. We show that these {\em clade-faces} are smaller dimensional BME
polytopes themselves. Finally, we show that for any order of joining nodes to
form a tree, there exists an associated distance matrix (i.e., dissimilarity
map) for which the NJ Algorithm returns the BME tree. More strongly, we show
that the BME cone and every NJ cone associated to a tree $T$ have an
intersection of positive measure.


  We propose a statistical method to test whether two phylogenetic trees with
given alignments are significantly incongruent. Our method compares the two
distributions of phylogenetic trees given by the input alignments, instead of
comparing point estimations of trees. This statistical approach can be applied
to gene tree analysis for example, detecting unusual events in genome evolution
such as horizontal gene transfer and reshuffling. Our method uses difference of
means to compare two distributions of trees, after embedding trees in a vector
space. Bootstrapping alignment columns can then be applied to obtain p-values.
To compute distances between means, we employ a "kernel trick" which speeds up
distance calculations when trees are embedded in a high-dimensional feature
space, e.g. splits or quartets feature space. In this pilot study, first we
test our statistical method's ability to distinguish between sets of gene trees
generated under coalescence models with species trees of varying dissimilarity.
We follow our simulation results with applications to various data sets of
gophers and lice, grasses and their endophytes, and different fungal genes from
the same genome. A companion toolkit, {\tt Phylotree}, is provided to
facilitate computational experiments.


  Quasi-stationary distributions, as discussed by Darroch & Seneta (1965), have
been used in biology to describe the steady state behaviour of population
models which, while eventually certain to become extinct, nevertheless maintain
an apparent stochastic equilibrium for long periods. These distributions have
some drawbacks: they need not exist, nor be unique, and their calculation can
present problems. In this paper, we give biologically plausible conditions
under which the quasi-stationary distribution is unique, and can be closely
approximated by distributions that are simple to compute.


  The ratio of males to females in a population is a meaningful characteristic
of sexual species. The reason for this biological property to be available to
the observers of nature seems to be a question never asked. Introducing the
notion of historically adapted populations as global minimizers of maintenance
cost functions, we propose a theoretical explanation for the reported stability
of this feature. This mathematical formulation suggests that sex ratio could be
considered as an indirect result shaped by the antagonism between the size of
the population and the finiteness of resources.


  Theoretically, communities at or near their equilibrium species number resist
entry of new species. Such 'biotic resistance' recently has been questioned
because of successful entry of alien species into diverse natural communities.
Data on 10,409 naturalizations of 5350 plant species over 16 sites dispersed
globally show exponential distributions for both species over sites and sites
over number of species shared. These exponentials signal a statistical
mechanics of species distribution, assuming two conditions. First, species and
sites are equivalent, either identical ('neutral'), or so complex that the
chance a species is in the right place at the right time is vanishingly small
('idiosyncratic'); the range of species and sites in our data disallows a
neutral explanation. Secondly, the total number of naturalisations is fixed in
any era by a 'regulator'. Previous correlation of species naturalization rates
with net primary productivity over time suggests that regulator is related to
productivity. We conclude that biotic resistance is a moving ceiling, with
resistance controlled by productivity. The general observation that the
majority of species occur naturally at only a few sites but only a few at many
now has a quantitative [exponential] character, offering the study of species'
distributions a previously unavailable rigor.


  The cyanobacterium Synechococcus elongatus uses both a protein
phosphorylation cycle and a transcription-translation cycle to generate
circadian rhythms that are highly robust against biochemical noise. We use
stochastic simulations to analyze how these cycles interact to generate stable
rhythms in growing, dividing cells. We find that a protein phosphorylation
cycle by itself is robust when protein turnover is low. For high decay or
dilution rates (and co mpensating synthesis rate), however, the
phosphorylation-based oscillator loses its integrity. Circadian rhythms thus
cannot be generated with a phosphorylation cycle alone when the growth rate,
and consequently the rate of protein dilution, is high enough; in practice, a
purely post-translational clock ceases to function well when the cell doubling
time drops below the 24 hour clock period. At higher growth rates, a
transcription-translation cycle becomes essential for generating robust
circadian rhythms. Interestingly, while a transcription-translation cycle is
necessary to sustain a phosphorylation cycle at high growth rates, a
phosphorylation cycle can dramatically enhance the robustness of a
transcription-translation cycle at lower protein decay or dilution rates. Our
analysis thus predicts that both cycles are required to generate robust
circadian rhythms over the full range of growth conditions.


  Various molecular interaction networks have been claimed to follow power-law
decay for their global connectivity distribution. It has been proposed that
there may be underlying generative models that explain this heavy-tailed
behavior by self-reinforcement processes such as classical or hierarchical
scale-free network models. Here we analyze a comprehensive data set of
protein-protein and transcriptional regulatory interaction networks in yeast,
an E. coli metabolic network, and gene activity profiles for different
metabolic states in both organisms. We show that in all cases the networks have
a heavy-tailed distribution, but most of them present significant differences
from a power-law model according to a stringent statistical test. Those few
data sets that have a statistically significant fit with a power-law model
follow other distributions equally well. Thus, while our analysis supports that
both global connectivity interaction networks and activity distributions are
heavy-tailed, they are not generally described by any specific distribution
model, leaving space for further inferences on generative models.


  Cancer is a proliferation disease affecting a genetically unstable cell
population, in which molecular alterations can be somatically inherited by
genetic, epigenetic or extragenetic transmission processes, leading to a
cooperation of neoplastic cells within tumoral tissue. The efflux protein
P-glycoprotein (P gp) is overexpressed in many cancer cells and has known
capacity to confer multidrug resistance to cytotoxic therapies. Recently,
cell-to-cell P-gp transfers have been shown. Herein, we combine experimental
evidence and a mathematical model to examine the consequences of an
intercellular P-gp trafficking in the extragenetic transfer of multidrug
resistance from resistant to sensitive cell subpopulations. We report
cell-to-cell transfers of functional P-gp in co-cultures of a P-gp
overexpressing human breast cancer MCF-7 cell variant, selected for its
resistance towards doxorubicin, with the parental sensitive cell line. We found
that P-gp as well as efflux activity distribution are progressively reorganized
over time in co-cultures analyzed by flow cytometry. A mathematical model based
on a Boltzmann type integro-partial differential equation structured by a
continuum variable corresponding to P-gp activity describes the cell
populations in co-culture. The mathematical model elucidates the population
elements in the experimental data, specifically, the initial proportions, the
proliferative growth rates, and the transfer rates of P-gp in the sensitive and
resistant subpopulations. We confirmed cell-to-cell transfer of functional
P-gp. The transfer process depends on the gradient of P-gp expression in the
donor-recipient cell interactions, as they evolve over time. Extragenetically
acquired drug resistance is an additional aptitude of neoplastic cells which
has implications in the diagnostic value of P-gp expression and in the design
of chemotherapy regimens


  We investigated the role of peptide folding stability in peptide
immunogenicity. It was the aim of this thesis to implement a stability
criterion based on energy computations using an AMBER force field, and to test
the implementation with a large dataset.


  Sequence organizations are viewed from two points: one is from informational
redundancy or informational correlation (IC) and another is from k-mer
frequency statistics. Two problems are investigated. The first is how the ICs
exceed the fluctuation bound and the order emerges from fluctuation in a genome
when the sequence length attains some critical value. We demonstrated that the
transition from fluctuation to order takes place at about sequence length
200-300 thousands bases for human and E coli genome. It means that the life
emerges from a region between macroscopic and microscopic. The second is about
the statistical law of the k-mer organization in a genome under the
evolutionary pressure and functional selection. We deduced a sum rule Q(k,N) on
the k-mer frequency deviations from the randomness in a N-long sequence of
genome and deduced the relations of Q(k,N) with k and N. We found that Q(k,N)
increases with length N at a constant rate for most genome sequences and
demonstrated that when the functional selection of k-mers is accumulated to
some critical value the ordering takes place. An important finding is the sum
rule correlated with the evolutionary complexity of the genome.


  This article considers a model of genealogy corresponding to a regular
exchangeable coalescent (also known as Xi-coalescent) started from a large
finite configuration, and undergoing neutral mutations. Asymptotic expressions
for the number of active lineages were obtained by the author in a previous
work. Analogous results for the number of active mutation-free lineages and the
combined lineage lengths are derived using the same martingale-based technique.
They are given in terms of convergence in probability, while extensions to
convergence in moments and convergence almost surely are discussed. The above
mentioned results have direct consequences on the sampling theory in the
Xi-coalescent setting. In particular, the regular Xi-coalescents that come down
from infinity (i.e., with locally finite genealogies), have an asymptotically
equal number of families under the corresponding infinite alleles and infinite
sites models. In special cases, quantitative asymptotic formulae for the number
of families that contain a fixed number of individuals can be given.


  Periodicity in population dynamics is a fundamental issue. In addition to
current species-specific analyses, allometry facilitates understanding of limit
cycles amongst different species. So far, body-size regressions have been
derived for the oscillation period of the population densities of warm-blooded
species, in particular herbivores. Here, we extend the allometric analysis to
other clades, allowing for a comparison between the obtained slopes and
intercepts. The oscillation periods were derived from databases and original
studies to cover a broad range of conditions and species. Then, values were
related to specific body size by regression analysis. For different groups of
herbivorous species, the oscillation period increased as a function of
individual mass as a power law with exponents of 0.11-0.27. The intercepts of
the resulting linear regressions indicated that cycle times for equally-sized
species increased from homeotherms up to invertebrates. Overall, cycle times
for predators did not scale to body size. Implications for these differences
were addressed in the light of intra- and interspecific delays.


  The dynamics of tumour evolution are not well understood. In this paper we
provide a statistical framework for evaluating the molecular variation observed
in different parts of a colorectal tumour. A multi-sample version of the Ewens
Sampling Formula forms the basis for our modelling of the data, and we provide
a simulation procedure for use in obtaining reference distributions for the
statistics of interest. We also describe the large-sample asymptotics of the
joint distributions of the variation observed in different parts of the tumour.
While actual data should be evaluated with reference to the simulation
procedure, the asymptotics serve to provide theoretical guidelines, for
instance with reference to the choice of possible statistics.


  We synthesise traditional unstructured food webs, allometric body size
scaling, trait-based modelling, and physiologically structured modelling to
provide a novel and ecologically relevant tool for size-structured food webs.
The framework allows food web models to include ontogenetic growth and
life-history omnivory at the individual level by resolving the population
structure of each species as a size-spectrum. Each species is characterised by
the trait 'size at maturation', and all model parameters are made species
independent through scaling with individual body size and size at maturation.
Parameter values are determined from cross-species analysis of fish communities
as life-history omnivory is widespread in aquatic systems, but may be
reparameterised for other systems. An ensemble of food webs is generated and
the resulting communities are analysed at four levels of organisation:
community level, species level, trait level, and individual level. The model
may be solved analytically by assuming that the community spectrum follows a
power law. The analytical solution provides a baseline expectation of the
results of complex food web simulations, and agrees well with the predictions
of the full model on 1) biomass distribution as a function of individual size,
2) biomass distribution as a function of size at maturation, and 3) relation
between predator-prey mass ratio of preferred and eaten food. The full model
additionally predicts the diversity distribution as a function of size at
maturation.


  We derive an expression for the variation between parallel trajectories in
phenotypic evolution, extending the well known result that predicts the mean
evolutionary path in adaptive dynamics or quantitative genetics. We show how
this expression gives rise to the notion of fluctuation domains - parts of the
fitness landscape where the rate of evolution is very predictable (due to
fluctuation dissipation) and parts where it is highly variable (due to
fluctuation enhancement). These fluctuation domains are determined by the
curvature of the fitness landscape. Regions of the fitness landscape with
positive curvature, such as adaptive valleys or branching points, experience
enhancement. Regions with negative curvature, such as adaptive peaks,
experience dissipation. We explore these dynamics in the ecological scenarios
of implicit and explicit competition for a limiting resource.


  This paper investigates the optimal harvesting strategy for a single species
living in random environments whose growth is given by a regime-switching
diffusion. Harvesting acts as a (stochastic) control on the size of the
population. The objective is to find a harvesting strategy which maximizes the
expected total discounted income from harvesting {\em up to the time of
extinction} of the species; the income rate is allowed to be state- and
environment-dependent. This is a singular stochastic control problem with both
the extinction time and the optimal harvesting policy depending on the initial
condition. One aspect of receiving payments up to the random time of extinction
is that small changes in the initial population size may significantly alter
the extinction time when using the same harvesting policy. Consequently, one no
longer obtains continuity of the value function using standard arguments for
either regular or singular control problems having a fixed time horizon. This
paper introduces a new sufficient condition under which the continuity of the
value function for the regime-switching model is established. Further, it is
shown that the value function is a viscosity solution of a coupled system of
quasi-variational inequalities. The paper also establishes a verification
theorem and, based on this theorem, an $\varepsilon$-optimal harvesting
strategy is constructed under certain conditions on the model. Two examples are
analyzed in detail.


  KIF1A, a processive single headed kinesin superfamily motor, hydrolyzes
Adenosine triphosphate (ATP) to move along a filamentous track called
microtubule. The stochastic movement of KIF1A on the track is characterized by
an alternating sequence of pause and translocation. The sum of the durations of
pause and the following translocation defines the dwell time. Using the NOSC
model (Nishinari et. al. PRL, {\bf 95}, 118101 (2005)) of individual KIF1A, we
systematically derive an analytical expression for the dwell time distribution.
More detailed information is contained in the probability densities of the
"conditional dwell times" $\tau_{\pm\pm}$ in between two consecutive steps each
of which could be forward (+) or backward (-). We calculate the probability
densities $\Xi_{\pm\pm}$ of these four conditional dwell times. However, for
the convenience of comparison with experimental data, we also present the two
distributions $\Xi_{\pm}^{*}$ of the times of dwell before a forward (+) and a
backward (-) step. In principle, our theoretical prediction can be tested by
carrying out single-molecule experiments with adequate spatio-temporal
resolution.


  A report of the meeting "Challenges in experimental data integration within
genome-scale metabolic models", Institut Henri Poincar\'e, Paris, October 10-11
2009, organized by the CNRS-MPG joint program in Systems Biology.


  The mapping between genotype and phenotype is encoded in the complex web of
epistatic interaction between genetic loci. In this rugged fitness landscape,
recombination processes, which tend to increase variation in the population,
compete with selection processes that tend to reduce genetic variation. Here we
show that the Bose-Einstein distribution describe the multiple stationary
states of a diploid population under this multi-loci evolutionary dynamics.
Moreover, the evolutionary process might undergo an interesting condensation
phase transition in the universality class of a Bose-Einstein condensation when
a finite fraction of pairs of linked loci, is fixed into given allelic states.
Below this phase transition the genetic variation within a species is
significantly reduced and only maintained by the remaining polymorphic loci.


  The thermodynamical properties of heterogeneous DNA sequences are computed by
path integral techniques applied to a nonlinear model Hamiltonian. The base
pairs relative displacements are interpreted as time dependent paths whose
amplitudes are consistent with the model potential for the hydrogen bonds
between complementary strands. The portion of configuration space contributing
to the partition function is determined, at any temperature, by selecting the
ensemble of paths which fulfill the second law of thermodynamics. For a short
DNA fragment, the denaturation is signaled by a succession of peaks in the
specific heat plots while the entropy grows continuously versus $T$. Thus, the
opening of the double strand with bubble formation appears as a smooth
crossover due to base pair fluctuation effects which are accounted for by the
path integral method. The multistep transition is driven by the AT-rich regions
of the DNA fragment. The base pairs path ensemble shows an enhanced degree of
cooperativity at about the same temperatures for which the specific heat peaks
occur. These findings establish a link between microscopic and macroscopic
signatures of the transition. The fractions of mean base pair stretchings are
computed by varying the AT base pairs content and taking some threshold values
for the occurrence of the molecule denaturation.


  As early indicated by Charles Darwin, languages behave and change very much
like living species. They display high diversity, differentiate in space and
time, emerge and disappear. A large body of literature has explored the role of
information exchanges and communicative constraints in groups of agents under
selective scenarios. These models have been very helpful in providing a
rationale on how complex forms of communication emerge under evolutionary
pressures. However, other patterns of large-scale organization can be described
using mathematical methods ignoring communicative traits. These approaches
consider shorter time scales and have been developed by exploiting both
theoretical ecology and statistical physics methods. The models are reviewed
here and include extinction, invasion, origination, spatial organization,
coexistence and diversity as key concepts and are very simple in their defining
rules. Such simplicity is used in order to catch the most fundamental laws of
organization and those universal ingredients responsible for qualitative
traits. The similarities between observed and predicted patterns indicate that
an ecological theory of language is emerging, supporting (on a quantitative
basis) its ecological nature, although key differences are also present. Here
we critically review some recent advances lying and outline their implications
and limitations as well as open problems for future research.


  Employing Monte Carlo simulations of semiflexible polymer rings in weak
spherical confinement a conformational transition to figure eight shaped,
writhed configurations is discovered and quantified.


  Single-molecule mechanical manipulation has enabled the quantitative
understanding of the kinetics of bond ruptures as well as protein unfolding
mechanism. Single-molecule experiments with theoretical models have allowed one
to gain insight into free energy landscape for chemical bond and/or protein
folding. For mechanically induced bond rupture, the bond-rupture kinetics may
be governed by loading device. However, the role of loading device on the
kinetics of mechanical rupture has been rarely received much attention until
recently. In this work, we have theoretically and/or computationally studied
the effect of loading-device stiffness on the kinetics of mechanical unfolding.
Specifically, we have considered a one-dimensional model for a bond rupture
whose kinetics is depicted by Kramers' theory. It is elucidated that the
kinetics of bond rupture is determined by force constant of loading device. The
Brownian dynamics simulation of a bond rupture is considered in order to
validate our theory. It is illustrated that the mean rupture force is dependent
on the force constant of a loading device, such that increase in the
loading-device stiffness leads to the higher bond rupture force. Moreover, we
have taken into account the computational simulation of mechanical unfolding of
a small protein, i.e. {\beta}-hairpin. Our simulation shows that unfolding
force is highly correlated with the stiffness of a loading device. Our
numerical and theoretical studies highlight the significance of a loading
device on the kinetics of mechanical unfolding of a chemical bond and/or a
folded domain in the single-molecule mechanical experiments.


  We consider approaches for controlling the buildup of noise by design of
gates for chemical and biomolecular computing, in order to realize stable,
scalable networks for multi-step information processing. Solvable rate-equation
models are introduced and used to illustrate several recently developed
concepts and methodologies. We also outline future challenges and possible
research directions.


  Network analysis is currently used in a myriad of contexts: from identifying
potential drug targets to predicting the spread of epidemics and designing
vaccination strategies, and from finding friends to uncovering criminal
activity. Despite the promise of the network approach, the reliability of
network data is a source of great concern in all fields where complex networks
are studied. Here, we present a general mathematical and computational
framework to deal with the problem of data reliability in complex networks. In
particular, we are able to reliably identify both missing and spurious
interactions in noisy network observations. Remarkably, our approach also
enables us to obtain, from those noisy observations, network reconstructions
that yield estimates of the true network properties that are more accurate than
those provided by the observations themselves. Our approach has the potential
to guide experiments, to better characterize network data sets, and to drive
new discoveries.


  We discuss the appropriate techniques for modelling the geometry of open
ended elastic polymer molecules. The molecule is assumed to have fixed
endpoints on a boundary surface. In particular we discuss the concept of the
winding number, a directional measure of the linking of two curves, which can
be shown to be invariant to the set of continuous deformations vanishing at the
polymer's end-point and which forbid it from passing through itself. This
measure is shown to be the appropriate constraint required to evaluate the
geometrical properties of a constrained DNA molecule. Using the net winding
measure we define a model of an open ended constrained DNA molecule which
combines the necessary constraint of self-avoidance with being analytically
tractable. This model builds upon the local models of Bouchiat and Mezard
(2000). In particular, we present a new derivation of the polar writhe
expression, which detects both the local winding of the curve and non local
winding between different sections of the curve. We then show that this
expression correctly tracks the net twisting of a DNA molecule subject to
rotation at the endpoints, unlike other definitions used in the literature.


  A neutral ecology model is simulated on an island chain, in which
neighbouring islands can exchange individuals but only the first island is able
to receive immigrants from a metacommunity. It is found by several measures
that biodiversity decreases along the chain, being highest for the first
island. Subtle changes in taxon abundance distributions can be detected when
islands in the chain are compared to diversity-matched single islands. The
results potentially apply to human microbial diversity, but highlight the
difficulty of using static single-site taxon abundance distributions to
discriminate between dispersal limitation mechanisms.


  We predict analytically that diagonal correlations of amino acid positions
within protein sequences statistically enhance protein propensity for
nonspecific binding. We use the term 'promiscuity' to describe such nonspecific
binding. Diagonal correlations represent statistically significant repeats of
sequence patterns where amino acids of the same type are clustered together.
The predicted effect is qualitatively robust with respect to the form of the
microscopic interaction potentials and the average amino acid composition. Our
analytical results provide an explanation for the enhanced diagonal
correlations observed in hubs of eukaryotic organismal proteomes [J. Mol. Biol.
409, 439 (2011)]. We suggest experiments that will allow direct testing of the
predicted effect.


  "Evidence and Evolution: the Logic behind the Science" was published in 2008
by Elliott Sober. It examines the philosophical foundations of the statistical
arguments used to evaluate hypotheses in evolutionary biology, based on simple
examples and likelihood ratios. The difficulty with reading the book from a
statistician's perspective is the reluctance of the author to engage into model
building and even less into parameter estimation. The first chapter nonetheless
constitutes a splendid coverage of the most common statistical approaches to
testing and model comparison, even though the advocation of the Akaike
information criterion against Bayesian alternatives is rather forceful. The
book also covers an examination of the "intelligent design" arguments against
the Darwinian evolution theory, predictably if unnecessarily resorting to
Popperian arguments to correctly argue that the creationist perspective fails
to predict anything. The following chapters cover the more relevant issues of
assessing selection versus drift and of testing for the presence of a common
ancestor. While remaining a philosophy treatise, Evidence and Evolution is
written in a way that is accessible to laymen, if rather unusual from a
statistician viewpoint, and the insight about testing issues gained from
Evidence and Evolution makes it a worthwhile read.


  We study several variants of the stochastic four-state rock-paper-scissors
game or, equivalently, cyclic three-species predator-prey models with conserved
total particle density, by means of Monte Carlo simulations on one- and
two-dimensional lattices. Specifically, we investigate the influence of spatial
variability of the reaction rates and site occupancy restrictions on the
transient oscillations of the species densities and on spatial correlation
functions in the quasi-stationary coexistence state. For small systems, we also
numerically determine the dependence of typical extinction times on the number
of lattice sites. In stark contrast with two-species stochastic Lotka-Volterra
systems, we find that for our three-species models with cyclic competition
quenched disorder in the reaction rates has very little effect on the dynamics
and the long-time properties of the coexistence state. Similarly, we observe
that site restriction only has a minor influence on the system's dynamical
properties. Our results therefore demonstrate that the features of the spatial
rock-paper-scissors system are remarkably robust with respect to model
variations, and stochastic fluctuations as well as spatial correlations play a
comparatively minor role.


  Rooted phylogenetic networks are often constructed by combining trees,
clusters, triplets or characters into a single network that in some
well-defined sense simultaneously represents them all. We review these four
models and investigate how they are related. In general, the model chosen
influences the minimum number of reticulation events required. However, when
one obtains the input data from two binary trees, we show that the minimum
number of reticulations is independent of the model. The number of
reticulations necessary to represent the trees, triplets, clusters (in the
softwired sense) and characters (with unrestricted multiple crossover
recombination) are all equal. Furthermore, we show that these results also hold
when not the number of reticulations but the level of the constructed network
is minimised. We use these unification results to settle several complexity
questions that have been open in the field for some time. We also give explicit
examples to show that already for data obtained from three binary trees the
models begin to diverge.


  Empirical substitution matrices represent the average tendencies of
substitutions over various protein families by sacrificing gene-level
resolution. We develop a codon-based model, in which mutational tendencies of
codon, a genetic code, and the strength of selective constraints against amino
acid replacements can be tailored to a given gene. First, selective constraints
averaged over proteins are estimated by maximizing the likelihood of each 1-PAM
matrix of empirical amino acid (JTT, WAG, and LG) and codon (KHG) substitution
matrices. Then, selective constraints specific to given proteins are
approximated as a linear function of those estimated from the empirical
substitution matrices.
  Akaike information criterion (AIC) values indicate that a model allowing
multiple nucleotide changes fits the empirical substitution matrices
significantly better. Also, the ML estimates of transition-transversion bias
obtained from these empirical matrices are not so large as previously
estimated. The selective constraints are characteristic of proteins rather than
species. However, their relative strengths among amino acid pairs can be
approximated not to depend very much on protein families but amino acid pairs,
because the present model, in which selective constraints are approximated to
be a linear function of those estimated from the JTT/WAG/LG/KHG matrices, can
provide a good fit to other empirical substitution matrices including cpREV for
chloroplast proteins and mtREV for vertebrate mitochondrial proteins.
  The present codon-based model with the ML estimates of selective constraints
and with adjustable mutation rates of nucleotide would be useful as a simple
substitution model in ML and Bayesian inferences of molecular phylogenetic
trees, and enables us to obtain biologically meaningful information at both
nucleotide and amino acid levels from codon and protein sequences.


  Statistical-mechanical lattice models for protein-DNA binding are well
established as a method to describe complex ligand binding equilibriums
measured in vitro with purified DNA and protein components. Recently, a new
field of applications has opened up for this approach since it has become
possible to experimentally quantify genome-wide protein occupancies in relation
to the DNA sequence. In particular, the organization of the eukaryotic genome
by histone proteins into a nucleoprotein complex termed chromatin has been
recognized as a key parameter that controls the access of transcription factors
to the DNA sequence. New approaches have to be developed to derive statistical
mechanical lattice descriptions of chromatin-associated protein-DNA
interactions. Here, we present the theoretical framework for lattice models of
histone-DNA interactions in chromatin and investigate the (competitive) DNA
binding of other chromosomal proteins and transcription factors. The results
have a number of applications for quantitative models for the regulation of
gene expression.


  We study the evolutionary dynamics of a maladapted population of
self-replicating sequences on strongly correlated fitness landscapes. Each
sequence is assumed to be composed of blocks of equal length and its fitness is
given by a linear combination of four independent block fitnesses. A mutation
affects the fitness contribution of a single block leaving the other blocks
unchanged and hence inducing correlations between the parent and mutant
fitness. On such strongly correlated fitness landscapes, we calculate the
dynamical properties like the number of jumps in the most populated sequence
and the temporal distribution of the last jump which is shown to exhibit a
inverse square dependence as in evolution on uncorrelated fitness landscapes.
We also obtain exact results for the distribution of records and extremes for
correlated random variables.


  We use spatial and non spatial models to argue that competition alone may
explain why two influenza strains do not usually coexist. The more virulent
strain is likely to crowd out the less virulent one. This can be seen as a
consequence of the Exclusion Principle of Ecology. We exhibit, however, a
spatial model for which coexistence is possible.


  We propose a simple model to compute the probability of success under a
quorum sensing strategy. We show that a quorum sensing strategy has a higher
probability of success than an individualistic strategy when, for instance, the
probability of success for a single individual is low and the cost of building
a quorum is not too high. On the other hand if the cost of building a quorum is
too high then the probability of success under quorum sensing always decreases
as a function of the quorum.


  Background: We study the statistical properties of fragment coverage in
genome sequencing experiments. In an extension of the classic Lander-Waterman
model, we consider the effect of the length distribution of fragments. We also
introduce the notion of the shape of a coverage function, which can be used to
detect abberations in coverage. The probability theory underlying these
problems is essential for constructing models of current high-throughput
sequencing experiments, where both sample preparation protocols and sequencing
technology particulars can affect fragment length distributions.
  Results: We show that regardless of fragment length distribution and under
the mild assumption that fragment start sites are Poisson distributed, the
fragments produced in a sequencing experiment can be viewed as resulting from a
two-dimensional spatial Poisson process. We then study the jump skeleton of the
the coverage function, and show that the induced trees are Galton-Watson trees
whose parameters can be computed.
  Conclusions: Our results extend standard analyses of shotgun sequencing that
focus on coverage statistics at individual sites, and provide a null model for
detecting deviations from random coverage in high-throughput sequence census
based experiments. By focusing on fragments, we are also led to a new approach
for visualizing sequencing data that should be of independent interest.


  Near the beginning of the century, Wright and Fisher devised an elegant,
mathematically tractable model of gene reproduction and replacement that laid
the foundation for contemporary population genetics. The Wright-Fisher model
and its extensions have given biologists powerful tools of statistical
inference that enabled the quantification of genetic drift and selection. Given
the utility of these tools, we often forget that their model - for
mathematical, and not biological reasons - makes assumptions that are violated
in most real-world populations. In this paper, I consider an alternative
framework that merges P. A. P. Moran's continuous-time Markov chain model of
allele frequency with the density dependent models of ecological competition
proposed by Gause, Lotka and Volterra, that, unlike Moran's model allow for a
stochastically varying -- but bounded -- population size. I require that allele
numbers vary according to a density-dependent population process, for which the
limiting law of large numbers is a dissipative, irreducible, competitive
dynamical system. Under the assumption that this limiting system admits a
codimension one submanifold of attractive fixed points -- a condition that
naturally generalises the weak selection regime of classical population
dynamics -- it is shown that for an appropriate rescaling of time, the finite
dimensional distributions of the original process converge to those of a
diffusion process on the submanifold. Weak convergence results are also
obtained for a related process.


  Protein phosphorylation is a reversible post-translational modification
commonly used by cell signaling networks to transmit information about the
extracellular environment into intracellular organelles for the regulation of
the activity and sorting of proteins within the cell. For this study we
reconstructed a literature-based mammalian kinase-substrate network from
several online resources. The interactions within this directed graph network
connect kinases to their substrates, through specific phosphosites including
kinase-kinase regulatory interactions. However, the "signs" of links,
activation or inhibition of the substrate upon phosphorylation, within this
network are mostly unknown. Here we show how we can infer the "signs"
indirectly using data from quantitative phosphoproteomics experiments applied
to mammalian cells combined with the literature-based kinase-substrate network.
Our inference method was able to predict the sign for 321 links and 153
phosphosites on 120 kinases, resulting in signed and directed subnetwork of
mammalian kinase-kinase interactions. Such an approach can rapidly advance the
reconstruction of cell signaling pathways and networks regulating mammalian
cells.


  Many functions have been recently defined to assess the similarity among
networks as tools for quantitative comparison. They stem from very different
frameworks - and they are tuned for dealing with different situations. Here we
show an overview of the spectral distances, highlighting their behavior in some
basic cases of static and dynamic synthetic and real networks.


  We model the immune surveillance of a pathogen which passes through $n$
immunologically distinct stages. The biological parameters of this system
induce a partial order on the stages, and this, in turn, determines which
stages will be subject to immune regulation. This corresponds to the system's
unique asymptotically stable fixed point.


  Single molecule and NMR measurements of protein dynamics increasingly uncover
the complexity of binding scenarios. Here we describe an extended
conformational selection model which embraces a repertoire of selection and
adjustment processes. Induced fit can be viewed as a subset of this repertoire,
whose contribution is affected by the bond-types stabilizing the interaction
and the differences between the interacting partners. We argue that protein
segments whose dynamics are distinct from the rest of the protein ('discrete
breathers') can govern conformational transitions and allosteric propagation
that accompany binding processes, and as such may be more sensitive to
mutational events. Additionally, we highlight the dynamic complexity of binding
scenarios as they relate to events such as aggregation and signalling, and the
crowded cellular environment.


  We analyze several aspects of the transport dynamics in the LH1-RC core of
purple bacteria, which consists basically in a ring of antenna molecules that
transport the energy into a target molecule, the reaction center, placed in the
center of the ring. We show that the periodicity of the system plays an
important role to explain the relevance of the initial state in the transport
efficiency. This picture is modified, and the transport enhanced for any
initial state, when considering that molecules have different energies, and
when including their interaction with the environment. We study this last
situation by using stochastic Schr{\"o}dinger equations, both for Markovian and
non-Markovian type of interactions.


  We present a novel statistical mechanics formalism for the theoretical
description of the process of protein folding$\leftrightarrow$unfolding
transition in water environment. The formalism is based on the construction of
the partition function of a protein obeying two-stage-like folding kinetics.
Using the statistical mechanics model of solvation of hydrophobic hydrocarbons
we obtain the partition function of infinitely diluted solution of proteins in
water environment. The calculated dependencies of the protein heat capacities
upon temperature are compared with the corresponding results of experimental
measurements for staphylococcal nuclease and metmyoglobin.


  Our models for detecting the effect of adaptation on population genomic
diversity are often predicated on a single newly arisen mutation sweeping
rapidly to fixation. However, a population can also adapt to a new situation by
multiple mutations of similar phenotypic effect that arise in parallel. These
mutations can each quickly reach intermediate frequency, preventing any single
one from rapidly sweeping to fixation globally (a "soft" sweep). Here we study
models of parallel mutation in a geographically spread population adapting to a
global selection pressure. The slow geographic spread of a selected allele can
allow other selected alleles to arise and spread elsewhere in the species
range. When these different selected alleles meet, their spread can slow
dramatically, and so form a geographic patchwork which could be mistaken for a
signal of local adaptation. This random spatial tessellation will dissipate
over time due to mixing by migration, leaving a set of partial sweeps within
the global population. We show that the spatial tessellation initially formed
by mutational types is closely connected to Poisson process models of
crystallization, which we extend. We find that the probability of parallel
mutation and the spatial scale on which parallel mutation occurs is captured by
a single characteristic length that reflects the expected distance a spreading
allele travels before it encounters a different spreading allele. This
characteristic length depends on the mutation rate, the dispersal parameter,
the effective local density of individuals, and to a much lesser extent the
strength of selection. We argue that even in widely dispersing species, such
parallel geographic sweeps may be surprisingly common. Thus, we predict, as
more data becomes available, many more examples of intra-species parallel
adaptation will be uncovered.


  The advent of X-ray Free Electron Lasers promises the possibility to
determine the structure of individual particles such as microcrystallites,
viruses and biomolecules from single-shot diffraction snapshots obtained before
the particle is destroyed by the intense femtosecond pulse. This program
requires the ability to determine the orientation of the particle giving rise
to each snapshot at signal levels as low as ~10-2 photons/pixel. Two apparently
different approaches have recently demonstrated this capability. Here we show
they represent different implementations of the same fundamental approach, and
identify the primary factors limiting their performance.


  Bilayer lipid membranes [BLMs] are an essential component of all biological
systems, forming a functional barrier for cells and organelles from the
surrounding environment. The lipid molecules that form membranes contain both
permanent and induced dipoles, and an electric field can induce the formation
of pores when the transverse field is sufficiently strong (electroporation).
Here, a phenomenological free energy is constructed to model the response of a
BLM to a transverse static electric field. The model contains a continuum
description of the membrane dipoles and a coupling between the headgroup
dipoles and the membrane tilt. The membrane is found to become unstable through
buckling modes, which are weakly coupled to thickness fluctuations in the
membrane. The thickness fluctuations, along with the increase in interfacial
area produced by membrane buckling, increase the probability of localized
membrane breakdown, which may lead to pore formation. The instability is found
to depend strongly on the strength of the coupling between the dipolar
headgroups and the membrane tilt as well as the degree of dipolar ordering in
the membrane.


  We present a new algorithm for the identification of bound regions from
ChIP-seq experiments. Our method for identifying statistically significant
peaks from read coverage is inspired by the notion of persistence in
topological data analysis and provides a non-parametric approach that is robust
to noise in experiments. Specifically, our method reduces the peak calling
problem to the study of tree-based statistics derived from the data. We
demonstrate the accuracy of our method on existing datasets, and we show that
it can discover previously missed regions and can more clearly discriminate
between multiple binding events. The software T-PIC (Tree shape Peak
Identification for ChIP-Seq) is available at
http://math.berkeley.edu/~vhower/tpic.html


  Entropy serves as a central observable in equilibrium thermodynamics.
However, many biological and ecological systems operate far from thermal
equilibrium. Here we show that entropy production can characterize the behavior
of such nonequilibrium systems. To this end we calculate the entropy production
for a population model that displays nonequilibrium behavior resulting from
cyclic competition. At a critical point the dynamics exhibits a transition from
large, limit-cycle like oscillations to small, erratic oscillations. We show
that the entropy production peaks very close to the critical point and tends to
zero upon deviating from it. We further provide analytical methods for
computing the entropy production which agree excellently with numerical
simulations.


  Electrophoretic separations of proteins are widely used in proteomic
analyses, and rely heavily on SDS electrophoresis. This mode of separation is
almost exclusively used when a single dimension separation is performed, and
generally represents the second dimension of two-dimensional separations.
Electrophoretic separations for proteomics use robust, well-established
protocols. However, many variations in almost all possible parameters have been
described in the literature over the years, and they may bring a decisive
advantage when the limits of the classical protocols are reached. The purpose
of this article is to review the most important of these variations, so that
the readers can be aware of how they can improve or tune protein separations
according to their needs. The chemical variations reviewed in this paper
encompass gel structure, buffer systems and detergents for SDS electrophoresis,
two-dimensional electrophoresis based on isoelectric focusing and
two-dimensional electrophoresis based on cationic zone electrophoresis.


  Every algorithm which can be executed on a computer can at least in principle
be realized in hardware, i.e. by a discrete physical system. The problem is
that up to now there is no programming language by which physical systems can
constructively be described. Such tool, however, is essential for the compact
description and automatic production of complex systems. This paper introduces
a programming language, called Akton-Algebra, which provides the foundation for
the complete description of discrete physical systems. The approach originates
from the finding that every discrete physical system reduces to a
spatiotemporal topological network of nodes, if the functional and metric
properties are deleted. A next finding is that there exists a homeomorphism
between the topological network and a sequence of symbols representing a
program by which the original nodal network can be reconstructed. Providing
Akton-Algebra with functionality turns it into a flow-controlled general data
processing language, which by introducing clock control and addressing can be
further transformed into a classical programming language. Providing
Akton-Algebra with metrics, i.e. the shape and size of the components, turns it
into a novel hardware system construction language.


  Motivated by possible astrophysical and biological applications we calculate
visible and near UV spectral lines of proflavine (C13H11N3,
3,6-diaminoacridine) in vacuum, as well as its anion, cation, and dication. The
pseudopotential density functional and time-dependent density functional
methods are used. We find a good agreement in spectral line positions
calculated by two real-time propagation methods and the Lanczos chain method.
Spectra of proflavine and its ions show characteristic UV lines which are good
candidates for a detection of these molecules in interstellar space and various
biological processes.


  Preliminary results are presented for the kinetics of phase separation in
three distinct models of protein aggregation. The first is a model of the
formation of spherical microcrystals of insulin via an initial formation of
fractal clusters of insulin. The results of our Brownian dynamics study of this
model are in qualitative agreement with a recent experimental study (Biophys.
Jour. 89 (2005) 3424-3433) of microcrystal formation from aqueous mixtures of
insulin. A second work involves a theory for the formation of metastable
bundles of sickle hemoglobin from fibers, based on a recent generic theory of
bundle formation (Phy. Rev. Lett. 99 (2007) 098101). We also discuss a model
for the microscopic formation of these fibers. Finally, we discuss preliminary
results for the kinetics of cluster formation for a six patch model of protein
crystallization.


  Starting from the hypothesis that the tubulin dimer is a conformationally
bistable molecule - fluctuating between a curved and a straight configuration
at room temperature - we develop a model for polymorphic dynamics of the
microtubule lattice. We show that tubulin bistability consistently explains
unusual dynamic fluctuations, the apparent length-stiffness relation of grafted
microtubules and the curved-helical appearance of microtubules in general.
Analyzing experimental data we conclude that taxol stabilized microtubules
exist in highly cooperative yet strongly fluctuating helical states. When
clamped by the end the microtubule undergoes an unusual zero energy motion - in
its effect reminiscent of a limited rotational hinge.


  Linearized catalytic reaction equations modeling e.g. the dynamics of genetic
regulatory networks under the constraint that expression levels, i.e. molecular
concentrations of nucleic material are positive, exhibit nontrivial dynamical
properties, which depend on the average connectivity of the reaction network.
In these systems the inflation of the edge of chaos and multi-stability have
been demonstrated to exist. The positivity constraint introduces a nonlinearity
which makes chaotic dynamics possible. Despite the simplicity of such minimally
nonlinear systems, their basic properties allow to understand fundamental
dynamical properties of complex biological reaction networks. We analyze the
Lyapunov spectrum, determine the probability to find stationary oscillating
solutions, demonstrate the effect of the nonlinearity on the effective in- and
out-degree of the active interaction network and study how the frequency
distributions of oscillatory modes of such system depend on the average
connectivity.


  Ecological communities exhibit pervasive patterns and inter-relationships
between size, abundance, and the availability of resources. We use scaling
ideas to develop a unified, model-independent framework for understanding the
distribution of tree sizes, their energy use and spatial distribution in
tropical forests. We demonstrate that the scaling of the tree crown at the
individual level drives the forest structure when resources are fully used. Our
predictions match perfectly with the scaling behaviour of an exactly solvable
self-similar model of a forest and are in good accord with empirical data. The
range, over which pure power law behaviour is observed, depends on the
available amount of resources. The scaling framework can be used for assessing
the effects of natural and anthropogenic disturbances on ecosystem structure
and functionality.


  Adaptive networks have been recently introduced in the context of disease
propagation on complex networks. They account for the mutual interaction
between the network topology and the states of the nodes. Until now, existing
models have been analyzed using low-complexity analytic formalisms, revealing
nevertheless some novel dynamical features. However, current methods have
failed to reproduce with accuracy the simultaneous time evolution of the
disease and the underlying network topology. In the framework of the adaptive
SIS model of Gross et al. [Phys. Rev. Lett. 96, 208701 (2006)], we introduce an
improved compartmental formalism able to handle this coevolutionary task
successfully. With this approach, we analyze the interplay and outcomes of both
dynamical elements, process and structure, on adaptive networks featuring
different degree distributions at the initial stage.


  We formulate a mathematical model for daily activities of a cow (eating,
lying down, and standing) in terms of a piecewise affine dynamical system. We
analyze the properties of this bovine dynamical system representing the single
animal and develop an exact integrative form as a discrete-time mapping. We
then couple multiple cow "oscillators" together to study synchrony and
cooperation in cattle herds. We comment on the relevant biology and discuss
extensions of our model. With this abstract approach, we not only investigate
equations with interesting dynamics but also develop interesting biological
predictions. In particular, our model illustrates that it is possible for cows
to synchronize \emph{less} when the coupling is increased.


  Analytical description of propagation phenomena on random networks has
flourished in recent years, yet more complex systems have mainly been studied
through numerical means. In this paper, a mean-field description is used to
coherently couple the dynamics of the network elements (nodes, vertices,
individuals...) on the one hand and their recurrent topological patterns
(subgraphs, groups...) on the other hand. In a SIS model of epidemic spread on
social networks with community structure, this approach yields a set of ODEs
for the time evolution of the system, as well as analytical solutions for the
epidemic threshold and equilibria. The results obtained are in good agreement
with numerical simulations and reproduce random networks behavior in the
appropriate limits which highlights the influence of topology on the processes.
Finally, it is demonstrated that our model predicts higher epidemic thresholds
for clustered structures than for equivalent random topologies in the case of
networks with zero degree correlation.


  Transcriptional activity has been shown to relate to the organization of
chromosomes in the eukaryotic nucleus and in the bacterial nucleoid. In
particular, highly transcribed genes, RNA polymerases and transcription factors
gather into discrete spatial foci called transcription factories. However, the
mechanisms underlying the formation of these foci and the resulting topological
order of the chromosome remain to be elucidated. Here we consider a
thermodynamic framework based on a worm-like chain model of chromosomes where
sparse designated sites along the DNA are able to interact whenever they are
spatially close-by. This is motivated by recurrent evidence that there exists
physical interactions between genes that operate together. Three important
results come out of this simple framework. First, the resulting formation of
transcription foci can be viewed as a micro-phase separation of the interacting
sites from the rest of the DNA. In this respect, a thermodynamic analysis
suggests transcription factors to be appropriate candidates for mediating the
physical interactions between genes. Next, numerical simulations of the polymer
reveal a rich variety of phases that are associated with different topological
orderings, each providing a way to increase the local concentrations of the
interacting sites. Finally, the numerical results show that both
one-dimensional clustering and periodic location of the binding sites along the
DNA, which have been observed in several organisms, make the spatial
co-localization of multiple families of genes particularly efficient.


  Current computational methods for exon-intron structure prediction from a
cluster of transcript (EST, mRNA) data do not exhibit the time and space
efficiency necessary to process large clusters of over than 20,000 ESTs and
genes longer than 1Mb. Guaranteeing both accuracy and efficiency seems to be a
computational goal quite far to be achieved, since accuracy is strictly related
to exploiting the inherent redundancy of information present in a large
cluster. We propose a fast method for the problem that combines two ideas: a
novel algorithm of proved small time complexity for computing spliced
alignments of a transcript against a genome, and an efficient algorithm that
exploits the inherent redundancy of information in a cluster of transcripts to
select, among all possible factorizations of EST sequences, those allowing to
infer splice site junctions that are highly confirmed by the input data. The
EST alignment procedure is based on the construction of maximal embeddings that
are sequences obtained from paths of a graph structure, called Embedding Graph,
whose vertices are the maximal pairings of a genomic sequence T and an EST P.
The procedure runs in time linear in the size of P, T and of the output.
PIntron, the software tool implementing our methodology, is able to process in
a few seconds some critical genes that are not manageable by other gene
structure prediction tools. At the same time, PIntron exhibits high accuracy
(sensitivity and specificity) when compared with ENCODE data. Detailed
experimental data, additional results and PIntron software are available at
http://www.algolab.eu/PIntron.


  We use Small Angle Neutron Scattering (SANS), with an original analysis
method, to obtain both the characteristic sizes and the inner composition of
lysozyme-pectin complexes depending on the charge density. Lysozyme is a
globular protein and pectin a natural anionic semiflexible polysaccharide with
a degree of methylation (DM) 0, 43 and 74. For our experimental conditions
(buffer ionic strength I = 2.5 10-2 mol/L and pH between 3 and 7), the
electrostatic charge of lysozyme is always positive (from 8 to 17 depending on
pH). The pectin charge per elementary chain segment is negative and can be
varied from almost zero to one through the change of DM and pH. The weight
molar ratio of lysozyme on pectin monomers is kept constant. The ratio of
negative charge content per volume to positive charge content per volume, -/+,
is varied between 10 and 0.007. On a local scale, for all charged pectins, a
correlation peak appears at 0.2 {\AA}-1 due to proteins clustering inside the
complexes. On a large scale, the complexes appear as formed of spherical
globules with a well defined radius of 10 to 50 nm, containing a few thousands
proteins. The volume fraction Phi of organic matter within the globules derived
from SANS absolute cross-sections is around 0.1. The protein stacking, which
occurs inside the globules, is enhanced when pectin is more charged, due to pH
or DM. The linear charge density of the pectin determines the size of the
globules for pectin chains of comparable molecular weights whether it is
controlled by the pH or the DM. The radius of the globules varies between 10 nm
and 50 nm. In conclusion the structure is driven by electrostatic interactions
and not by hydrophobic interactions. The molecular weight also has a large
influence on the structure of the complexes since long chains tend to form
larger globules. This maybe one reason why DM and pH are not completely
equivalent in our system since DM 0 has a short mass, but this may not be the
only one. For very low pectin charge (-/+ = 0.07), globules do not appear and
the scattering signals a gel-like structure. We did not observe any
beads-on-a-string structure.


  Binding, lateral diffusion and exchange are fundamental dynamic processes
involved in protein association with cellular membranes. In this study, we
developed numerical simulations of lateral diffusion and exchange of
fluorophores in membranes with arbitrary bleach geometry and exchange of the
membrane localized fluorophore with the cytosol during Fluorescence Recovery
after Photobleaching (FRAP) experiments. The model simulations were used to
design FRAP experiments with varying bleach region sizes on plasma-membrane
localized wild type GFP-Ras2 with a dual lipid anchor and mutant GFP-Ras2C318S
with a single lipid anchor in live yeast cells to investigate diffusional
mobility and the presence of any exchange processes operating in the time scale
of our experiments. Model parameters estimated using data from FRAP experiments
with a 1 micron x 1 micron bleach region-of-interest (ROI) and a 0.5 micron x
0.5 micron bleach ROI showed that GFP-Ras2, single or dual lipid modified,
diffuses as single species with no evidence of exchange with a cytoplasmic
pool. This is the first report of Ras2 mobility in yeast plasma membrane. The
methods developed in this study are generally applicable for studying diffusion
and exchange of membrane associated fluorophores using FRAP on commercial
confocal laser scanning microscopes.


  Using modern technology, it is now common to survey microbial communities by
sequencing DNA or RNA extracted in bulk from a given environment. Comparative
methods are needed that indicate the extent to which two communities differ
given data sets of this type. UniFrac, a method built around a somewhat ad hoc
phylogenetics-based distance between two communities, is one of the most
commonly used tools for these analyses. We provide a foundation for such
methods by establishing that if one equates a metagenomic sample with its
empirical distribution on a reference phylogenetic tree, then the weighted
UniFrac distance between two samples is just the classical
Kantorovich-Rubinstein (KR) distance between the corresponding empirical
distributions. We demonstrate that this KR distance and extensions of it that
arise from incorporating uncertainty in the location of sample points can be
written as a readily computable integral over the tree, we develop $L^p$
Zolotarev-type generalizations of the metric, and we show how the p-value of
the resulting natural permutation test of the null hypothesis "no difference
between the two communities" can be approximated using a functional of a
Gaussian process indexed by the tree. We relate the $L^2$ case to an ANOVA-type
decomposition and find that the distribution of its associated Gaussian
functional is that of a computable linear combination of independent $\chi_1^2$
random variables.


  Viral quasispecies can be regarded as a swarm of genetically related mutants
or a quasispecies (QS). A common formalism to approach QS is the
replicator-mutator equation (RME). However, a problem with the RME is how to
quantify the interaction coefficients between viral variants. Here, this is
addressed by adopting an ecological perspective and resorting to the niche
theory of competing communities, which assumes that the utilization of
resources primarily determines ecological segregation between competing
individuals (the different viral variants that constitute the QS). Using this
novel combination of RME plus the ecological concept of niche overlapping, for
describing QS, we explore the population distributions of viral variants that
emerge, as well as the corresponding dynamics. We observe that the population
distribution requires very long transients both to A) reach equilibrium and B)
to show a clear dominating master sequence. Based on different independent and
recent experimental evidence, we find that when some cooperation or
facilitation between variants is included in appropriate doses we can solve
both A) and B). We show that a useful quantity to calibrate the degree of
cooperation is the Shannon entropy. Therefore, in order to get a typical
quasispecies profile, it seems that pure competition is not enough. Rather,
some degree of cooperation among viral variants is needed. This has several
biological implications that might contribute to shed light on the mechanisms
operating in QS dynamics and to understand the QS as a whole entity.


  We study the evolution of cooperation in spatial public goods games with four
competing strategies: cooperators, defectors, punishing cooperators, and
punishing defectors. To explore the robustness of the cooperation-promoting
effect of costly punishment, besides the usual strategy adoption dynamics we
also apply strategy mutations. As expected, frequent mutations create kind of
well-mixed conditions, which support the spreading of defectors. However, when
the mutation rate is small, the final stationary state does not significantly
differ from the state of the mutation-free model, independently of the values
of the punishment fine and cost. Nevertheless, the mutation rate affects the
relaxation dynamics. Rare mutations can largely accelerate the spreading of
costly punishment. This is due to the fact that the presence of defectors
breaks the balance of power between both cooperative strategies, which leads to
a different kind of dynamics.


  Since Darwin, species trees have been used as a simplified description of the
relationships which summarize the complicated network $N$ of reality. Recent
evidence of hybridization and lateral gene transfer, however, suggest that
there are situations where trees are inadequate. Consequently it is important
to determine properties that characterize networks closely related to $N$ and
possibly more complicated than trees but lacking the full complexity of $N$.
  A connected surjective digraph map (CSD) is a map $f$ from one network $N$ to
another network $M$ such that every arc is either collapsed to a single vertex
or is taken to an arc, such that $f$ is surjective, and such that the inverse
image of a vertex is always connected. CSD maps are shown to behave well under
composition. It is proved that if there is a CSD map from $N$ to $M$, then
there is a way to lift an undirected version of $M$ into $N$, often with added
resolution. A CSD map from $N$ to $M$ puts strong constraints on $N$.
  In general, it may be useful to study classes of networks such that, for any
$N$, there exists a CSD map from $N$ to some standard member of that class.


  The concept of a protein diffusing in its free energy folding landscape has
been fruitful for both theory and experiment. Yet the choice of the reaction
coordinate (RC) introduces an undesirable degree of arbitrariness into the
problem. We analyze extensive simulation data of an alpha-helix in explicit
water solvent as it stochastically folds and unfolds. The free energy profiles
for different RCs exhibit significant variation, some having an activation
barrier, others not. We show that this variation has little effect on the
predicted folding kinetics if the diffusivity profiles are properly taken into
account. This kinetic quasi-universality is rationalized by an RC rescaling,
which, due to the reparameterization invariance of the Fokker-Planck equation,
allows the combination of free energy and diffusivity effects into a single
function, the rescaled free energy profile. This rescaled free energy indeed
shows less variation among different RCs than the bare free energy and
diffusivity profiles separately do, if we properly distinguish between RCs that
contain knowledge of the native state and those that are purely geometric in
nature. Our method for extracting diffusivity profiles is easily applied to
experimental single molecule time series data and might help to reconcile
conflicts that arise when comparing results from different experimental probes
for the same protein.


  We study the stochastic susceptible-infected-recovered (SIR) model with
time-dependent forcing using analytic techniques which allow us to disentangle
the interaction of stochasticity and external forcing. The model is formulated
as a continuous time Markov process, which is decomposed into a deterministic
dynamics together with stochastic corrections, by using an expansion in inverse
system size. The forcing induces a limit cycle in the deterministic dynamics,
and a complete analysis of the fluctuations about this time-dependent solution
is given. This analysis is applied when the limit cycle is annual, and after a
period-doubling when it is biennial. The comprehensive nature of our approach
allows us to give a coherent picture of the dynamics which unifies past work,
but which also provides a systematic method for predicting the periods of
oscillations seen in whooping cough and measles epidemics.


  Understanding under what conditions interacting populations, whether they be
plants, animals, or viral particles, coexist is a question of theoretical and
practical importance in population biology. Both biotic interactions and
environmental fluctuations are key factors that can facilitate or disrupt
coexistence. To better understand this interplay between these deterministic
and stochastic forces, we develop a mathematical theory extending the nonlinear
theory of permanence for deterministic systems to stochastic difference and
differential equations. Our condition for coexistence requires that there is a
fixed set of weights associated with the interacting populations and this
weighted combination of populations' invasion rates is positive for any
(ergodic) stationary distribution associated with a subcollection of
populations. Here, an invasion rate corresponds to an average per-capita growth
rate along a stationary distribution. When this condition holds and there is
sufficient noise in the system, we show that the populations approach a unique
positive stationary distribution. Moreover, we show that our coexistence
criterion is robust to small perturbations of the model functions. Using this
theory, we illustrate that (i) environmental noise enhances or inhibits
coexistence in communities with rock-paper-scissor dynamics depending on
correlations between interspecific demographic rates, (ii) stochastic variation
in mortality rates has no effect on the coexistence criteria for discrete-time
Lotka-Volterra communities, and (iii) random forcing can promote genetic
diversity in the presence of exploitative interactions.


  The past decade has seen a revived interest in the unavoidable or intrinsic
noise in biochemical and genetic networks arising from the finite copy number
of the participating species. That is, rather than modeling regulatory networks
in terms of the deterministic dynamics of concentrations, we model the dynamics
of the probability of a given copy number of the reactants in single cells.
Most of the modeling activity of the last decade has centered on stochastic
simulation of individual realizations, i.e., Monte-Carlo methods for generating
stochastic time series. Here we review the mathematical description in terms of
probability distributions, introducing the relevant derivations and
illustrating several cases for which analytic progress can be made either
instead of or before turning to numerical computation.


  We introduce a theory of sequential causal inference in which learners in a
chain estimate a structural model from their upstream teacher and then pass
samples from the model to their downstream student. It extends the population
dynamics of genetic drift, recasting Kimura's selectively neutral theory as a
special case of a generalized drift process using structured populations with
memory. We examine the diffusion and fixation properties of several drift
processes and propose applications to learning, inference, and evolution. We
also demonstrate how the organization of drift process space controls fidelity,
facilitates innovations, and leads to information loss in sequential learning
with and without memory.


  The importance of stochasticity within biological systems has been shown
repeatedly during the last years and has raised the need for efficient
stochastic tools. We present SABRE, a tool for stochastic analysis of
biochemical reaction networks. SABRE implements fast adaptive uniformization
(FAU), a direct numerical approximation algorithm for computing transient
solutions of biochemical reaction networks. Biochemical reactions networks
represent biological systems studied at a molecular level and these reactions
can be modeled as transitions of a Markov chain. SABRE accepts as input the
formalism of guarded commands, which it interprets either as continuous-time or
as discrete-time Markov chains. Besides operating in a stochastic mode, SABRE
may also perform a deterministic analysis by directly computing a mean-field
approximation of the system under study. We illustrate the different
functionalities of SABRE by means of biological case studies.


  This paper presents new results about the optimization based generation of
chemical reaction networks (CRNs) of higher deficiency. Firstly, it is shown
that the graph structure of the realization containing the maximal number of
reactions is unique if the set of possible complexes is fixed. Secondly, a
mixed integer programming based numerical procedure is given for computing a
realization containing the minimal/maximal number of complexes. Moreover, the
linear inequalities corresponding to full reversibility of the CRN realization
are also described. The theoretical results are illustrated on meaningful
examples.


  In evolutionary dynamics, the probability that a mutation spreads through the
whole population, having arisen in a single individual, is known as the
fixation probability. In general, it is not possible to find the fixation
probability analytically given the mutant's fitness and the topological
constraints that govern the spread of the mutation, so one resorts to
simulations instead. Depending on the topology in use, a great number of
evolutionary steps may be needed in each of the simulation events, particularly
in those that end with the population containing mutants only. We introduce two
techniques to accelerate the determination of the fixation probability. The
first one skips all evolutionary steps in which the number of mutants does not
change and thereby reduces the number of steps per simulation event
considerably. This technique is computationally advantageous for some of the
so-called layered networks. The second technique, which is not restricted to
layered networks, consists of aborting any simulation event in which the number
of mutants has grown beyond a certain threshold value, and counting that event
as having led to a total spread of the mutation. For large populations, and
regardless of the network's topology, we demonstrate, both analytically and by
means of simulations, that using a threshold of about 100 mutants leads to an
estimate of the fixation probability that deviates in no significant way from
that obtained from the full-fledged simulations. We have observed speedups of
two orders of magnitude for layered networks with 10000 nodes.


  We present a physical mechanism that can cause the mitotic spindle to
oscillate. The driving force for this mechanism emerges from the polymerization
of astral microtubules interacting with the cell cortex. We show that Brownian
ratchet model for growing microtubules reaching the cell cortex, mediate an
effective mass to the spindle body and therefore force it to oscillate. We
compare the predictions of this mechanism with the previous mechanisms which
were based on the effects of motor proteins. Finally we combine the effects of
microtubules polymerization and motor proteins, and present the detailed phase
diagram for possible oscillating states.


  The Guard-Guardee model for plant immunity describes how resistance proteins
(guards) in host cells monitor host target proteins (guardees) that are
manipulated by pathogen effector proteins. A recently suggested extension of
this model includes decoys, which are duplicated copies of guardee proteins,
and which have the sole function to attract the effector and, when modified by
the effector, trigger the plant immune response. Here we present a
proof-of-principle model for the functioning of decoys in plant immunity,
quantitatively developing this experimentally-derived concept. Our model links
the basic cellular chemistry to the outcomes of pathogen infection and
resulting fitness costs for the host. In particular, the model allows
identification of conditions under which it is optimal for decoys to act as
triggers for the plant immune response, and of conditions under which it is
optimal for decoys to act as sinks that bind the pathogen effectors but do not
trigger an immune response.


  A variety of genome transformations can occur as a microbial population
adapts to a large environmental change. In particular, genomic surveys indicate
that, following the transition to an obligate, host-dependent symbiont, the
density of transposons first rises, then subsequently declines over
evolutionary time. Here, we show that these observations can be accounted for
by a class of generic stochastic models for the evolution of genomes in the
presence of continuous selection and gene duplication. The models use a fitness
function that allows for partial contributions from multiple gene copies, is an
increasing but bounded function of copy number, and is optimal for one fully
adapted gene copy. We use Monte Carlo simulation to show that the dynamics
result in an initial rise in gene copy number followed by a subsequent fall due
to adaptation to the new environmental parameters. These results are robust for
reasonable gene duplication and mutation parameters when adapting to a novel
target sequence. Our model provides a generic explanation for the dynamics of
microbial transposon density following a large environmental changes such as
host restriction.