  Social animals have to take into consideration the behaviour of conspecifics
when making decisions to go by their daily lives. These decisions affect their
fitness and there is therefore an evolutionary pressure to try making the right
choices. In many instances individuals will make their own choices and the
behaviour of the group will be a democratic integration of all decisions.
However, in some instances it can be advantageous to follow the choice of a few
individuals in the group if they have more information regarding the situation
that has arisen. Here I provide early evidence that decisions about shifts in
activity states in a population of bottlenose dolphin follow such a decision
making process. This unshared consensus is mediated by a non-vocal signal which
can be communicated globally within the dolphin school. These signals are
emitted by individuals that tend to have more information about the behaviour
of potential competitors because of their position in the social network. I
hypothesise that this decision making process emerged from the social structure
of the population and the need to maintain mixed-sex schools.


  We present a mathematical model, based on the compilation and statistical
processing of radiocarbon dates, of the transition from the Mesolithic to the
Neolithic, from about 7,000 to 4,000 BC in Europe. The arrival of the Neolithic
is traditionally associated with the establishment of farming-based economies;
yet in considerable areas of north-eastern Europe it is linked with the
beginning of pottery-making in the context of foraging-type communities.
Archaeological evidence, radiocarbon dates and genetic markers are consistent
with the spread of farming from a source in the Near East. However, farming was
less important in the East; the Eastern and Western Neolithic have distinct
signatures. We use a population dynamics model to suggest that this distinction
can be attributed to the presence of two waves of advance, one from the Near
East, and another through Eastern Europe. Thus, we provide a quantitative
framework in which a unified interpretation of the Western and Eastern
Neolithic can be developed.


  In this work we develop a theory of interaction of randomly patterned
surfaces as a generic prototype model of protein-protein interactions. The
theory predicts that pairs of randomly superimposed identical (homodimeric)
random patterns have always twice as large magnitude of the energy fluctuations
with respect to their mutual orientation, as compared with pairs of different
(heterodimeric) random patterns. The amplitude of the energy fluctuations is
proportional to the square of the average pattern density, to the square of the
amplitude of the potential and its characteristic length, and scales linearly
with the area of surfaces. The greater dispersion of interaction energies in
the ensemble of homodimers implies that strongly attractive complexes of random
surfaces are much more likely to be homodimers, rather than heterodimers. Our
findings suggest a plausible physical reason for the anomalously high fraction
of homodimers observed in real protein interaction networks.


  We study front propagation when an invading species competes with a resident;
we assume nearest-neighbor preemptive competition for resources in an
individual-based, two-dimensional lattice model. The asymptotic front velocity
exhibits power-law dependence on the difference between the two species' clonal
propagation rates (key ecological parameters). The mean-field approximation
behaves similarly, but the power law's exponent slightly differs from the
individual-based model's result. We also study roughening of the front, using
the framework of non-equilibrium interface growth. Our analysis indicates that
initially flat, linear invading fronts exhibit Kardar-Parisi-Zhang (KPZ)
roughening in one transverse dimension. Further, this finding implies, and is
also confirmed by simulations, that the temporal correction to the asymptotic
front velocity is of ${\cal O}(t^{-2/3})$.


  Stochastic gene expression has been implicated in a variety of cellular
processes, including cell differentiation and disease. In this issue of Cell,
Weinberger et al. (2005) take an integrated computational-experimental approach
to study the Tat transactivation feedback loop in HIV-1 and show that
fluctuations in a key regulator, Tat, can result in a phenotypic bifurcation.
This phenomenon is observed in an isogenic population where individual cells
display two distinct expression states corresponding to latent and productive
infection by HIV-1. These findings demonstrate the importance of stochastic
gene expression in molecular "decision-making."


  Supplemental data for "Stochastic Gene Expression in a Lentiviral Positive
Feedback Loop: HIV-1 Tat Fluctuations Drive Phenotypic Diversity"
[q-bio.MN/0608002, Cell. 2005 Jul 29;122(2):169-82].


  We use a generalised version of the individual-based Tangled Nature model of
evolutionary ecology to study the relationship between ecosystem structure and
evolutionary history. Our evolved model ecosystems typically exhibit
interaction networks with exponential degree distributions and an inverse
dependence between connectance and species richness. We use a simplified
network evolution model to demonstrate that the observed degree distributions
can occur as a consequence of partial correlations in the inheritance process.
Futher to this, in the limit of low connectance and maximal correlation,
distributions of power law form, $P(k){\propto}1/k$, can be achieved. We also
show that a hyperbolic relationship between connectance and species richness,
$C{\sim}1/D$ can arise as a consequence of probabilistic constraints on the
evolutionary search process.


  Adhesion-dependent cells actively sense the mechanical properties of their
environment through mechanotransductory processes at focal adhesions, which are
integrin-based contacts connecting the extracellular matrix to the
cytoskeleton. Here we present first steps towards a quantitative understanding
of focal adhesions as mechanosensors. It has been shown experimentally that
high levels of force are related to growth of and signaling at focal adhesions.
In particular, activation of the small GTPase Rho through focal adhesions leads
to the formation of stress fibers. Here we discuss one way in which force might
regulate the internal state of focal adhesions, namely by modulating the
internal rupture dynamics of focal adhesions. A simple two-spring model shows
that the stiffer the environment, the more efficient cellular force is built up
at focal adhesions by molecular motors interacting with the actin filaments.


  Bistability is a major mechanism for cellular decision making and usually
results from positive feedback in biochemical control systems. Here we show
theoretically that bistability between unbound and bound states of adhesion
clusters results from positive feedback mediated by structural rather than
biochemical processes, namely by receptor-ligand dissociation and association
dynamics which depend non-linearly on mechanical force and receptor-ligand
separation. For small cell-matrix adhesions, we find rapid switching between
unbound and bound states, which in the initial stages of adhesion allows the
cell to explore its environment through many transient adhesions.


  A fissioning organism may purge unrepairable damage by bequeathing it
preferentially to one of its daughters. Using the mathematical formalism of
superprocesses, we propose a flexible class of analytically tractable models
that allow quite general effects of damage on death rates and splitting rates
and similarly general damage segregation mechanisms. We show that, in a
suitable regime, the effects of randomness in damage segregation at fissioning
are indistinguishable from those of randomness in the mechanism of damage
accumulation during the organism's lifetime. Moreover, the optimal population
growth is achieved for a particular finite, non-zero level of combined
randomness from these two sources. In particular, when damage accumulates
deterministically, optimal population growth is achieved by a moderately
unequal division of damage between the daughters. Too little or too much
division is sub-optimal. Connections are drawn both to recent experimental
results on inheritance of damage in protozoans, to theories of the evolution of
aging, and to models of resource division between siblings.


  We explore the minimal conditions for sustainable cooperation on a spatially
distributed population of memoryless, unconditional strategies (cooperators and
defectors) in presence of unbiased, non contingent mobility in the context of
the Prisoner's Dilemma game. We find that cooperative behavior is not only
possible but may even be enhanced by such an ``always-move'' rule, when
compared with the strongly viscous (``never-move'') case. In addition, mobility
also increases the capability of cooperation to emerge and invade a population
of defectors, what may have a fundamental role in the problem of the onset of
cooperation.


  The study of synchronization in biological systems is essential for the
understanding of the rhythmic phenomena of living organisms at both molecular
and cellular levels. In this paper, by using simple dynamical systems theory,
we present a novel mechanism, named transient resetting, for the
synchronization of uncoupled biological oscillators with stimuli. This
mechanism not only can unify and extend many existing results on (deterministic
and stochastic) stimulus-induced synchrony, but also may actually play an
important role in biological rhythms. We argue that transient resetting is a
possible mechanism for the synchronization in many biological organisms, which
might also be further used in medical therapy of rhythmic disorders. Examples
on the synchronization of neural and circadian oscillators are presented to
verify our hypothesis.


  Given a weighted undirected graph $G=(V,E,d)$, the Molecular Distance
Geometry Problem (MDGP) is that of finding a function $x:G\to \mathbb{R}^{3}$,
where $||x(u)-x(v)||=d(u,v)$ for each $\{u,v\}\in E$. We show that under a few
assumptions usually satisfied in proteins, the MDGP can be formulated as a
search in a discrete space. We call this MDGP subclass the Discretizable MDGP
(DMDGP). We show that the DMDGP is \textbf{NP}-complete and we propose an
algorithm, called Branch-and-Prune (BP), which solves the DMDGP exactly. The BP
algorithm performs exceptionally well in terms of solution accuracy and can
find all solutions to any DMDGP instance. We successfully test the BP algorithm
on several randomly generated instances.


  We propose a model for describing the growth on an untreated tumor, which is
characterized in a simple way by a minimal number of parameters with a
well-defined physical interpretation. The model is motivated by invoking the
Master Equation and the Principle of Detailed Balance in the present context,
and it is easily generalizable to include the effects of various types of
therapies. In the simplest version that we consider here, it leads to a linear
equation that describes the population growth in a dynamic environment, for
which a complete solution can be given in terms of the integral of the growth
rate. The essential features of the general solution for this case are
illustrated with a few examples.


  We study the influence of stochastic effects due to finite population size in
the evolutionary dynamics of populations interacting in the multi-person
Prisoner's Dilemma game. This paper is an extension of the investigation
presented in a recent paper [Eriksson and Lindgren (2005), J. Theor. Biol.
232(3), 399]. One of the main results of the previous study is that there are
modes of dynamic behaviour, such as limit cycles and fixed points, that are
maintained due to a non-zero mutation level, resulting in a significantly
higher level of cooperation than was reported in earlier studies. In the
present study, we investigate two mechanisms in the evolutionary dynamics for
finite populations: (i) a stochastic model of the mutation process, and (ii) a
stochastic model of the selection process. The most evident effect comes from
the second extension, where we find that a previously stable limit cycle is
replaced by a trajectory that to a large extent is close to a fixed point that
is stable in the deterministic model. The effect is strong even when population
size is as large as 10,000. The effect of the first mechanism is less
pronounced, and an argument for this difference is given.


  In repeated interactions between individuals, we do not expect that exactly
