  We establish a connection between finite fields and finite dynamical systems.
We show how this connection can be used to shed light on some problems in
finite dynamical systems and in particular, in linear systems.


  The present work is devoted to describe a set of rules explaining the
discriminating versus non-discriminating behavior of the di-basic stages and to
characterize the role of each base in determining such a behavior. Bases are
analyze as dual entities characterized by its chemical type and the number of H
bonds involved in the codon-anticodon interaction. A codon is characterized as
an asymmetric informative entity whose global informative capacity results from
the spatially organized combinatory of the 6 proper-ties assigned by the 3
bases.


  There are many different models--both continuous and discrete--used to
describe gene mutation fixation. In particular, the Moran process, the Kimura
equation and the replicator dynamics are all well known models, that might lead
to different conclusions. We present a discussion of a unified framework to
embrace all these models, in the large population regime.


  We propose a method for detecting differential gene expression that exploits
the correlation between genes. Our proposal averages the univariate scores of
each feature with the scores in correlation neighborhoods. In a number of real
and simulated examples, the new method often exhibits lower false discovery
rates than simple t-statistic thresholding. We also provide some analysis of
the asymptotic behavior of our proposal. The general idea of
correlation-sharing can be applied to other prediction problems involving a
large number of correlated features. We give an example in protein mass
spectrometry.


  A pedigree is a directed graph in which each vertex (except the founder
vertices) has two parents. The main result in this paper is a construction of
an infinite family of counter examples to a reconstruction problem on
pedigrees, thus negatively answering a question of Steel and Hein. Some
positive reconstruction results are also presented. The problem of counting
distinct (mutually non-isomorphic) pedigrees is considered. The known lower and
upper bounds on the number of pedigrees are improved upon, and their relevance
to pedigree reconstruction from DNA sequence data is discussed. It is shown
that the information theoretic bound on the number of segregating sites in the
sequence data that is minimally essential for reconstructing pedigrees would
not significantly change with improved enumerative estimates.


  A substantial focus of research in molecular biology are gene regulatory
networks: the set of transcription factors and target genes which control the
involvement of different biological processes in living cells. Previous
statistical approaches for identifying gene regulatory networks have used gene
expression data, ChIP binding data or promoter sequence data, but each of these
resources provides only partial information. We present a Bayesian hierarchical
model that integrates all three data types in a principled variable selection
framework. The gene expression data are modeled as a function of the unknown
gene regulatory network which has an informed prior distribution based upon
both ChIP binding and promoter sequence data. We also present a variable
weighting methodology for the principled balancing of multiple sources of prior
information. We apply our procedure to the discovery of gene regulatory
relationships in Saccharomyces cerevisiae (Yeast) for which we can use several
external sources of information to validate our results. Our inferred
relationships show greater biological relevance on the external validation
measures than previous data integration methods. Our model also estimates
synergistic and antagonistic interactions between transcription factors, many
of which are validated by previous studies. We also evaluate the results from
our procedure for the weighting for multiple sources of prior information.
Finally, we discuss our methodology in the context of previous approaches to
data integration and Bayesian variable selection.


  Genes are often regulated in living cells by proteins called transcription
factors (TFs) that bind directly to short segments of DNA in close proximity to
specific genes. These binding sites have a conserved nucleotide appearance,
which is called a motif. Several recent studies of transcriptional regulation
require the reduction of a large collection of motifs into clusters based on
the similarity of their nucleotide composition. We present a principled
approach to this clustering problem based upon a Bayesian hierarchical model
that accounts for both within- and between-motif variability. We use a
Dirichlet process prior distribution that allows the number of clusters to vary
and we also present a novel generalization that allows the core width of each
motif to vary. This clustering model is implemented, using a Gibbs sampling
strategy, on several collections of transcription factor motif matrices. Our
clusters provide a means by which to organize transcription factors based on
binding motif similarities, which can be used to reduce motif redundancy within
large databases such as JASPAR and TRANSFAC. Finally, our clustering procedure
has been used in combination with discovery of evolutionarily-conserved motifs
to predict co-regulated genes. An alternative to our Dirichlet process prior
distribution is explored but shows no substantive difference in the clustering
results for our datasets. Our Bayesian clustering model based on the Dirichlet
process has several advantages over traditional clustering methods that could
make our procedure appropriate and useful for many clustering applications.


  This paper discusses the problem of identifying differentially expressed
groups of genes from a microarray experiment. The groups of genes are
externally defined, for example, sets of gene pathways derived from biological
databases. Our starting point is the interesting Gene Set Enrichment Analysis
(GSEA) procedure of Subramanian et al. [Proc. Natl. Acad. Sci. USA 102 (2005)
15545--15550]. We study the problem in some generality and propose two
potential improvements to GSEA: the maxmean statistic for summarizing
gene-sets, and restandardization for more accurate inferences. We discuss a
variety of examples and extensions, including the use of gene-set scores for
class predictions. We also describe a new R language package GSA that
implements our ideas.


  The reconstruction of large phylogenetic trees from data that violates
clocklike evolution (or as a supertree constructed from any m input trees)
raises a difficult question for biologists - how can one assign relative dates
to the vertices of the tree? In this paper we investigate this problem,
assuming a uniform distribution on the order of the inner vertices of the tree
(which includes, but is more general than, the popular Yule distribution on
trees). We derive fast algorithms for computing the probability that (i) any
given vertex in the tree was the j--th speciation event (for each j), and (ii)
any one given vertex is earlier in the tree than a second given vertex. We show
how the first algorithm can be used to calculate the expected length of any
given interior edge in any given tree that has been generated under either a
constant-rate speciation model, or the coalescent model.


  In a phylogenetic tree, we often don't have information about the time a
speciation event (inner node) occured. Under a neutral model for speciation, I
develop fast algorithms for calculating the probability that an inner node i is
the k-th speciation event. For the Yule and the coalescent model, I develop an
edge length estimation as well. Various properties of the Yule model are
discussed throughout the thesis.


  We consider a random fitness landscape on the space of haploid diallelic
genotypes with n genetic loci, where each genotype is considered either
inviable or viable depending on whether or not there are any incompatibilities
among its allele pairs. We suppose that each allele pair in the set of all
possible allele pairs on the n loci is independently incompatible with
probability p=c/(2n). We examine the connectivity of the viable genotypes under
single locus mutations and show that, for 0<c<1, the number of clusters of
viable genotypes in this landscape converges weakly (in n) to N=2^{Psi} where
Psi is Poisson distributed; while for c>1, there are no viable genotypes with
probability converging to one. The genotype space is equivalent to the
n-dimensional hypercube and the viable genotypes are solutions to a random
2-SAT problem, so the same result holds for the connectivity of solutions in
the hypercube to a random 2-SAT problem.


  We consider a branching model introduced by Kimmel for cell division with
parasite infection. Cells contain proliferating parasites which are shared
randomly between the two daughter cells when they divide. We determine the
probability that the organism recovers, meaning that the asymptotic proportion
of contaminated cells vanishes. We study the tree of contaminated cells, give
the asymptotic number of contaminated cells and the asymptotic proportions of
contaminated cells with a given number of parasites. This depends on domains
inherited from the behavior of branching processes in random environment (BPRE)
and given by the bivariate value of the means of parasite offsprings. In one of
these domains, the convergence of proportions holds in probability, the limit
is deterministic and given by the Yaglom quasistationary distribution.
Moreover, we get an interpretation of the limit of the Q-process as the
size-biased quasistationary distribution.


  Bertoin and Le Gall (2003) introduced a certain probability measure valued
