  One of the main problems of drug design is that of optimizing the
drug--target interaction. In the case in which the target is a viral protein
displaying a high mutation rate, a second problem arises, namely the eventual
development of resistance. We wish to suggest a scheme for the design of
non--conventional drugs which do not face any of these problems and apply it to
the case of HIV--1 protease. It is based on the knowledge that the folding of
single--domain proteins, like e.g. each of the monomers forming the HIV--1--PR
homodimer, is controlled by local elementary structures (LES), stabilized by
local contacts among hydrophobic, strongly interacting and highly conserved
amino acids which play a central role in the folding process. Because LES have
evolved over myriads of generations to recognize and strongly interact with
each other so as to make the protein fold fast as well as to avoid aggregation
with other proteins, highly specific (and thus little toxic) as well as
effective folding--inhibitor drugs suggest themselves: short peptides (or
eventually their mimetic molecules), displaying the same amino acid sequence of
that of LES (p--LES). Aside from being specific and efficient, these inhibitors
are expected not to induce resistance: in fact, mutations which successfully
avoid their action imply the destabilization of one or more LES and thus should
lead to protein denaturation. Making use of Monte Carlo simulations within the
framework of a simple although not oversimplified model, which is able to
reproduce the main thermodynamic as well as dynamic properties of monoglobular
proteins, we first identify the LES of the HIV--1--PR and then show that the
corresponding p--LES peptides act as effective inhibitors of the folding of the
protease which do not create resistance.


  The initial analysis of the recently sequenced genome of Acanthamoeba
polyphaga Mimivirus, the largest known double-stranded DNA virus, predicted a
proteome of size and complexity more akin to small parasitic bacteria than to
other nucleo-cytoplasmic large DNA viruses, and identified numerous functions
never before described in a virus. It has been proposed that the Mimivirus
lineage could have emerged before the individualization of cellular organisms
from the 3 domains of life. An exhaustive in silico analysis of the non-coding
moiety of all known viral genomes, now uncovers the unprecedented perfect
conservation of a AAAATTGA motif in close to 50% of the Mimivirus genes. This
motif preferentially occurs in genes transcribed from the predicted leading
strand and is associated with functions required early in the viral infectious
cycle, such as transcription and protein translation. A comparison with the
known promoter of unicellular eukaryotes, in particular amoebal protists,
strongly suggests that the AAAATTGA motif is the structural equivalent of the
TATA box core promoter element. This element is specific to the Mimivirus
lineage, and may correspond to an ancestral promoter structure predating the
radiation of the eukaryotic kingdoms. This unprecedented conservation of core
promoter regions is another exceptional features of Mimivirus, that again
raises the question of its evolutionary origin.


  In this paper, we aim to give a tutorial for undergraduate students studying
statistical methods and/or bioinformatics. The students will learn how data
visualization can help in genomic sequence analysis. Students start with a
fragment of genetic text of a bacterial genome and analyze its structure. By
means of principal component analysis they ``discover'' that the information in
the genome is encoded by non-overlapping triplets. Next, they learn how to find
gene positions. This exercise on PCA and K-Means clustering enables active
study of the basic bioinformatics notions. Appendix 1 contains program listings
that go along with this exercise. Appendix 2 includes 2D PCA plots of triplet
usage in moving frame for a series of bacterial genomes from GC-poor to GC-rich
ones. Animated 3D PCA plots are attached as separate gif files. Topology
(cluster structure) and geometry (mutual positions of clusters) of these plots
depends clearly on GC-content.


  The discovery and genome analysis of Acanthamoeba polyphaga Mimivirus, the
largest known DNA virus, challenged much of the accepted dogma regarding
viruses. Its particle size (>400 nm), genome length (1.2 million bp) and huge
gene repertoire (911 protein coding genes) all contribute to blur the
established boundaries between viruses and the smallest parasitic cellular
organisms. Phylogenetic analyses also suggested that the Mimivirus lineage
could have emerged prior to the individualization of cellular organisms from
the three established domains, triggering a debate that can only be resolved by
generating and analyzing more data. The next step is then to seek some evidence
that Mimivirus is not the only representative of its kind and determine where
to look for new Mimiviridae. An exhaustive similarity search of all Mimivirus
predicted proteins against all publicly available sequences identified many of
their closest homologues among the Sargasso Sea environmental sequences.
Subsequent phylogenetic analyses suggested that unknown large viruses
evolutionarily closer to Mimivirus than to any presently characterized species
exist in abundance in the Sargasso Sea. Their isolation and genome sequencing
could prove invaluable in understanding the origin and diversity of large DNA
viruses, and shed some light on the role they eventually played in the
emergence of eukaryotes.


  Metabolic networks are known to be scale free but the evolutionary origin of