the same situation will occur from one time to another. Contrary to what is
common in models of repeated games in the literature, most real situations may
differ a lot and they are seldom completely symmetric. The purpose of this
paper is to discuss a simple model of cognitive processing in the context of a
repeated interaction with varying payoffs. The interaction between players is
modelled by a repeated game with random observable payoffs. Cooperation is not
simply associated with a certain action but needs to be understood as a
phenomenon of the behaviour in the repeated game. The players are thus faced
with a more complex situation, compared to the Prisoner's Dilemma that has been
widely used for investigating the conditions for cooperation in evolving
populations. Still, there are robust cooperating strategies that usually evolve
in a population of players. In the cooperative mode, these strategies select an
action that allows for maximizing the sum of the payoff of the two players in
each round, regardless of the own payoff. Two such players maximise the
expected total long-term payoff. If the opponent deviates from this scheme, the
strategy invokes a punishment action, which aims at lowering the opponent's
score for the rest of the (possibly infinitely) repeated game. The introduction
of mistakes to the game actually pushes evolution towards more cooperative
strategies even though the game becomes more difficult.


  We define and analyze an inductive reasoning game of voluntary yearly
vaccination in order to establish whether or not a population of individuals
acting in their own self-interest would be able to prevent influenza epidemics.
We find that epidemics are rarely prevented. We also find that severe epidemics
may occur without the introduction of pandemic strains. We further address the
situation where market incentives are introduced to help ameliorating
epidemics. Surprisingly, we find that vaccinating families exacerbates
epidemics. However, a public health program requesting prepayment of
vaccinations may significantly ameliorate influenza epidemics.


  We consider character sequences evolving on a phylogenetic tree under the
TKF91 model. We show that as the sequence lengths tend to infinity the the
topology of the phylogenetic tree and the edge lengths are determined by any
one of (a) the alignment of sequences (b) the collection of sequence lengths.
We also show that the probability of any homology structure on a collection of
sequences related by a TKF91 process on a tree is independent of the root
location.
  Keywords: phylogenetics, DNA sequence evolution models, identifiability,
alignment


  We build networks of genetic similarity in which the nodes are organisms
sampled from biological populations. The procedure is illustrated by
constructing networks from genetic data of a marine clonal plant. An important
feature in the networks is the presence of clone subgraphs, i.e. sets of
organisms with identical genotype forming clones. As a first step to understand
the dynamics that has shaped these networks, we point up a relationship between
a particular degree distribution and the clone size distribution in the
populations. We construct a dynamical model for the population dynamics,
focussing on the dynamics of the clones, and solve it for the required
distributions. Scale free and exponentially decaying forms are obtained
depending on parameter values, the first type being obtained when clonal growth
is the dominant process. Average distributions are dominated by the power law
behavior presented by the fastest replicating populations.


  Persistency is the property, for differential equations in $\R^n$, that
solutions starting in the positive orthant do not approach the boundary. For
chemical reactions and population models, this translates into the
non-extinction property: provided that every species is present at the start of
the reaction, no species will tend to be eliminated in the course of the
reaction. This paper provides checkable conditions for persistence of chemical
species in reaction networks, using concepts and tools from Petri net theory.
Nontrivial examples are provided to illustrate the theory.


  RNA molecules form a sequence-specific self-pairing pattern at low
temperatures. We analyze this problem using a random pairing energy model as
well as a random sequence model that includes a base stacking energy in favor
of helix propagation. The free energy cost for separating a chain into two
equal halves offers a quantitative measure of sequence specific pairing. In the
low temperature glass phase, this quantity grows quadratically with the
logarithm of the chain length, but it switches to a linear behavior of entropic
origin in the high temperature molten phase. Transition between the two phases
is continuous, with characteristics that resemble those of a disordered elastic
manifold in two dimensions. For designed sequences, however, a power-law
distribution of pairing energies on a coarse-grained level may be more
appropriate. Extreme value statistics arguments then predict a power-law growth
of the free energy cost to break a chain, in agreement with numerical
simulations. Interestingly, the distribution of pairing distances in the ground
state secondary structure follows a remarkable power-law with an exponent -4/3,
independent of the specific assumptions for the base pairing energies.


  When bacteria are grown on a mixture of two growth-limiting substrates, they
exhibit a rich spectrum of substrate consumption patterns including diauxic
growth, simultaneous consumption, and bistable growth. In previous work, we
showed that a minimal model accounting only for enzyme induction and dilution
captures all the substrate consumption patterns. Here, we construct the
bifurcation diagram of the minimal model. The bifurcation diagram explains
several general properties of mixed-substrate growth. (1) In almost all cases
of diauxic growth, the "preferred" substrate is the one that, by itself,
supports a higher specific growth rate. In the literature, this property is
often attributed to optimality of regulatory mechanisms. Here, we show that the
minimal model, which contains only induction, displays the property under
fairly general conditions. This suggests that the higher growth rate of the
preferred substrate is an intrinsic property of the induction and dilution
kinetics.(2) The model explains the phenotypes of various mutants containing
lesions in the regions encoding for the operator, repressor, and peripheral
enzymes. A particularly striking phenotype is the "reversal of the diauxie" in
which the wild-type and mutant strains consume the very same two substrates in
opposite order. This phenotype is difficult to explain in terms of molecular
mechanisms, but it turns out to be a natural consequence of the model. We show
furthermore that the model is robust. The key property of the model, namely,
the competitive dynamics of the enzymes, is preserved even if the model is
modified to account for various regulatory mechanisms. Finally, the model has
important implications for size regulation in development, since it suggests
that protein dilution is one mechanism for coupling patterning and growth.


  We examine microscopic mechanisms for coupling stochastic oscillators so that
they display similar and correlated temporal variations. Unlike oscillatory
motion in deterministic dynamical systems, complete synchronization of
stochastic oscillators does not occur, but appropriately defined oscillator
phase variables coincide. This is illustrated in model chemical systems and
genetic networks that produce oscillations in the dynamical variables, and we
show that suitable coupling of different networks can result in their phase
synchronization.


  Many complex systems display self-organized critical states characterized by
1/f frequency scaling of power spectra. Global variables such as the
electroencephalogram, scale as 1/f, which could be the sign of self-organized
critical states in neuronal activity. By analyzing simultaneous recordings of
global and neuronal activities, we confirm the 1/f scaling of global variables
for selected frequency bands, but show that neuronal activity is not consistent
with critical states. We propose a model of 1/f scaling which does not rely on
critical states, and which is testable experimentally.


  With the launch of second line anti-retroviral therapy for HIV infected
individuals, there has been an increased expectation on surviving period of
people with HIV. We consider previously well-known models in HIV epidemiology
where the parameter for incubation period is used as one of the important
components to explain the dynamics of the variables. Such models are extended
here to explain the dynamics with respect to a given therapy that prolongs life
of an HIV infected individual. A deconvolution method is demonstrated for
estimation of parameters in the situations when no-therapy and multiple
therapies are given to the infected population. The models and deconvolution
method are extended in order to study the impact of therapy in age-structured
populations. A generalization for a situation when n-types of therapies are
available is given. Models are demonstrated using hypothetical data and
sensitivity of the parameters are also computed.


  We use a path integral representation to solve the Eigen and Crow-Kimura
molecular evolution models for the case of multiple fitness peaks with
arbitrary fitness and degradation functions. In the general case, we find that
the solution to these molecular evolution models can be written as the optimum
of a fitness function, with constraints enforced by Lagrange multipliers and
with a term accounting for the entropy of the spreading population in sequence
space. The results for the Eigen model are applied to consider virus or cancer
proliferation under the control of drugs or the immune system.


  In this paper we use a simple toy model to explore the function of the gene
Osteosarcoma-9. We are in particular interested in understanding the role of
this gene as a potent anti-apoptotic factor. The theoretical description is
constrained by experimental data from induction of apoptosis in cells where
OS-9 is overexpressed. The data available suggest that OS-9 promotes cell
viability and confers resistance to apoptosis, potentially implicating OS-9 in
the survival of cancer cells. Three different apoptosis inducing mechanisms
were tested and are here modelled. More complex and realistic models are also
discussed.


  We study the connection between the evolutionary replicator dynamics and the
number of Nash equilibria in large random bi-matrix games. Using techniques of
disordered systems theory we compute the statistical properties of both, the
fixed points of the dynamics and the Nash equilibria. Except for the special
case of zero-sum games one finds a transition as a function of the so-called
co-operation pressure between a phase in which there is a unique stable fixed
point of the dynamics coinciding with a unique Nash equilibrium, and an
unstable phase in which there are exponentially many Nash equilibria with
statistical properties different from the stationary state of the replicator
equations. Our analytical results are confirmed by numerical simulations of the
replicator dynamics, and by explicit enumeration of Nash equilibria.


  Longevity of a taxonomic group is an important issue in understanding the
dynamics of evolution. In this respect a key observation is that genera,
families or orders can each be assigned a characteristic average lifetime [Van
Valen, L., (1973) Evolutionary Theory 1, 1-30]. Using the fossil marine animal
genera database [Sepkoski, J.J.Jr. (2002) A Compendium of Fossil Marine Animal
Genera, Bull. Am. Paleontol. 363, 563 pp.] we here examine key determinants for
robustness of a higher taxonomic group in terms of the characteristics of its
constituents. We find insignificant correlation between the size of an order
and its stability against extinction, whereas we observe amazingly large
correlation between the lifetime of an order and the lifetime of its
constituent genera.


  We study the fluctuation dynamics of localized denaturation bubbles in