Markov process that describes the evolution of a population, such that a sample
from this population would exhibit a genealogy given by the so-called
$\Lambda$-coalescent, or coalescent with multiple collisions, introduced
independently by Pitman (1999) and Sagitov (1999). We show how this process can
be extended to the case where lineages can experience mutations. Regenerative
compositions enter naturally into this model, which is somewhat surprising,
considering a negative result by M{\"o}hle (2007).


  We address the problem of studying the toric ideals of phylogenetic
invariants for a general group-based model on an arbitrary claw tree. We focus
on the group $\mathbb Z_2$ and choose a natural recursive approach that extends
to other groups. The study of the lattice associated with each phylogenetic
ideal produces a list of circuits that generate the corresponding lattice
ideal. In addition, we describe explicitly a quadratic lexicographic
Gr\"{o}bner basis of the toric ideal of invariants for the claw tree on an
arbitrary number of leaves. Combined with a result of Sturmfels and Sullivant,
this implies that the phylogenetic ideal of every tree for the group $\mathbb
Z_2$ has a quadratic Gr\"{o}bner basis. Hence, the coordinate ring of the toric
variety is a Koszul algebra.


  The Kimura 3-parameter model on a tree of n leaves is one of the most used in
phylogenetics. The affine algebraic variety W associated to it is a toric
variety. We study its geometry and we prove that it is isomorphic to a
geometric quotient of the affine space by a finite group acting on it. As a
consequence, we are able to study the singularities of W and prove that the
biologically meaningful points are smooth points. Then we give an algorithm for
constructing a set of minimal generators of the localized ideal at these
points, for an arbitrary number of leaves n. This leads to a major improvement
of phylogenetic reconstruction methods based on algebraic geometry.


  In this paper, we will analyze the behavior of the Neighbor Joining algorithm
on five taxa and we will show that the partition of the sample (data) space for
estimation of a tree topology with five taxa into subspaces, within each of
which the Neighbor Joining algorithm returns the same tree topology. A key of
our method to partition the sample space is the action of the symmetric group
$S_5$ on the set of distance matrices by changing the labels of leaves. The
method described in this paper can be generalized to trees with more than five
taxa.


  This work deals with two problems arising in mathematical ecology. The first
problem is concerned with diploid branching particle models and its behavior
when rapid stirring is added to the interaction. The particle models involve
two types of particles, male and female, and branching can only occur when both
types of particles are present. We show that if the branching rate is
sufficiently large, this particle model has a nontrivial stationary
distribution, i.e. one that does not concentrate all weight on the all-0 state,
using a comparison argument due to R. Durrett. We also show extinction for
small branching rates, thereby establishing the existence of a phase
transition. We then add two different rapid stirring mechanisms to the
interactions and show that for the particle models with rapid stirring, there
also exist nontrivial stationary distribution(s); for this, we analyze the
limiting PDE and establish a condition on the PDE that guarantees existence of
nontrivial stationary distributions for sufficient fast stirring.
  The second problem deals with a model of sympatric speciation, i.e.
speciation in the absence of geographical separation, originally proposed by U.
Dieckmann and M. Doebeli in 1999. We modify their original model to obtain
several constant-population particle models. We concentrate on a
continuous-time model that converges to a deterministic dynamical system as the
number of particles becomes large. We establish various results regarding
whether speciation occurs by studying the existence of bimodal stationary
distributions for the limiting dynamical system.


  Recently, Adami and coworkers have been able to measure the information
content of digital organisms living in their {\em Avida} artificial life
system. They show that over time, the organisms behave like Maxwell's demon,
accreting information (or complexity) as they evolve. In {\em Avida} the
organisms don't interact with each other, merely reproduce at a particular rate
(their fitness), and attempt to evaluate an externally given arithmetic
function in order win bonus fitness points. Measuring the information content
of a digital organism is essentially a process of counting the number of
genotypes that give rise to the same phenotype.
  Whilst Avidan organisms have a particularly simple phenotype, Tierran
organisms interact with each other, giving rise to an ecology of phenotypes. In
this paper, I discuss techniques for comparing pairs of Tierran organisms to
determine if they are phenotypically equivalent. I then discuss a method for
computing an estimate of the number of phenotypically equivalent genotypes that
is more accurate than the ``hot site'' estimate used by Adami's group. Finally,
I report on an experimental analysis of a Tierra run.


  At Alife VI, Mark Bedau proposed some evolutionary statistics as a means of
classifying different evolutionary systems. Ecolab, whilst not an artificial
life system, is a model of an evolving ecology that has advantages of
mathematical tractability and computational simplicity. The Bedau statistics
are well defined for Ecolab, and this paper reports statistics measured for
typical Ecolab runs, as a function of mutation rate. The behaviour ranges from
class 1 (when mutation is switched off), through class 3 at intermediate
mutation rates (corresponding to scale free dynamics) to class 2 at high
mutation rates. The class 3/class 2 transition corresponds to an error
threshold. Class 4 behaviour, which is typified by the Biosphere, is
characterised by unbounded growth in diversity. It turns out that Ecolab is
governed by an inverse relationship between diversity and connectivity, which
also seems likely of the Biosphere. In Ecolab, the mutation operator is
conservative with respect to connectivity, which explains the boundedness of
diversity. The only way to get class 4 behaviour in Ecolab is to develop an
evolutionary dynamics that reduces connectivity of time.


  AIDS is the pandemic of our era. A disease that scares us not only because it
is fatal but also because its insidious time course makes us all potential
carriers long before it hands us our heads in a basket. The strange three stage
dynamics of aids is also one of the major puzzles in describing the disease
theoretically (1). Aids starts, like most diseases, in a peak of virus
expression (2,3), which is practically wiped out by the immune system. However
it then remains in the body at a low level of expression until later (some time
years later) when there is an outbreak of the disease which terminally cripples
the immune system causing death from various common pathogens. In this paper we
show, using a microscopic simulation, that the time course of AIDS is
determined by the interactions of the virus and the immune cells in the shape
space of antigens and that it is the virus's ability to move more rapidly in
this space (it's high mutability) that causes the time course and eventual
'victory' of the disease. These results open the way for further experimental
and therapeutic conclusions in the ongoing battle with the HIV epidemic.


  A dynamical model of an ecological community is analyzed within a "mean-field
approximation" in which one of the species interacts with the combination of
all of the other species in the community. Within this approximation the model
may be formulated as a master equation describing a one-step stochastic
process. The stationary distribution is obtained in closed form and is shown to
reduce to a logseries or lognormal distribution, depending on the values that
the parameters describing the model take on. A hyperbolic relationship between
the connectance of the matrix of interspecies interactions and the average
number of species, exists for a range of parameter values. The time evolution
of the model at short and intermediate times is analyzed using van Kampen's
approximation, which is valid when the number of individuals in the community
is large. Good agreement with numerical simulations is found. The large time
behavior, and the approach to the stationary state, is obtained by solving the
equation for the generating function of the probability distribution. The
analytical results which follow from the analysis are also in good agreement
with direct simulations of the model.


  Simulations of the coevolution of many interacting species are performed
using the Webworld model. The model has a realistic set of predator-prey
equations that describe the population dynamics of the species for any
structure of the food web. The equations account for competition between
species for the same resources, and for the diet choice of predators between
alternative prey according to an evolutionarily stable strategy. The set of
species present undergoes long-term evolution due to speciation and extinction
events. We summarize results obtained on the macro-evolutionary dynamics of
speciations and extinctions, and on the statistical properties of the food webs
that are generated by the model. Simulations begin from small numbers of
species and build up to larger webs with relatively constant species number on
average. The rate of origination and extinction of species are relatively high,
but remain roughly balanced throughout the simulations. When a `parent' species
undergoes speciation, the `child' species usually adds to the same trophic
level as the parent. The chance of the child species surviving is significantly
higher if the parent is on the second or third trophic level than if it is on
the first level, most likely due to a wider choice of possible prey for species
on higher levels. Addition of a new species sometimes causes extinction of
existing species. The parent species has a high probability of extinction
because it has strong competition with the new species. Non-parental
competitors of the new species also have a significantly higher extinction
probability than average, as do prey of the new species. Predators of the new
species are less likely than average to become extinct.


  A mechanism of sympatric speciation is presented based on the
interaction-induced developmental plasticity of phenotypes. First, phenotypes
of individuals with identical genotypes split into a few groups, according to
instability in the developmental dynamics that are triggered with the
competitive interaction among individuals. Then, through mutational change of
genes, the phenotypic differences are fixed to genes, until the groups are
completely separated in genes as well as phenotypes. It is also demonstrated
that the proposed theory leads to hybrid sterility under sexual recombination,
and thus speciation is completed in the sense of reproductive isolation. As a
result of this post-mating isolation, the mating preference evolves later. When
there are two alleles, the correlation between alleles is formed, to
consolidate the speciation. When individuals are located in space, different
species are later segregated spatially, implying that the speciation so far
regarded to be allopatric may be a result of sympatric speciation.
Relationships with previous theories, frequency-dependent selection,
reinforcement, Baldwin's effect, phenotypic plasticity, and resource
competition are briefly discussed. Relevance of the results to natural
evolution are discussed, including punctuated equilibrium, incomplete
penetrance in mutants, and the change in flexibility in genotype-phenotype
correspondence. Finally, it is discussed how our theory is confirmed both in
field and in laboratory (in an experiment with the use of E coli.).


  The origin of heredity is studied as a recursive state in a replicating
proto-cell consisting of many molecule species in mutually catalyzing reaction
networks. Protocells divide when the number of molecules, increasing due to
replication, exceeds a certain threshold. We study how the chemicals in a
catalytic network can form recursive production states in the presence of
errors in the replication process. Depending on the balance between the total
number of molecules in a cell and the number of molecule species, we have found
three phases; a phase without a recursive production state, a phase with
itinerancy over a few recursive states, and a phase with fixed recursive
production states. Heredity is realized in the latter two phases where molecule
species that are population-wise in the minority are preserved and control the
phenotype of the cell. It is shown that evolvability is realized in the
itinerancy phase, where a change in the number of minority molecules controls a
change of the chemical state.


  A generalization of the standard susceptible-infectious-removed (SIR)
stochastic model for epidemics in sparse random networks is introduced which
incorporates contact tracing in addition to random screening. We propose a
deterministic mean-field description which yields quantitative agreement with
stochastic simulations on random graphs. We also analyze the role of contact
tracing in epidemics control in small-world networks and show that its
effectiveness grows as the rewiring probability is reduced.


  We study the ABC model (A + B --> 2B, B + C --> 2C, C + A --> 2A), and its
counterpart: the three--component neutral drift model (A + B --> 2A or 2B, B +
C --> 2B or 2C, C + A --> 2C or 2A.) In the former case, the mean field
approximation exhibits cyclic behaviour with an amplitude determined by the
initial condition. When stochastic phenomena are taken into account the
amplitude of oscillations will drift and eventually one and then two of the
three species will become extinct. The second model remains stationary for all
initial conditions in the mean field approximation, and drifts when stochastic
phenomena are considered. We analyzed the distribution of first extinction
times of both models by simulations and from the point of view of the
Fokker-Planck equation. Survival probability vs. time plots suggest an
exponential decay. For the neutral model the extinction rate is inversely
proportional to the system size, while the cyclic model exhibits anomalous
behaviour for small system sizes. In the large system size limit the extinction
times for both models will be the same. This result is compatible with the
smallest eigenvalue obtained from the numerical solution of the Fokker-Planck
equation. We also studied the long--time behaviour of the probability
distribution. The exponential decay is found to be robust against certain
changes, such as the three reactions having different rates.


  We address the problem of the relative importance of the intrinsic chaos and
the external noise in determining the complexity of population dynamics. We use
a recently proposed method for studying the complexity of nonlinear random
dynamical systems. The new measure of complexity is defined in terms of the
average number of bits per time-unit necessary to specify the sequence
generated by the system. This measure coincides with the rate of divergence of
nearby trajectories under two different realizations of the noise. In
particular, we show that the complexity of a nonlinear time-series model
constructed from sheep populations comes completely from the environmental
variations. However, in other situations, intrinsic chaos can be the crucial
factor. This method can be applied to many other systems in biology and
physics.


  Bedau has developed a general set of evolutionary statistics that quantify
the adaptive component of evolutionary processes. On the basis of these
measures, he has proposed a set of 4 classes of evolutionary system. All
artificial life sytems so far looked at fall into the first 3 classes, whereas
the biosphere, and possibly the human economy belongs to the 4th class. The
challenge to the artificial life community is to identify exactly what is
difference between these natural evolutionary systems, and existing artificial
life systems. At ALife VII, I presented a study using an artificial
evolutionary ecology called \EcoLab. Bedau's statistics captured the
qualitative behaviour of the model. \EcoLab{} exhibited behaviour from the
first 3 classes, but not class 4, which is characterised by unbounded growth in
diversity. \EcoLab{} exhibits a critical surface given by an inverse
relationship between connectivity and diversity, above which the model cannot
tarry long. Thus in order to get unbounded diversity increase, there needs to
be a corresponding connectivity reducing (or food web pruning) process. This
paper reexamines this question in light of two possible processes that reduce
ecosystem connectivity: a tendency for specialisation and increase in
biogeographic zones through continental drift.


  Of all the issues discussed at {\em Alife VII: Looking Forward,
  Looking Backward}, the issue of whether it was possible to create an
artificial life system that exhibits {\em open-ended evolution} of novelty is
by far the biggest. Of the 14 open problems settled on as a result of debate at
the conference, some 6 are directly, or indirectly related to this issue. Most
people equate open-ended evolution with complexity growth, although a priori
these seem to be different things. In this paper I report on experiments to
measure the complexity of Tierran organisms, and show the results for a {\em
size-neutral} run of Tierra. In this run, no increase in organismal complexity
was observed, although organism size did increase through the run. This result
is discussed, offering some signposts on path to solving the issue of open
ended evolution.


  We investigate the formation process of nonlinear vibrational modes
representing broad H-bridge multi--site breathers in a DNA--shaped double
strand.
  Within a network model of the double helix we take individual motions of the
bases within the base pair plane into account. The resulting H-bridge
deformations may be asymmetric with respect to the helix axis. Furthermore the
covalent bonds may be deformed distinctly in the two backbone strands.
  Unlike other authors that add different extra terms we limit the interaction
to the hydrogen bonds within each base pair and the covalent bonds along each
strand. In this way we intend to make apparent the effect of the characteristic
helicoidal structure of DNA. We study the energy exchange processes related
with the relaxation dynamics from a non-equilibrium conformation. It is
demonstrated that the twist-opening relaxation dynamics of a radially distorted
double helix attains an equilibrium regime characterized by a multi-site
H-bridge breather.


  We present long Monte Carlo simulations of a simple model of biological
macroevolution in which births, deaths, and mutational changes in the genome
take place at the level of individual organisms. The model displays punctuated
equilibria and flicker noise with a 1/f-like power spectrum, consistent with
some current theories of evolutionary dynamics.


  We study numerically the mechanical stability and elasticity properties of
duplex DNA molecules within the frame of a network model incorporating
microscopic degrees of freedom related with the arrangement of the base pairs.
We pay special attention to the opening-closing dynamics of double-stranded DNA
molecules which are forced into non-equilibrium conformations. Mechanical
stress imposed at one terminal end of the DNA molecule brings it into a
partially opened configuration. We examine the subsequent relaxation dynamics
connected with energy exchange processes between the various degrees of freedom
and structural rearrangements leading to complete recombination to the
double-stranded conformation. The similarities and differences between the
relaxation dynamics for a planar ladder-like DNA molecule and a twisted one are
discussed in detail. In this way we show that the attainment of a
quasi-equilibrium regime proceeds faster in the case of the twisted DNA form
than for its thus less flexible ladder counterpart. Furthermore we find that
the velocity of the complete recombination of the DNA molecule is lower than
the velocity imposed by the forcing unit which is in compliance with the
experimental observations for the opening-closing cycle of DNA molecules.


  We present a study by linear stability analysis and large-scale Monte Carlo
simulations of a simple model of biological coevolution. Selection is provided
through a reproduction probability that contains quenched, random interspecies
interactions, while genetic variation is provided through a low mutation rate.
Both selection and mutation act on individual organisms. Consistent with some
current theories of macroevolutionary dynamics, the model displays
intermittent, statistically self-similar behavior with punctuated equilibria.
The probability density for the lifetimes of ecological communities is well
approximated by a power law with exponent near -2, and the corresponding power
spectral densities show 1/f noise (flicker noise) over several decades. The
long-lived communities (quasi-steady states) consist of a relatively small
number of mutualistically interacting species, and they are surrounded by a
``protection zone'' of closely related genotypes that have a very low
probability of invading the resident community. The extent of the protection
zone affects the stability of the community in a way analogous to the height of
the free-energy barrier surrounding a metastable state in a physical system.
Measures of biological diversity are on average stationary with no discernible
trends, even over our very long simulation runs of approximately 3.4x10^7
generations.


  We investigate a class of continuum models for the motion of a
two-dimensional biological group under the influence of nonlocal social
interactions. The dynamics may be uniquely decomposed into incompressible
motion and potential motion. When the motion is purely incompressible, the
model possesses solutions which have constant population density and sharp
boundaries for all time. Numerical simulations of these "swarm patches'' reveal
rotating mill-like swarms with circular cores and spiral arms. The sign of the
social interaction term determines the direction of the rotation, and the
interaction length scale affects the degree of spiral formation. When the
motion is purely potential, the social interaction term has the meaning of
repulsion or attraction depending on its sign. For the repulsive case, the
population spreads and the density profile is smoothed. With increasing
interaction length scale, the motion becomes more convective and experiences
slower diffusive smoothing. For the attractive case, the population
self-organizes into regions of high and low density. The characteristic length
scale of the density pattern is predicted and confirmed by numerical
simulations.


  We study the ABC model in the cyclic competition and neutral drift versions,
with mutations and migrations introduced into the model. When stochastic
phenomena are taken into account, there are three distinct regimes in the
model. (i) In the "fixation" regime, the first extinction time scales with the
system size N and has an exponential distribution, with an exponent that
depends on the mutation/migration probability per particle. (ii) In the
"diversity" regime, the order parameter remains nonzero for very long times,
and becomes zero only rarely, almost never for large system sizes. (iii) In the
critical regime, the first extinction time has a power-law distribution with
exponent -1. The transition corresponds to a crossover from diffusive behaviour
to Gaussian fluctuations about a stable solution. The analytical results are
checked against computer simulations of the model.


  Some recent results for a three--dimensional, semi--classical, tight--binding
model for DNA show that there are two types of polarons, named radial and twist
polarons, that can transport charge along the DNA molecule. However, the
existence of two types of base pairs in real DNA, makes it crucial to find out
if charge transport also exist in DNA chains with different base pairs. In this
paper we address this problem in its simple case, an homogeneous chain except
for a single different base pair, what we call a base-pair inhomogeneity, and
its effect on charge transport. Radial polarons experience either reflection or
trapping. However, twist polarons are good candidates for charge transport
along real DNA. This transport is also very robust with respect to weak
parametric and diagonal disorder.


  We investigate the charge transport in synthetic DNA polymers built up from
single types of base pairs. In the context of a polaron-like model, for which
an electronic tight-binding system and bond vibrations of the double helix are
coupled, we present estimates for the electron-vibration coupling strengths
utilizing a quantum-chemical procedure. Subsequent studies concerning the
mobility of polaron solutions, representing the state of a localized charge in
unison with its associated helix deformation, show that the system for
poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers, respectively possess
quantitatively distinct transport properties. While the former supports
unidirectionally moving electron breathers attributed to highly efficient
long-range conductivity the breather mobility in the latter case is
comparatively restrained inhibiting charge transport. Our results are in
agreement with recent experimental results demonstrating that poly(dG)-poly(dC)
DNA molecules acts as a semiconducting nanowire and exhibits better conductance
than poly(dA)-poly(dT) ones.


  Motivation: Microarray experiments result in large scale data sets that
require extensive mining and refining to extract useful information. We have
been developing an efficient novel algorithm for nonmetric multidimensional
scaling (nMDS) analysis for very large data sets as a maximally unsupervised
data mining device. We wish to demonstrate its usefulness in the context of
bioinformatics. In our motivation is also an aim to demonstrate that
intrinsically nonlinear methods are generally advantageous in data mining.
Results: The Pearson correlation distance measure is used to indicate the
dissimilarity of the gene activities in transcriptional response of cell
cycle-synchronized human fibroblasts to serum [Iyer et al., Science vol. 283,
p83 (1999)]. These dissimilarity data have been analyzed with our nMDS
algorithm to produce an almost circular arrangement of the genes. The temporal
expression patterns of the genes rotate along this circular arrangement. If an
appropriate preparation procedure may be applied to the original data set,
linear methods such as the principal component analysis (PCA) could achieve
reasonable results, but without data preprocessing linear methods such as PCA
cannot achieve a useful picture. Furthermore, even with an appropriate data
preprocessing, the outcomes of linear procedures are not as clearcut as those
by nMDS without preprocessing.


  We find a Lax pair for the geometrically exact discrete Hamiltonian equations
for the discrete elastic rod. This is paper III of a series.


  We demonstrate how direct simulation of stochastic, individual-based models
can be combined with continuum numerical analysis techniques to study the
dynamics of evolving diseases. % Sidestepping the necessity of obtaining
explicit population-level models, the approach analyzes the (unavailable in
closed form) `coarse' macroscopic equations, estimating the necessary
quantities through appropriately initialized, short `bursts' of
individual-based dynamic simulation. % We illustrate this approach by analyzing
a stochastic and discrete model for the evolution of disease agents caused by
point mutations within individual hosts. % Building up from classical SIR and
SIRS models, our example uses a one-dimensional lattice for variant space, and
assumes a finite number of individuals. % Macroscopic computational tasks
enabled through this approach include stationary state computation, coarse
projective integration, parametric continuation and stability analysis.


  We consider a predator-prey model of planktonic population dynamics, of
excitable character, living in an open and chaotic fluid flow, i.e., a state of
fluid motion in which fluid trajectories are unbounded but a chaotic region
exists that is restricted to a localized area. Despite that excitability is a
transient phenomenon and that fluid trajectories are continuously leaving the
system, there is a regime of parameters where the excitation remains
permanently in the system, given rise to a persistent plankton bloom. This
regime is reached when the time scales associated to fluid stirring become
slower than the ones associated to biological growth.


  Anomalous kinetics of infective (e.g., autocatalytic) reactions in open,
nonhyperbolic chaotic flows are important for many applications in biological,
chemical, and environmental sciences. We present a scaling theory for the
singular enhancement of the production caused by the universal, underlying
fractal patterns. The key dynamical invariant quantities are the effective
fractal dimension and effective escape rate, which are primarily determined by
the hyperbolic components of the underlying dynamical invariant sets. The
theory is general as it includes all previously studied hyperbolic reactive
dynamics as a special case. We introduce a class of dissipative embedding maps
for numerical verification.


  We present a theory of cultural evolution based upon a renormalization group
scheme. We consider rational but cognitively limited agents who optimize their
decision making process by iteratively updating and refining the mental
representation of their natural and social environment. These representations
are built around the most important degrees of freedom of their world. Cultural
coherence among agents is defined as the overlap of mental representations and
is characterized using an adequate order parameter. As the importance of social
interactions increases or agents become more intelligent, we observe and
quantify a series of dynamic phase transitions by which cultural coherence
advances in the society. A similar phase transition may explain the so-called
"cultural explosion" in human evolution some 50,000 years ago.


  The concept of neutral evolutionary networks being a significant factor in
evolutionary dynamics was first proposed by Huynen {\em et al.} about 7 years
ago. In one sense, the principle is easy to state -- because most mutations to
an organism are deleterious, one would expect that neutral mutations that don't
affect the phenotype will have disproportionately greater representation
amongst successor organisms than one would expect if each mutation was equally
likely. So it was with great surprise that I noted neutral mutations being very
rare in a visualisation of phylogenetic trees generated in {\em Tierra}, since
I already knew that there was a significant amount of neutrality in the Tierra
genotype-phenotype map. The paper reports on an investigation into this
mystery.


  A symmetrical cubic discrete coupled logistic equation is proposed to model
the symbiotic interaction of two isolated species. The coupling depends on the
population size of both species and on a positive constant $\lambda$, named the
mutual benefit. Different dynamical regimes are obtained when the mutual
benefit is modified. For small $\lambda$, the species become extinct. For
increasing $\lambda$, the system stabilizes in a synchronized state or
oscillates in a 2 periodic orbit. For the greatest permitted values of
$\lambda$, the dynamics evolves into a quasiperiodic, into a chaotic scenario
or into extinction. The basins for these regimes are visualized as coloured
figures on the plane. These patterns suffer different change as consequence of
basins' bifurcations. The use of the critical curves let us to determine the
influence of the zones with different number of first rank preimages in those
bifurcation mechanisms.


  In this paper we try to end the debate concerning the suitability of
different updating schemes in random Boolean networks (RBNs). We quantify for
the first time loose attractors in asyncrhonous RBNs, which allows us to
analyze the complexity reduction related to different updating schemes. We also
report that all updating schemes yield very similar critical stability values,
meaning that the "edge of chaos" does not depend much on the updating scheme.
After discussion, we conclude that synchonous RBNs are justifiable theoretical
models of biological networks.


  Ecolab and Webworld are both models of evolution produced by adding evolution
to ecological equations. They differ primarily in the form of the ecological
equations. Both models are self-organised to a state where extinctions balance
speciations. However, Ecolab shows evidence of this self-organised state being
critical, whereas Webworld does not. This paper examines the self-organised
states of these two models and suggest the likely cause of the difference. Also
the lifetime distribution for a mean field version of Ecolab is computed,
showing that the fat tail of the distribution is due to coevolutionary adaption
of the species.


  The concept of neutral evolutionary networks being a significant factor in
evolutionary dynamics was first proposed by Huynen {\em et al.} about 7 years
ago. In one sense, the principle is easy to state -- because most mutations to
an organism are deleterious, one would expect that neutral mutations that don't
affect the phenotype will have disproportionately greater representation
amongst successor organisms than one would expect if each mutation was equally
likely. So it was with great surprise that I noted neutral mutations being very
rare in a visualisation of phylogenetic trees generated in {\em Tierra}, since
I already knew that there was a significant amount of neutrality in the Tierra
genotype-phenotype map. It turns out that competition for resources between
host and parasite inhibits neutral evolution.


  The concept of cross diffusion is applied to some biological systems. The
conditions for persistence and Turing instability in the presence of cross
diffusion are derived. Many examples including: predator-prey, epidemics (with
and without delay), hawk-dove-retaliate and prisoner's dilemma games are given.


  Recent papers have considered moving breathers (MBs) in DNA models including
long range interaction due to the dipole moments of the hydrogen bonds. We have
recalculated the value of the charge transfer when hydrogen bonds stretch using
quantum chemical methods which takes into account the whole nucleoside pairs.
We explore the consequences of this value on the properties of MBs, including
the range of frequencies for which they exist and their effective masses. They
are able to travel through bending points with fairly large curvatures provided
that their kinetic energy is larger than a minimum energy which depends on the
curvature. These energies and the corresponding velocities are also calculated
in function of the curvature.


  We present the first evolved solutions to a computational task within the
Neuronal Organism Evolution model (Norgev) of artificial neural network
development. These networks display a remarkable robustness to external noise
sources, and can regrow to functionality when severely damaged. In this
framework, we evolved a doubling of network functionality (double-NAND
circuit). The network structure of these evolved solutions does not follow the
logic of human coding, and instead more resembles the decentralized dendritic
connection pattern of more biological networks such as the 'C. elegans' brain.


  We consider a spatially distributed population dynamics model with excitable
predator-prey dynamics, where species propagate in space due to their taxis
with respect to each other's gradient in addition to, or instead of, their
diffusive spread. Earlier, we have described new phenomena in this model in one
spatial dimension, not found in analogous systems without taxis: reflecting and
self-splitting waves. Here we identify new phenomena in two spatial dimensions:
unusual patterns of meander of spirals, partial reflection of waves, swelling
wavetips, attachment of free wave ends to wave backs, and as a result, a novel
mechanism of self-supporting complicated spatio-temporal activity, unknown in
reaction-diffusion population models.


  A cubic discrete coupled logistic equation is proposed to model the
predator-prey problem. The coupling depends on the population size of both
species and on a positive constant $\lambda$, which could depend on the prey
reproduction rate and on the predator hunting strategy. Different dynamical
regimes are obtained when $\lambda$ is modified. For small $\lambda$, the
species become extinct. For a bigger $\lambda$, the preys survive but the
predators extinguish. Only when the prey population reaches a critical value
then predators can coexist with preys. For increasing $\lambda$, a bistable
regime appears where the populations apart of being stabilized in fixed
quantities can present periodic, quasiperiodic and chaotic oscillations.
Finally, bistability is lost and the system settles down in a steady state, or,
for the biggest permitted $\lambda$, in an invariant curve. We also present the
basins for the different regimes. The use of the critical curves lets us
determine the influence of the zones with different number of first rank
preimages in the bifurcation mechanisms of those basins.


  If one isolated species is supposed to evolve following the logistic mapping,
then we are tempted to think that the dynamics of two species can be expressed
by a coupled system of two discrete logistic equations. As three basic
relationships between two species are present in Nature, namely symbiosis,
predator-prey and competition, three different models are obtained. Each model
is a cubic two-dimensional discrete logistic-type equation with its own
dynamical properties: stationarity, periodicity, quasi-periodicity and chaos.
Furthermore, these models could be considered as the basic ingredients to
construct more complex interactions in the ecological networks.


  The eigenvalues and eigenvectors of the connectivity matrix of complex
networks contain information about its topology and its collective behavior. In
particular, the spectral density $\rho(\lambda)$ of this matrix reveals
important network characteristics: random networks follow Wigner's semicircular
law whereas scale-free networks exhibit a triangular distribution. In this
paper we show that the spectral density of hierarchical networks follow a very
different pattern, which can be used as a fingerprint of modularity. Of
particular importance is the value $\rho(0)$, related to the homeostatic
response of the network: it is maximum for random and scale free networks but
very small for hierarchical modular networks. It is also large for an actual
biological protein-protein interaction network, demonstrating that the current
leading model for such networks is not adequate.


  The goal of this tutorial is to promote interest in the study of random
Boolean networks (RBNs). These can be very interesting models, since one does
not have to assume any functionality or particular connectivity of the networks
to study their generic properties. Like this, RBNs have been used for exploring
the configurations where life could emerge. The fact that RBNs are a
generalization of cellular automata makes their research a very important
topic. The tutorial, intended for a broad audience, presents the state of the
art in RBNs, spanning over several lines of research carried out by different
groups. We focus on research done within artificial life, as we cannot exhaust
the abundant research done over the decades related to RBNs.


  Coherence evolution of two food web models can be obtained under the stirring
effect of chaotic advection. Each food web model sustains a three--level
trophic system composed of interacting predators, consumers and vegetation.
These populations compete for a common limiting resource in open flows with
chaotic advection dynamics. Here we show that two species (the top--predators)
of different colonies chaotically advected by a jet--like flow can synchronize
their evolution even without migration interaction. The evolution is
charaterized as a phase synchronization. The phase differences (determined
through the Hilbert transform) of the variables representing those species show
a coherent evolution.


  An analytic theory of species abundance patterns (SAPs) in biological
networks is presented. The theory is based on multispecies replicator dynamics
equivalent to the Lotka-Volterra equation, with diverse interspecies
interactions. Various SAPs observed in nature are derived from a single
parameter. The abundance distribution is formed like a widely observed
left-skewed lognormal distribution. As the model has a general form, the result
can be applied to similar patterns in other complex biological networks, e.g.
gene expression.


  Self-Organized Criticality (SOC) phenomena could have a significant effect on
the dynamics of ecosystems. The Bak-Sneppen (BS) model is a simple and robust
model of biological evolution that exhibits punctuated equilibrium behavior.
Here we will introduce random version of BS model. Also we generalize the
single objective BS model to a multiobjective one.


  Selection theory requires multiple, distinct, simultaneously-actualized
states. In cognition, each thought or cognitive state changes the 'selection
pressure' against which the next is evaluated; they are not simultaneously
selected amongst. Creative thought is more a matter of honing in a vague idea
through redescribing successive iterations of it from different real or
imagined perspectives; in other words, actualizing potential through exposure
to different contexts. It has been proven that the mathematical description of
contextual change of state introduces a non-Kolmogorovian probability
distribution, and a classical formalism such as selection theory cannot be
used. This paper argues that creative thought evolves not through a Darwinian
process, but a process of context-driven actualization of potential.


  A route to chaos is studied in 3-dimensional maps of logistic type.
Mechanisms of period doubling for invariant closed curves (ICC) are found for
specific 3-dimensional maps. These bifurcations cannot be observed for ICC in
the 2-dimensional case. When the parameter of the system is modified, localized
oscillations occur on the ICC that give rise to weakly chaotic rings, then to
chaotic attractors, which finally disappear by contact bifurcations. These maps
can be considered as models for the symbiotic interaction of three species.


  Two symmetrically coupled logistic equations are proposed to mimic the
competitive interaction between two species. The phenomena of coexistence,
oscillations and chaos are present in this cubic discrete system. This work,
together with two other similar ones recently published by the authors,
completes a triptych dedicated to the two-species relationships present in
Nature, namely the symbiosis, the predator-prey and the competition. These
models can be used as basic ingredients to build up more complex interactions
in the ecological networks.


  An array system of coupled maps is proposed as a model for economy evolution.
The local dynamics of each map or agent is controlled by two parameters. One of
them represents the growth capacity of the agent and the other one is a control
term representing the local environmental pressure which avoids an exponential
growth. The asymptotic state of the system evolution displays a complex
behavior. The distribution of the maps values in this final regime is of power
law type. In the model, inequality emerges as a result of the dynamical
processes taking place in the microscopic scales.


  We study sustained oscillations in two-dimensional oscillator systems driven
by Rayleigh-type negative friction. In particular we investigate the influence
of mismatch of the two frequencies. Further we study the influence of external
noise and nonlinearity of the conservative forces. Our consideration is
restricted to the case that the driving is rather weak and that the forces show
only weak deviations from radial symmetry. For this case we provide results for
the attractors and the bifurcations of the system. We show that for rational
relations of the frequencies the system develops several rotational excitations
with right/left symmetry, corresponding to limit cycles in the four-dimensional
phase space. The corresponding noisy distributions have the form of hoops or
tires in the four-dimensional space. For irrational frequency relations, as
well as for increasing strength of driving or noise the periodic excitations
are replaced by chaotic oscillations.


  Evolution depends on the possibility of successfully exploring fitness
landscapes via mutation and recombination. With these search procedures,
exploration is difficult in "rugged" fitness landscapes, where small mutations
can drastically change functionalities in an organism. Random Boolean networks
(RBNs), being general models, can be used to explore theories of how evolution
can take place in rugged landscapes; or even change the landscapes.
  In this paper, we study the effect that redundant nodes have on the
robustness of RBNs. Using computer simulations, we have found that the addition
of redundant nodes to RBNs increases their robustness. We conjecture that
redundancy is a way of "smoothening" fitness landscapes. Therefore, redundancy
can facilitate evolutionary searches. However, too much redundancy could reduce
the rate of adaptation of an evolutionary process. Our results also provide
supporting evidence in favour of Kauffman's conjecture (Kauffman, 2000, p.195).


  The improbability of a spontaneously generated self-assembling molecule has
suggested that life began with a set of simpler, collectively replicating
elements, such as an enclosed autocatalytic set of polymers (or protocell).
Since replication occurs without a self-assembly code, acquired characteristics
are inherited. Moreover, there is no strict distinction between alive and dead;
one can only infer that a protocell was alive if it replicates. These features
of early life render natural selection inapplicable to the description of its
change-of-state because they defy its underlying assumptions. Moreover, natural
selection describes only randomly generated novelty; it cannot describe the
emergence of form at the interface between organism and environment.
Self-organization is also inadequate because it is restricted to interactions
amongst parts; it too cannot account for context-driven change. A modified
version of selection theory or self-organization would not work because the
description of change-of-state through interaction with an incompletely
specified context has a completely different mathematical structure, i.e.
entails a non-Kolmogorovian probability model. It is proposed that the
evolution of early life is appropriately described as lineage transformation
through context-driven actualization of potential (CAP), with self-organized
change-of-state being a special case of no contextual influence, and
competitive exclusion of less fit individuals through a selection-like process
possibly (but not necessarily) playing a secondary role. It is argued that
natural selection played an important role in evolution only after genetically
mediated replication was established.


  Aggregated phenomena in social sciences and economics are highly dependent on
the way individuals interact. To help understanding the interplay between
socio-economic activities and underlying social networks, this paper studies a
sequential prisoner's dilemma with binary choice. It proposes an analytical and
computational insight about the role of endogenous networks in emergence and
sustainability of cooperation and exhibits an alternative to the choice and
refusal mechanism that is often proposed to explain cooperation. The study
focuses on heterogeneous equilibriums and emergence of cooperation from an
all-defector state that are the two stylized facts that this model successfully
reconstructs.


  The study on the collective dynamics of synchronization among genetic
oscillators is essential for the understanding of the rhythmic phenomena of
living organisms at both molecular and cellular levels. Genetic oscillators are
biochemical networks, which can generally be modelled as nonlinear dynamic
systems. We show in this paper that many genetic oscillators can be transformed
into Lur'e form by exploiting the special structure of biological systems. By
using control theory approach, we provide a theoretical method for analyzing
the synchronization of coupled nonidentical genetic oscillators. Sufficient
conditions for the synchronization as well as the estimation of the bound of
the synchronization error are also obtained. To demonstrate the effectiveness
of our theoretical results, a population of genetic oscillators based on the
Goodwin model are adopted as numerical examples.


  The evolution of complexity has been a central theme for Biology and
Artificial Life (Bonner, 1988; Bedau et al., 2000). Complexification has been
interpreted in different ways: as a process of diversification between evolving
units (Bonner, 1988) or as a scaling process that is related to the idea of
transitions between different levels of complexity (Smith and Szathmary, 1995).
There have been previous workshops on this topic, e.g. (Heylighen et al., 1999;
Lenaerts et al., 2002), but many open questions still remain. Trying to answer
these questions from a general perspective might not immediately produce
concrete answers for the different fields that are interested in this topic.
Consequently, this workshop is organised with a particular focus in mind: The
emergence of complexity through evolutionary mechanisms. Primarily we want to
have a discussion on evolutionary and related dynamics as mechanisms for
producing complexity. Furthermore, we want to bring together historical and
novel research in this context.


  We investigate epidemic models with spatial structure based on the cellular
automata method. The construction of the cellular automata is from the study by
Weimar and Boon about the reaction-diffusion equations [Phys. Rev. E 49, 1749
(1994)]. Our results show that the spatial epidemic models exhibit the
spontaneous formation of irregular spiral waves at large scales within the
domain of chaos. Moreover, the irregular spiral waves grow stably. The system
also shows a spatial period-2 structure at one dimension outside the domain of
chaos. It is interesting that the spatial period-2 structure will break and
transform into a spatial synchronous configuration in the domain of chaos. Our
results confirm that populations embed and disperse more stably in space than
they do in nonspatial counterparts.


  We apply random matrix theory to complex networks. We show that nearest
neighbor spacing distribution of the eigenvalues of the adjacency matrices of
various model networks, namely scale-free, small-world and random networks
follow universal Gaussian orthogonal ensemble statistics of random matrix
theory. Secondly we show an analogy between the onset of small-world behavior,
quantified by the structural properties of networks, and the transition from
Poisson to Gaussian orthogonal ensemble statistics, quantified by Brody
parameter characterizing a spectral property. We also present our analysis for
a protein-protein interaction network in budding yeast.


  Chaotic itinerancy is a universal dynamical concept that describes itinerant
motion among many different ordered states through chaotic transition in
dynamical systems. Unlike the expectation of the prevalence of chaotic
itinerancy in high-dimensional systems, we identify chaotic itinerant behavior
from a relatively simple ecological system, which consists only of two coupled
consumer-resource pairs. The system exhibits chaotic bursting activity, in
which the explosion and the shrinkage of the population alternate indefinitely,
while the explosion of one pair co-occurs with the shrinkage of the other pair.
We analyze successfully the bursting activity in the framework of chaotic
itinerancy, and find that large duration times of bursts tend to cluster in
time, allowing the effective burst prognosis. We also investigate the control
schemes on the bursting activity, and demonstrate that invoking the competitive
rise of the consumer in one pair can even elongate the burst of the other pair
rather than shorten it.


  We investigate a model of evolutionary dynamics on a smooth landscape which
features a ``mutator'' allele whose effect is to increase the mutation rate. We
show that the expected proportion of mutators far from equilibrium, when the
fitness is steadily increasing in time, is governed solely by the transition
rates into and out of the mutator state. This results is a much faster rate of
fitness increase than would be the case without the mutator allele. Near the
fitness equilibrium, however, the mutators are severely suppressed, due to the
detrimental effects of a large mutation rate near the fitness maximum. We
discuss the results of a recent experiment on natural selection of E. coli in
the light of our model.


  We propose a simple discrete stochastic model for calcium dynamics in living
cells. Specifically, the calcium concentration distribution is assumed to give
rise to a set of probabilities for the opening/closing of channels which
release calcium thereby changing those probabilities. We study this model in
one dimension, analytically in the mean-field limit of large number of channels
per site N, and numerically for small N. As the number of channels per site is
increased, the transition from a non-propagating region of activity to a
propagating one changes in nature from one described by directed percolation to
that of deterministic depinning in a spatially discrete system. Also, for a
small number of channels a propagating calcium wave can leave behind a novel
fluctuation-driven state, in a parameter range where the limiting deterministic
model exhibits only single pulse propagation.


  For taxonomic levels higher than species, the abundance distributions of
number of subtaxa per taxon tend to approximate power laws, but often show
strong deviationns from such a law. Previously, these deviations were
attributed to finite-time effects in a continuous time branching process at the
generic level. Instead, we describe here a simple discrete branching process
which generates the observed distributions and find that the distribution's
deviation from power-law form is not caused by disequilibration, but rather
that it is time-independent and determined by the evolutionary properties of
the taxa of interest. Our model predicts-with no free parameters-the
rank-frequency distribution of number of families in fossil marine animal
orders obtained from the fossil record. We find that near power-law
distributions are statistically almost inevitable for taxa higher than species.
The branching model also sheds light on species abundance patterns, as well as
on links between evolutionary processes, self-organized criticality and
fractals.


  The macroscopic curvature induced in double helical B-DNA by regularly
repeated adenine tracts (A-tracts) is a long known, but still unexplained
phenomenon. This effect plays a key role in DNA studies because it is unique in
the amount and the variety of the available experimental information and,
therefore, is likely to serve as a gate to the unknown general mechanisms of
recognition and regulation of genome sequences. The dominating idea in the
recent years was that, in general, macroscopic bends in DNA are caused by long
range electrostatic repulsion between phosphate groups when some of them are
neutralized by proximal external charges. In the case of A-tracts this may be
specifically bound solvent counterions. Here we report about molecular dynamics
simulations where a correct static curvature in a DNA fragment with phased
adenine tracts emerges spontaneously in conditions where any role of
counterions or long range electrostatic effects can be excluded.


  We investigate the evolutionary processes behind the development and
optimization of multiple threads of execution in digital organisms using the
avida platform, a software package that implements Darwinian evolution on
populations of self-replicating computer programs. The system is seeded with a
linearly executed ancestor capable only of reproducing its own genome, whereas
its underlying language has the capacity for multiple threads of execution
(i.e., simultaneous expression of sections of the genome.) We witness the
evolution to multi-threaded organisms and track the development of distinct
expression patterns. Additionally, we examine both the evolvability of
multi-threaded organisms and the level of thread differentiation as a function
of environmental complexity, and find that differentiation is more pronounced
in complex environments.


  The macroscopic curvature induced in the double helical B-DNA by regularly
repeated adenine tracts (A-tracts) is a long known, but still unexplained
phenomenon. This effect plays a key role in DNA studies because it is unique in
the amount and the variety of the available experimental information and,
therefore, is likely to serve as a gate to the unknown general mechanisms of
recognition and regulation of genome sequences. We report the results of
molecular dynamics simulations of a 25-mer B-DNA fragment with a sequence
including three A-tract phased with the helical screw. It represents the first
model system where properly directed static curvature emerges spontaneously in
conditions excluding any initial bias except the base pair sequence. The effect
has been reproduced in three independent MD trajectories of 10-20 ns with
important qualitative details suggesting that the final bent state is a strong
attractor of trajectories form a broad domain of the conformational space. The
ensemble of curved conformations, however, reveals significant microscopic
heterogeneity in contradiction to all existing theoretical models of bending.
Analysis of these unexpected observations leads to a new, significantly
different hypothesis of the possible mechanism of intrinsic bends in the double
helical DNA.


  The macroscopic curvature induced in the double helical B-DNA by regularly
repeated adenine tracts (A-tracts) plays an exceptional role in structural
studies of DNA because this effect presents the most well-documented example of
sequence specific conformational modulations. Recently, a new hypothesis of its
physical origin has been put forward, based upon the results of molecular
dynamics simulations of a 25-mer fragment with three A-tracts phased with the
helical screw. Its sequence, however, had never been encountered in
experimental studies, but was constructed empirically so as to maximize the
magnitude of bending in specific computational conditions. Here we report the
results of a similar investigation of another 25-mer B-DNA fragment now with a
natural base pair sequence found in a bent locus of a minicircle DNA. It is
shown that the static curvature of a considerable magnitude and stable
direction towards the minor grooves of A-tracts emerges spontaneously in
conditions excluding any initial bias except the base pair sequence. Comparison
of the bending dynamics of these two DNA fragments reveals both qualitative
similarities and interesting differences. The results suggest that the A-tract
induced bending obtained in simulations reproduces the natural phenomenon and
validates the earlier conclusions concerning its possible mechanism.


  In order to make a case for or against a trend in the evolution of complexity
in biological evolution, complexity needs to be both rigorously defined and
measurable. A recent information-theoretic (but intuitively evident) definition
identifies genomic complexity with the amount of information a sequence stores
about its environment. We investigate the evolution of genomic complexity in
populations of digital organisms and monitor in detail the evolutionary
transitions that increase complexity. We show that because natural selection
forces genomes to behave as a natural ``Maxwell Demon'', within a fixed
environment genomic complexity is forced to increase.


  A three-dimensional model of the genetic vocabulary is proposed in terms of
determinative degree, an introduced characteristic of nucleotide, which
reflects absolute difference between purin and pyrimidin bases in one DNA
strand. In the fremework of this model manifest symmetry and group properties
are described. It is proposed to use the determinative degree of codons to
analyze genetic texts, and to explain the biological sense of various
functional regions of DNA.


  As a response to ligand binding, transmembrane cell receptors often enhance
their clustering, or oligomerization, during the signalling process. Here we
present a statistical mechanical model which combines the aspects of clustering
and signalling. In this model, receptors float on the surface, while for two
neighboring receptors, there is an interaction energy dependent on their
conformational states. On the other hand, ligand binding of a receptor shifts
the energy difference between the two conformational states. Due to thermal
fluctuation, the effects of clustering and signalling are statistical average
quantities. This model reduces to a floating Ising model with a random field.
We calculate the signalling in a grand canonical ensemble mean field approach,
using Hubbard-Stratonovich transformation and replica method. Monte Carlo
simulations are also performed. Essential biological features are obtained in
our model.


  We investigate the relationship between the average fitness decay due to
single mutations and the strength of epistatic interactions in genetic
sequences. We observe that epistatic interactions between mutations are
correlated to the average fitness decay, both in RNA secondary structure
prediction as well as in digital organisms replicating in silico. This
correlation implies that during adaptation, epistasis and average mutational
effect cannot be optimized independently. In experiments with RNA sequences
evolving on a neutral network, the selective pressure to decrease the
mutational load then leads to a reduction of the amount of sequences with
strong antagonistic interactions between deleterious mutations in the
population.


  Various symmetries connected with purine-pyrimidine content of DNA sequences
are studied in terms of the intruduced determinative degree, a new
characteristics of nucleotide which is connected with codon usage. A
numerological explanation of CG pressure is proposed. A classification of DNA
sequences is given. Calculations with real sequences show that
purine-pyrimidine symmetry increases with growing of organization. A new small
parameter which characterizes the purine-pyrimidine symmetry breaking is
proposed for the DNA theory.


  Autocatalytic reaction system with a small number of molecules is studied
numerically by stochastic particle simulations. A novel state due to
fluctuation and discreteness in molecular numbers is found, characterized as
extinction of molecule species alternately in the autocatalytic reaction loop.
Phase transition to this state with the change of the system size and flow is
studied, while a single-molecule switch of the molecule distributions is
reported. Relevance of the results to intracellular processes are briefly
discussed.


  We study the evolution of asexual microorganisms with small mutation rate in
fluctuating environments, and develop techniques that allow us to expand the
formal solution of the evolution equations to first order in the mutation rate.
Our method can be applied to both discrete time and continuous time systems.
While the behavior of continuous time systems is dominated by the average
fitness landscape for small mutation rates, in discrete time systems it is
instead the geometric mean fitness that determines the system's properties. In
both cases, we find that in situations in which the arithmetic (resp.
geometric) mean of the fitness landscape is degenerate, regions in which the
fitness fluctuates around the mean value present a selective advantage over
regions in which the fitness stays at the mean. This effect is caused by the
vanishing genetic diffusion at low mutation rates. In the absence of strong
diffusion, a population can stay close to a fluctuating peak when the peak's
height is below average, and take advantage of the peak when its height is
above average.


  It is proposed that cancer results from the breakdown of universal control
mechanisms which developed in mutual association as part of the historical
process that brought individual cells together into multi-cellular communities.
By systematically comparing the genomes of uni-celled with multi-celled
organisms, one might be able to identify the most promising sites for
intervention aimed at restoring the damaged control mechanisms and thereby
arresting the cancer.


  A conceptual difficulty in the Hooke's-law description of ideal Gaussian
polymer-chain elasticity is sometimes apparent in analyses of experimental data
or in physical models designed to simulate the behavior of biopolymers. The
problem, the tendency of a chain to collapse in the absence of external forces,
is examined in the following examples: DNA-stretching experiments, gel
electrophoresis, and protein folding. We demonstrate that the application of a
statistical-mechanically derived repulsive force, acting between the chain
ends, whose magnitude is proportional to the absolute temperature and inversely
proportional to the scalar end separation removes this difficulty.


  We study the evolution of large but finite asexual populations evolving in
fitness landscapes in which all mutations are either neutral or strongly
deleterious. We demonstrate that despite the absence of higher fitness
genotypes, adaptation takes place as regions with more advantageous
distributions of neutral genotypes are discovered. Since these discoveries are
typically rare events, the population dynamics can be subdivided into separate
epochs, with rapid transitions between them. Within one epoch, the average
fitness in the population is approximately constant. The transitions between
epochs, however, are generally accompanied by a significant increase in the
average fitness. We verify our theoretical considerations with two analytically
tractable bitstring models.


  Microtubule dynamic instability arises from the hydrolysis of GTP bound to
the beta-monomer of the tubulin dimer. The conformational change induced by
hydrolysis is unknown, but microtubules disassemble into protofilaments of
GDP-bound tubulin that curve away from the microtubule axis. This paper
presents the unfolding of a portion of the tubulin molecule into the
microtubule interior as a plausible, unifying explanation for diverse
structural and kinetic features of microtubules. This is the first specific
structural hypothesis for the hydrolysis induced conformational change of
tubulin that simultaneously explains weakening of lateral bonds, bending about
longitudinal bonds, changes in protofilament supertwist associated with GTP
hydrolysis, structural features of GDP-tubulin double rings, faster disassembly
at higher temperatures and slower disassembly in the presence of glycerol and
deuterium oxide. The hypothesis suggests further theoretical investigations and
direct experimental tests.


  A self-organizing approach is proposed for gene finding based on the model of
codon usage for coding regions and positional preference for noncoding regions.
The symmetry between the direct and reverse coding regions is adopted for
reducing the number of parameters. Without requiring prior training, parameters
are estimated by iteration. By employing the window sliding technique and
likelihood ratio, a very accurate segmentation is obtained.


  We investigate the competition between two quasispecies residing on two
disparate neutral networks. Under the assumption that the two neutral networks
have different topologies and fitness levels, it is the mutation rate that
determines which quasispecies will eventually be driven to extinction. For
small mutation rates, we find that the quasispecies residing on the neutral
network with the lower replication rate will disappear. For higher mutation
rates, however, the faster replicating sequences may be outcompeted by the
slower replicating ones in case the connection density on the second neutral
network is sufficiently high. Our analytical results are in excellent agreement
with flow-reactor simulations of replicating RNA sequences.


  Bacterial chemotaxis is controlled by the conformational changes of the
receptors, in response to the change of the ambient chemical concentration. In
a statistical mechanical approach, the signalling due to the conformational
changes is a thermodynamic average quantity, dependent on the temperature and
the total energy of the system, including both ligand-receptor interaction and
receptor-receptor interaction. This physical theory suggests to biology a new
understanding of cooperation in ligand binding and receptor signalling
problems. How much experimental support of this approach can be obtained from
the currently available data? What are the parameter values? What is the
practical information for experiments? Here we make comparisons between the
theory and recent experimental results. Although currently comparisons can only
be semi-quantitative or qualitative, consistency is clearly shown. The theory
also helps to sort a variety of data.


  We propose a solution on the stopping criterion in segmenting inhomogeneous
DNA sequences with complex statistical patterns. This new stopping criterion is
based on Bayesian Information Criterion (BIC) in the model selection framework.
When this stopping criterion is applied to a left telomere sequence of yeast
Saccharomyces cerevisiae and the complete genome sequence of bacterium
Escherichia coli, borders of biologically meaningful units were identified
(e.g. subtelomeric units, replication origin, and replication terminus), and a
more reasonable number of domains was obtained. We also introduce a measure
called segmentation strength which can be used to control the delineation of
large domains. The relationship between the average domain size and the
threshold of segmentation strength is determined for several genome sequences.


  Previous divide-and-conquer segmentation analyses of DNA sequences do not
provide a satisfactory stopping criterion for the recursion. This paper
proposes that segmentation be considered as a model selection process. Using
the tools in model selection, a limit for the stopping criterion on the relaxed
end can be determined. The Bayesian information criterion, in particular,
provides a much more stringent stopping criterion than what is currently used.
Such a stringent criterion can be used to delineate larger DNA domains. A
relationship between the stopping criterion and the average domain size is
empirically determined, which may aid in the determination of isochore borders.


  A comprehensive, semi-quantitative model for the thermodynamics of
hydrophobic solvation is presented. The model is based on a very simple premise
suggested by the scaled particle theory and treats both solute and solvent
molecules as hard spheres. A connection between the peculiarly large
heat-capacity change for hydrophobic solvation and the large temperature
dependence of the thermal expansivity of water is found. Analysis reveals a
possible physical origin for the converging behavior of solvation entropies for
a series of homologous hydrophobic compounds. The model suggests that the low
solubility and the large heat-capacity change of hydrophobic solvation stem
from two distinct aspects of water molecules: the static geometry of the
molecule and the dynamic hydrogen bonding network, respectively.


  We introduce a model of molecular evolution in which the fitness of an
individual depends both on its own and on the parent's genotype. The model can
be solved by means of a nonlinear mapping onto the standard quasispecies model.
The dependency on the parental genotypes cancels from the mean fitness, but not
from the individual sequence concentrations. For finite populations, the
position of the error threshold is very sensitive to the influence from parent
genotypes. In addition to biological applications, our model is important for
understanding the dynamics of self-replicating computer programs.


  An approach based on using the idea of distinguished coding phase in explicit
form for identification of protein-coding regions (exons) in whole genome has
been proposed. For several genomes an optimal window length for averaging
GC-content function and calculating codon frequencies has been found.
Self-training procedure based on clustering in multidimensional space of
triplet frequencies is proposed. For visualization of data in the space of
triplet requiencies method of elastic maps was applied.


  The macroscopic curvature of double helical DNA induced by regularly repeated
adenine tracts is well-known but still puzzling. Its physical origin remains
controversial even though it is perhaps the best-documented sequence modulation
of DNA structure. We report here the results of comparative theoretical and
experimental studies of bending dynamics in 35-mer DNA fragments. This length
appears large enough for the curvature to be distinguished by gel
electrophoresis. Two DNA fragments, with identical base pair composition, but
different sequences are compared. In the first one, a single A-tract motif was
four times repeated in phase with the helical screw whereas the second sequence
was "random". Both calculations and experiments indicate that the A-tract DNA
is distinguished by the large static curvature and characteristic bending
dynamics, suggesting that the computed effect corresponds to the experimental
phenomenon. The results poorly agree with earlier views that attributed a
decisive role in DNA bending to sequence specific base pair stacking or binding
of solvent counterions, but lend additional support to the hypothesis of a
compressed frustrated state of the backbone as the principal physical cause of
the static curvature. We discuss the possible ways of experimental verification
of this hypothesis.


  T-cell activation plays a critical role in the initiation and propagation of
HIV-1 infection and yet transient activation of the immune system is a normal
response to immunization. In this study we extend a simple of model of HIV
infection to include T cell activation by vaccination. We then show that the
model can reproduce many but not all of the features of the post-tetanus
immunization rise in viral load observed and reported on by Stanley et al., [N.
Engl. J. Med., 334:1222--1230 (1996)]. Amplitudes and approximate timing of
post-immunization maximum viremia were matched in nine of twelve cases; in
patients with double peaks of nearly equal amplitude the later peaks were
matched. Patients with single early peaks in their viremia were best fit with
relatively high antigen clearance rate constants, while patients with early and
late peaks in viremia of equal or nearly equal amplitude were best-fit with
much lower clearance rates. The rate of clearance may depend on the level of
pre-existing antibody, which in turn would be expected to depend on the time
since the individual under study received a previous tetanus vaccination.


  We investigate the evolutionary dynamics of a finite population of sequences
adapting to NK fitness landscapes. We find that, unlike in the case of an
infinite population, the average fitness in a finite population is maximized at
a small but finite, rather than vanishing, mutation rate. The highest local
maxima in the landscape are visited for even larger mutation rates, close to a
transition point at which the population delocalizes (i.e., leaves the fitness
peak at which it was localized) and starts traversing the sequence space. If
the mutation rate is increased even further, the population undergoes a second
transition and loses all sensitivity to fitness peaks. This second transition
corresponds to the standard error threshold transition first described by
Eigen. We discuss the implications of our results for biological evolution and
for evolutionary optimization techniques.


  We study an autocatalytic system consisting of several interacting chemical
species. We observe a strong dependence of the concentrations of the chemicals
on the size of the system. This dependence is caused by the discrete nature of
the molecular concentrations. Two basic mechanisms responsible for them are
identified and elucidated. The relevance of the transitions to processes in
biochemical systems and in micro-reactors is briefly discussed.


  It is natural to ask how bacteria coordinate metabolism when depletion of an
essential nutrient limits their growth, and they must slow their entire rate of
biosynthesis. A major nutrient with a fluctuating abundance is nitrogen. The
growth rate of enteric bacteria under nitrogen-limiting conditions is known to
correlate with the internal concentration of free glutamine, the glutamine
pool. Here we compare the patterns of utilization of L-glutamine and
L-glutamate, the two central intermediates of nitrogen metabolism. Monomeric
precursors of all of the cell's macromolecules -- proteins, nucleic acids, and
surface polymers -- require the amide group of glutamine at the first dedicated
step of biosynthesis. This is the case even though only a minority (~12%) of
total cell nitrogen derives from glutamine. In contrast, the amino group of
glutamate, which provides the remainder of cell nitrogen, is generally required
late in biosynthetic pathways, e.g. in transaminase reactions for amino acid
synthesis. We propose that the pattern of glutamine dependence coordinates the
decrease in biosynthesis under conditions of nitrogen limitation. Hence, the
glutamine pool plays a global regulatory role in the cell.


  Random models of evolution are instrumental in extracting rates of
microscopic evolutionary mechanisms from empirical observations on genetic
variation in genome sequences. In this context it is necessary to know the
statistical properties of empirical observables (such as the local homozygosity
for instance). Previous work relies on numerical results or assumes Gaussian
approximations for the corresponding distributions. In this paper we give an
analytical derivation of the statistical properties of the local homozygosity
and other empirical observables assuming selective neutrality. We find that
such distributions can be very non-Gaussian.


  We introduce a probabilistic model for protein sliding motion along DNA
during the search of a target sequence. The model accounts for possible effects
due to sequence-dependent interaction between the nonspecific DNA and the
protein. As an example, we focus on T7 RNA-polymerase and exploit the available
information about its interaction at the promoter site in order to investigate
the influence of bacteriophage T7 DNA sequence on the dynamics of the sliding
process. Hydrogen bonds in the major groove are used as the main
sequence-dependent interaction between RNA-polymerase and DNA. The resulting
dynamical properties and the possibility of an experimental verification are
discussed in details. We show that, while at large times the process reaches a
pure diffusive regime, it initially displays a sub-diffusive behavior. The
subdiffusive regime can lasts sufficiently long to be of biological interest.


  We introduce a model of DNA sequence evolution which can account for biases
in mutation rates that depend on the identity of the neighboring bases. An
analytic solution for this class of non-equilibrium models is developed by
adopting well-known methods of nonlinear dynamics. Results are presented for
the CpG-methylation-deamination process which dominates point substitutions in
vertebrates. The dinucleotide frequencies generated by the model (using
empirically obtained mutation rates) match the overall pattern observed in
non-coding DNA. A web-based tool has been constructed to compute single- and
dinucleotide frequencies for arbitrary neighbor-dependent mutation rates.
Alsoprovided is the backward procedure to infer the mutation rates using
maximum likelihood analysis given the observed single- and dinucleotide
frequencies. Reasonable estimates of the mutation rates can be obtained very
efficiently, using generic non-coding DNA sequences as input, after masking
outlong homonucleotide subsequences. Our method is much more convenient and
versatile to use than the traditional method of deducing mutation rates by
counting mutation events in carefully chosen sequences. More generally, our
approach provides a more realistic but still tractable description of
non-coding genomic DNA, and may be used as a null model for various sequence
analysis applications.


  The mutation and selection of regulatory DNA sequences is presented as an
ideal model system of molecular evolution where genotype, phenotype, and
fitness can be explicitly and independently characterized. In this theoretical
study, we construct an explicit model for the evolution of regulatory
sequences, making use of the known biophysics of the binding of regulatory
proteins to DNA sequences, under the assumption that fitness of a sequence
depends only on its binding affinity to the regulatory protein. The model is
confined to the mean field (i.e., infinite population size) limit. Using
realistic values for all parameters, we determine the minimum fitness advantage
needed to maintain a binding sequence, demonstrating explicitly the ``error
threshold'' below which a binding sequence cannot survive the accumulated
effect of mutation over long time. The commonly observed ``fuzziness'' in
binding motifs arises naturally as a consequence of the balance between
selection and mutation in our model. In addition, we devise a simple model for
the evolution of multiple binding sequences in a given regulatory region. We
find the number of evolutionarily stable binding sequences to increase in a
step-like fashion with increasing fitness advantage, if multiple regulatory
proteins can synergistically enhance gene transcription. We discuss possible
experimental approaches to resolve open questions raised by our study.


  We study theoretical ``design principles'' for transcription factor-DNA
interaction in bacteria, focusing particularly on the statistical interaction
of the transcription factors (TF's) with the genomic background (i.e., the
genome without the target sites). We introduce and motivate the concept of
`programmability', i.e. the ability to set the threshold concentration for TF
binding over a wide range merely by mutating the binding sequence of a target
site. This functional demand, together with physical constraints arising from
the thermodynamics and kinetics of TF-DNA interaction, leads us to a narrow
range of ``optimal'' interaction parameters. We find that this parameter set
agrees well with experimental data for the interaction parameters of a few
exemplary prokaryotic TF's. This indicates that TF-DNA interaction is indeed
programmable. We suggest further experiments to test whether this is a general
feature for a large class of TF's.


  We discuss two new approaches to extract relevant biological information on
the Transcription Factors (and in particular to identify their binding
sequences) from the statistical distribution of oligonucleotides in the
upstream region of the genes. Both the methods are based on the notion of a
``regulatory network'' responsible for the various expression patterns of the
genes. In particular we concentrate on families of coregulated genes and look
for the simultaneous presence in the upstream regions of these genes of the
same set of transcription factor binding sites. We discuss two instances which
well exemplify the features of the two methods: the coregulation of glycolysis
in Drosophila melanogaster and the diauxic shift in Saccharomyces cerevisiae.


  We present a model of the infection of Hantavirus in deer mouse, Peromyscus
maniculatus, based on biological observations of the system in the North
American Southwest. The results of the analysis shed light on relevant
observations of the biological system, such as the sporadical disappearance of
the infection, and the existence of foci or ``refugia'' that perform as
reservoirs of the virus when environmental conditions are less than optimal.


  Gene regulation in eukaryotes is mainly effected through transcription
factors binding to rather short recognition motifs generally located upstream
of the coding region. We present a novel computational method to identify
regulatory elements in the upstream region of eukaryotic genes. The genes are
grouped in sets sharing an overrepresented short motif in their upstream
sequence. For each set, the average expression level from a microarray
experiment is determined: If this level is significantly higher or lower than
the average taken over the whole genome, then the overerpresented motif shared
by the genes in the set is likely to play a role in their regulation. The
method was tested by applying it to the genome of Saccharomyces cerevisiae,
using the publicly available results of a DNA microarray experiment, in which
expression levels for virtually all the genes were measured during the diauxic
shift from fermentation to respiration. Several known motifs were correctly
identified, and a new candidate regulatory sequence was determined.


  Traveling waves are analyzed in a model of the Hantavirus infection of deer
mice. The existence of two kinds of wave phenomena is predicted. An
environmental parameter governs a transition between two regimes of
propagation. In one of them the front of infection lags behind at a constant
rate. In the other, fronts of susceptible and infected mice travel at the same
speed, separated by a constant delay. The dependence of the delay on system
parameters is analyzed numerically and through a piecewise linearization.


  Due to its polyionic character the DNA double helix is stable and
biologically active only in salty aqueous media where its charge is compensated
by solvent counterions. Monovalent metal ions are ubiquitous in DNA environment
and they are usually considered as the possible driving force of
sequence-dependent modulations of DNA structure that make it recognizable by
proteins. In an effort to directly examine this hypothesis, MD simulations of
DNA in a water drop surrounded by vacuum were carried out, which relieves the
requirement of charge neutrality. Surprisingly, with zero concentration of
counterions a dodecamer DNA duplex appears metastable and its structure remains
similar to that observed in experiments.


  A lattice based model of a protein is used to study the dimerization
equilibrium of the GCN4 leucine zipper. Replica exchange Monte Carlo is used to
determine the free energy of both the monomeric and dimeric forms as a function
of temperature. The method of coincidences is then introduced to explicitly
calculate the entropy loss associated with dimerization, and from it the free
energy difference between monomer and dimer, as well as the corresponding
equilibrium reaction constant. We find that the entropy loss of dimerization is
a strong function of energy (or temperature), and that it is much larger than
previously estimated, especially for high energy states. The results confirm
that it is possible to study the dimerization equilibrium of GCN4 at
physiological concentrations within the reduced representation of the protein
employed.