this structural property is not clearly understood. One way of studying the
dynamical process is to compare the metabolic networks of species that have
arisen at different points in evolution and hence are related to each other to
varying extents. We have compared the reaction sets of each metabolite across
and within 15 groups of species. For a given pair of species and a given
metabolite, the number $\Delta k$ of reactions of the metabolite that appear in
the metabolic network of only one species and not the other is a measure of the
distance between the two networks. While $\Delta k$ is small within groups of
related species and large across groups, we find its probability distribution
to be $\sim (\Delta k)^{-\gamma'}$ where $\gamma'$ is a universal exponent that
is the same within and across groups. This exponent equals, upto statistical
uncertainties, the exponent $\gamma$ in the scale free degree distribution
$\sim k^{-\gamma}$. We argue that this, as well as our finding that $\Delta k$
is approximately linearly correlated with the degree $k$ of the metabolite, is
evidence of a `proportionate change' process in evolution. We also discuss some
molecular mechanisms that might be responsible for such an evolutionary
process.


  To unveil the logic of cell from a level of chemical reaction dynamics, we
need to clarify how ensemble of chemicals can autonomously produce the set of
chemical, without assuming a specific external control echanism. A cell
consists of a huge number of chemical species that catalyze each other. Often
the number of each molecule species is not so large, and accordingly the number
fluctuations in each molecule speciescan be large. In the amidst of such
diversity and large fluctuations, how can a cell make recursive production? On
the other hand, a cell can change its state to evolve to a different type over
a longer time span. How are reproduction and evolution compatible? We address
these questions, based on several model studies with catalytic reaction
network. In the present survey paper, we first formulate basic questions on the
recursiveness and evolvability of a cell, and then state the standpoint of our
research to answer the questions, that is termed as 'constructive biology'.
Based on this standpoint, we present general strategy of modeling a cell as a
chemical reaction network.


  During division it is of primary importance for a cell to correctly determine
the site of cleavage. The bacterium Escherichia coli divides in the center,
producing two daughter cells of equal size. Selection of the center as the
correct division site is in part achieved by the Min-proteins. They oscillate
between the two cell poles and thereby prevent division at these locations.
Here, a phenomenological description for these oscillations is presented, where
lateral interactions between proteins on the cell membrane play a key role.
Solutions to the dynamic equations are compared to experimental findings. In
particular, the temporal period of the oscillations is measured as a function
of the cell length and found to be compatible with the theoretical prediction.


  The theory of island biogeography[1] asserts that an island or a local
community approaches an equilibrium species richness as a result of the
interplay between the immigration of species from the much larger metacommunity
source area and local extinction of species on the island (local community).
Hubbell[2] generalized this neutral theory to explore the expected steady-state
distribution of relative species abundance (RSA) in the local community under
restricted immigration. Here we present a theoretical framework for the unified
neutral theory of biodiversity[2] and an analytical solution for the
distribution of the RSA both in the metacommunity (Fisher's logseries) and in
the local community, where there are fewer rare species. Rare species are more
extinction-prone, and once they go locally extinct, they take longer to
re-immigrate than do common species. Contrary to recent assertions[3], we show
that the analytical solution provides a better fit, with fewer free parameters,
to the RSA distribution of tree species on Barro Colorado Island (BCI)[4] than
the lognormal distribution[5,6].


  We investigate the scaling properties of the Penna model, which has become a
popular tool for the study of population dynamics and evolutionary problems in
recent years. We find that the model generates a normalised age distribution
for which a simple scaling rule is proposed, that is able to reproduce
qualitative features for all genome sizes.


  Motivated by the results of recent laboratory experiments (Yoshida et al.
Nature, 424, 303-306 (2003)) as well as many earlier field observations that
evolutionary changes can take place in ecosystems over relatively short
ecological time scales, several ``unified'' mathematical models of evolutionary
ecology have been developed over the last few years with the aim of describing
the statistical properties of data related to the evolution of ecosystems.
Moreover, because of the availability of sufficiently fast computers, it has
become possible to carry out detailed computer simulations of these models. For
the sake of completeness and to put these recent developments in the proper
perspective, we begin with a brief summary of some older models of ecological
phenomena and evolutionary processes. However, the main aim of this article is
to review critically these ``unified'' models, particularly those published in
the physics literature, in simple language that makes the new theories
accessible to wider audience.


  We use a microscopic model of population dynamics, a modified version of the
well known Penna model, to study some aspects of microevolution. This research
is motivated by recent reports on the effect of selective hunting on the gene
pool of bighorn sheep living in the Ram Mountain region, in Canada. Our model
finds a sharp transition in the structure of the gene pool as some threshold
for the number of animals hunted is reached.


  The growth function of populations is central in biomathematics. The main
dogma is the existence of density dependence mechanisms, which can be modelled
with distinct functional forms that depend on the size of the population. One
important class of regulatory functions is the $\theta$-logistic, which
generalises the logistic equation. Using this model as a motivation, this paper
introduces a simple dynamical reformulation that generalises many growth
functions. The reformulation consists of two equations, one for population
size, and one for the growth rate. Furthermore, the model shows that although
population is density-dependent, the dynamics of the growth rate does not
depend either on population size, nor on the carrying capacity. Actually, the
growth equation is uncoupled from the population size equation, and the model
has only two parameters, a Malthusian parameter $\rho$ and a competition
coefficient $\theta$. Distinct sign combinations of these parameters reproduce
not only the family of $\theta$-logistics, but also the van Bertalanffy,
Gompertz and Potential Growth equations, among other possibilities. It is also
shown that, except for two critical points, there is a general size-scaling
relation that includes those appearing in the most important allometric
theories, including the recently proposed Metabolic Theory of Ecology. With
this model, several issues of general interest are discussed such as the growth
of animal population, extinctions, cell growth and allometry, and the effect of
environment over a population.


  We address the controversial hot question concerning the validity of the
loose coupling versus the lever-arm theories in the actomyosin dynamics by
re-interpreting and extending the phenomenological washboard potential model
proposed by some of us in a previous paper. In this new model a Brownian motion
harnessing thermal energy is assumed to co-exist with the deterministic swing
of the lever-arm, to yield an excellent fit of the set of data obtained by some
of us on the sliding of Myosin II heads on immobilized actin filaments under
various load conditions. Our theoretical arguments are complemented by accurate
numerical simulations, and the robustness of the model is tested via different
choices of parameters and potential profiles.


  Network motifs, the recurring regulatory structural patterns in networks, are
able to self-organize to produce networks. Three major motifs, feedforward
loop, single input modules and bi-fan are found in gene regulatory networks.
The large ratio of genes to transcription factors (TFs) in genomes leads to a
sharing of TFs by motifs and is sufficient to result in network
self-organization. We find a common design principle of these motifs: short
transcript's half-life (THL) TFs are significantly enriched in motifs and hubs.
This enrichment becomes one of the driving forces for the emergence of the
network scale-free topology and allows the network to quickly adapt to
environmental changes. Most feedforward loops and bi-fans contain at least one
short THL TF, which can be seen as a criterion for self-assembling these
motifs. We have classified the motifs according to their short THL TF content.
We show that the percentage of the different motif subtypes varies in different
cellular conditions.


  Recently there has been a lot of interest in identifying modules at the level
of genetic and metabolic networks of organisms, as well as in identifying
single genes and reactions that are essential for the organism. A goal of
computational and systems biology is to go beyond identification towards an
explanation of specific modules and essential genes and reactions in terms of
specific structural or evolutionary constraints. In the metabolic networks of
E. coli, S. cerevisiae and S. aureus, we identified metabolites with a low
degree of connectivity, particularly those that are produced and/or consumed in
just a single reaction. Using FBA we also determined reactions essential for
growth in these metabolic networks. We find that most reactions identified as
essential in these networks turn out to be those involving the production or
consumption of low degree metabolites. Applying graph theoretic methods to
these metabolic networks, we identified connected clusters of these low degree
metabolites. The genes involved in several operons in E. coli are correctly
predicted as those of enzymes catalyzing the reactions of these clusters. We
independently identified clusters of reactions whose fluxes are perfectly
correlated. We find that the composition of the latter `functional clusters' is
also largely explained in terms of clusters of low degree metabolites in each
of these organisms. Our findings mean that most metabolic reactions that are
essential can be tagged by one or more low degree metabolites. Those reactions
are essential because they are the only ways of producing or consuming their
respective tagged metabolites. Furthermore, reactions whose fluxes are strongly
correlated can be thought of as `glued together' by these low degree
metabolites.


  In the study of the conformational behavior of complex systems, such as
proteins, several related statistical measures are commonly used to compare two
different potential energy functions. Among them, the Pearson's correlation
coefficient r has no units and allows only semi-quantitative statements to be
made. Those that do have units of energy and whose value may be compared to a
physically relevant scale, such as the root mean square deviation (RMSD), the
mean error of the energies (ER), the standard deviation of the error (SDER) or
the mean absolute error (AER), overestimate the distance between potentials.
Moreover, their precise statistical meaning is far from clear. In this article,
a new measure of the distance between potential energy functions is defined
which overcomes the aforementioned difficulties. In addition, its precise
physical meaning is discussed, the important issue of its additivity is
investigated and some possible applications are proposed. Finally, two of these
applications are illustrated with practical examples: the study of the van der
Waals energy, as implemented in CHARMM, in the Trp-Cage protein (PDB code 1L2Y)
and the comparison of different levels of the theory in the ab initio study of
the Ramachandran map of the model peptide HCO-L-Ala-NH2.


  We discuss recent developments in the modeling of negative autoregulated
genetic networks. In particular, we consider the temporal evolution of the
population of mRNA and proteins in simple networks using rate equations. In the
limit of low copy numbers, fluctuation effects become significant and more
adequate modeling is then achieved using the master equation formalism. The
analogy between regulatory gene networks and chemical reaction networks on dust
grains in the interstellar medium is discussed. The analysis and simulation of
complex reaction networks are also considered.


  Pilot studies of structured treatment interruptions (STI) in HIV therapy have
shown that patients can maintain low viral loads whilst benefiting from reduced
treatment toxicity. However, a recent STI clinical trial reported a high degree
of virologic failure. Here we present a novel hypothesis that could explain
virologic failure to STI and provides new insights of great clinical relevance.
We analyze a classic mathematical model of HIV within-host viral dynamics and
find that nonlinear parametric resonance occurs when STI are added to the
model; resonance is observed as virologic failure. We use the model to simulate
clinical trial data and to calculate patient-specific resonant spectra. We gain
two important insights. Firstly, within an STI trial, we determine that
patients who begin with similar viral loads can be expected to show extremely
different virologic responses as a result of resonance. Thus, high
heterogeneity of patient response within a STI clinical trial is to be
expected. Secondly and more importantly, we determine that virologic failure is
not simply due to STI or patient characteristics; rather it is the result of a
complex dynamic interaction between STI and patient viral dynamics. Hence, our
analyses demonstrate that no universal regimen with periodic interruptions will
be effective for all patients. On the basis of our results, we suggest that
immunologic and virologic parameters should be used to design patient-specific
STI regimens.


  The evolutionary origin of universal statistics in biochemical reaction
network is studied, to explain the power-law distribution of reaction links and
the power-law distributions of chemical abundances. Using cell models with
catalytic reaction network, we find evidence that the power-law distribution in
abundances of chemicals emerges by the selection of cells with higher growth
speeds. Through the further evolution, this inhomogeneity in chemical
abundances is shown to be embedded in the distribution of links, leading to the
power-law distribution. These findings provide novel insights into the nature
of network evolution in living cells.


  Protein structure is generally conceptualized as the global arrangement or of
smaller, local motifs of helices, sheets, and loops. These regular, recurring
secondary structural elements have well-understood and standardized definitions
in terms of amino acid backbone geometry and the manner in which hydrogen
bonding requirements are satisfied. Recently, "tube" models have been proposed
to explain protein secondary structure in terms of the geometrically optimal
packing of a featureless cylinder. However, atomically detailed simulations
demonstrate that such packing considerations alone are insufficient for
defining secondary structure; both excluded volume and hydrogen bonding must be
explicitly modeled for helix formation. These results have fundamental
implications for the construction and interpretation of realistic and
meaningful biomacromolecular models.


  The recent algorithm for five split tRNA-genes in N.equitans is new . It
locates missing tRNA-trp, tRNA-imet, tRNA-glu and tRNA-his . But the split
tRNA-trp(CCA) solution is anomalous ; the tRNA-imet lacks cognition elements
for aminoacylation . In view therefore we present here alternate non-split
composite solutions for tRNA-trp, tRNA-imet, tRNA-glu and tRNA-his .


  Convective counterparts of variants of the nonlinear Fisher equation which
describes reaction diffusion systems in population dynamics are studied with
the help of an analytic prescription and shown to lead to interesting
consequences for the evolution of population densities. The initial value
problem is solved explicitly for some cases and for others it is shown how to
find traveling wave solutions analytically. The effect of adding diffusion to
the convective equations is first studied through exact analysis through a
piecewise linear representation of the nonlinearity. Using an appropriate small
parameter suggested by that analysis, a perturbative treatment is developed to
treat the case in which the convective evolution is augmented by a small amount
of diffusion.


  Markov chain Monte Carlo algorithms play a key role in the Bayesian approach
to phylogenetic inference. In this paper, we present the first theoretical work
analyzing the rate of convergence of several Markov chains widely used in
phylogenetic inference. We analyze simple, realistic examples where these
Markov chains fail to converge quickly. In particular, the data studied are
generated from a pair of trees, under a standard evolutionary model. We prove
that many of the popular Markov chains take exponentially long to reach their
stationary distribution. Our construction is pertinent since it is well known
that phylogenetic trees for genes may differ within a single organism. Our
results shed a cautionary light on phylogenetic analysis using Bayesian
inference and highlight future directions for potential theoretical work.


  We discuss a simple and powerful approach for the ab initio identification of
cis-regulatory motifs involved in transcriptional regulation. The method we
present integrates several elements: human-mouse comparison, statistical
analysis of genomic sequences and the concept of coregulation. We apply it to a
complete scan of the human genome. By using the catalogue of conserved upstream
sequences collected in the CORG database we construct sets of genes sharing the
same overrepresented motif (short DNA sequence) in their upstream regions both
in human and in mouse. We perform this construction for all possible motifs
from 5 to 8 nucleotides in length and then filter the resulting sets looking
for two types of evidence of coregulation: first, we analyze the Gene Ontology
annotation of the genes in the set, searching for statistically significant
common annotations; second, we analyze the expression profiles of the genes in
the set as measured by microarray experiments, searching for evidence of
coexpression. The sets which pass one or both filters are conjectured to
contain a significant fraction of coregulated genes, and the upstream motifs
characterizing the sets are thus good candidates to be the binding sites of the
TF's involved in such regulation. In this way we find various known motifs and
also some new candidate binding sites.


  The architecture of the network of protein-protein physical interactions in
Saccharomyces cerevisiae is exposed through the combination of two
complementary theoretical network measures, betweenness centrality and
`Q-modularity'. The yeast interactome is characterized by well-defined
topological modules connected via a small number of inter-module protein
interactions. Should such topological inter-module connections turn out to
constitute a form of functional coordination between the modules, we speculate
that this coordination is occurring typically in a pair-wise fashion, rather
than by way of high-degree hub proteins responsible for coordinating multiple
modules. The unique non-hub-centric hierarchical organization of the
interactome is not reproduced by gene duplication-and-divergence stochastic
growth models that disregard global selective pressures.


  Topological analysis of large scale protein-protein interaction networks
(PINs) is important for understanding the organisational and functional
principles of individual proteins. The number of interactions that a protein
has in a PIN has been observed to be correlated with its indispensability.
Essential proteins generally have more interactions than the non-essential
ones. We show here that the lethality associated with removal of a protein from
the yeast proteome correlates with different centrality measures of the nodes
in the PIN, such as the closeness of a protein to many other proteins, or the
number of pairs of proteins which need a specific protein as an intermediary in
their communications, or the participation of a protein in different protein
clusters in the PIN. These measures are significantly better than random
selection in identifying essential proteins in a PIN. Centrality measures based
on graph spectral properties of the network, in particular the subgraph
centrality, show the best performance in identifying essential proteins in the
yeast PIN. Subgraph centrality gives important structural information about the
role of individual proteins and permits the selection of possible targets for
rational drug discovery through the identification of essential proteins in the
PIN.


  The causes and implications of the regional variations in the spread of the
incipient agriculture in Europe remain poorly understood. We apply population
dynamics models to study the dispersal of the Neolithic in Europe from a
localized area in the Near East, solving the two-dimensional reaction-diffusion
equation on a spherical surface. We focus on the role of major river paths and
coastlines in the advance of farming to model the rapid advances of the Linear
Pottery (LBK) and the Impressed Ware traditions along the Danube-Rhine corridor
and the Mediterranean coastline respectively. We argue that the random walk of
individuals, which results in diffusion of the population, can be anisotropic
in those areas. The standard reaction-diffusion equation is thus supplemented
with advection-like terms confined to the proximity of major rivers and
coastlines. The model allows for the spatial variation in both the human
mobility (diffusivity) and the carrying capacity of landscapes, reflecting the
local altitude and latitude. This approach can easily be generalised to include
other environmental factors, such as the bioproductivity of landscapes. Our
model successfully accounts for the regional variations in the spread of the
Neolithic, consistent with the radiocarbon dated data, and reproduces a time
delay in the spread of farming to the Eastern Europe and Scandinavia.


  Structure predictions of helical membrane proteins have been designed to take
advantage of the structural autonomy of secondary structure elements, as
postulated by the two-stage model of Engelman and Popot. In this context, we
investigate structure calculation strategies for two membrane proteins with
different functions, sizes, aminoacid compositions, and topologies: the
glycophorin A homodimer (a paradigm for close inter-helical packing in membrane
proteins) and aquaporin (a channel protein). Our structure calculations are
based on two alternative folding schemes: a one-step simulated annealing from
an extended chain conformation, and a two-step procedure inspired by the
grid-search methods traditionally used in membrane protein predictions. In this
framework, we investigate rationales for the utilization of sparse NMR data
such as distance-based restraints and residual dipolar couplings in structure
calculations of helical membrane proteins.


  By means of the nonlinear modeling technique (NM technique), we find the
nonlinear deterministic structures in the promoter regions (PPRs) of DNA
sequences, called deterministic valleys (DVs) in this paper. These DVs prefer
to occur much more outside of a special region around TSS. The number,
positions and shapes of the DVs are basically different for different PPRs.
Generally, these DVs do not occur in the CpG islands, which tells us that they
should be special structures with new biological functions rather than the CpG
islands. PACS numbers: 87.14.Gg, 87.10.1e, 05.45.2a


  A computer simulation based on individual births and deaths gives a
biodiversity increasing from cold to warm climates, in agreement with reality.
Complexity of foodwebs increases with time and at a higher rate at low
latitudes, and there is a higher rate of species creation at low latitudes.
Keeping many niches empty makes the results correspond more closely to natural
gradients.


  A microscopic agent dynamical model for diploid age-structured populations is
used to study evolution of polymorphism and sympatric speciation. The
underlying ecology is represented by a unimodal distribution of resources of
some width. Competition among individuals is also described by a similar
distribution, and its strength is maximum for individuals with the same
phenotype and decreases with distance in phenotype space as a gaussian, with
some width. These two widths define the model's phase space, in which we
identify the regions where an autonomous emergence of stable polymorphism or
speciation is more likely.


  Proteins created by combinatorial methods in vitro are an important source of
information for understanding sequence-structure-function relationships.
Alignments of folded proteins from combinatorial libraries can be analyzed
using methods developed for naturally occurring proteins, but this neglects the
information contained in the unfolded sequences of the library. We introduce
two algorithms, logistic regression and excess information analysis, that use
both the folded and unfolded sequences and compare them against contingency
table and statistical coupling analysis, which only use the former. The test
set for this benchmark study is a library of fictitious proteins that fold
according to a hypothetical energy model. Of the four methods studied, only
logistic regression is able to correctly recapitulate the energy model from the
sequence alignment. The other algorithms predict spurious interactions between
alignment positions with strong but individual influences on protein stability.
When present in the same protein, stabilizing amino acids tend to lower the
energy below the threshold needed for folding. As a result, their frequencies
in the alignment can be correlated even if the positions do not interact. We
believe any algorithm that neglects the nonlinear relationship between folding
and energy is susceptible to this error.


  We demonstrate diversification rather than optimisation for highly
interacting organisms in a well mixed biological system by means of a simple
model and reference to experiment, and find the cause to be the complex network
of interactions formed, allowing species less well adapted to an environment to
flourish by co-interaction over the `best' species. This diversification can be
considered as the construction of many co-evolutionary niches by the network of
interactions between species. Evidence for this comes from work with the
bacteria Escherichia coli, which may coexist with their own mutants under
certain conditions. Diversification only occurs above a certain threshold
interaction strength, below which competitive exclusion occurs.


  Yet in spite of advances in hospital treatment, hospitals continue to be a
breeding ground for several airborne diseases and for diseases that are
transmitted through close contacts like SARS, methicillin-resistant
Staphylococcus aureus (MRSA), norovirus infections and tuberculosis (TB). Here
we extract contact networks for up to 295,108 inpatients for durations up to
two years from a database used for administrating a local public healthcare
system serving a population of 1.9 million individuals. Structural and
dynamical properties of the network of importance for the transmission of
contagious diseases are then analyzed by methods from network epidemiology. The
contact networks are found to be very much determined by an extreme (age
independent) variation in duration of hospital stays and the hospital
structure. We find that that the structure of contacts between in-patients
exhibit structural properties, such as a high level of transitivity,
assortativity and variation in number of contacts, that are likely to be of
importance for the transmission of less contagious diseases. If these
properties are considered when designing prevention programs the risk for and
the effect of epidemic outbreaks may be decreased.


  In this paper we present a collection of results pertaining to haplotyping.
The first set of results concerns the combinatorial problem of reconstructing
haplotypes from incomplete and/or imperfectly sequenced haplotype data. More
specifically, we show that an interesting, restricted case of Minimum Error
Correction (MEC) is NP-hard, point out problems in earlier claims about a
related problem, and present a polynomial-time algorithm for the ungapped case
of Longest Haplotype Reconstruction (LHR). Secondly, we present a polynomial
time algorithm for the problem of resolving genotype data using as few
haplotypes as possible (the Pure Parsimony Haplotyping Problem, PPH) where each
genotype has at most two ambiguous positions, thus solving an open problem
posed by Lancia et al in "Haplotyping Populations by Pure Parsimony: Complexity
of Exact and Approximation Algorithms."


  The transcription of DNA into mRNA is initiated and aided by a number of
transcription factors (TFs), proteins with DNA-binding regions that attach
themselves to binding sites in the DNA (transcription factor binding sites,
TFBSs). As it has become apparent that both TFs and TFBSs are highly variable,
tools are needed to quantify the strength of the interaction resulting from a
certain TF variant binding to a certain TFBS. We used a simple way to predict
interactions between protein and DNA: given experimental cases from the
literature where the interaction strength between two sequences has been
quantified, we created training vectors for rule induction by regarding each
amino acid and nucleotide position as a single feature in the example vector.
The resulting interaction strength was used as the target class or value. These
training vectors were then used to build a rule induction model. We applied the
rule induction method to three protein families - transcription factors from
the NF-kappaB, the early-growth-response (EGR), and the paired domain groups -
and their corresponding DNA targets. These three prediction problems increase
in complexity with regard to model building, and thus gave us a good range for
validation. The main focus of the study was the most complex problem: paired
domain-DNA binding. For this problem, we also found sequence/binding strength
correlations using measures from information theory. Prediction results were
uniformly good: the rule induction approach was able to correctly order all of
nine unseen examples for one NF-kappaB protein, achieved a correlation
coefficient of 0.52 on unseen (and noisy) examples in the EGR case, and reached
a classification accuracy of 69.7% for the Paired domain as evaluated by
cross-validation.


  This article adopts game theory to build a model for explaining the predation
behavior of animals.We assume that both the prey and the preydator have two
stratigies in this game,the active one and the passive one.By calculating the
outcome and the income of energy in different stratigies, we find the solution
to analyze the different evolution path of both the prey and the predator.A
simulation result approximately represents the correctness of our model.


  In order to extend the results obtained with minimal lattice models to more
realistic systems, we study a model where proteins are described as a chain of
20 kinds of structureless amino acids moving in a continuum space and
interacting through a contact potential controlled by a 20x20 quenched random
matrix. The goal of the present work is to design and characterize amino acid
sequences folding to the SH3 conformation, a 60-residues recognition domain
common to many regulatory proteins. We show that a number of sequences can
fold, starting from a random conformation, to within a distance root mean
square deviation (dRMSD) of 2.6A from the native state. Good folders are those
sequences displaying in the native conformation an energy lower than a
sequence--independent threshold energy.


  We recently introduced a physical model [Hoang et al., P. Natl. Acad. Sci.
USA (2004), Banavar et al., Phys. Rev. E (2004)] for proteins which
incorporates, in an approximate manner, several key features such as the
inherent anisotropy of a chain molecule, the geometrical and energetic
constraints placed by the hydrogen bonds and sterics, and the role played by
hydrophobicity. Within this framework, marginally compact conformations
resembling the native state folds of proteins emerge as broad competing minima
in the free energy landscape even for a homopolymer. Here we show how the
introduction of sequence heterogeneity using a simple scheme of just two types
of amino acids, hydrophobic (H) and polar (P), and sequence design allows a
selected putative native fold to become the free energy minimum at low
temperature. The folding transition exhibits thermodynamic cooperativity, if
one neglects the degeneracy between two different low energy conformations
sharing the same fold topology.


  In eukaryote nucleosome, DNA wraps around a histone octamer in a left-handed
way. We study the process of chirality formation of nucleosome with Brownian
dynamics simulation. We model the histone octamer with a quantitatively
adjustable chirality: left-handed, right-handed or non-chiral, and simulate the
dynamical wrapping process of a DNA molecule on it. We find that the chirality
of a nucleosome formed is strongly dependent on that of the histone octamer,
and different chiralities of the histone octamer induce its different rotation
directions in the wrapping process of DNA. In addition, a very weak chirality
of the histone octamer is quite enough for sustaining the correct chirality of
the nucleosome formed. We also show that the chirality of a nucleosome may be
broken at elevated temperature.


  We introduce and discuss a minimal individual-based model for influenza
dynamics. The model takes into account the effects of specific immunization
against viral strains, but also infectivity randomness and the presence of a
short-lived strain transcending immunity recently suggested in the literature.
We show by simulations that the resulting model exhibits substitution of viral
strains along the years, but that their divergence remains bounded. We also
show that dropping any of these features results in a drastically different
behavior, leading either to the extinction of the disease, to the proliferation
of the viral strains, or to their divergence.


  The pentaradial organisation of echinoderms is postulated to have evolved as
the result of the reorganisation of the internal U-shaped mesentery of the
intestinal tract during inflation of the trunk of a pterobranch-like ancestor.
Under this scenario, loops of the mesentery developed between five hydraulic
bulges by three different mechanisms: (1) by the early formation of two
additional loops resulting directly in five hydraulic bulges, (2) by the
subsequent formation of two additional loops resulting in the formation of
initially three and then five hydraulic bulges or (3) by the inflation of the
body without formation of loops. Accordingly there are at least three main
evolutionary pathways within the echinoderms. The anatomical structures, such
as the ambulacral system or the skeletal capsule, which characterize the
echinoderms, developed independently in both lineages. From all the three
pathways various body structures can be derived matching those found in fossil
and extant echinoderms. The eleutherozoans and the earliest pentaradial
echinoderms most likely evolved from the direct-pentaradial pathway whereas
triradiate echinoderms and those which show a pentaradial organisation
superimposed on triradial symmetry, likely evolved along the indirect pathway.
The asymmetric or bilateral symmetric echinoderms evolved from the third
paythway. The most important morphological transformations leading to the
directly and indirectly pentaradial echinoderms are discussed and described.


  Trypsin and chymotrypsin are both serine proteases with high sequence and
structural similarities, but with different substrate specificity. Previous
experiments have demonstrated the critical role of the two loops outside the
binding pocket in controlling the specificity of the two enzymes. To understand
the mechanism of such a control of specificity by distant loops, we have used
the Gaussian Network Model to study the dynamic properties of trypsin and
chymotrypsin and the roles played by the two loops. A clustering method was
introduced to analyze the correlated motions of residues. We have found that
trypsin and chymotrypsin have distinct dynamic signatures in the two loop
regions which are in turn highly correlated with motions of certain residues in
the binding pockets. Interestingly, replacing the two loops of trypsin with
those of chymotrypsin changes the motion style of trypsin to chymotrypsin-like,
whereas the same experimental replacement was shown necessary to make trypsin
have chymotrypsin's enzyme specificity and activity. These results suggest that
the cooperative motions of the two loops and the substrate-binding sites
contribute to the activity and substrate specificity of trypsin and
chymotrypsin.


  When echinoderms are conceptualized as hydraulic entities, the early
evolution of this group can be presented in a scenario which describes how a
bilateral ancestor (an enteropneust-like organism) gradually evolved into a
pentaradial echinoderm. According to this scenario, the arms are outgrowths
from the anterior/posterior body axis of the bilateral pterobranchia-like
intermediate. These outgrowths developed when the originally U-shaped mesentery
of the intestinal tract formed loops, and correspondingly, the tensile chords
of the mesentery were attached to the body wall in five loops. The wall faces
between these regions of tensile chords could bulge out under the hydraulic
pressure of the body coelom. The originally more or less round body cavity was
deformed into a pneu with five bulges. The loops of the gut forced a roughly
symmetric arrangement, which was enhanced by a physical fact: five pneus as
well as one pneu with five internal tethers, naturally adopt a pentaradial
pattern of "minimum contact surfaces", as the most economic arrangement. These
evolutionary transformations were accompanied by certain histological
modifications, such as the development of mutable connective tissues and
skeletal elements that fused to ossicles and provided shape stabilization in
the form of a calcareous skeleton in the tissues of the body wall. The
resultant organism was an ancestral eleutherozoan echinoderm (Ur-Echinoderm).


  The dynamic theory of inhomogeneous populations developed during the last
decade predicts several essential new dynamic regimes applicable even to the
well-known, simple population models. We show that, in an inhomogeneous
population with a distributed reproduction coefficient, the entire initial
distribution of the coefficient should be used to investigate real population
dynamics. In the general case, neither the average rate of growth nor the
variance or any finite number of moments of the initial distribution is
sufficient to predict the overall population growth. We developed methods for
solving the heterogeneous models and explored the dynamics of the total
population size together with the reproduction coefficient distribution. We
show that, typically, there exists a phase of hyper-exponential growth that
precedes the well-known exponential phase of population growth in a free
regime. The developed formalism is applied to models of global demography and
the problem of population explosion predicted by the known hyperbolic formula
of world population growth. We prove here that the hyperbolic formula presents
an exact solution to the Malthus model with an exponentially distributed
reproduction coefficient and that population explosion is a corollary of
certain implicit unrealistic assumptions. Alternative models of world
population growth are derived; they show a notable phenomenon, a transition
from protracted hyperbolical growth (the phase of hyper-exponential
development) to the brief transitional phase of exponential growth and,
subsequently, to stabilization. The model solutions are consistent with real
data and produce relatively accurate forecasts.


  The emergence and spreading of chirality on the early Earth is considered by
studying a set of reaction-diffusion equations based on a polymerization model.
It is found that effective mixing of the early oceans is necessary to reach the
present homochiral state. The possibility of introducing mass extinctions and
modifying the emergence rate of life is discussed.


  The differences between uni-directional and bi-directional polymerization are
considered. The uni-directional case is discussed in the framework of the RNA
world. Similar to earlier models of this type, where polymerization was assumed
to proceed in a bi-directional fashion (presumed to be relevant to peptide
nucleic acids), left-handed and right-handed monomers are produced via an
autocatalysis from an achiral substrate. The details of the bifurcation from a
racemic solution to a homochiral state of either handedness is shown to be
remarkably independent of whether the polymerization in uni-directional or
bi-directional. Slightly larger differences are seen when dissociation is
allowed and the dissociation fragments are being recycled into the achiral
substrate.


  A variety of viruses tightly pack their genetic material into protein capsids
that are barely large enough to enclose the genome. In particular, in
bacteriophages, forces as high as 60 pN are encountered during packaging and
ejection, produced by DNA bending elasticity and self-interactions. The high
forces are believed to be important for the ejection process, though the extent
of their involvement is not yet clear. As a result, there is a need for
quantitative models and experiments that reveal the nature of the forces
relevant to DNA ejection. Here we report measurements of the ejection forces
for two different mutants of bacteriophage lambda, lambda b221cI26 and lambda
cI60, which differ in genome length by ~30%. As expected for a force-driven
ejection mechanism, the osmotic pressure at which DNA release is completely
inhibited varies with the genome length: we find inhibition pressures of 15 atm
and 25 atm, respectively, values that are in agreement with our theoretical
calculations.


  A model of the dynamics of natural rotifer populations is described as a
discrete nonlinear map depending on three parameters, which reflect
characteristics of the population and environment. Model dynamics and their
change by variation of these parameters were investigated by methods of
bifurcation theory. A phase-parametric portrait of the model was constructed
and domains of population persistence (stable equilibrium, periodic and
a-periodic oscillations of population size) as well as population extinction
were identified and investigated. The criteria for population persistence and
approaches to determining critical parameter values are described. The results
identify parameter values that lead to population extinction under various
environmental conditions. They further illustrate that the likelihood of
extinction can be substantially increased by small changes in environmental
quality, which shifts populations into new dynamical regimes.


  The correlations of primary and secondary structures were analyzed using
proteins with known structure from Protein Data Bank. The correlation values of
amino acid type and the eight secondary structure types at distant position
were calculated for distances between -25 and 25. Shapes of the diagrams
indicate that amino acids polarity and capability for hydrogen bonding have
influence on the secondary structure at some distances. Clear preference of
most of the amino acids towards certain secondary structure type classifies
amino acids into four groups: alpha-helix admirers, strand admirers, turn and
bend admirers and the others. Group four consists of His and Cis, the amino
acids that do not show clear preference for any secondary structure. Amino
acids from a group have similar physicochemical properties, and the same
structural characteristics. The results suggest that amino acid preference for
secondary structure type is based on the structural characteristics at Cb and
Cg atoms of amino acid. alpha-helix admirers do not have polar heteroatoms on
Cb and Cg atoms, nor branching or aromatic group on Cb atom. Amino acids that
have aromatic groups or branching on Cb atom are strand admirers. Turn and bend
admirers have polar heteroatom on Cb or Cg atoms or do not have Cb atom at all.
Our results indicate that polarity and capability for hydrogen bonding have
influence on the secondary structure at some distance, and that amino acid
preference for secondary structure is caused by structural properties at Cb or
Cg atoms.


  We apply the computational methodology of phase retrieval to the problem of
folding heteropolymers. The ground state fold of the polymer is defined by the
intersection of two sets in the configuration space of its constituent
monomers: a geometrical chain constraint and a threshold constraint on the
contact energy. A dynamical system is then defined in terms of the projections
to these constraint sets, such that its fixed points solve the set intersection
problem. We present results for two off-lattice HP models: one with only
rotameric degrees of freedom, and one proposed by Stillinger et al. with
flexible bond angles. Our phase retrieval inspired algorithm is competitive
with more established algorithms and even finds lower energy folds for one of
the longer polymer chains.


  The coexistence of distinct templates is a common feature of the diverse
proposals advanced to resolve the information crisis of prebiotic evolution.
However, achieving robust template coexistence turned out to be such a
difficult demand that only a class of models, the so-called package models,
seems to have met it so far. Here we apply Wright's Island formulation of group
selection to study the conditions for the coexistence of two distinct template
types confined in packages (vesicles) of finite capacity. In particular, we
show how selection acting at the level of the vesicles can neutralize the
pressures towards the fixation of any one of the template types (random drift)
and of the type with higher replication rate (deterministic competition). We
give emphasis to the role of the distinct generation times of templates and
vesicles as yet another obstacle to coexistence.


  Gene duplication is key to molecular evolution in all three domains of life
and may be the first step in the emergence of new gene function. It is a well
recognized feature in large DNA viruses, but has not been studied extensively
in the largest known virus to date, the recently discovered Acanthamoeba
Polyphaga Mimivirus. Here we present a systematic analysis of gene and genome
duplication events in the Mimivirus genome. We find that one third of the
Mimivirus genes are related to at least one other gene in the Mimivirus genome,
either through a large segmental genome duplication event that occurred in the
more remote past, either through more recent gene duplication events, which
often occur in tandem. This shows that gene and genome duplication played a
major role in shaping the Mimivirus genome. Using multiple alignments together
with remote homology detection methods based on Hidden Markov Model comparison,
we assign putative functions to some of the paralogous gene families. We
suggest that a large part of the duplicated Mimivirus gene families are likely
to interfere with important host cell processes, such as transcription control,
protein degradation, and cell regulatory processes. Our findings support the
view that large DNA viruses are complex evolving organisms, possibly deeply
rooted within the tree of life, and oppose the paradigm that viral evolution is
dominated by lateral gene acquisition, at least in what concerns large DNA
viruses.


  We introduce a three-dimensional stochastic reaction-diffusion model to
describe MinD/MinE dynamical structures in Escherichia coli. This model
spontaneously generates pole-to-pole oscillations of the membrane-associated
MinD proteins, MinE ring, as well as filaments of the membrane-associated MinD
proteins. Experimental data suggest MinD filaments are two-stranded. In order
to model them we assume that each membrane-associated MinD protein can form up
to three bonds with adjacent membrane associated MinD molecules and that MinE
induced hydrolysis strongly depends on the number of bonds MinD has
established.


  The functionality of proteins is related to their structure in the native
state. Protein structures are made up of emergent building blocks of helices
and almost planar sheets. A simple coarse-grained geometrical model of a
flexible tube barely subject to compaction provides a unified framework for
understanding the common character of globular proteins.We argue that a recent
critique of the tube idea is not well founded.


  Complex systems of polynomial equations have to be set up and solved
algebraically in order to obtain analytic solutions for maximum likelihood on
phylogenetic trees. This has restricted the types of systems previously
resolved to the simplest models - three and four taxa under a molecular clock,
with just two state characters. In this work we give, for the first time,
analytic solutions for a family of trees with four state characters, like
normal DNA or RNA. The model of substitution we use is the Jukes-Cantor model,
and the trees are on three taxa under molecular clock, namely rooted triplets.
  We employ a number of approaches and tools to solve this system: Spectral
methods (Hadamard conjugation), a new representation of variables (the path-set
spectrum), and algebraic geometry tools (the resultant of two polynomials). All
these, combined with heavy application of computer algebra packages (Maple),
let us derive the desired solution.


  In most stochastic models of molecular sequence evolution the probability of
each possible pattern of homologous characters at a site is estimated
numerically. However in the case of Kimura's three-substitution-types (K3ST)
model, these probabilities can be expressed analytically by Hadamard
conjugation as a function of the phylogeny T and the substitution probabilities
on each edge of T, together with an analytic inverse function. In this paper we
produce a direct proof of these results, using pathset distances which
generalise pairwise distances between sequences. This interpretation allows us
to apply Hadamard conjugation to a number of topical problems in the
mathematical analysis of sequence evolution.


  A deterministic mutation-selection model in the sequence space approach is
investigated. Genotypes are identified with two-letter sequences. Mutation is
modelled as a Markov process, fitness functions are of Hopfield type, where the
fitness of a sequence is determined by the Hamming distances to a number of
predefined patterns. Using a maximum principle for the population mean fitness
in equilibrium, the error threshold phenomenon is studied for quadratic
Hopfield-type fitness functions with small numbers of patterns. Different from
previous investigations of the Hopfield model, the system shows error threshold
behaviour not for all fitness functions, but only for certain parameter values.


  To gain a deeper insight into cellular processes such as transcription and
translation, one needs to uncover the mechanisms controlling the
configurational changes of nucleic acids. As a step toward this aim, we present
here a novel mesoscopic-level computational model that provides a new window
into nucleic acid dynamics. We model a single-stranded nucleic as a polymer
chain whose monomers are the nucleosides. Each monomer comprises a bead
representing the sugar molecule and a pin representing the base. The bead-pin
complex can rotate about the backbone of the chain. We consider pairwise
stacking and hydrogen-bonding interactions. We use a modified Monte Carlo
dynamics that splits the dynamics into translational bead motion and rotational
pin motion. By performing a number of tests we first show that our model is
physically sound. We then focus on the study of a the kinetics of a DNA
hairpin--a single-stranded molecule comprising two complementary segments
joined by a non-complementary loop--studied experimentally. We find that
results from our simulations agree with experimental observations,
demonstrating that our model is a suitable tool for the investigation of the
hybridization of single strands.


  Much recent work has explored molecular and population-genetic constraints on
the rate of protein sequence evolution. The best predictor of evolutionary rate
is expression level, for reasons which have remained unexplained. Here, we
hypothesize that selection to reduce the burden of protein misfolding will
favor protein sequences with increased robustness to translational missense
errors. Pressure for translational robustness increases with expression level
and constrains sequence evolution. Using several sequenced yeast genomes,
global expression and protein abundance data, and sets of paralogs traceable to
an ancient whole-genome duplication in yeast, we rule out several confounding
effects and show that expression level explains roughly half the variation in
Saccharomyces cerevisiae protein evolutionary rates. We examine causes for
expression's dominant role and find that genome-wide tests favor the
translational robustness explanation over existing hypotheses that invoke
constraints on function or translational efficiency. Our results suggest that
proteins evolve at rates largely unrelated to their functions, and can explain
why highly expressed proteins evolve slowly across the tree of life.


  We analyse forty seven chloroplastid genes of the large subunit of RuBisCO,
from the Algal order Ectocarpales, sourced from GenBank. Codon-usage weighted
by the nucleotide base bias defines our score called the
Codon-Impact-Parameter. This score is used to obtain phylogenetic relations
amongst the 47 Ectocarpales. We compare our classification with the ones done
earlier.


  We present a model that considers the maturation of the antibody population
following primary antigen presentation as a global optimization problem. The
trade-off that emerges from our model describes the balance between the safety
of mutations that lead to local improvements in affinity and the necessity of
the system to undergo global reconfigurations in the antibody's shape in order
to achieve its goals, in this example of fast-paced evolution. The parameter p
which quantifies this trade-off appears to be itself both robust and evolvable.
This parallels the rapidity and consistency of the optimization operating
during the biologic response. In this paper, we explore the robust qualities
and evolvability of this tunable control parameter, p.


  We formulate the dynamical model for the anti-tumor immune response based on
intercellular cytokine-mediated interactions with the interleukin-2 (IL-2)
taken into account. The analysis shows that the expression level of tumor
antigens on antigen presenting cells has a distinct influence on the tumor
dynamics. At low antigen presentation a progressive tumor growth takes place to
the highest possible value. At high antigen presentation there is a decrease in
tumor size to some value when the dynamical equilibrium between the tumor and
the immune system is reached. In the case of the medium antigen presentation
both these regimes can be realized depending on the initial tumor size and the
condition of the immune system. A pronounced immunomodulating effect (the
suppression of tumor growth and the normalization of IL-2 concentration) is
established by considering the influence of low-intensity electromagnetic
microwaves as a parametric perturbation of the dynamical system. This finding
is in qualitative agreement with the recent experimental results on
immunocorrective effects of centimeter electromagnetic waves in tumor-bearing
mice.


  The recently discovered Acanthamoeba polyphaga Mimivirus is the largest known
DNA virus. Its particle size (>400 nm), genome length (1.2 million bp) and
large gene repertoire (911 protein coding genes) blur the established
boundaries between viruses and parasitic cellular organisms. In addition, the
analysis of its genome sequence identified new types of genes not expected to
be seen in a virus, such as aminoacyl-tRNA synthetases and other central
components of the translation machinery. In this article, we examine how the
finding of a giant virus for the first time overlapping with the world of
cellular organisms in terms of size and genome complexity might durably
influence the way we look at microbial biodiversity, and force us to
fundamentally revise our classification of life forms. We propose to introduce
the word "girus" to recognize the intermediate status of these giant DNA
viruses, the genome complexity of which make them closer to small parasitic
prokaryotes than to regular viruses.


  Giant double-stranded DNA viruses (such as record breaking Acanthamoeba
polyphaga Mimivirus), with particle sizes of 0.2 to 0.6 micron, genomes of 300
kbp to 1.200 kbp, and commensurate complex gene contents, constitute an
evolutionary mystery. They challenge the common vision of viruses,
traditionally seen as highly streamlined genomes optimally fitted to the
smallest possible -filterable- package. Such giant viruses are now discovered
in increasing numbers through the systematic sampling of ocean waters as well
as freshwater aquatic environments, where they play a significant role in
controlling phyto- and bacterio- plankton populations. The 4th algal virus
workshop showed that the study of these ecologically important viruses is now
massively entering the genomic era, promising a better understanding of their
diversity and, hopefully, some insights on their origin and the evolutionary
forces that shaped their genomes.


  A nucleotides sequence is identified, in the two (four) letters alphabet, by
the the labels of a vector state of an irreducible representation of U_q(sl(2))
(U_q(sl(2) + sl(2))), in the limit q -> 0. A master equation for the
distribution function is written, where the intensity of the one-spin flip is
assumed to depend from the variation of the labels of the state. In the two
letters approximation, the numerically computed equilibrium distribution for
short sequences is nicely fitted by a Yule distribution, which is the observed
distribution of the ranked short oligonucleotides frequency in DNA. The four
letter alphabet description, applied to the codons, is able to reproduce the
form of the fitted rank ordered usage frequencies distribution.


  A gene's rate of sequence evolution is among the most fundamental
evolutionary quantities in common use, but what determines evolutionary rates
has remained unclear. Here, we show that the two most commonly used methods to
disentangle the determinants of evolutionary rate, partial correlation analysis
and ordinary multivariate regression, produce misleading or spurious results
when applied to noisy biological data. To overcome these difficulties, we
employ an alternative method, principal component regression, which is a
multivariate regression of evolutionary rate against the principal components
of the predictor variables. We carry out the first combined analysis of seven
predictors (gene expression level, dispensability, protein abundance, codon
adaptation index, gene length, number of protein-protein interactions, and the
gene's centrality in the interaction network). Strikingly, our analysis reveals
a single dominant component which explains 40-fold more variation in
evolutionary rate than any other, suggesting that protein evolutionary rate has
a single determinant among the seven predictors. The dominant component
explains nearly half the variation in the rate of synonymous and protein
evolution. Our results support the hypothesis that selection against the cost
of translation-error-induced protein misfolding governs the rate of synonymous
and protein sequence evolution in yeast.


  We investigate the folding behavior of protein sequences by numerically
studying all sequences with maximally compact lattice model through exhaustive
enumeration. We get the prion-like behavior of protein folding. Individual
proteins remaining stable in the isolated native state may change their
conformations when they aggregate. We observe the folding properties as the
interfacial interaction strength changes, and find that the strength must be
strong enough before the propagation of the most stable structures happens.


  Many biological networks have been labelled scale-free as their degree
distribution can be approximately described by a powerlaw distribution. While
the degree distribution does not summarize all aspects of a network it has
often been suggested that its functional form contains important clues as to
underlying evolutionary processes that have shaped the network. Generally
determining the appropriate functional form for the degree distribution has
been fitted in an ad-hoc fashion.
  Here we apply formal statistical model selection methods to determine which
functional form best describes degree distributions of protein interaction and
metabolic networks. We interpret the degree distribution as belonging to a
class of probability models and determine which of these models provides the
best description for the empirical data using maximum likelihood inference,
composite likelihood methods, the Akaike information criterion and
goodness-of-fit tests. The whole data is used in order to determine the
parameter that best explains the data under a given model (e.g. scale-free or
random graph). As we will show, present protein interaction and metabolic
network data from different organisms suggests that simple scale-free models do
not provide an adequate description of real network data.


  We review some of our recent results obtained within the scope of simple
lattice models and Monte Carlo simulations that illustrate the role of native
geometry in the folding kinetics of two state folders.


  Ion channels are proteins with a hole down the middle embedded in cell
membranes. Membranes form insulating structures and the channels through them
allow and control the movement of charged particles, spherical ions, mostly
Na+, K+, Ca++, and Cl-. Membranes contain hundreds or thousands of types of
channels, fluctuating between open conducting, and closed insulating states.
Channels control an enormous range of biological function by opening and
closing in response to specific stimuli using mechanisms that are not yet
understood in physical language. Open channels conduct current of charged
particles following laws of Brownian movement of charged spheres rather like
the laws of electrodiffusion of quasi-particles in semiconductors. Open
channels select between similar ions using a combination of electrostatic and
'crowded charge' (Lennard-Jones) forces. The specific location of atoms and the
exact atomic structure of the channel protein seems much less important than
certain properties of the structure, namely the volume accessible to ions and
the effective density of fixed and polarization charge. There is no sign of
other chemical effects like delocalization of electron orbitals between ions
and the channel protein. Channels play a role in biology as important as
transistors in computers, and they use rather similar physics to perform part
of that role. Understanding their fluctuations awaits physical insight into the
source of the variance and mathematical analysis of the coupling of the
fluctuations to the other components and forces of the system.


  We study the probabilities of evolution based on random mutations and natural
selection. We conclude that evolution to multicellular eukaryots, or even
prokaryots, is unlikely to be the result of only random mutations. Complex
organisms have evolved through several mechanisms besides random mutations,
namely DNA recombination, adaptive mutations, and acquisition of foreign DNA.
We conclude that all living organisms, in addition to being self-organizing and
reproducing (autopoyetic), have built-in mechanisms of evolution, some of which
respond in very specific ways to environmental stress.


  We present an individual-based model for two interacting populations
diffusing on lattices in which a strong natural selection develops
spontaneously. The models combine traditional local predator-prey dynamics with
random walks. Individual's mobility is considered as an inherited trait. Small
variations upon inheritance, mimicking mutations, provide variability on which
natural selection may act. Although the dynamic rules defining the models do
not explicitly favor any mobility values, we found that the average mobility of
both populations tend to be maximized in various situations. In some situations
there is evidence of polymorphism, indicated by an adaptive landscape with many
local maxima. We provide evidence relating selective pressure for high mobility
with pattern formation.


  Features unique to a transfer-RNA are recognized by the corresponding
tRNA-synthetase. Keeping this in view we isolate the discriminating features of
all archaeal tRNA. These are our identity elements. Further, we investigate
tRNA-characteristics that delineate the different orders of archaea.


  Background: Recent studies in a growing number of organisms have yielded
accumulating evidence that a significant portion of the non-coding region in
the genome is transcribed. We address this issue in the yeast Saccharomyces
cerevisiae.
  Results: Taking into account the absence of a significantly large yeast EST
database, we use microarray expression data collected for genomic regions
erroneously believed to be coding to study the expression pattern of non-coding
regions in the Saccharomyces cerevisiae genome. We find that at least 164 out
of 589 (28%) such regions are expressed under specific biological conditions.
In particular, looking at the probes that are located opposing other known
genes at the same genomic locus, we find that 88 out of 341 (26%) of these
genes support antisense transcription. The expression patterns of these
antisense genes are positively correlated. We validate these results using
RT-PCR on a sample of 6 non-coding transcripts.
  Conclusions: 1. The yeast genome is transcribed on a scale larger than
previously assumed. 2. Correlated transcription of antisense genes is abundant
in the yeast genome. 3. Antisense genes in yeast are non-coding.


  A number of recent studies have hypothesized that sense-antisense RNA
transcript pairs create dsRNA duplexes that undergo extensive A-to-I RNA
editing. Here we studied human and mouse genomic antisense regions, and found
that the editing level in these areas is negligible. This observation puts in
question the scope of sense-antisense duplexes formation in-vivo, which is the
basis for a number of proposed regulatory mechanisms.


  We model the dynamical states of the C-termini of tubulin dimers that
comprise neuronal microtubules. We use molecular dynamics and other
computational tools to explore the time-dependent behavior of conformational
states of a C-terminus of tubulin within a microtubule and assume that each
C-terminus interacts via screened Coulomb forces with the surface of a tubulin
dimer, with neighboring C-termini and also with any adjacent
microtubule-associated protein 2 (MAP2). Each C-terminus can either bind to the
tubulin surface via one of several positively charged regions or can be allowed
to explore the space available in the solution surrounding the dimer. We find
that the preferential orientation of each C-terminus is away from the tubulin
surface but binding to the surface may also take place, albeit at a lower
probability. The results of our model suggest that perturbations generated by
the C-termini interactions with counter-ions surrounding a MAP2 may propagate
over distances greater than those between adjacent microtubules. Thus, the MAP2
structure is able to act as a kind of biological wire (or a cable) transmitting
local electrostatic perturbations resulting in ionic concentration gradients
from one microtubule to another. We briefly discuss implications the current
dynamic modeling may have on synaptic activation and potentiation.


  We propose a stochastic Gillespie scheme to describe the temporal
fluctuations of local denaturation zones in double-stranded DNA as a single
molecule time series. It is demonstrated that the model recovers the
equilibrium properties. We also study measurable dynamical quantities such as
the bubble size autocorrelation function. This efficient computational approach
will be useful to analyse in detail recent single molecule experiments on
clamped homopolymer breathing domains, to probe the parameter values of the
underlying Poland-Scheraga model, as well as to design experimental conditions
for similar setups.


  The two strands of a DNA molecule with a repetitive sequence can pair into
many different basepairing patterns. For perfectly periodic sequences, early
bulk experiments of Poerschke indicate the existence of a sliding process,
permitting the rapid transition between different relative strand positions
[Biophys. Chem. 2 (1974) 83]. Here, we use a detailed theoretical model to
study the basepairing dynamics of periodic and nearly periodic DNA. As
suggested by Poerschke, DNA sliding is mediated by basepairing defects (bulge
loops), which can diffuse along the DNA. Moreover, a shear force f on opposite
ends of the two strands yields a characteristic dynamic response: An outward
average sliding velocity v~1/N is induced in a double strand of length N,
provided f is larger than a threshold f_c. Conversely, if the strands are
initially misaligned, they realign even against an external force less than
f_c. These dynamics effectively result in a viscoelastic behavior of DNA under
shear forces, with properties that are programmable through the choice of the
DNA sequence. We find that a small number of mutations in periodic sequences
does not prevent DNA sliding, but introduces a time delay in the dynamic
response. We clarify the mechanism for the time delay and describe it
quantitatively within a phenomenological model. Based on our findings, we
suggest new dynamical roles for DNA in artificial nanoscale devices. The
basepairing dynamics described here is also relevant for the extension of
repetitive sequences inside genomic DNA.


  We construct a statistical ensemble of games, where in each independent
subensemble we have two players playing the same game.
  We derive the mean payoffs per move of the representative players of the
game, and we evaluate all the deterministic policies with finite memory. In
particular, we show that if one of the players has a generalized tit-for-tat
policy, the mean payoff per move of both players is the same, forcing the
equalization of the mean payoffs per move of both players. In the case of
symmetric, non-cooperative and dilemmatic games, we show that generalized
tit-for-tat policies together with the condition of not being the first to
defect, leads to the highest mean payoffs per move for the players.


  We report the results of an experimental and theoretical study of the
electronic and structural properties of a key eumelanin precursor -
5,6,-dihydroxyindole-2-carboxylic acid (DHICA) and its dimeric forms. We have
used optical spectroscopy to follow the oxidative polymerization of DHICA to
eumelanin, and observe red shifting and broadening of the absorption spectrum
as the reaction proceeds. First principles density functional theory
calculations indicate that DHICA oligomers (possible reaction products of
oxidative polymerization) have red shifted HOMO-LUMO gaps with respect to the
monomer. Furthermore, different bonding configurations (leading to oligomers
with different structures) produce a range of gaps. These experimental and
theoretical results lend support to the chemical disorder model where the broad
band monotonic absorption characteristic of all melanins is a consequence of
the superposition of a large number of inhomogeneously broadened Gaussian
transitions associated with each of the components of a melanin ensemble. These
results suggest that the traditional model of eumelanin as an amorphous organic
semiconductor is not required to explain its optical properties, and should be
thoroughly re-examined. These results have significant implications for our
understanding of the physics, chemistry and biological function of these
important biological macromolecules. Indeed, one may speculate that the robust
functionality of melanins in vitro is a direct consequence of its
heterogeneity, i.e. chemical disorder is a "low cost" natural resource in these
systems.


  We propose a mechanism for mechanical regulation at the membrane of living
cells, based on the exchange of membrane area between the cell membrane and a
membrane reservoir. The reservoir is composed of invaginated membrane
microdomains which are liable to flatten upon increase of membrane strain,
effectively controlling membrane tension. We show that the domain shape
transition is first order, allowing for coexistence between flat and
invaginated domains. During coexistence, the membrane tension is controlled by
the domains elasticity and by the kinetics of the shape transition. We show
that the tension of the plasma membrane of living cells is inherently transient
and dynamical, and that valuable insights into the organization of the cell
membrane can be obtained by studying the variation of the cell membrane tension
upon mechanical perturbation.


  We propose a criterion for optimal parameter selection in coarse-grained
models of proteins, and develop a refined elastic network model (ENM) of bovine
trypsinogen. The unimodal density-of-states distribution of the trypsinogen ENM
disagrees with the bimodal distribution obtained from an all-atom model;
however, the bimodal distribution is recovered by strengthening interactions
between atoms that are backbone neighbors. We use the backbone-enhanced model
to analyze allosteric mechanisms of trypsinogen, and find relatively strong
communication between the regulatory and active sites.


  We report here a new entropic mechanism of protein thermostability due to
residual dynamics of rotamer isomerization in native state. All-atom
simulations show that Lysines have much greater number of accessible rotamers
than Arginines in folded states of proteins. This finding suggests that Lysines
would preferentially entropically stabilize the native state. Indeed we show in
computational experiments that Arginine-to-Lysine amino acid substitutions
result in noticeable stabilization of proteins. We then hypothesize that if
evolution uses this physical mechanisms in its strategies of thermophilic
adaptation then hyperthermostable organisms would have much greater content of
Lysines in their proteomes than of comparable in size and similarly charged
Arginines.. Consistent with that, high-throughput comparative analysis of
complete proteomes shows extremely strong bias towards Arginine-to-Lysine
replacement in hyperthermophilic organisms and overall much greater content of
Lysines than Arginines in hyperthermophiles. This finding cannot be explained
by GC compositional biases. Our study provides an example of how analysis of a
delicate physical mechanism of thermostability helps to resolve a puzzle in
comparative genomics as to why aminoacid compositions of hyperthermophilic
proteomes are significantly biased towards Lysines but not Arginines


  Every weighted tree corresponds naturally to a cooperative game that we call
a "tree game"; it assigns to each subset of leaves the sum of the weights of
the minimal subtree spanned by those leaves. In the context of phylogenetic
trees, the leaves are species and this assignment captures the diversity
present in the coalition of species considered. We consider the Shapley value
of tree games and suggest a biological interpretation. We determine the linear
transformation M that shows the dependence of the Shapley value on the edge
weights of the tree, and we also compute a null space basis of M. Both depend
on the "split counts" of the tree. Finally, we characterize the Shapley value
on tree games by four axioms, a counterpart to Shapley's original theorem on
the larger class of cooperative games.


  Our paper computationally explores the extinction dynamics of an animal
species effected by a sudden spike in mortality due to an extreme event. In our
study, the animal species has a 2-year life cycle and is endowed with a high
survival probability under normal circumstances. Our proposed approach does not
involve any restraining assumptions concerning environmental variables or
predator-prey relationships. Rather it is based on the simple premise that if
observed on an year-to-year basis, the population size will be noted to either
have gone up or come down as compared to last year. The conceptualization is
borrowed from the theory of asset pricing in stochastic finance. Our results
indicate that an extreme event with a maximum shock size (i.e. the maximum
number of immediate mortalities that may be caused by an extreme event)
exceeding two-thirds the size of the pristine population can potentially drive
any animal species with a 2-year life cycle to extinction for any fecundity
level.


  A model describing the dynamics related to the spreading of non-lethal
infectious diseases in a fixed-size population is proposed. The model consists
of a non-linear delay-differential equation describing the time evolution of
the increment in the number of infectious individuals and depends upon a
limited number of parameters. Predictions are in good qualitative agreement
with data on influenza.


  It has been claimed that protein-protein interaction (PPI) networks are
scale-free based on the observation that the node degree sequence follows a
power law. Here we argue that these claims are likely to be based on erroneous
statistical analysis. Typically, the supporting data are presented using
frequency-degree plots. We show that such plots can be misleading, and should
correctly be replaced by rank-degree plots. We provide two PPI network examples
in which the frequency-degree plots appear linear on a log-log scale, but the
rank-degree plots demonstrate that the node degree sequence is far from a power
law. We conclude that at least these PPI networks are not scale-free.


  We introduce a symbolic representation like set theory to consider ecologic
interactions between species (ECOSET). The ecologic exclusion principle (Gause)
is put in a symbolic way and used as operational tool to consider more complex
cases like interaction with sterile species (SIT technique), two species with
two superposed sources (niche differentiation) and N+P species competing by N
resources, etc. Displacement (regional or characters) is also considered by
using this basic tool. Our symbolic notation gives us an operative and easy way
to consider elementary process in ecology. Some experimental data (laboratory
or field) for ecologic process are re-considered under the optic of this
set-theory.


  We report the discovery of methanogenic archaea that account for abrupt
factor 10 increases in methane concentration found by E. Brook at depths of
2954 and 3036 m in the GISP2 (Greenland Ice Sheet Project 2) ice core. The
total microbial concentration we measured with direct cell counts tracks the
excesses of methanogens that we identified by their F420 fluorescence. The
highly localized (<1 m thick) layers of methanogens suggest flow induced mixing
of layers of microbe laden anaerobic basal ice with glacial ice. The metabolic
rate we found for microbes at 2954 and 3036 m lies roughly on the Arrhenius
line for microbes imprisoned in rock, sediment, and basal ice. Equating the
loss rate of methane recently discovered in the Martian atmosphere to the
production rate by possible methanogens, we estimate that their Martian habitat
would be at a temperature of ~0 degrees C and that the concentration, if
uniformly distributed in a 100 m thick layer, would be 0.04 cell per cubic
centimeter.


  The deep-water fauna of the Mediterranean is characterized by an absence of
distinctive characteristics and by a relative impoverishment. Both are a result
of events after the Messinian salinity crisis (Late Miocene). The three main
classes of phenomena involved in producing or recording these effects are
analysed and discussed: - Historical: Sequential faunal changes during the
Pliocene and thereafter in particular those during the Quaternary glaciations
and still in progress. - Bathymetric: Changes in the vertical aspects of the
Bathyal and Abyssal zones that took place under peculiar conditions, i.e.
homothermy, a relative oligotrophy, the barrier of the Gibraltar sill, and
water mass movement. The deeper the habitat of a species in the Mediterranean,
the more extensive is its distribution elsewhere. - Geographical: There are
strong affinities and relationships between Mediterranean and Atlantic faunas.
Endemic species remain a biogeographical problem. Species always become smaller
in size eastward where they occupy a progressively deeper habitat. Thus, the
existing deep Mediterranean Sea appears to be younger than any other deep-sea
constituent of the World Ocean.


  Some microbes appear to be able to metabolize in glacial ice or permafrost.
The rate depends on temperature, nutrient level, and bioelement availability,
among other factors. I have developed a plausible argument that they do this
while confined in veins filled with acidic or saline solution that provides
nutrients and elements necessary for growth. Here I develop this scenario
further and discuss some of its implications for ice-covered planetary bodies
and for the the origin of life. An accompanying paper in the conference
proceedings (Bay et al.) discusses plans to test this hypothesis using
epifluorescence microscopy of pristine, unmelted ice samples and an optical
biospectrologging tool to assay living and dead microbes in boreholes in
glacial ice.


  The degree distribution of many biological and technological networks has
been described as a power-law distribution. While the degree distribution does
not capture all aspects of a network, it has often been suggested that its
functional form contains important clues as to underlying evolutionary
processes that have shaped the network. Generally, the functional form for the
degree distribution has been determined in an ad-hoc fashion, with clear
power-law like behaviour often only extending over a limited range of
connectivities. Here we apply formal model selection techniques to decide which
probability distribution best describes the degree distributions of protein
interaction networks. Contrary to previous studies this well defined approach
suggests that the degree distribution of many molecular networks is often
better described by distributions other than the popular power-law
distribution. This, in turn, suggests that simple, if elegant, models may not
necessarily help in the quantitative understanding of complex biological
processes.\


  As cell proliferation is limited due to the loss of telomere repeats in DNA
of normal somatic cells during division, telomere attrition can possibly play
an important role in determining the maximum life span of organisms as well as
contribute to the process of biological ageing. With computer simulations of
cell culture development in organisms, which consist of tissues of normal
somatic cells with finite growth, we otain an increase of life span and life
expectancy for longer telomeric DNA in the zygote. By additionally considering
a two-mutation model for carcinogenesis and indefinite proliferation by the
activation of telomerase, we demonstrate that the risk of dying due to cancer
can outweigh the positive effect of longer telomeres on the longevity.


  The Bacillus subtilis Spo0J/Soj proteins, implicated in chromosome
segregation and transcriptional regulation, show striking dynamics: Soj
undergoes irregular relocations from pole to pole or nucleoid to nucleoid. Here
we report on a mathematical model of the Soj dynamics. Our model, which is
closely based on the available experimental data, readily generates dynamic Soj
relocations. We show that the irregularity of the relocations may be due to the
stochastic nature of the underlying Spo0J/Soj interactions and diffusion. We
propose explanations for the behavior of several Spo0J/Soj mutants including
the "freezing" of the Soj dynamics observed in filamentous cells. Our approach
underlines the importance of incorporating stochastic effects when modelling
spatiotemporal protein dynamics inside cells.


This paper has been withdrawn by the author.


  Integrins mediate cell adhesion and are essential receptors for the
development and functioning of multicellular organisms. Integrin activation is
known to require both ligand and talin binding and to correlate with cluster
formation but the activation mechanism and precise roles of these processes are
not yet resolved. Here mathematical modeling, with known experimental
parameters, is used to show that the binding of a stabilizing factor, such as
talin, is alone insufficient to enable ligand-dependent integrin activation for
all observed conditions; an additional positive feedback is required.


  The protonation of N2 bound to the active center of nitrogenase has been
investigated using state-of-the-art DFT calculations. Dinitrogen in the
bridging mode is activated by forming two bonds to Fe sites, which results in a
reduction of the energy for the first hydrogen transfer by 123 kJ/mol. The
axial binding mode with open sulfur bridge is less reactive by 30 kJ/mol and
the energetic ordering of the axial and bridged binding mode is reversed in
favor of the bridging dinitrogen during the first protonation. Protonation of
the central ligand is thermodynamically favorable but kinetically hindered. If
the central ligand is protonated, the proton is transferred to dinitrogen
following the second protonation. Protonation of dinitrogen at the Mo site does
not lead to low-energy intermediates.


  We investigate the dynamics of a simple epidemiological model for the
invasion by a pathogen strain of a population where another strain circulates.
We assume that reinfection by the same strain is possible but occurs at a
reduced rate due to acquired immunity. The rate of reinfection by a distinct
strain is also reduced due to cross-immunity. Individual based simulations of
this model on a `small-world' network show that the host contact network
structure significantly affects the outcome of such an invasion, and as a
consequence will affect the patterns of pathogen evolution. In particular, host
populations interacting through a 'small-world' network of contacts support
lower prevalence of infection than well-mixed populations, and the region in
parameter space for which an invading strain can become endemic and coexist
with the circulating strain is smaller, reducing the potential to accommodate
pathogen diversity. We discuss the underlying mechanisms for the reported
effects, and we propose an effective mean-field model to account for the
contact structure of the host population in 'small-world' networks.


  Sporulation in Bacillus subtilis serves as a paradigm for the development of
two different cell types (mother cell and prespore) from a single cell. The
mechanism by which the two different developmental programs are initiated has
been much studied but is not well understood. With the help of existing and new
experimental results, a mathematical model has been developed that reproduces
all published in vitro experiments and makes new predictions about the
properties of the system in vivo.


  New characters based on the arrangement and morphology of dermal denticles
show that sawfishes can be divided into two distinctive groups. The first
group, comprising the knifetooth sawfish Anoxypristis cuspidata, is
characterized by tricuspid denticles variably located on both dorsal and
ventral parts of the body. The second group is represented by species of the
genus Pristis, showing an uniform and homogenous dermal covering of
monocuspidate denticles on both dorsal and ventral sides of the body and within
the buccopharyngeal cavity. Pristis is further divided into two subgroups: the
first comprises species with denticles lacking any keels and furrows (the
smalltooth sawfish Pristis pectinata, the green sawfish Pristis zijsron and the
dwarf sawfish Pristis clavata); the second comprises species with denticles
presenting keels and furrows well differentiated on their anterior part (the
common sawfish Pristis pristis, the largetooth sawfish Pristis perotteti and
the greattooth sawfish Pristis microdon). This investigation of the dermal
covering provides results which agree with studies that separate the same two
species groups of Pristis on the basis of other morphological data.


  This paper describes a Bayesian statistical method for determining the
genetic basis of a complex genetic trait. The method uses a sample of unrelated
individuals classified into two groups, for example cases and controls. Each
group is assumed to have been genotyped at a battery of marker loci using a
laboratory effort efficient technique called DNA pooling. The aim is to detect
and map a quantitative trait locus (QTL) that is not one of the typed markers.
The method works by conducting an exact Bayesian analysis under a number of
simplifying population genetic assumptions that are somewhat unrealistic.
Despite this, the method is shown to perform acceptably on datasets simulated
under a more realistic model, and furthermore is shown to outperform classical
single point methods.


  In our previous studies, we developed discrete-space Birth, Death and
Innovation Models (BDIM) of genome evolution. These models explain the origin
of the characteristic Pareto distribution of paralogous gene family sizes in
genomes, and model parameters that provide for the evolution of these
distributions within a realistic timeframe have been identified. Here we
develop the diffusion version of BDIM whose dynamics is described by the
Fokker-Plank equation and the stationary solution could be any specified Pareto
function. The diffusion models have time-dependent solutions of a special kind,
namely, the generalized self-similar solutions, which describe the transition
from one stationary distribution of the system to another; this provides for
the possibility of examining the temporal dynamics of genome evolution.
Analysis of the generalized self-similar solutions of the diffusion BDIM
reveals a biphasic curve of genome growth in which the initial, relatively
short, self-accelerating phase is followed by a prolonged phase of slow
deceleration. In biological terms, this regime of evolution can be tentatively
interpreted as a punctuated-equilibrium-like phenomenon such that whereby
evolutionary transitions are accompanied by rapid gene amplification and
innovation, followed by slow relaxation to a new stationary state.


  We show that simple stochastic models of genome evolution lead to power law
asymptotics of protein domain family size distribution. These models, called
Birth, Death and Innovation Models (BDIM), represent a special class of
balanced birth-and-death processes, in which domain duplication and deletion
rates are asymptotically equal up to the second order. The simplest, linear
BDIM shows an excellent fit to the observed distributions of domain family size
in diverse prokaryotic and eukaryotic genomes. However, the stochastic version
of the linear BDIM explored here predicts that the actual size of large
paralogous families is reached on an unrealistically long timescale. We show
that introduction of non-linearity, which might be interpreted as interaction
of a particular order between individual family members, allows the model to
achieve genome evolution rates that are much better compatible with the current
estimates of the rates of individual duplication/loss events.