heteropolymer DNA with a master equation and complementary stochastic
simulation based on novel DNA stability data. A significant dependence of
opening probability and waiting time between bubble events on the local DNA
sequence is revealed and quantified for a biological sequence of the T7
bacteriophage. Quantitative agreement with data from fluorescence correlation
spectroscopy (FCS) is demonstrated.


  We propose a novel network growth model coupled with the competition
interaction to simulate macroevolution. Our work shows that the competition
plays an important role in macroevolution and it is more rational to describe
the interaction between species by network structures. Our model presents a
complete picture of the development of phyla and the splitting process. It is
found that periodic mass extinction occurred in our networks without any
extraterrestrial factors and the lifetime distribution of species is very close
to fossil record. We also perturb networks with two scenarios of mass
extinctions on different hierarchic levels in order to study their recovery.


  We present a maskless microscope projection lithography system (MPLS), in
which photomasks have been replaced by a Digital Micromirror Device type
spatial light modulator (DMD, Texas Instruments). Employing video projector
technology high resolution patterns, designed as bitmap images on the computer,
are displayed using a micromirror array consisting of about 786000 tiny
individually addressable tilting mirrors. The DMD, which is located in the
image plane of an infinity corrected microscope, is projected onto a substrate
placed in the focal plane of the microscope objective. With a 5x(0.25 NA) Fluar
microscope objective, a fivefold reduction of the image to a total size of 9
mm2 and a minimum feature size of 3.5 microns is achieved. Our system can be
used in the visible range as well as in the near UV (with a light intensity of
up to 76 mW/cm2 around the 365 nm Hg-line). We developed an inexpensive and
simple method to enable exact focusing and controlling of the image quality of
the projected patterns. Our MPLS has originally been designed for the
light-directed in situ synthesis of DNA microarrays. One requirement is a high
UV intensity to keep the fabrication process reasonably short. Another demand
is a sufficient contrast ratio over small distances (of about 5 microns). This
is necessary to achieve a high density of features (i.e. separated sites on the
substrate at which different DNA sequences are synthesized in parallel fashion)
while at the same time the number of stray light induced DNA sequence errors is
kept reasonably small. We demonstrate the performance of the apparatus in
light-directed DNA chip synthesis and discuss its advantages and limitations.


  We extend our previously developed general approach (1) to study a
phenomenological model in which the simulated packing of hard, attractive
spheres on a prolate spheroid surface with convexity constraints produces
structures identical to those of prolate virus capsid structures. Our
simulation approach combines the traditional Monte Carlo method with the method
of random sampling on an ellipsoidal surface and a convex hull searching
algorithm. Using this approach we study the assembly and structural origin of
non-icosahedral, elongated virus capsids, such as two aberrant flock house
virus (FHV) particles and the prolate prohead of bacteriophage phi29, and
discuss the implication of our simulation results in the context of recent
experimental findings.


  The need to understand the assembly kinetics of fibril formation has become
urgent because of the realization that soluble oligomers of amyloidogenic
peptides may be even more neurotoxic than the end product, namely, the amyloid
fibrils. In order to fully understand the routes to fibril formation one has to
characterize the major species in the assembly pathways. The characterization
of the energetics and dynamics of oligomers (dimers, trimers etc) is difficult
using experiments alone because they undergo large conformational fluctuations.
In this context, carefully planned molecular dynamics simulation studies,
computations using coarse-grained models, and bioinformatic analysis have given
considerable insights into the early events in the route to fibril formation.
Here, we describe progress along this direction using examples taken largely
from our own work. In this chapter, we focus on aspects of protein aggregation
using Abeta-peptides and prion proteins as examples.


  We show that the nearest neighbour distribution of distances between basis
pairs of some intron-less and intron-containing coding regions are the same
when a procedure, called {\em unfolding}, is applied. Such a procedure consists
in separating the secular variations from the oscillatory terms. The form of
the distribution obtained is quite similar to that of a random, i.e.
Poissonian, sequence. This is done for the HUMBMYH7CD, DROMYONMA, HUMBMYH7 and
DROMHC sequences. The first two correspond to highly coding regions while the
last two correspond to non-coding regions. We also show that the distributions
before the unfolding procedure depend on the secular part but, after the
unfolding procedure we obtain an striking result: all distributions are similar
to each other. The result becomes independent of the content of introns or the
species we have chosen. This is in contradiction with the results obtained with
the detrended fluctuation analysis in which the correlations yield different
results for intron-less and intron-containing regions.


  Large-scale structural patterns commonly occur in network models of complex
systems including a skewed node degree distribution and small-world topology.
These patterns suggest common organizational constraints and similar functional
consequences. Here, we investigate a structural pattern termed pathway
proliferation. Previous research enumerating pathways that link species
determined that as pathway length increases, the number of pathways tends to
increase without bound. We hypothesize that this pathway proliferation
influences the flow of energy, matter, and information in ecosystems. In this
paper, we clarify the pathway proliferation concept, introduce a measure of the
node--node proliferation rate, describe factors influencing the rate, and
characterize it in 17 large empirical food-webs. During this investigation, we
uncovered a modular organization within these systems. Over half of the
food-webs were composed of one or more subgroups that were strongly connected
internally, but weakly connected to the rest of the system. Further, these
modules had distinct proliferation rates. We conclude that pathway
proliferation in ecological networks reveals subgroups of species that will be
functionally integrated through cyclic indirect effects.


  Genomes and genes diversify during evolution; however, it is unclear to what
extent genes still retain the relationship among species. Model species for
molecular phylogenetic studies include yeasts and viruses whose genomes were
sequenced as well as plants that have the fossil-supported true phylogenetic
trees available. In this study, we generated single gene trees of seven yeast
species as well as single gene trees of nine baculovirus species using all the
orthologous genes among the species compared. Homologous genes among seven
known plants were used for validation of the fi nding. Four algorithms: maximum
parsimony, minimum evolution, maximum likelihood, and neighbor-joining, were
used. Trees were reconstructed before and after weighting the DNA and protein
sequence lengths among genes. Rarely a gene can always generate the "true tree"
by all the four algorithms. However, the most frequent gene tree, termed
"maximum gene-support tree" (MGS tree, or WMGS tree for the weighted one), in
yeasts, baculoviruses, or plants was consistently found to be the "true tree"
among the species. The results provide insights into the overall degree of
divergence of orthologous genes of the genomes analyzed and suggest the
following: 1) The true tree relationship among the species studied is still
maintained by the largest group of orthologous genes; 2) There are usually more
orthologous genes with higher similarities between genetically closer species
than between genetically more distant ones; and 3) The maximum gene-support
tree refl ects the phylogenetic relationship among species in comparison.
  Keywords: genome, gene evolution, molecular phylogeny, true tree


  36 single genes of six plants inferred 18 unique trees using maximum
parsimony. Such incongruence is an important issue and how to reconstruct the
congruent tree still is one of the most challenges in molecular phylogenetics.
For resolving this problem, a genome-wide EST data mining approach was
systematically investigated by retrieving a large size of EST data of 144
shared genes of six green plants from GenBank. The results show that the
concatenated alignments approach overcame incongruence among single-gene
phylogenies and successfully reconstructed the congruent tree of six species
with 100% jackknife support across each branch when 144 genes was used.
Jackknife supports of correct branches increased with number of genes linearly,
but those of wrong branches also increased linearly. For inferring the
congruent tree, the minimum 30 genes were required. This approach may provide
potential power in resolving conflictions of phylogenies.
  Keywords: Genome-wide; Data mining; EST; Phylogeny; Congruent tree; Jackknife
support; Plants.


  We investigate the relationship between the nested organization of
mutualistic systems and their robustness against the extinction of species. We
establish that a nested pattern of contacts is the best possible one as far as
robustness is concerned, but only when the least linked species have the
greater probability of becoming extinct. We introduce a coefficient that
provides a quantitative measure of the robustness of a mutualistic system.


  The process of `Evolutionary Diffusion', i.e. reproduction with local
mutation but without selection in a biological population, resembles standard
Diffusion in many ways. However, Evolutionary Diffusion allows the formation of
local peaks with a characteristic width that undergo drift, even in the
infinite population limit. We analytically calculate the mean peak width and
the effective random walk step size, and obtain the distribution of the peak
width which has a power law tail. We find that independent local mutations act
as a diffusion of interacting particles with increased stepsize.


  We analyse the classical model of competition between three species studied
by May and Leonard ({\it SIAM J Appl Math} \textbf{29} (1975) 243-256) with the
approaches of singularity analysis and symmetry analysis to identify values of
the parameters for which the system is integrable. We observe some striking
relations between critical values arising from the approach of dynamical
systems and the singularity and symmetry analyses.


  We studied the mathematical relations between species abundance distributions
(SADs) and species-area relationships (SARs) and found that a power-law SAR can
be generally derived from a power-law SAD without a special assumption such as
the ``canonical hypothesis''. In the present analysis, an SAR-exponent is
obtained as a function of an SAD-exponent for a finite number of species. We
also studied the inverse problem, from SARs to SADs, and found that a power-SAD
can be derived from a power-SAR under the condition that the functional form of
the corresponding SAD is invariant for changes in the number of species. We
also discuss general relationships among lognormal SADs, the broken-stick model
(exponential SADs), linear SARs and logarithmic SARs. These results suggest the
existence of a common mechanism for SADs and SARs, which could prove a useful
tool for theoretical and experimental studies on biodiversity and species
coexistence.


  We explore the complex dynamical behavior of simple predator-prey models of
biological coevolution that account for interspecific and intraspecific
competition for resources, as well as adaptive foraging behavior. In long
kinetic Monte Carlo simulations of these models we find quite robust 1/f-like
noise in species diversity and population sizes, as well as power-law
distributions for the lifetimes of individual species and the durations of
quiet periods of relative evolutionary stasis. In one model, based on the
Holling Type II functional response, adaptive foraging produces a metastable
low-diversity phase and a stable high-diversity phase.


  We are looking at local protein interaction networks from the perspective of
directed, labeled graphs with quantitative values for monotonic changes in
concentrations. These systems can be used to perform stability analysis for a
stable attractor, given initial values. They can also show re-configuration of
whole system states by dynamic insertion of links, given specific patterns of
input. The latter issue seems particularly relevant for the concept of
multistability in cellular memory. We attempt to show that this level of
analysis is well-suited for a number of relevant biological subsystems, such as
transactivation in cardiac myocytes or G-protein coupling to adrenergic
receptors. In particular, we analyse the 'motif' of an "overflow gate" as a
concentration-dependent system reconfiguration.


  The mating pathway in \emph{Saccharomyces cerevisiae} is one of the best
understood signal transduction pathways in eukaryotes. It transmits the mating
signal from plasma membrane into the nucleus through the G-protein coupled
receptor and the mitogen-activated protein kinase (MAPK) cascade. According to
the current understandings of the mating pathway, we construct a system of
ordinary differential equations to describe the process. Our model is
consistent with a wide range of experiments, indicating that it captures some
main characteristics of the signal transduction along the pathway.
Investigation with the model reveals that the shuttling of the scaffold protein
and the dephosphorylation of kinases involved in the MAPK cascade cooperate to
regulate the response upon pheromone induction and to help preserving the
fidelity of the mating signaling. We explored factors affecting the
dose-response curves of this pathway and found that both negative feedback and
concentrations of the proteins involved in the MAPK cascade play crucial role.
Contrary to some other MAPK systems where signaling sensitivity is being
amplified successively along the cascade, here the mating signal is transmitted
through the cascade in an almost linear fashion.


  Innate and adaptive immunity are connected via antigen processing and
presentation (APP), which results in the presentation of antigenic peptides to
T cells in the complex with the major histocompatibility (MHC) determinants.
MHC class II (MHC II) determinants present antigens to CD4+ T cells, which are
the main regulators of the immune response. Their genes are transcribed from
compact promoters that form first the MHC II enhanceosome, which contains
DNA-bound activators and then the MHC II transcriptosome with the addition of
the class II transactivator (CIITA). CIITA is the master regulator of MHC II
transcription. It is expressed constitutively in dendritic cells (DC) and
mature B cells and is inducible in most other cell types. Three isoforms of
CIITA exist, depending on cell type and inducing signals. CIITA is regulated at
the levels of transcription and post-translational modifications, which are
still not very clear. Inappropriate immune responses are found in several
diseases, including cancer and autoimmunity. Since CIITA regulates the
expression of MHC II genes, it is involved directly in the regulation of the
immune response. The knowledge of CIITA will facilitate the manipulation of the
immune response and might contribute to the treatment of these diseases.


  The formation mechanism of operons remains controversial despite the proposal
of many models. Although acquisition of genes from other species, horizontal
gene transfer, is considered to occur, definitive concrete cases have been
unavailable. It is desirable to select horizontally transferred genes reliably
and examine their relationship to operons. We here developed a method to
identify candidates of horizontally transferred genes based on minimization of
gene cluster insertions/deletions. To select a benchmark set of positively
horizontally transferred genes against which the candidate set can be
appraised, we devised another procedure using intergenetic alignments.
Comparison with the benchmark set of horizontally transferred genes
demonstrated the absence of a significant number of false positives in the
candidates, showing that the method identifies horizontally transferred genes
with a high degree of confidence. Horizontally transferred genes constitute at
least 5.5% of the genes in Escherichia, Shigella, and Salmonella and ~46% of
which originate from other gamma-proteobacteria. Not only informational genes,
but also operational genes (those involved in housekeeping) are horizontally
transferred less frequently than expected. A gene-cluster analysis of
Escherichia coli K-12 operons revealed that horizontal transfer produced four
entire operons and expanded two operons, but deletion of intervening genes
accounts for the formation of no operons. We propose that operons generally
form by horizontal gene transfer. We further suggest that genes with related
essential functions tend to reside in conserved operons, while genes in
nonconserved operons generally confer slight advantage to the organisms and
frequently undergo horizontal transfer and decay.


  We study evolutionary game dynamics in finite populations. We analyze an
evolutionary process, which we call pairwise comparison, for which we adopt the
ubiquitous Fermi distribution function from statistical mechanics. The inverse
temperature in this process controls the intensity of selection, leading to a
unified framework for evolutionary dynamics at all intensities of selection,
from random drift to imitation dynamics. We derive, for the first time, a
simple closed formula which determines the feasibility of cooperation in finite
populations, whenever cooperation is modeled in terms of any symmetric
two-person game. In contrast with previous results, the present formula is
valid at all intensities of selection and for any initial condition. We
investigate the evolutionary dynamics of cooperators in finite populations, and
study the interplay between intensity of selection and the remnants of interior
fixed points in infinite populations, as a function of a given initial number
of cooperators, showing how this interplay strongly affects the approach to
fixation of a given trait in finite populations, leading to counter-intuitive
results at different intensities of selection.


  In stochastic dynamical systems, different concepts of stability can be
obtained in different limits. A particularly interesting example is
evolutionary game theory, which is traditionally based on infinite populations,
where strict Nash equilibria correspond to stable fixed points that are always
evolutionarily stable. However, in finite populations stochastic effects can
drive the system away from strict Nash equilibria, which gives rise to a new
concept for evolutionary stability. The conventional and the new stability
concepts may apparently contradict each other leading to conflicting
predictions in large yet finite populations. We show that the two concepts can
be derived from the frequency dependent Moran process in different limits. Our
results help to determine the appropriate stability concept in large finite
populations. The general validity of our findings is demonstrated showing that
the same results are valid employing vastly different co-evolutionary
processes.


  We propose that the distribution of DNA words in genomic sequences can be
primarily characterized by a double Pareto-lognormal distribution, which
explains lognormal and power-law features found across all known genomes. Such
a distribution may be the result of completely random sequence evolution by
duplication processes. The parametrization of genomic word frequencies allows
for an assessment of significance for frequent or rare sequence motifs.


  Excitability is an attribute of life, and is a driving force in the descent
of complexity. Cellular electrical activity as realized by membrane proteins
that act as either channels or transporters is the basis of excitability.
Electrical signaling is mediated by a wave of action potentials, which consist
of synchronous redistribution of ionic gradients down ion channels. Ion
channels select for the passage of a particular ion species. Potassium ion
channels are gated by a variety of stimuli, including membrane voltage. Sodium
and calcium channels are gated only by membrane voltage, suggesting the
conservative argument that voltage-gated potassium channels are the founding
members of the voltage-gated ion channel superfamily. The principal focus of
this work is the investigation of the complement of potassium ion channels in
our genome and its generalizabilty. An array of issues relevant to excitability
is addressed, and a range of engagement in questions regarding the unity of
life is proffered.


  Proteases regulate various aspects of the life cycle in all organisms by
cleaving specific peptide bonds. Their action is so central for biochemical
processes that at least 2% of any known genome encodes for proteolytic enzymes.
Here we show that selected proteases pairs, despite differences in oligomeric
state, catalytic residues and fold, share a common structural organization of
functionally relevant regions which are further shown to undergo similar
concerted movements. The structural and dynamical similarities found
pervasively across evolutionarily distant clans point to common mechanisms for
peptide hydrolysis.


  The clustering of B cell receptor (BCR) molecules and the formation of the
protein segregation structure known as the immunological synapse appears to
precede antigen (Ag) uptake by B cells. The mature B cell synapse is
characterized by a central cluster of BCR/Ag molecular complexes surrounded by
a ring of LFA-1/ICAM-1 complexes. Recent experimental evidence shows receptor
clustering in B cells can occur via mechanical or signaling-driven processes.
An alternative mechanism of diffusion and affinity-dependent binding has been
proposed to explain synapse formation in the absence of signaling-driven
processes. In this work, we investigated the biophysical mechanisms that drive
immunological synapse formation in B cells across the physiological range of
BCR affinity (KA~10^6-10^10 M-1) through computational modeling. Our
computational approach is based on stochastic simulation of diffusion and
reaction events with a clearly defined mapping between probabilistic parameters
of our model and their physical equivalents. We show that a
diffusion-and-binding mechanism is sufficient to drive synapse formation only
at low BCR affinity and for a relatively stiff B cell membrane that undergoes
little deformation. We thus predict the need for alternative mechanisms: a
difference in the mechanical properties of BCR/Ag and LFA-1/ICAM-1 bonds and/or
signaling driven processes.


  The formal consideration of the concept of interaction in thermodynamic
analysis makes it possible to deduce, in the broadest terms, new results
related to the coevolution of interacting systems, irrespective of their
distance from thermodynamic equilibrium. In this paper I prove the existence of
privileged coevolution trajectories characterized by the minimum joint
production of internal entropy, a conclusion that extends Prigogine theorem to
systems evolving far from thermodynamic equilibrium. Along these trajectories
of minimum internal entropy production one of the system goes always ahead of
the other with respect to equilibrium.


  I propose here a new concept of information based on two relevant aspects of
its expression. The first related to the undeniable fact that the expression of
information modifies the physical state of its receiver. The second to the
arbitrariness of such physical changes. In fact, these changes are not
deducible from physical laws but from a code established arbitrarily. Thus,
physical information is proposed here as the capacity of producing arbitrary
changes. Once defined physical information from this physical point of view, I
deduce some basic properties of informed systems. These properties (renewal,
self-reproducing, evolution, diversification) are immediately recognized as the
attributes most characteristic of living beings, the only natural informed
systems we know. I also propose here a double evaluation of information. The
former is an absolute measure of the physical effects of its expression based
on Einstein's probability. The second is a functional measure based on the
probability that an informed systems attain a given objective as consequence of
the expression of its information.


  Mechanical unfolding of RNA structures, ranging from hairpins to ribozymes,
using laser optical tweezer (LOT) experiments have begun to reveal the features
of the energy landscape that cannot be easily explored using conventional
experiments. Upon application of constant force ($f$), RNA hairpins undergo
cooperative transitions from folded to unfolded states whereas subdomains of
ribozymes unravel one at a time. Here, we use a self-organized polymer (SOP)
model and Brownian dynamics simulations to probe mechanical unfolding at
constant force and constant-loading rate of four RNA structures of varying
complexity. Our work shows (i) the response of RNA to force is largely
determined by the native structure; (ii) only by probing mechanical unfolding
over a wide range of forces can the underlying energy landscape be fully
explored.


  Bacteria and their bacteriophages are the most abundant, widespread and
diverse groups of biological entities on the planet. In an attempt to
understand how the interactions between bacteria, virulent phages and temperate
phages might affect the diversity of these groups, we developed a novel
stochastic network model for examining the co-evolution of these ecologies. In
our approach, nodes represent whole species or strains of bacteria or phages,
rather than individuals, with "speciation" and extinction modelled by
duplication and removal of nodes. Phage-bacteria links represent host-parasite
relationships and temperate-virulent phage links denote prophage-encoded
resistance. The effect of horizontal transfer of genetic information between
strains was also included in the dynamical rules. The observed networks evolved
in a highly dynamic fashion but the ecosystems were prone to collapse (one or
more entire groups going extinct). Diversity could be stably maintained in the
model only if the probability of speciation was independent of the diversity.
Such an effect could be achieved in real ecosystems if the speciation rate is
primarily set by the availability of ecological niches.


  We study the qualitative properties of population cycles in a predator-prey
system where genetic variability allows contemporary rapid evolution of the
prey. Previous numerical studies have found that prey evolution in response to
changing predation risk can have major quantitative and qualitative effects on
predator-prey cycles, including: (i) large increases in cycle period, (ii)
changes in phase relations (so that predator and prey are cycling exactly out
of phase, rather than the classical quarter-period phase lag), and (iii)
"cryptic" cycles in which total prey density remains nearly constant while
predator density and prey traits cycle. Here we focus on a chemostat model
motivated by our experimental system [Fussmann et al. 2000,Yoshida et al. 2003]
with algae (prey) and rotifers (predators), in which the prey exhibit rapid
evolution in their level of defense against predation. We show that the effects
of rapid prey evolution are robust and general, and furthermore that they occur
in a specific but biologically relevant region of parameter space: when traits
that greatly reduce predation risk are relatively cheap (in terms of reductions
in other fitness components), when there is coexistence between the two prey
types and the predator, and when the interaction between predators and
undefended prey alone would produce cycles. Because defense has been shown to
be inexpensive, even cost-free, in a number of systems [Andersson and Levin
1999, Gagneux et al. 2006,Yoshida et al. 2004], our discoveries may well be
reproduced in other model systems, and in nature. Finally, some of our key
results are extended to a general model in which functional forms for the
predation rate and prey birth rate are not specified.


  Dynamical properties of the transcriptional regulatory network of {\it
Escherichia coli} and {\it Saccharomyces cerevisiae} are studied within the
framework of random Boolean functions. The dynamical response of these networks
to a single point mutation is characterized by the number of mutated elements
as a function of time and the distribution of the relaxation time to a new
stationary state, which turn out to be different in both networks. Comparison
with the behavior of randomized networks reveals relevant structural
characteristics other than the mean connectivity, namely the organization of
circuits and the functional form of the in-degree distribution. The abundance
of single-element circuits in {\it E. coli} and the power-law in-degree
distribution of {\it S. cerevisiae} shift their dynamics towards marginal
stability overcoming the restrictions imposed by their mean connectivities,
which is argued to be related to the simultaneous presence of robustness and
adaptivity in living organisms.


  Dynamics of a social population is analyzed taking into account some physical
constraints on individual behavior and decision making abilities. The model,
based on Evolutionary Game Theory, predicts that a population has to pass
through a series of different games, e.g as a consequence of environmental
fluctuations, in order to develop social cooperation and communication skills.
It differs from the general assumption that evolution of cooperation, the so
called Cooperation Paradox, can be explained by a single set of rules for
intra-population competitions. The developed methods, potentially, have a
practical value for some learning optimization problems in multiagent, e.g.
corporate, environment.


  Determination of sizes and flexibilities of RNA molecules is important in
understanding the nature of packing in folded structures and in elucidating
interactions between RNA and DNA or proteins. Using the coordinates of the
structures of RNA in the Protein Data Bank we find that the size of the folded
RNA structures, measured using the radius of gyration, $R_G$, follows the Flory
scaling law, namely, $R_G =5.5 N^{1/3}$ \AA where N is the number of
nucleotides. The shape of RNA molecules is characterized by the asphericity
$\Delta$ and the shape $S$ parameters that are computed using the eigenvalues
of the moment of inertia tensor. From the distribution of $\Delta$, we find
that a large fraction of folded RNA structures are aspherical and the
distribution of $S$ values shows that RNA molecules are prolate ($S>0$). The
flexibility of folded structures is characterized by the persistence length
$l_p$. By fitting the distance distribution function $P(r)$ to the worm-like
chain model we extracted the persistence length $l_p$. We find that $l_p\approx
1.5 N^{0.33}$ \AA. The dependence of $l_p$ on $N$ implies the average length of
helices should increases as the size of RNA grows. We also analyze packing in
the structures of ribosomes (30S, 50S, and 70S) in terms of $R_G$, $\Delta$,
$S$, and $l_p$. The 70S and the 50S subunits are more spherical compared to
most RNA molecules. The globularity in 50S is due to the presence of an
unusually large number (compared to 30S subunit) of small helices that are
stitched together by bulges and loops. Comparison of the shapes of the intact
70S ribosome and the constituent particles suggests that folding of the
individual molecules might occur prior to assembly.


  Complex networks are universal, arising in fields as disparate as sociology,
physics, and biology. In the past decade, extensive research into the
properties and behaviors of complex systems has uncovered surprising
commonalities among the topologies of different systems. Attempts to explain
these similarities have led to the ongoing development and refinement of
network models and graph-theoretical analysis techniques with which to
characterize and understand complexity. In this tutorial, we demonstrate
through illustrative examples, how network measures and models have contributed
to the elucidation of the organization of complex systems.


  Background:Prediction of protein three-dimensional structures from amino acid
sequences is a long-standing goal in computational/molecular biology. The
successful discrimination of protein folds would help to improve the accuracy
of protein 3D structure prediction. Results: In this work, we propose a method
based on linear discriminant analysis (LDA) for recognizing proteins belonging
to 30 different folds using the occurrence of amino acid residues in a set of
1612 proteins. The present method could discriminate the globular proteins from
30 major folding types with the sensitivity of 37%, which is comparable to or
better than other methods in the literature. A web server has been developed
for predicting the folding type of the protein from amino acid sequence and it
is available at http://granular.com/PROLDA/. Conclusions:Linear discriminant
analysis based on amino acid occurrence could successfully recognize protein
folds. The present method has several advantages such as, (i) it directly
predicts the folding type of a protein without performing pair-wise
comparisons, (ii) it can discriminate folds among large number of proteins and
(iii) it is very fast to obtain the results. This is a simple method, which can
be easily incorporated in any other structure prediction algorithms.


  We address phylogenetic reconstruction when the data is generated from a
mixture distribution. Such topics have gained considerable attention in the
biological community with the clear evidence of heterogeneity of mutation
rates. In our work, we consider data coming from a mixture of trees which share
a common topology, but differ in their edge weights (i.e., branch lengths). We
first show the pitfalls of popular methods, including maximum likelihood and
Markov chain Monte Carlo algorithms. We then determine in which evolutionary
models, reconstructing the tree topology, under a mixture distribution, is
(im)possible. We prove that every model whose transition matrices can be
parameterized by an open set of multi-linear polynomials, either has
non-identifiable mixture distributions, in which case reconstruction is
impossible in general, or there exist linear tests which identify the topology.
This duality theorem, relies on our notion of linear tests and uses ideas from
convex programming duality. Linear tests are closely related to linear
invariants, which were first introduced by Lake, and are natural from an
algebraic geometry perspective.


  We consider a broad class of stochastic lattice predator-prey models, whose
main features are overviewed. In particular, this article aims at drawing a
picture of the influence of spatial fluctuations, which are not accounted for
by the deterministic rate equations, on the properties of the stochastic
models. Here, we outline the robust scenario obeyed by most of the lattice
predator-prey models with an interaction "a' la Lotka-Volterra". We also show
how a drastically different behavior can emerge as the result of a subtle
interplay between long-range interactions and a nearest-neighbor exchange
process.


  The native three dimensional structure of a single protein is determined by
the physico chemical nature of its constituent amino acids. The twenty
different types of amino acids, depending on their physico chemical properties,
can be grouped into three major classes - hydrophobic, hydrophilic and charged.
We have studied the anatomy of the weighted and unweighted networks of
hydrophobic, hydrophilic and charged residues separately for a large number of
proteins. Our results show that the average degree of the hydrophobic networks
has significantly larger value than that of hydrophilic and charged networks.
The average degree of the hydrophilic networks is slightly higher than that of
charged networks. The average strength of the nodes of hydrophobic networks is
nearly equal to that of the charged network; whereas that of hydrophilic
networks has smaller value than that of hydrophobic and charged networks. The
average strength for each of the three types of networks varies with its
degree. The average strength of a node in charged networks increases more
sharply than that of the hydrophobic and hydrophilic networks. Each of the
three types of networks exhibits the 'small-world' property. Our results
further indicate that the all amino acids' networks and hydrophobic networks
are of assortative type. While maximum of the hydrophilic and charged networks
are of assortative type, few others have the characteristics of disassortative
mixing of the nodes. We have further observed that all amino acids' networks
and hydrophobic networks bear the signature of hierarchy; whereas the
hydrophilic and charged networks do not have any hierarchical signature.


  Iron is an essential trace-element for most organisms. However, because high
concentration of free intracellular iron is cytotoxic, cells have developed
complex regulatory networks that keep free intracellular iron concentration at
optimal range, allowing the incorporation of the metal into iron-using enzymes
and minimizing damage to the cell. We built a mathematical model of the network
that controls iron uptake and usage in the bacterium Escherichia coli to
explore the dynamics of iron flow. We simulate the effect of sudden decrease or
increase in the extracellular iron level on intracellular iron distribution.
Based on the results of simulations we discuss the possible roles of the small
RNA RyhB and the Fe-S cluster assembly systems in the optimal redistribution of
iron flows. We suggest that Fe-S cluster assembly is crucial to prevent the
accumulation of toxic levels of free intracellular iron when the environment
suddenly becomes iron rich.


  The formation and regulation of macromolecular complexes provides the
backbone of most cellular processes, including gene regulation and signal
transduction. The inherent complexity of assembling macromolecular structures
makes current computational methods strongly limited for understanding how the
physical interactions between cellular components give rise to systemic
properties of cells. Here we present a stochastic approach to study the
dynamics of networks formed by macromolecular complexes in terms of the
molecular interactions of their components. Exploiting key thermodynamic
concepts, this approach makes it possible to both estimate reaction rates and
incorporate the resulting assembly dynamics into the stochastic kinetics of
cellular networks. As prototype systems, we consider the lac operon and phage
lambda induction switches, which rely on the formation of DNA loops by proteins
and on the integration of these protein-DNA complexes into intracellular
networks. This cross-scale approach offers an effective starting point to move
forward from network diagrams, such as those of protein-protein and DNA-protein
interaction networks, to the actual dynamics of cellular processes.


  The formation of DNA loops by proteins and protein complexes is ubiquitous to
many fundamental cellular processes, including transcription, recombination,
and replication. Here we review recent advances in understanding the properties
of DNA looping in its natural context and how they propagate to the cellular
behavior through gene regulation. The results of connecting the molecular
properties with cellular physiology indicate that looping of DNA in vivo is
much more complex and easier than predicted from current models and reveals a
wealth of previously unappreciated details.


  We investigate a continuous time, probability measure-valued dynamical system
that describes the process of mutation-selection balance in a context where the
population is infinite, there may be infinitely many loci, and there are weak
assumptions on selective costs. Our model arises when we incorporate very
general recombination mechanisms into a previous model of mutation and
selection from Steinsaltz, Evans and Wachter (2005) and take the relative
strength of mutation and selection to be sufficiently small. The resulting
dynamical system is a flow of measures on the space of loci. Each such measure
is the intensity measure of a Poisson random measure on the space of loci: the
points of a realization of the random measure record the set of loci at which
the genotype of a uniformly chosen individual differs from a reference wild
type due to an accumulation of ancestral mutations. Our motivation for working
in such a general setting is to provide a basis for understanding
mutation-driven changes in age-specific demographic schedules that arise from
the complex interaction of many genes, and hence to develop a framework for
understanding the evolution of aging.
  We establish the existence and uniqueness of the dynamical system, provide
conditions for the existence and stability of equilibrium states, and prove
that our continuous-time dynamical system is the limit of a sequence of
discrete-time infinite population mutation-selection-recombination models in
the standard asymptotic regime where selection and mutation are weak relative
to recombination and both scale at the same infinitesimal rate in the limit.


  Results are reported concerning the formation of spatial patterns in the
two-species ratio-dependent predator-prey model driven by spatial
colored-noise. The results show that there is a critical value with respect to
the intensity of spatial noise for this system when the parameters are in the
Turing space, above which the regular spatial patterns appear in two
dimensions, but under which there are not regular spatial patterns produced. In
particular, we investigate in two-dimensional space the formation of regular
spatial patterns with the spatial noise added in the side and the center of the
simulation domain, respectively.


  We study statistical properties of interacting protein-like surfaces and
predict two strong, related effects: (i) statistically enhanced self-attraction
of proteins; (ii) statistically enhanced attraction of proteins with similar
structures. The effects originate in the fact that the probability to find a
pattern self-match between two identical, even randomly organized interacting
protein surfaces is always higher compared with the probability for a pattern
match between two different, promiscuous protein surfaces. This theoretical
finding explains statistical prevalence of homodimers in protein-protein
interaction networks reported earlier. Further, our findings are confirmed by
the analysis of curated database of protein complexes that showed highly
statistically significant overrepresentation of dimers formed by structurally
similar proteins with highly divergent sequences (superfamily heterodimers). We
predict that significant fraction of heterodimers evolved from homodimers with
the negative design evolutionary pressure applied against promiscuous homodimer
formation. This is achieved through the formation of highly specific contacts
formed by charged residues as demonstrated both in model and real superfamily
heterodimers


  An important part of the analysis of bio-molecular networks is to detect
different functional units. Different functions are reflected in a different
evolutionary dynamics, and hence in different statistical characteristics of
network parts. In this sense, the {\em global statistics} of a biological
network, e.g., its connectivity distribution, provides a background, and {\em
local deviations} from this background signal functional units. In the
computational analysis of biological networks, we thus typically have to
discriminate between different statistical models governing different parts of
the dataset. The nature of these models depends on the biological question
asked. We illustrate this rationale here with three examples: identification of
functional parts as highly connected \textit{network clusters}, finding
\textit{network motifs}, which occur in a similar form at different places in
the network, and the analysis of \textit{cross-species network correlations},
which reflect evolutionary dynamics between species.


  No influence was seen when in two models with memory effects the populations
were drastically decreased after equilibrium was established, and then allowed
to increase again.


  Why sex is maintained in nature is a fundamental question in biology. Natural
genetic transformation (NGT) is a sexual process by which bacteria actively
take up exogenous DNA and use it to replace homologous chromosomal sequences.
As it has been demonstrated, the role of NGT in repairing deleterious mutations
under constant selection is insufficient for its survival, and the lack of
other viable explanations have left no alternative except that DNA uptake
provides nucleotides for food. Here we develop a novel simulation approach for
the long-term dynamics of genome organization (involving the loss and
acquisition of genes) in a bacterial species consisting of a large number of
spatially distinct populations subject to independently fluctuating ecological
conditions. Our results show that in the presence of weak inter-population
migration NGT is able to subsist as a mechanism to reload locally lost,
intermittently selected genes from the collective gene pool of the species
through DNA uptake from migrants. Reloading genes and combining them with those
in locally adapted genomes allow individual cells to re-adapt faster to
environmental changes. The machinery of transformation survives under a wide
range of model parameters readily encompassing real-world biological
conditions. These findings imply that the primary role of NGT is not to serve
the cell with food, but to provide homologous sequences for restoring genes
that have disappeared from or become degraded in the local population.


  Evolutionary Prisoner's Dilemma games with quenched inhomogeneities in the
spatial dynamical rules are considered. The players following one of the two
pure strategies (cooperation or defection) are distributed on a two-dimensional
lattice. The rate of strategy adoption from a randomly chosen neighbors are
controlled by the payoff difference and a two-value pre-factor $w$
characterizing the players whom the strategy learned from. The reduced teaching
activity of players is distributed randomly with concentrations $\nu$ at the
beginning and fixed further on. Numerical and analytical calculations are
performed to study the concentration of cooperators as a function of $w$ and
$\nu$ for different noise levels and connectivity structures. Significant
increase of cooperation is found within a wide range of parameters for this
dynamics. The results highlight the importance of asymmetry characterizing the
exchange of master-follower role during the strategy adoptions.


  The "Machiavellian intelligence" hypothesis (or the "social brain"
hypothesis) posits that large brains and distinctive cognitive abilities of
humans have evolved via intense social competition in which social competitors
developed increasingly sophisticated "Machiavellian" strategies as a means to
achieve higher social and reproductive success. Here we build a mathematical
model aiming to explore this hypothesis. In the model, genes control brains
which invent and learn strategies (memes) which are used by males to gain
advantage in competition for mates. We show that the dynamics of intelligence
has three distinct phases. During the dormant phase only newly invented memes
are present in the population. During the cognitive explosion phase the
population's meme count and the learning ability, cerebral capacity
(controlling the number of different memes that the brain can learn and use),
and Machiavellian fitness of individuals increase in a runaway fashion. During
the saturation phase natural selection resulting from the costs of having large
brains checks further increases in cognitive abilities. Overall, our results
suggest that the mechanisms underlying the "Machiavellian intelligence"
hypothesis can indeed result in the evolution of significant cognitive
abilities on the time scale of 10 to 20 thousand generations. We show that
cerebral capacity evolves faster and to a larger degree than learning ability.
Our model suggests that there may be a tendency toward a reduction in cognitive
abilities (driven by the costs of having a large brain) as the reproductive
advantage of having a large brain decreases and the exposure to memes increases
in modern societies.


  Sea urchin feeding fronts are a striking example of spatial pattern formation
in an ecological system. If it is assumed that urchins are asocial, and that
they move randomly, then the formation of these dense fronts is an apparent
paradox. The key lies in observations that urchins move further in areas where
their algal food is less plentiful. This naturally leads to the accumulation of
urchins in areas with abundant algae. If urchin movement is represented as a
random walk, with a step size that depends on algal concentration, then their
movement may be described by a Fokker-Planck diffusion equation. For certain
combinations of algal growth and urchin grazing, travelling wave solutions are
obtained. Two dimensional simulations of urchin algal dynamics show that an
initially uniformly distributed urchin population, grazing on an alga with a
smoothly varying density, may form a propagating front separating two sharply
delineated regions. On one side of the front algal density is uniformly low,
and on the other side of the front algal density is uniformly high. Bounds on
when stable fronts will form are obtained in terms of urchin density and
grazing, and algal growth.


  This paper addresses the question of how population diffusion affects the
formation of the spatial patterns in the spatial epidemic model by Turing
mechanisms. In particular, we present theoretical analysis to results of the
numerical simulations in two dimensions. Moreover, there is a critical value
for the system within the linear regime. Below the critical value the spatial
patterns are impermanent, whereas above it stationary spot and stripe patterns
can coexist over time. We have observed the striking formation of spatial
patterns during the evolution, but the isolated ordered spot patterns don't
emerge in the space.


  We describe a simple numerical simulation, suitable for an undergraduate
project (or graduate problem set), of the Brownian motion of a particle in a
Hooke-law potential well. Understanding this physical situation is a practical
necessity in many experimental contexts, for instance in single molecule
biophysics; and its simulation helps the student to appreciate the dynamical
character of thermal equilibrium. We show that the simulation succeeds in
capturing behavior seen in experimental data on tethered particle motion.


  This research explores the role and representation of network structure for
LTI Systems. We demonstrate that transfer functions contain no structural
information without more assumptions being made about the system, assumptions
that we believe are unreasonable when dealing with truly complex systems. We
then introduce Dynamical Structure Functions as an alternative, graphical-model
based representation of LTI systems that contain both dynamical and structural
information of the system. We use Dynamical Structure to prove necessary and
sufficient conditions for estimating structure from data, and demonstrate, for
example, the danger of attempting to use steady-state information to estimate
network structure.


  Viruses self-assemble from identical capsid proteins and their genome
consisting, for example, of a long single stranded (ss) RNA. For a big class of
T = 3 viruses capsid proteins have long positive N-terminal tails. We explore
the role played by the Coulomb interaction between the brush of positive
N-terminal tails rooted at the inner surface of the capsid and the negative ss
RNA molecule. We show that viruses are most stable when the total contour
length of ss RNA is close to the total length of the tails. For such a
structure the absolute value of the total RNA charge is approximately twice
larger than the charge of the capsid. This conclusion agrees with structural
data.


  Plants both lose water and take in carbon dioxide through microscopic
stomatal pores, each of which is regulated by a surrounding pair of guard
cells. During drought, the plant hormone abscisic acid (ABA) inhibits stomatal
opening and promotes stomatal closure, thereby promoting water conservation.
Here we synthesize experimental results into a consistent guard cell signal
transduction network for ABA-induced stomatal closure, and develop a dynamic
model of this process. Our model captures the regulation of more than forty
identified network components, and accords well with previous experimental
results at both the pathway and whole cell physiological level. Our analysis
reveals the novel predictions that the disruption of membrane depolarizability,
anion efflux, actin cytoskeleton reorganization, cytosolic pH increase, the
phosphatidic acid pathway or of K+ efflux through slowly activating K+ channels
at the plasma membrane lead to the strongest reduction in ABA responsiveness.
Initial experimental analysis assessing ABA-induced stomatal closure in the
presence of cytosolic pH clamp imposed by the weak acid butyrate is consistent
with model prediction. Our method can be readily applied to other biological
signaling networks to identify key regulatory components in systems where
quantitative information is limited.


  Finite-size fluctuations arising in the dynamics of competing populations may
have dramatic influence on their fate. As an example, in this article, we
investigate a model of three species which dominate each other in a cyclic
manner. Although the deterministic approach predicts (neutrally) stable
coexistence of all species, for any finite population size, the intrinsic
stochasticity unavoidably causes the eventual extinction of two of them.


  This paper proposes a new mathematical approach to characterize native
protein structures based on the discrete differential geometry of tetrahedron
tiles. In the approach, local structure of proteins is classified into finite
types according to shape. And one would obtain a number sequence representation
of protein structures automatically. As a result, it would become possible to
quantify structural preference of amino-acids objectively. And one could use
the wide variety of sequence alignment programs to study protein structures
since the number sequence has no internal structure.
  The programs and this paper with clear figures are available from
http://www.genocript.com.


  Competition within the immune system may degrade immune control of viral
infections. We formalize the evolution that occurs in both HIV-1 and the immune
system quasispecies. Inclusion of competition in the immune system leads to a
novel balance between the immune response and HIV-1, in which the eventual
outcome is HIV-1 escape rather than control. The analytical model reproduces
the three stages of HIV-1 infection. We propose a vaccine regimen that may be
able to reduce competition between T cells, potentially eliminating the third
stage of HIV-1.


  Among the distance based algorithms in phylogenetic tree reconstruction, the
neighbor-joining algorithm has been a widely used and effective method. We
propose a new algorithm which counts the number of consistent quartets for
cherry picking with tie breaking. We show that the success rate of the new
algorithm is almost equal to that of neighbor-joining. This gives an
explanation of the qualitative nature of neighbor-joining and that of
dissimilarity maps from DNA sequence data. Moreover, the new algorithm always
reconstructs correct trees from quartet consistent dissimilarity maps.


  Cell wall proteins are essential constituents of plant cell walls; they are
involved in modifications of cell wall components, wall structure, signaling
and interactions with plasma membrane proteins at the cell surface. The
application of proteomic approaches to the cell wall compartment raises
important questions: are there technical problems specific to cell wall
proteomics? What kinds of proteins can be found in Arabidopsis walls? Are some
of them unexpected? What sort of post-translational modifications have been
characterized in cell wall proteins to date? The purpose of this review is to
discuss the experimental results obtained to date using proteomics, as well as
some of the new questions challenging future research.


  The ultimate goal of proteomic analysis of a cell compartment should be the
exhaustive identification of resident proteins; excluding proteins from other
cell compartments. Plant cell walls possess specific difficulties. Several
reported procedures to isolate cell walls for proteomic analyses led to the
isolation of a high proportion (more than 50%) of predicted intracellular
proteins. The rationales of several published procedures to isolate cell walls
for proteomics were analyzed, with regard to the bioinformatic-predicted
subcellular localization of the identified proteins. A new procedure was
developed to prepare cell walls from etiolated hypocotyls of Arabidopsis
thaliana. After salt extraction, a high proportion of proteins predicted to be
secreted was released (73%), belonging to the same functional classes as
proteins identified using previously described protocols. The new cell wall
preparation described in this paper gives the lowest proportion of proteins
predicted to be intracellular when compared to available protocols. The
application of its principles should lead to a more realistic view of the cell
wall proteome, at least for the weakly bound CWP extractable by salts. In
addition, it offers a clean cell wall preparation for subsequent extraction of
strongly bound CWP.


  Smallpox is a deadly disease that can be intentionally reintroduced into the
human population as a bioweapon. While host gene expression microarray
profiling can be used to detect infection, the analysis of this information
using unsupervised and supervised classification techniques can produce
contradictory results. Here, we present a novel computational approach to
incorporate molecular genome annotation features that are key for identifying
early infection biomarkers (EIB). Our analysis identified 58 EIBs expressed in
peripheral blood mononuclear cells (PBMCs) collected from 21 cynomolgus
macaques (Macaca fascicularis) infected with two variola strains via aerosol
and intravenous exposure. The level of expression of these EIBs was correlated
with disease progression and severity. No overlap between the EIBs
co-expression and protein interaction data reported in public databases was
found. This suggests that a pathogen-specific re-organization of the gene
expression and protein interaction networks occurs during infection. To
identify potential genome-wide protein interactions between variola and humans,
we performed a protein domain analysis of all smallpox and human proteins. We
found that only 55 of the 161 protein domains in smallpox are also present in
the human genome. These co-occurring domains are mostly represented in proteins
involved in blood coagulation, complement activation, angiogenesis,
inflammation, and hormone transport. Several of these proteins are within the
EIBs category and suggest potential new targets for the development of
therapeutic countermeasures.


  We apply the methods of optimal experimental design to a differential
equation model for epidermal growth factor receptor (EGFR) signaling,
trafficking, and down-regulation. The model incorporates the role of a recently
discovered protein complex made up of the E3 ubiquitin ligase, Cbl, the guanine
exchange factor (GEF), Cool-1 (Beta-Pix), and the Rho family G protein Cdc42.
The complex has been suggested to be important in disrupting receptor
down-regulation. We demonstrate that the model interactions can accurately
reproduce the experimental observations, that they can be used to make
predictions with accompanying uncertainties, and that we can apply ideas of
optimal experimental design to suggest new experiments that reduce the
uncertainty on unmeasurable components of the system.


  We develop a general method to explore how the function performed by a
biological network can constrain both its structural and dynamical network
properties. This approach is orthogonal to prior studies which examine the
functional consequences of a given structural feature, for example a scale free
architecture. A key step is to construct an algorithm that allows us to
efficiently sample from a maximum entropy distribution on the space of boolean
dynamical networks constrained to perform a specific function, or cascade of
gene expression. Such a distribution can act as a "functional null model" to
test the significance of any given network feature, and can aid in revealing
underlying evolutionary selection pressures on various network properties.
Although our methods are general, we illustrate them in an analysis of the
yeast cell cycle cascade. This analysis uncovers strong constraints on the
architecture of the cell cycle regulatory network as well as significant
selection pressures on this network to maintain ordered and convergent
dynamics, possibly at the expense of sacrificing robustness to structural
perturbations.


  Biomolecular networks have to perform their functions robustly. A robust
function may have preferences in the topological structures of the underlying
network. We carried out an exhaustive computational analysis on network
topologies in relation to a patterning function in Drosophila embryogenesis. We
found that while the vast majority of topologies can either not perform the
required function or only do so very fragilely, a small fraction of topologies
emerges as particularly robust for the function. The topology adopted by
Drosophila, that of the segment polarity network, is a top ranking one among
all topologies with no direct autoregulation. Furthermore, we found that all
robust topologies are modular--each being a combination of three kinds of
modules. These modules can be traced back to three sub-functions of the
patterning function and their combinations provide a combinatorial variability
for the robust topologies. Our results suggest that the requirement of
functional robustness drastically reduces the choices of viable topology to a
limited set of modular combinations among which nature optimizes its choice
under evolutionary and other biological constraints.


  An attempt to use phylogenetic invariants for tree reconstruction was made at
the end of the 80s and the beginning of the 90s by several authors (the initial
idea due to Lake and Cavender and Felsenstein in 1987. However, the efficiency
of methods based on invariants is still in doubt, probably because these
methods only used few generators of the set of phylogenetic invariants. The
method studied in this paper was first introduced by Casanellas, Garcia and
Sullivant and it is the first method based on invariants that uses the whole
set of generators for DNA data. The simulation studies performed in this paper
prove that it is a very competitive and highly efficient phylogenetic
reconstruction method, especially for non-homogeneous models on phylogenetic
trees.


  We model a plant leaf by using two-dimensional domain walls with internal
structures. Such domain walls can be found as soliton solutions in field theory
describing magnetic materials. The radiation scattered by such domain walls
behaves quite similar to the spectral reflectance of plant leaves. The model
nicely simulates the spectral reflectance of a plant leaf as a function of the
wavelength.


  Inspired by Minority Games, we constructed a novel individual-level game of
adaptive decision-making based on the dilemma of deciding whether to
participate in voluntary influenza vaccination programs. The proportion of the
population vaccinated (i.e., the vaccination coverage) determines epidemic
severity. Above a critical vaccination coverage, epidemics are prevented; hence
individuals find it unnecessary to vaccinate. The adaptive dynamics of the
decisions directly affect influenza epidemiology and, conversely, influenza
epidemiology strongly influences decision-making. This feedback mechanism
creates a unique self-organized state where epidemics are prevented. This state
is attracting, but unstable; thus epidemics are rarely prevented. This result
implies that vaccination will have to be mandatory if the public health
objective is to prevent influenza epidemics. We investigated how collective
behavior changes when public health programs are implemented. Surprisingly,
programs requiring advance payment for several years of vaccination prevents
severe epidemics, even with voluntary vaccination. Prevention is determined by
the individuals' adaptability, memory, and number of pre-paid vaccinations.
Notably, vaccinating families exacerbates and increases the frequency of severe
epidemics.


  We study the dynamics of a population subject to selective pressures,
evolving either on RNA neutral networks or in toy fitness landscapes. We
discuss the spread and the neutrality of the population in the steady state.
Different limits arise depending on whether selection or random drift are
dominant. In the presence of strong drift we show that observables depend
mainly on $M \mu$, $M$ being the population size and $\mu$ the mutation rate,
while corrections to this scaling go as 1/M: such corrections can be quite
large in the presence of selection if there are barriers in the fitness
landscape. Also we find that the convergence to the large $M \mu$ limit is
linear in $1/M \mu$. Finally we introduce a protocol that minimizes drift; then
observables scale like 1/M rather than $1/(M\mu)$, allowing one to determine
the large $M$ limit faster when $\mu$ is small; furthermore the genotypic
diversity increases from $O(\ln M)$ to $O(M)$.


  Proteins synthesized on membrane-bound ribosomes are sorted at the Golgi
apparatus level for delivery to various cellular destinations: the plasma
membrane or the extracellular space, and the lytic vacuole or lysosome. Sorting
involves the assembly of vesicles, which preferentially package soluble
proteins with a common destination. The selection of proteins for a particular
vesicle type involves the recognition of proteins by specific receptors, such
as the vacuolar sorting receptors for vacuolar targeting. Most eukaryotic
organisms have one or two receptors to target proteins to the lytic vacuole.
Surprisingly, plants have several members of the same family, seven in
Arabidopsis thaliana. Why do plants have so many proteins to sort soluble
proteins to their respective destinations? The presence of at least two types
of vacuoles, lytic and storage, seems to be a partial answer. In this review we
analyze the last experimental evidence supporting the presence of different
subfamilies of plant vacuolar sorting receptors.


  This submission has been withdrawn by arXiv administration due to complaints
of misuse of institutional affiliation.


  We introduce a novel concept, the minimal molecular surface (MMS), as a new
paradigm for the theoretical modeling of biomolecule-solvent interfaces. When a
less polar macromolecule is immersed in a polar environment, the surface free
energy minimization occurs naturally to stabilizes the system, and leads to an
MMS separating the macromolecule from the solvent. For a given set of atomic
constraints (as obstacles), the MMS is defined as one whose mean curvature
vanishes away from the obstacles. An iterative procedure is proposed to compute
the MMS. Extensive examples are given to validate the proposed algorithm and
illustrate the new concept. We show that the MMS provides an indication to
DNA-binding specificity. The proposed algorithm represents a major step forward
in minimal surface generation.


  Asterias (\url{http://www.asterias.info}) is an integrated collection of
freely-accessible web tools for the analysis of gene expression and aCGH data.
Most of the tools use parallel computing (via MPI). Most of our applications
allow the user to obtain additional information for user-selected genes by
using clickable links in tables and/or figures. Our tools include:
normalization of expression and aCGH data; converting between different types
of gene/clone and protein identifiers; filtering and imputation; finding
differentially expressed genes related to patient class and survival data;
searching for models of class prediction; using random forests to search for
minimal models for class prediction or for large subsets of genes with
predictive capacity; searching for molecular signatures and predictive genes
with survival data; detecting regions of genomic DNA gain or loss. The
capability to send results between different applications, access to additional
functional information, and parallelized computation make our suite unique and
exploit features only available to web-based applications.


  We present a model for a biological reactor in which bacteria tend to
aggregate in flocs, as encountered in wastewater treatment plants. The
influence of this flocculation on the growth dynamics of the bacteria is
studied. We argue that a description in terms of a specific growth rate is
possible when the flocculation dynamics is much faster than the other processes
in the system. An analytical computation shows that in this case, the growth
rate is density-dependent, i.e., depends both on the substrate and the biomass
density. When the flocculation time scale overlaps with the other time scales
present in the system, the notion of specific growth rate becomes problematic.
However, we show numerically that a density-dependent growth rate can still
accurately describe the system response to certain perturbations.


  Kinesin and related motor proteins utilize ATP fuel to propel themselves
along the external surface of microtubules in a processive and directional
fashion. We show that the observed step-like motion is possible through time
varying charge distributions furnished by the ATP hydrolysis circle while the
static charge configuration on the microtuble provides the guide for motion.
Thus, while the chemical hydrolysis energy induces appropriate local
conformational changes, the motor translational energy is fundamentally
electrostatic. Numerical simulations of the mechanical equations of motion show
that processivity and directionality are direct consequences of the
ATP-dependent electrostatic interaction between the different charge
distributions of kinesin and microtubule.


  Stretching of a protein by a fluid flow is compared to that in a force-clamp
apparatus. The comparison is made within a simple topology-based dynamical
model of a protein in which the effects of the flow are implemented using
Langevin dynamics. We demonstrate that unfolding induced by a uniform flow
shows a richer behavior than that in the force clamp. The dynamics of unfolding
is found to depend strongly on the selection of the amino acid, usually one of
the termini, which is anchored. These features offer potentially wider
diagnostic tools to investigate structure of proteins compared to experiments
based on the atomic force microscopy.


  Common fragile sites (cfs) are specific regions in the human genome that are
particularly prone to genomic instability under conditions of replicative
stress. Several investigations support the view that common fragile sites play
a role in carcinogenesis. We discuss a genome-wide approach based on graph
theory and Gene Ontology vocabulary for the functional characterization of
common fragile sites and for the identification of genes that contribute to
tumour cell biology. CFS were assembled in a network based on a simple measure
of correlation among common fragile site patterns of expression. By applying
robust measurements to capture in quantitative terms the non triviality of the
network, we identified several topological features clearly indicating
departure from the Erdos-Renyi random graph model. The most important outcome
was the presence of an unexpected large connected component far below the
percolation threshold. Most of the best characterized common fragile sites
belonged to this connected component. By filtering this connected component
with Gene Ontology, statistically significant shared functional features were
detected. Common fragile sites were found to be enriched for genes associated
to the immune response and to mechanisms involved in tumour progression such as
extracellular space remodeling and angiogenesis. Our results support the
hypothesis that fragile sites serve a function; we propose that fragility is
linked to a coordinated regulation of fragile genes expression.