  There is an intrinsic relationship between the molecular evolution in
primordial period and the properties of genomes and proteomes of contemporary
species. The genomic data may help us understand the driving force of evolution
of life at molecular level. In absence of evidence, numerous problems in
molecular evolution had to fall into a twilight zone of speculation and
controversy in the past. Here we show that delicate structures of variations of
genomic base compositions and amino acid frequencies resulted from the genetic
code evolution. And the driving force of evolution of life also originated in
the genetic code evolution. The theoretical results on the variations of amino
acid frequencies and genomic base compositions agree with the experimental
observations very well, not only in the variation trends but also in some fine
structures. Inversely, the genomic data of contemporary species can help
reconstruct the genetic code chronology and amino acid chronology in primordial
period. Our results may shed light on the intrinsic mechanism of molecular
evolution and the genetic code evolution.


  We numerically investigate the six-species predator-prey game in complex
networks as well as in $d$-dimensional hypercubic lattices with $d=1,2,..., 6$.
The interaction topology of the six species contains two loops, each of which
is composed of cyclically predating three species. As the mutation rate $P$ is
lowered below the well-defined phase transition point, the $Z_2$ symmetry
related with the interchange of the two loops is spontaneously broken, and it
has been known that the system develops the defensive alliance in which three
cyclically predating species defend each other against the invasion of other
species. In the small-world network structure characterized by the rewiring
probability $\alpha$, the phase diagram shows the reentrant behavior as
$\alpha$ is varied, indicating a twofold role of the shortcuts. In
$d$-dimensional regular hypercubic lattices, the system also exhibits the
reentrant phase transition as $d$ is increased. We identify universality class
of the phase transition and discuss the proper mean-field limit of the system.


  We propose a principle of consistency between different hierarchical levels
of biological systems. Given a consistency between molecule replication and
cell reproduction, universal statistical laws on cellular chemical abundances
are derived and confirmed experimentally. They include a power law distribution
of gene expressions, a lognormal distribution of cellular chemical abundances
over cells, and embedding of the power law into the network connectivity
distribution. Second, given a consistency between genotype and phenotype, a
general relationship between phenotype fluctuations by genetic variation and
isogenic phenotypic fluctuation by developmental noise is derived. Third, we
discuss the chaos mechanism for stem cell differentiation with autonomous
regulation, resulting from a consistency between cell reproduction and growth
of the cell ensemble.


  Natively unfolded proteins exist as an ensemble of flexible conformations
lacking a well defined tertiary structure along a large portion of their
polypeptide chain. Despite the absence of a stable configuration, they are
involved in important cellular processes. In this work we used from three
indicators of folding status, derived from the analysis of mean packing and
mean contact energy of a protein sequence as well as from VSL2, a disorder
predictor, and we combined them into a consensus score to identify natively
unfolded proteins in several genomes from Archaea, Bacteria and Eukarya. We
found a high correlation among the number of predicted natively unfolded
proteins and the number of proteins in the genomes. More specifically, the
number of natively unfolded proteins scaled with the number of proteins in the
genomes, with exponent 1.81 +- 0.10. This scaling law may be important to
understand the relation between the number of natively unfolded proteins and
their roles in cellular processes.


  Protein domains are found on genomes with notable statistical distributions,
which bear a high degree of similarity. Previous work has shown how these
distributions can be accounted for by simple models, where the main ingredients
are probabilities of duplication, innovation, and loss of domains. However, no
one so far has addressed the issue that these distributions follow definite
trends depending on protein-coding genome size only. We present a stochastic
duplication/innovation model, falling in the class of so-called Chinese
Restaurant Processes, able to explain this feature of the data. Using only two
universal parameters, related to a minimal number of domains and to the
relative weight of innovation to duplication, the model reproduces two
important aspects: (a) the populations of domain classes (the sets, related to
homology classes, containing realizations of the same domain in different
proteins) follow common power-laws whose cutoff is dictated by genome size, and
(b) the number of domain families is universal and markedly sublinear in genome
size. An important ingredient of the model is that the innovation probability
decreases with genome size. We propose the possibility to interpret this as a
global constraint given by the cost of expanding an increasingly complex
interactome.


  The tolerance of bacterial populations to biocidal or antibiotic treatment
has been well documented in both biofilm and planktonic settings. However,
there is still very little known about the mechanisms that produce this
tolerance. Evidence that small, non-mutant subpopulations of bacteria are not
affected by antibiotic challenge has been accumulating and provides an
attractive explanation for the failure of typical dosing protocols. Although a
dosing challenge can kill all the susceptible bacteria, the remaining persister
cells can serve as a source of population regrowth. We give a robust condition
for the failure of a periodic dosing protocol for a general chemostat model,
which supports the mathematical conclusions and simulations of an earlier, more
specialized batch model. Our condition implies that the treatment protocol
fails globally, in the sense that a mixed bacterial population will ultimately
persist above a level that is independent of the initial composition of the
population. We also give a sufficient condition for treatment success, at least
for initial population compositions near the steady state of interest,
corresponding to bacterial washout. Finally, we investigate how the speed at
which the bacteria are wiped out depends on the duration of administration of
the antibiotic. We find that this dependence is not necessarily monotone,
implying that optimal dosing does not necessarily correspond to continuous
administration of the antibiotic. Thus, genuine periodic protocols can be more
advantageous in treating a wide variety of bacterial infections.


  The Dynamic Template model of biological cell membranes and the cytoplasm as
spatially organised fluid layers is extended to plant cells, and is shown to
offer a feasible shear driven mechanism for the co-alignment of internal and
external fibres observed during growth and tropic responses


  The Dynamic Template is applied to muscle cells, and is shown to provide a
feedback mechanism that can co-ordinate the rearrangements of cytoskeletal
proteins that accompany contraction


  Scaffolding proteins that direct the assembly of multiple kinases into a
spatially localized signaling complex are often essential for the maintenance
of an appropriate biological response. Although scaffolds are widely believed
to have dramatic effects on the dynamics of signal propagation, the mechanisms
that underlie these consequences are not well understood. Here, Monte Carlo
simulations of a model kinase cascade are used to investigate how the temporal
characteristics of signaling cascades can be influenced by the presence of
scaffold proteins. Specifically, we examine the effects of spatially localizing
kinase components on a scaffold on signaling dynamics. The simulations indicate
that a major effect that scaffolds exert on the dynamics of cell signaling is
to control how the activation of protein kinases is distributed over time.
Scaffolds can influence the timing of kinase activation by allowing for kinases
to become activated over a broad range of times, thus allowing for signaling at
both early and late times. Scaffold concentrations that result in optimal
signal amplitude also result in the broadest distributions of times over which
kinases are activated. These calculations provide insights into one mechanism
that describes how the duration of a signal can potentially be regulated in a
scaffold mediated protein kinase cascade. Our results illustrate another
complexity in the broad array of control properties that emerge from the
physical effects of spatially localizing components of kinase cascades on
scaffold proteins.


  We seek a route to the equilibrium where all the agents cooperate in the
iterated prisoner's dilemma game on a two-dimensional plane, focusing on the
role of tit-for-tat strategy. When a time horizon, within which a strategy can
recall the past, is one time step, an equilibrium can be achieved as
cooperating strategies dominate the whole population via proliferation of
tit-for-tat. Extending the time horizon, we filter out poor strategies by
simplified replicator dynamics and observe a similar evolutionary pattern to
reach the cooperating equilibrium. In particular, the rise of a modified
tit-for-tat strategy plays a central role, which implies how a robust strategy
is adopted when provided with an enhanced memory capacity.


  The interaction cutoff contribution to the ruggedness of protein-protein
energy landscape (the artificial ruggedness) is studied in terms of relative
energy fluctuations for 1/r^n potentials based on a simplistic model of a
protein complex. Contradicting the principle of minimal frustration, the
artificial ruggedness exists for short cutoffs and gradually disappears with
the cutoff increase. The critical values of the cutoff were calculated for each
of eleven popular power-type potentials with n=0-9, 12 and for two thresholds
of 5% and 10%. The artificial ruggedness decreases to tolerable thresholds for
cutoffs longer than the critical ones. The results showed that for both
thresholds the critical cutoff is a non-monotonic function of the potential
power n. The functions reach the maximum at n=3-4 and then decrease with the
increase of the potential power. The difference between two cutoffs for 5% and
10% artificial ruggedness becomes negligible for potentials decreasing faster
than 1/r^12. The results suggest that cutoffs longer than critical ones can be
recommended for protein-protein potentials.


  An investigation was conducted to study the robustness of the results
obtained from the cellular automata model which describes the spread of the HIV
infection within lymphoid tissues [R. M. Zorzenon dos Santos and S. Coutinho,
Phys. Rev. Lett. 87, 168102 (2001)]. The analysis focussed on the dynamic
behavior of the model when defined in lattices with different symmetries and
dimensionalities. The results illustrated that the three-phase dynamics of the
planar models suffered minor changes in relation to lattice symmetry variations
and, while differences were observed regarding dimensionality changes,
qualitative behavior was preserved. A further investigation was conducted into
primary infection and sensitiveness of the latency period to variations of the
model's stochastic parameters over wide ranging values. The variables
characterizing primary infection and the latency period exhibited power-law
behavior when the stochastic parameters varied over a few orders of magnitude.
The power-law exponents were approximately the same when lattice symmetry
varied, but there was a significant variation when dimensionality changed from
two to three. The dynamics of the three-dimensional model was also shown to be
insensitive to variations of the deterministic parameters related to cell
resistance to the infection, and the necessary time lag to mount the specific
immune response to HIV variants. The robustness of the model demonstrated in
this work reinforce that its basic hypothesis are consistent with the
three-stage dynamic of the HIV infection observed in patients.


  Myosin V was the first discovered processive motor from the myosin family. It
has therefore been subject of a number of mechanical, kinetic, optical and
structural studies and now belongs to the best characterised motor proteins.
This effort has been accompanied by a number of different theoretical models.
In this article we give an overview of them and discuss what they have in
common and where the open questions are. The latter include the existence of
sub-steps, the process that limits the run length, the nature of backward
steps, the flexibility of the lever arm and the state of the lead head.


  The origin of diversification and coexistence of genes and species have been
traditionally studied in isolated biological levels. Ecological and
evolutionary views have focused on the mechanisms that enable or constrain
species coexistence, genetic variation and the genetics of speciation, but a
unified theory linking those approaches is still missing. Here we introduce
evolutionary graphs in the context of neutral theories of molecular evolution
and biodiversity to provide a framework that simultaneously addresses
speciation rate and joint genetic and species diversities. To illuminate this
question we also study two models of evolution on graphs with fitness
differences, which provide insights on how genetic and ecological dynamics
drive the speed of diversification. Neutral evolution generates the highest
speed of speciation, species richness (i.e. five times and twice as many
species as compared to genetic and ecological graphs, respectively) and
genetic--species diversity (i.e., twice as many as genetic and ecological
graphs, respectively). Thus the speed of speciation, the genetic--species
diversity and coexistence can differ dramatically depending on whether genetic
factors versus ecological factors drive the evolution of the system. By linking
molecular, sexual and trophic behavior at ecological and evolutionary scales,
interacting graphs can illuminate the origin and evolution of diversity and
organismal coexistence.


  Mechanical characterization of protein molecules has played a role on gaining
insight into the biological functions of proteins, since some proteins perform
the mechanical function. Here, we present the mesoscopic model of biological
protein materials composed of protein crystals prescribed by Go potential for
characterization of elastic behavior of protein materials. Specifically, we
consider the representative volume element (RVE) containing the protein
crystals represented by alpha-carbon atoms, prescribed by Go potential, with
application of constant normal strain to RVE. The stress-strain relationship
computed from virial stress theory provides the nonlinear elastic behavior of
protein materials and their mechanical properties such as Young's modulus,
quantitatively and/or qualitatively comparable to mechanical properties of
biological protein materials obtained from experiments and/or atomistic
simulations. Further, we discuss the role of native topology on the mechanical
properties of protein crystals. It is shown that parallel strands (hydrogen
bonds in parallel) enhances the mechanical resilience of protein materials.


  Cell is the structural and functional unit of life. This Resource Letter
serves as a guide to the literature on nano-machines which drive not only
intracellular movements, but also motility of the cell. These machines are
usually proteins or macromolecular assemblies which require appropriate fuel
for their operations. Although, traditionally, these machines were subjects of
investigation in biology and biochemistry, increasing use of the concepts and
techniques of physics in recent years have contributed to the quantitative
understanding of the fundamental principles underlying their operational
mechanisms. The possibility of exploiting these principles for the design and
control of artificial nano-machines has opened up a new frontier in the
bottom-up approach to nano-technology.


  The nature of epistasis has important consequences for the evolutionary
significance of sex and recombination. Recent efforts to find negative
epistasis as source of negative linkage disequilibrium and associated long-term
sex advantage have yielded little support. Sign epistasis, where the sign of
the fitness effects of alleles varies across genetic backgrounds, is
responsible for ruggedness of the fitness landscape with implications for the
evolution of sex that have been largely unexplored. Here, we describe fitness
landscapes for two sets of strains of the asexual fungus \emph{Aspergillus
niger} involving all combinations of five mutations. We find that $\sim 30$% of
the single-mutation fitness effects are positive despite their negative effect
in the wild-type strain, and that several local fitness maxima and minima are
present. We then compare adaptation of sexual and asexual populations on these
empirical fitness landscapes using simulations. The results show a general
disadvantage of sex on these rugged landscapes, caused by the break down by
recombination of genotypes escaping from local peaks. Sex facilitates escape
from a local peak only for some parameter values on one landscape, indicating
its dependence on the landscape's topography. We discuss possible reasons for
the discrepancy between our results and the reports of faster adaptation of
sexual populations.


  Background: A partially random target selection method was developed to
design and produce affinity reagents (target) to any protein query. It is based
on the recent concept of Proteomic Code (for review see Biro, 2007 [1]) which
suggests that significant number of amino acids in specifically interacting
proteins are coded by partially complementary codons. It means that the 1st and
3rd residues of codons coding many co-locating amino acids are complementary
but the 2nd may but not necessarily complementary: like 5'-AXG-3'/3'-CXT-5'
codon pair, where X is any nucleotide.
  Results: A mixture of 45 residue long, reverse, partially complementary
oligonucleotide sequences (target pool) were synthesized to selected epitopes
of query mRNA sequences. The 2nd codon residues were randomized. The target
oligonucleotide pool was inserted into vectors, expressed and the protein
products were screened for affinity to the query in Bacterial Two-Hybrid
System. The best clones were used for larger-scale protein syntheses and
characterization. It was possible to design and produce specific and with high
affinity reacting (Kd:100 nM) oligopeptide reagents to GAL4 query
oligopeptides.
  Conclusions: Second codon residue randomization is a promising method to
design and produce affinity peptides to any protein sequences. The method has
the potential to be a rapid, inexpensive, high throughput, non-immunoglobulin
based alternative to recent in vivo antibody generating procedures.


  Questions on possible relationship between phenotypic plasticity and
evolvability, as well as that between robustness and evolution have been
addressed over decades in the field of evolution-development. By introducing an
evolutionary stability assumption on the distribution of phenotype and
genotype, we establish quantitative relationships on plasticity, phenotypic
fluctuations, and evolvability. Derived are proportionality among plasticity as
a responsiveness of phenotype against environmental change, variances of
phenotype fluctuations of genetic and developmental origins, and evolution
speed. Confirmation of the relationships is given by numerical experiments of a
gene expression dynamics model with an evolving transcription network, whereas
verifications by laboratory evolution experiments are also discussed. These
results provide quantitative formulation on canalization and genetic
assimilation, in terms of fluctuations of gene expression levels.


  The cytoplasm of a living cell is crowded with several macromolecules of
different shapes and sizes. Molecular diffusion in such a medium becomes
anomalous due to the presence of macromolecules and diffusivity is expected to
decrease with increase in macromolecular crowding. Moreover, many cellular
processes are dependent on molecular diffusion in the cell cytosol. The
enzymatic reaction rate has been shown to be affected by the presence of such
macromolecules. A simple numerical model is proposed here based on percolation
and diffusion in disordered systems to study the effect of macromolecular
crowding on the enzymatic reaction rates. The model explains qualitatively some
of the experimental observations.


  Social networks exhibit strikingly systematic patterns across a wide range of
human contexts. While genetic variation accounts for a significant portion of
the variation in many complex social behaviors, the heritability of egocentric
social network attributes is unknown. Here we show that three of these
attributes (in-degree, transitivity, and centrality) are heritable. We then
develop a "mirror network" method to test extant network models and show that
none accounts for observed genetic variation in human social networks. We
propose an alternative "Attract and Introduce" model with two simple forms of
heterogeneity that generates significant heritability as well as other
important network features. We show that the model is well suited to real
social networks in humans. These results suggest that natural selection may
have played a role in the evolution of social networks. They also suggest that
modeling intrinsic variation in network attributes may be important for
understanding the way genes affect human behaviors and the way these behaviors
spread from person to person.


  Heat-induced mobility of nucleosomes along DNA is an experimentally
well-studied phenomenon. A recent experiment shows that the repositioning is
modified in the presence of minor-groove binding DNA ligands. We present here a
stochastic three-state model for the diffusion of a nucleosome along DNA in the
presence of such ligands. It allows us to describe the dynamics and the steady
state of such a motion analytically. The analytical results are in excellent
agreement with numerical simulations of this stochastic process.With this
model, we study the response of a nucleosome to an external force and how it is
affected by the presence of ligands.


  Evolutionary mechanism in a self-organized system cause some functional
changes that force to adapt new conformation of the interaction pattern between
the components of that system. Measuring the structural differences one can
retrace the evolutionary relation between two systems. We present a method to
quantify the topological distance between two networks of different sizes,
finding that the architectures of the networks are more similar within the same
class than the outside of their class. With 43 cellular networks of different
species, we show that the evolutionary relationship can be elucidated from the
structural distances.


  In this paper, we explore the feasibility of using coarse-grained models to
simulate the self-assembly of DNA nanostructures. We introduce a simple model
of DNA where each nucleotide is represented by two interaction sites
corresponding to the phosphate-sugar backbone and the base. Using this model,
we are able to simulate the self-assembly of both DNA duplexes and Holliday
junctions from single-stranded DNA. We find that assembly is most successful in
the temperature window below the melting temperatures of the target structure
and above the melting temperature of misbonded aggregates. Furthermore, in the
case of the Holliday junction, we show how a hierarchical assembly mechanism
reduces the possibility of becoming trapped in misbonded configurations. The
model is also able to reproduce the relative melting temperatures of different
structures accurately, and allows strand displacement to occur.


  Our aim is to build a set of rules, such that reasoning over temporal
dependencies within gene regulatory networks is possible. The underlying
transitions may be obtained by discretizing observed time series, or they are
generated based on existing knowledge, e.g. by Boolean networks or their
nondeterministic generalization. We use the mathematical discipline of formal
concept analysis (FCA), which has been applied successfully in domains as
knowledge representation, data mining or software engineering. By the attribute
exploration algorithm, an expert or a supporting computer program is enabled to
decide about the validity of a minimal set of implications and thus to
construct a sound and complete knowledge base. From this all valid implications
are derivable that relate to the selected properties of a set of genes. We
present results of our method for the initiation of sporulation in Bacillus
subtilis. However the formal structures are exhibited in a most general manner.
Therefore the approach may be adapted to signal transduction or metabolic
networks, as well as to discrete temporal transitions in many biological and
nonbiological areas.


  Background: The secondary structure and complexity of mRNA influences its
accessibility to regulatory molecules (proteins, micro-RNAs), its stability and
its level of expression. The mobile elements of the RNA sequence, the wobble
bases, are expected to regulate the formation of structures encompassing coding
sequences.
  Results: The sequence/folding energy (FE) relationship was studied by
statistical, bioinformatic methods in 90 CDS containing 26,370 codons. I found
that the FE (dG) associated with coding sequences is significant and negative
(407 kcal/1000 bases, mean +/- S.E.M.) indicating that these sequences are able
to form structures. However, the FE has only a small free component, less than
10% of the total. The contribution of the 1st and 3rd codon bases to the FE is
larger than the contribution of the 2nd (central) bases. It is possible to
achieve a ~ 4-fold change in FE by altering the wobble bases in synonymous
codons. The sequence/FE relationship can be described with a simple algorithm,
and the total FE can be predicted solely from the sequence composition of the
nucleic acid. The contributions of different synonymous codons to the FE are
additive and one codon cannot replace another. The accumulated contributions of
synonymous codons of an amino acid to the total folding energy of an mRNA is
strongly correlated to the relative amount of that amino acid in the translated
protein.
  Conclusion: Synonymous codons are not interchangable with regard to their
role in determining the mRNA FE and the relative amounts of amino acids in the
translated protein, even if they are indistinguishable in respect of amino acid
coding.


  The analysis of eight molecular datasets involving human and teleost examples
along with morphological samples from several groups of Neotropical electric
fish (Order: Gymnotiformes) were used in this thesis to test the dynamics of
both intraspecific variation and interspecific diversity. In terms of
investigating molecular interspecific diversity among humans, two experimental
exercises were performed. A cladistic exchange experiment tested for the extent
of discontinuity and interbreeding between H. sapiens and neanderthal
populations. As part of the same question, another experimental exercise tested
the amount of molecular variance resulting from simulations which treated
neanderthals as being either a local population of modern humans or as a
distinct subspecies. Finally, comparisons of hominid populations over time with
fish species helped to define what constitutes taxonomically relevant
differences between morphological populations as expressed among both trait
size ranges and through growth patterns that begin during ontogeny. Compared to
the subdivision found within selected teleost species, H. sapiens molecular
data exhibited little variation and discontinuity between geographical regions.
Results of the two experimental exercises concluded that neanderthals exhibit
taxonomic distance from modern H. sapiens. However, this distance was not so
great as to exclude the possibility of interbreeding between the two
subspecific groups. Finally, a series of characters were analyzed among species
of Neotropical electric fish. These analyses were compared with hominid
examples to determine what constituted taxonomically relevant differences
between populations as expressed among specific morphometric traits that
develop during the juvenile phase.


  The phenomenon of gene conservation is an interesting evolutionary problem
related to speciation and adaptation. Conserved genes are acted upon in
evolution in a way that preserves their function despite other structural and
functional changes going on around them. The recent availability of
whole-genomic data from closely related species allows us to test the
hypothesis that a core genome present in a hypothetical common ancestor is
inherited by all sister taxa. Furthermore, this core genome should serve
essential functions such as genetic regulation and cellular repair.
Whole-genome sequences from three strains of bacteria (Shewanella sp.) were
used in this analysis. The open reading frames (ORFs) for each identified and
putative gene were used for each genome. Reciprocal Blast searches were
conducted on all three genomes, which distilled a list of thousands of genes to
68 genes that were identical across taxa. Based on functional annotation, these
genes were identified as housekeeping genes, which confirmed the original
hypothesis. This method could be used in eukaryotes as well, in particular the
relationship between humans, chimps, and macaques.


  We present an alternative method for calculating likelihoods in molecular
phylogenetics. Our method is based on partial likelihood tensors, which are
generalizations of partial likelihood vectors, as used in Felsenstein's
approach. Exploiting a lexicographic sorting and partial likelihood tensors, it
is possible to obtain significant computational savings. We show this on a
range of simulated data by enumerating all numerical calculations that are
required by our method and the standard approach.


  We propose two simple probability models to compute the probability of an
influenza pandemic. Under a random walk model the probability that all
pandemics between times 0 and 300 occur by time 150 is 1/2. Under a Poisson
model with mean inter arrival time of 30 years the probability that no pandemic
occurs during at least 60 years is 14%. These probabilities are much higher
than generally perceived. So yes the next influenza pandemic will happen but
maybe much later than generally thought.


  Cancer progression is an evolutionary process that is driven by mutation and
selection in a population of tumor cells. We discuss mathematical models of
cancer progression, starting from traditional multistage theory. Each stage is
associated with the occurrence of genetic alterations and their fixation in the
population. We describe the accumulation of mutations using conjunctive
Bayesian networks, an exponential family of waiting time models in which the
occurrence of mutations is constrained to a partial temporal order. Two
opposing limit cases arise if mutations either follow a linear order or occur
independently. We derive exact analytical expressions for the waiting time
until a specific number of mutations have accumulated in these limit cases as
well as for the general conjunctive Bayesian network. Finally, we analyze a
stochastic population genetics model that explicitly accounts for mutation and
selection. In this model, waves of clonal expansions sweep through the
population at equidistant intervals. We present an approximate analytical
expression for the waiting time in this model and compare it to the results
obtained for the conjunctive Bayesian networks.


  We study a model for a protein searching for a target, using facilitated
diffusion, on a DNA molecule confined in a finite volume. The model includes
three distinct pathways for facilitated diffusion: (a) sliding - in which the
protein diffuses along the contour of the DNA (b) jumping - where the protein
travels between two sites along the DNA by three-dimensional diffusion, and
finally (c) intersegmental transfer - which allows the protein to move from one
site to another by transiently binding both at the same time. The typical
search time is calculated using scaling arguments which are verified
numerically. Our results suggest that the inclusion of intersegmental transfer
(i) decreases the search time considerably (ii) makes the search time much more
robust to variations in the parameters of the model and (iii) that the optimal
search time occurs in a regime very different than that found for models which
ignore intersegmental transfers. The behavior we find is rich and shows
surprising dependencies, for example, on the DNA length.


  Background: There is a 3-fold redundancy in the Genetic Code; most amino
acids are encoded by more than one codon. These synonymous codons are not used
equally; there is a Codon Usage Bias (CUB). This article will provide novel
information about the origin and evolution of this bias.
  Results: Codon Usage Bias (CUB, defined here as deviation from equal usage of
synonymous codons) was studied in 113 species. The average CUB was 29.3 +/-
1.1% (S.E.M, n=113) of the theoretical maximum and declined progressively with
evolution and increasing genome complexity. A Pan-Genomic Codon Usage Frequency
(CUF) Table was constructed to describe genome-wide relationships among codons.
Significant correlations were found between the number of synonymous codons and
(i) the frequency of the respective amino acids (ii) the size of CUB. Numerous,
statistically highly significant, internal correlations were found among codons
and the nucleic acids they comprise. These strong correlations made it possible
to predict missing synonymous codons (wobble bases) reliably from the remaining
codons or codon residues.
  Conclusions: The results put the concept of "codon bias" into a novel
perspective. The internal connectivity of codons indicates that all synonymous
codons might be integrated parts of the Genetic Code with equal importance in
maintaining its functional integrity.


  The lac operon in Escherichia coli has been studied extensively and is one of
the earliest gene systems found to undergo both positive and negative control.
The lac operon is known to exhibit bistability, in the sense that the operon is
either induced or uninduced. Many dynamical models have been proposed to
capture this phenomenon. While most are based on complex mathematical
formulations, it has been suggested that for other gene systems network
topology is sufficient to produce the desired dynamical behavior.
  We present a Boolean network as a discrete model for the lac operon. We
include the two main glucose control mechanisms of catabolite repression and
inducer exclusion in the model and show that it exhibits bistability. Further
we present a reduced model which shows that lac mRNA and lactose form the core
of the lac operon, and that this reduced model also exhibits the same dynamics.
This work corroborates the claim that the key to dynamical properties is the
topology of the network and signs of interactions.


  This paper proposes rigorous geometrical treatment of bioelectrodynamics,
underpinning two fast-growing biomedical research fields: bioelectromagnetism,
which deals with the ability of life to produce its own electromagnetism, and
bioelectromagnetics, which deals with the effect on life from external
electromagnetism.
  Keywords: Bioelectrodynamics, exterior geometrical machinery, Dirac-Feynman
quantum electrodynamics, functional electrical stimulation


  Survival probability within a certain time horizon T is a common measure of
population viability. The choice of T implicitly involves a time preference,
similar to economic discounting: Conservation success is evaluated at the time
horizon T, while all effects that occur later than T are not considered.
Despite the obvious relevance of the time horizon, ecological studies seldom
analyze its impact on the evaluation of conservation options. In this paper, we
show that, while the choice of T does not change the ranking of conservation
options for single species under stationary conditions, it may substantially
change conservation decisions for multiple species. We conclude that it is of
crucial importance to investigate the sensitivity of model results to the
choice of the time horizon or other measures of time preference when
prioritizing biodiversity conservation efforts.


  Understanding the patterns and processes of diversification of life in the
planet is a key challenge of science. The Tree of Life represents such
diversification processes through the evolutionary relationships among the
different taxa, and can be extended down to intra-specific relationships. Here
we examine the topological properties of a large set of interspecific and
intraspecific phylogenies and show that the branching patterns follow
allometric rules conserved across the different levels in the Tree of Life, all
significantly departing from those expected from the standard null models. The
finding of non-random universal patterns of phylogenetic differentiation
suggests that similar evolutionary forces drive diversification across the
broad range of scales, from macro-evolutionary to micro-evolutionary processes,
shaping the diversity of life on the planet.


  We study the dynamics of a knot in a semiflexible polymer confined to a
narrow channel of width comparable to the polymers' persistence length. Using a
combination of Brownian dynamics simulations and a coarse-grained stochastic
model, we characterize the coupled dynamics of knot size variation and knot
diffusion along the polymer, which ultimately leads to spontaneous unknotting.
We find that the knot grows to macroscopic size before disappearing.
Interestingly, an external force applied to the ends of the confined polymer
speeds up spontaneous unknotting.


  We simulate the assembly dynamics of icosahedral capsids from subunits that
interconvert between different conformations (or quasi-equivalent states). The
simulations identify mechanisms by which subunits form empty capsids with only
one morphology but adaptively assemble into different icosahedral morphologies
around nanoparticle cargoes with varying sizes, as seen in recent experiments
with brome mosaic virus (BMV) capsid proteins. Adaptive cargo encapsidation
requires moderate cargo-subunit interaction strengths; stronger interactions
frustrate assembly by stabilizing intermediates with incommensurate curvature.
We compare simulation results to experiments with cowpea chlorotic mottle virus
empty capsids and BMV capsids assembled on functionalized nanoparticles and
suggest new cargo encapsidation experiments. Finally, we find that both empty
and templated capsids maintain the precise spatial ordering of subunit
conformations seen in the crystal structure even if interactions that preserve
this arrangement are favored by as little as the thermal energy, consistent
with experimental observations that different subunit conformations are highly
similar.


  Based on a multi-agent model, we investigate how target waves emerge from a
population dynamics with cyclical interactions among three species. We show
that the periodically injecting source in a small central area can generate
target waves in a two-dimensional lattice system. By detecting the temporal
period of species' concentration at the central area, three modes of target
waves can be distinguished. Those different modes result from the competition
between local and global oscillations induced by cyclical interactions: Mode A
corresponds to a synchronization of local and global oscillations, Mode B
results from an intermittent synchronization, and Mode C corresponds to the
case when the frequency of the local oscillation is much higher than that of
the global oscillation. This work provides insights into pattern formation in
biologic and ecologic systems that are totally different from the extensively
studied diffusion systems driven by chemical reactions.


  The secretion of vesicles for intracellular transport often rely on the
aggregation of specialized membrane-bound proteins into a coat able to curve
cell membranes. The nucleation and growth of a protein coat is a kinetic
process that competes with the energy-consuming turnover of coat components
between the membrane and the cytosol. We propose a generic kinetic description
of coat assembly and the formation of coated vesicles, and discuss its
implication to the dynamics of COP vesicles that traffic within the Golgi and
with the Endoplasmic Reticulum. We show that stationary coats of fixed area
emerge from the competition between coat growth and the recycling of coat
components, in a fashion resembling the treadmilling of cytoskeletal filaments.
We further show that the turnover of coat components allows for a highly
sensitive switching mechanism between a quiescent and a vesicle producing
membrane, upon a slowing down of the exchange kinetics. We claim that the
existence of this switching behaviour, also triggered by factors such as the
presence of cargo and variation of the membrane mechanical tension, allows for
efficient regulation of vesicle secretion. We propose a model, supported by
different experimental observations, in which vesiculation of secretory
membranes is impaired by the energy consuming desorption of coat proteins,
until the presence of cargo or other factors triggers a dynamical switch into a
vesicle producing state.


  Time series analysis of fossil biodiversity of marine invertebrates in the
Paleobiology Database (PBDB) shows a significant periodicity at approximately
63 My, in agreement with previous analyses based on the Sepkoski database. I
discuss how this result did not appear in a previous analysis of the PBDB. The
existence of the 63 My periodicity, despite very different treatment of
systematic error in both PBDB and Sepkoski databases strongly argues for
consideration of its reality in the fossil record. Cross-spectral analysis of
the two datasets finds that a 62 My periodicity coincides in phase by 1.6 My,
equivalent to better than the errors in either measurement. Consequently, the
two data sets not only contain the same strong periodicity, but its peaks and
valleys closely correspond in time. Two other spectral peaks appear in the PBDB
analysis, but appear to be artifacts associated with detrending and with the
increased interval length. Sampling-standardization procedures implemented by
the PBDB collaboration suggest that the signal is not an artifact of sampling
bias. Further work should focus on finding the cause of the 62 My periodicity.


  The genetic code is nearly universal, and the arrangement of the codons in
the standard codon table is highly non-random. The three main concepts on
origin and evolution of the code are the stereochemical theory; the coevolution
theory; and the error minimization theory. These theories are not mutually
exclusive and are also compatible with the frozen accident hypothesis.
Mathematical analysis of the structure and possible evolutionary trajectories
of the code shows that it is highly robust to translational error but there is
a huge number of more robust codes, so that the standard code potentially could
evolve from a random code via a short sequence of codon series reassignments.
Thus, much of the evolution that led to the standard code can be interpreted as
a combination of frozen accident with selection for translational error
minimization although contributions from coevolution of the code with metabolic
pathways and/or weak affinities between amino acids and nucleotide triplets
cannot be ruled out. However, such scenarios for the code evolution are based
on formal schemes whose relevance to the actual primordial evolution is
uncertain, so much caution in interpretation is necessary. A real understanding
of the code's origin and evolution is likely to be attainable only in
conjunction with a credible scenario for the evolution of the coding principle
itself and the translation system.


  Unfolded proteins may contain native or non-native residual structure, which
has important implications for the thermodynamics and kinetics of folding as
well as for misfolding and aggregation diseases. However, it has been
universally accepted that residual structure should not affect the global size
scaling of the denatured chain, which obeys the statistics of random coil
polymers. Here we use a single-molecule optical technique, fluorescence
correlation spectroscopy, to probe the denatured state of set of repeat
proteins containing an increasing number of identical domains, from two to
twenty. The availability of this set allows us to obtain the scaling law for
the unfolded state of these proteins, which turns out to be unusually compact,
strongly deviating from random-coil statistics. The origin of this unexpected
behavior is traced to the presence of extensive non-native polyproline II
helical structure, which we localize to specific segments of the polypeptide
chain. We show that the experimentally observed effects of PPII on the size
scaling of the denatured state can be well-described by simple polymer models.
Our findings suggest an hitherto unforeseen potential of non-native structure
to induce significant compaction of denatured proteins, affecting significantly
folding pathways and kinetics.


  We present an explicit unified stochastic model of fluctuations in population
size due to random birth, death, density-dependent competition and
environmental fluctuations. Stochastic dynamics provide insight into small
populations, including processes such as extinction, that cannot be correctly
treated by deterministic methods. We present exact analytical and
simulation-based results for extinction times of our stochastic model and
compare the different effects of environmental stochasticity and intrinsic
demographic stochasticity. We use both the discrete master equation approach
and an exact mapping to a Fokker-Planck equation (the Poisson method) and
stochastic equation, showing they are precisely equivalent. We also calculate
approximate extinction times using a steepest descent method. This model can
readily be extended to accommodate metapopulation structure and genetic
variation in the population and thus represents a step towards a
microscopically explicit synthesis of population dynamics and population
genetics.


  Segmentation in arthropod embryogenesis represents a well-known example of
body plan diversity. Striped patterns of gene expression that lead to the
future body segments appear simultaneously or sequentially in long and short
germ-band development, respectively. Regulatory genes relevant for stripe
formation are evolutionarily conserved among arthropods, therefore the
differences in the observed traits are thought to have originated from how the
genes are wired. To reveal the basic differences in the network structure, we
have numerically evolved hundreds of gene regulatory networks that produce
striped patterns of gene expression. By analyzing the topologies of the
generated networks, we show that the characteristics of stripe formation in
long and short germ-band development are determined by Feed-Forward Loops
(FFLs) and negative Feed-Back Loops (FBLs) respectively. Network architectures,
gene expression patterns and knockout responses exhibited by the artificially
evolved networks agree with those reported in the fly Drosophila melanogaster
and the beetle Tribolium castaneum. For other arthropod species, principal
network architectures that remain largely unknown are predicted.


  The extinction time of an isolated population can be exponentially reduced by
a periodic modulation of its environment. We investigate this effect using, as
an example, a stochastic branching-annihilation process with a time-dependent
branching rate. The population extinction is treated in eikonal approximation,
where it is described as an instanton trajectory of a proper reaction
Hamiltonian. The modulation of the environment perturbs this trajectory and
synchronizes it with the modulation phase. We calculate the corresponding
change in the action along the instanton using perturbation techniques
supported by numerical calculations. The techniques include a first-order
theory with respect to the modulation amplitude, a second-order theory in the
spirit of the Kapitsa pendulum effect, and adiabatic theory valid for low
modulation frequencies.


  Recently, a number of drug-therapy, disease, drug, and drug-target networks
have been introduced. Here we suggest novel methods for network-based
prediction of novel drug targets and for improvement of drug efficiency by
analysing the effects of drugs on the robustness of cellular networks.


  The depletion theory of nanoparticles immersed in a semidilute polymer
solution is reinterpreted in terms of depleted chains of polymer segments.
Limitations and extensions of mean-field theory are discussed. An explicit
expression for the interaction between two small spheres is derived. The
depletion free energy for a particle of general shape is given in terms of the
capacitance or effective Stokes radius. This affords a close to quantitative
explanation for the effect of polymer on protein precipitation.


  Transcription is regulated through interplay between transcription factors,
an RNA polymerase(RNAP), and a promoter. Even for a simple repressive
transcription factor that disturbs promoter activity at the initial binding of
RNAP, its repression level is not determined solely by the dissociation
constant of transcription factor but is sensitive to the time scales of
processes in RNAP. We first analyse the promoter activity under strong
repression by a slow binding repressor, in which case transcriptions occur in a
burst, followed by a long quiescent period while a repressor binds to the
operator; the number of transcriptions, the bursting and the quiescent times
are estimated by reaction rates. We then examine interference effect from an
opposing promoter, using the correlation function of transcription initiations
for a single promoter. The interference is shown to de-repress the promoter
because RNAP's from the opposing promoter most likely encounter the repressor
and remove it in case of strong repression. This de-repression mechanism should
be especially prominent for the promoters that facilitate fast formation of
open complex with the repressor whose binding rate is slower than about 1/sec.
Finally, we discuss possibility of this mechanism for high activity of promoter
PR in the hyp-mutant of lambda phage.


  We propose, in this article, an analysis of the Darwin's approach to
sociality. Sociality is perfectly integrated into the selective model, and is
caused by the same process as struggle for existence. Thus, the selective
process does not prohibit demonstrating cooperation, but on the contrary,
explains it. We will see that the unification of sociality, in a more adequate
representation to mammals, is not an expression of a naive anthropomorphism,
but of a methodological anthropomorphism; normal consequence of his gradualist
approach of phylogeny.


  Background The theoretical requirements for a genetic code were well defined
and modeled by George Gamow and Francis Crick in the 50-es. Their models
failed. However the valid Genetic Code, provided by Nirenberg and Matthaei in
1961, ignores many theoretical requirements for a perfect Code. Something is
simply missing from the canonical Code.
  Results The 3x redundancy of the Genetic code is usually explained as a
necessity to increase the resistance of the mutation resistance of the genetic
information. However it has many additional roles. 1.) It has a periodical
structure which corresponds to the physico-chemical and structural properties
of amino acids. 2.) It provides physico-chemical definition of codon
boundaries. 3.) It defines a code for amino acid co-locations (interactions) in
the coded proteins. 4.) It regulates, through wobble bases the free folding
energy (and structure) of mRNAs. I shortly review the history of the Genetic
Code as well as my own published observations to provide a novel, original
explanation of its redundancy.
  Conclusions The redundant Genetic Code contains biological information which
is additional to the 64/20 definition of amino acids. This additional
information is used to define the 3D structure of coding nucleic acids as well
as the coded proteins and it is called the Proteomic Code and mRNA Assisted
Protein Folding.


  This is the first report, to our knowledge, on a systematic method for
constructing a large scale kinetic metabolic model with incomplete information
on kinetic parametersr, and its initial application to the modeling of central
metabolism of Methylobacterium extorquens AM1, a methylotrophic and
environmental important bacterium, with all necessary constraints. Through a
systematic and consistent procedure of finding a set of parameters in the
physiological range we overcome an outstanding difficulty in large scale
kinetic modeling: the requirement for a massive number of enzymatic reaction
parameters. We are able to construct the kinetic model based on general
biological considerations and incomplete experimental kinetic parameters. The
success of our approach with incompletely input information is guaranteed by
two known principles in biology, the robustness of the system and the
cooperation among its various parts. (Will be pleased to be informed on other
methodologies dealing with same type of problems: aoping@u.washington.edu)


  We compare three basic kinds of discrete mathematical models used to portray
phylogenetic relationships among species and higher taxa: phylogenetic trees,
Hennig trees and Nelson cladograms. All three models are trees, as that term is
commonly used in mathematics; the difference between them lies in the
biological interpretation of their vertices and edges. Phylogenetic trees and
Hennig trees carry exactly the same information, and translation between these
two kinds of trees can be accomplished by a simple algorithm. On the other
hand, evolutionary concepts such as monophyly are represented as different
mathematical substructures are represented differently in the two models. For
each phylogenetic or Hennig tree, there is a Nelson cladogram carrying the same
information, but the requirement that all taxa be represented by leaves
necessarily makes the representation less efficient. Moreover, we claim that it
is necessary to give some interpretation to the edges and internal vertices of
a Nelson cladogram in order to make it useful as a biological model. One
possibility is to interpret internal vertices as sets of characters and the
edges as statements of inclusion; however, this interpretation carries little
more than incomplete phenetic information. We assert that from the standpoint
of phylogenetics, one is forced to regard each internal vertex of a Nelson
cladogram as an actual (albeit unsampled) species simply to justify the use of
synapomorphies rather than symplesiomorphies.


  A metabolic model can be represented as bipartite graph comprising linked
reaction and metabolite nodes. Here it is shown how a network of conserved
fluxes can be assigned to the edges of such a graph by combining the reaction
fluxes with a conserved metabolite property such as molecular weight. A similar
flux network can be constructed by combining the primal and dual solutions to
the linear programming problem that typically arises in constraint-based
modelling. Such constructions may help with the visualisation of flux
distributions in complex metabolic networks. The analysis also explains the
strong correlation observed between metabolite shadow prices (the dual linear
programming variables) and conserved metabolite properties. The methods were
applied to recent metabolic models for Escherichia coli, Saccharomyces
cerevisiae, and Methanosarcina barkeri. Detailed results are reported for E.
coli; similar results were found for the other organisms.


  We consider an infinite asexual population with a mutator allele which can
elevate mutation rates. With probability $f$, a transition from nonmutator to
mutator state occurs but the reverse transition is forbidden. We find that at
$f=0$, the population is in the state with minimum mutation rate and at
$f=f_c$, a phase transition occurs between a mixed phase with both nonmutators
and mutators and a pure mutator phase. We calculate the critical probability
$f_c$ and the total mutator fraction $Q$ in the mixed phase exactly. Our
predictions for $Q$ are in agreement with those seen in microbial populations
in static environments.


  This paper deals with ratio-dependent predator-prey systems with delay. We
will investigate under what conditions delay cannot cause instability in higher
dimension. We give an example when delay causes instability.


  In this paper we extend the recent work of C.A. Braumann \cite{B2007} to the
case of stochastic differential equation with random coefficients. Furthermore,
the relationship of the It\^o-Stratonovich stochastic calculus to studies of
random population growth is also explained.


  We study the statistical properties of the SIR epidemics in heterogeneous
networks, when an epidemic is defined as only those SIR propagations that reach
or exceed a minimum size s_c. Using percolation theory to calculate the average
fractional size <M_SIR> of an epidemic, we find that the strength of the
spanning link percolation cluster $P_{\infty}$ is an upper bound to <M_SIR>.
For small values of s_c, $P_{\infty}$ is no longer a good approximation, and
the average fractional size has to be computed directly. The value of s_c for
which $P_{\infty}$ is a good approximation is found to depend on the
transmissibility T of the SIR. We also study Q, the probability that an SIR
propagation reaches the epidemic mass s_c, and find that it is well
characterized by percolation theory. We apply our results to real networks
(DIMES and Tracerouter) to measure the consequences of the choice s_c on
predictions of average outcome sizes of computer failure epidemics.


  Understanding the melting of short DNA sequences probes DNA at the scale of
the genetic code and raises questions which are very different from those posed
by very long sequences, which have been extensively studied. We investigate
this problem by combining experiments and theory. A new experimental method
allows us to make a mapping of the opening of the guanines along the sequence
as a function of temperature. The results indicate that non-local effects may
be important in DNA because an AT-rich region is able to influence the opening
of a base pair which is about 10 base pairs away. An earlier mesoscopic model
of DNA is modified to correctly describe the time scales associated to the
opening of individual base pairs well below melting, and to properly take into
account the sequence. Using this model to analyze some characteristic sequences
for which detailed experimental data on the melting is available [Montrichok et
al. 2003 Europhys. Lett. {\bf 62} 452], we show that we have to introduce
non-local effects of AT-rich regions to get acceptable results. This brings a
second indication that the influence of these highly fluctuating regions of DNA
on their neighborhood can extend to some distance.


  When it is viewed at the scale of a base pair, DNA appears as a nonlinear
lattice. Modelling its properties is a fascinating goal. The detailed
experiments that can be performed on this system impose constraints on the
models and can be used as a guide to improve them. There are nevertheless many
open problems, particularly to describe DNA at the scale of a few tens of base
pairs, which is relevant for many biological phenomena.


  Chemotaxis receptors in E. coli form clusters at the cell poles and also
laterally along the cell body, and this clustering plays an important role in
signal transduction. Recently, experiments using flourrescence imaging have
shown that, during cell growth, lateral clusters form at positions
approximately periodically spaced along the cell body. In this paper, we
demonstrate within a lattice model that such spatial organization could arise
spontaneously from a stochastic nucleation mechanism. The same mechanism may
explain the recent observation of periodic aggregates of misfolded proteins in
E. coli.


  It has been suggested that galactic shock asymmetry induced by our galaxy's
infall toward the Virgo Cluster may be a source of periodicity in cosmic ray
exposure as the solar system oscillates perpendicular to the galactic plane.
Here we investigate a mechanism by which cosmic rays might affect terrestrial
biodiversity, ionization and dissociation in the atmosphere, resulting in
depletion of ozone and a resulting increase in the dangerous solar UVB flux on
the ground, with an improved ionization background computation averaged over a
massive ensemble (about 7 x 10^5) shower simulations. We study minimal and full
exposure to the postulated extragalactic background. The atmospheric effects
are greater than with our earlier, simplified ionization model. At the lower
end of the range effects are too small to be of serious consequence. At the
upper end of the range, ~6 % global average loss of ozone column density
exceeds that currently experienced due to effects such as accumulated
chlorofluorocarbons. The intensity is less than a nearby supernova or galactic
gamma-ray burst, but the duration would be about 10^6 times longer. Present UVB
enhancement from current ozone depletion ~3% is a documented stress on the
biosphere, but a depletion of the magnitude found at the upper end of our range
would double the global average UVB flux. For estimates at the upper end of the
range of the cosmic ray variability over geologic time, the mechanism of
atmospheric ozone depletion may provide a major biological stress, which could
easily bring about major loss of biodiversity. Future high energy astrophysical
observations will resolve the question of whether such depletion is likely.


  We investigate the protein expression pattern of the lamB gene in Escherichia
coli LE392. The gene product LamB is an important membrane protein for maltose
transport into cells but it is also exploited by bacteriophage lambda for
infection. Although our bacterial population is clonal, stochastic gene
expression leads to a majority population with a large receptor number and a
minority population with a small receptor number. We find that the LamB
receptor distribution p(n) of the small-n population is scale invariant with
the exponent depending on growth conditions. A heuristic model is proposed that
relates the observed exponent to the protein production rate.


  Recent field investigations of photosynthetic bacteria living in geothermal
hot spring environments have revealed surprisingly complex ecosystems, with an
unexpected level of genetic diversity. One case of particular interest involves
the distribution along hot spring thermal gradients of genetically distinct
bacterial strains that differ in their preferred temperatures for reproduction
and photosynthesis. In such systems, a single variable, temperature, defines
the relevant environmental variation. In spite of this, each region along the
thermal gradient exhibits multiple strains of photosynthetic bacteria adapted
to several distinct thermal optima, rather than the expected single thermal
strain adapted to the local environmental temperature. Here we analyze
microbiology data from several ecological studies to show that the thermal
distribution field data exhibit several universal features independent of
location and specific bacterial strain. These include the distribution of
optimal temperatures of different thermal strains and the functional dependence
of the net population density on temperature. Further, we present a simple
population dynamics model of these systems that is highly constrained by
biophysical data and by physical features of the environment. This model can
explain in detail the observed diversity of different strains of the
photosynthetic bacteria. It also reproduces the observed thermal population
distributions, as well as certain features of population dynamics observed in
laboratory studies of the same organisms.


  Renaturation and hybridization reactions lead to the pairing of complementary
single-stranded nucleic acids. We present here a theoretical investigation of
the mechanism of these reactions in vitro under thermal conditions (dilute
solutions of single-stranded chains, in the presence of molar concentrations of
monovalent salts and at elevated temperatures). The mechanism follows a
Kramers' process, whereby the complementary chains overcome a potential barrier
through Brownian motion. The barrier originates from a single rate-limiting
nucleation event in which the first complementary base pairs are formed. The
reaction then proceeds through a fast growth of the double helix. For the DNA
of bacteriophages T7, T4 and $\phi$X174 as well as for Escherichia coli DNA,
the bimolecular rate $k_2$ of the reaction increases as a power law of the
average degree of polymerization $<N>$ of the reacting single- strands: $k_2
\prop <N>^\alpha$. This relationship holds for $100 \leq <N> \leq 50 000$ with
an experimentally determined exponent $\alpha = 0.51 \pm 0.01$. The length
dependence results from a thermodynamic excluded-volume effect. The reacting
single-stranded chains are predicted to be in universal good solvent
conditions, and the scaling law is determined by the relevant equilibrium
monomer contact probability. The value theoretically predicted for the exponent
is $\alpha = 1-\nu \theta_2$, where $\nu$ is Flory's swelling exponent ($nu
approx 0.588$) and $\theta_2$ is a critical exponent introduced by des
Cloizeaux ($\theta_2 \approx 0.82$), yielding $\alpha = 0.52 \pm 0.01$, in
agreement with the experimental results.


  A recent application of the peptide strings concept has yielded novel
perceptions on cell growth regulation, for instance that of oncoprotein
metastasis. Here, this interdisciplinary approach at the boundary between
physics and biology has been applied to gain a more profound insight into
rheumatoid arthritis. As a result of the present investigation, this disease
could be viewed as due to a metabolic dysregulation/syndrome-associated
breakdown in the immunoglobulin A-based surveillance of the potentially
pathogenic fungus Candida albicans that subsequently engenders a widespread
self-destruction through cross-reactive auto-epitopes, ultimately amounting to
the systemic predominance of a pro-inflammatory peptide string. Its therapeutic
counterpart equally proposed in this report might serve as a model for future
strategies against autoimmunity.


  We propose a theoretical model for species that is focused on intrinsic
organization and processes. Specifically, we regard species to be networks of
organisms integrated by reproductive mechanisms (e.g., conjugation, meiosis,
syngamy) regardless of degree or type of divergence, longevity, size, fate, or
other criteria. Ramifications of viewing units of diversity in terms of
inherent organization include the following: (1) Species, subspecies,
population, deme, Evolutionary Significant Unit (ESU), and related terms cannot
be distinguished on the basis of inherent organization and, consequently, do
not reflect a hierarchy of organization. In other words, only one model is
required to explain systems composed of organisms. (2) Temporarily isolated
networks of organisms possess the same intrinsic organization and, therefore,
the same ontological status as permanently isolated systems. (3) Units of
diversity are recognized regardless of kind or degree of divergence. (4) Fate
is not an intrinsic property of systems and is an inappropriate consideration
in species models. (5) Networks are recognized at the moment they become
isolated; therefore, resultant classifications reflect the causal events that
generate diversity (e.g., vicariance) rather than post-vicariance events such
as character evolution and, consequently, have greater historical relevance
(e.g., biogeography). (6) Because the proposed model is based on organization,
it is similar to models in other disciplines (e.g., astronomy, chemistry, and
physics) and should lead to greater unification of thought within biology and
among scientific disciplines. (7) Lineage concepts are problematic because they
do not describe a unique level of biological organization and because
ancestor-descendant relationships are not unique to living systems.


  We study the role of microtubule movement in bidirectional organelle
transport in Drosophila S2 cells and show that EGFP-tagged peroxisomes in cells
serve as sensitive probes of motor induced, noisy cytoskeletal motions.
Multiple peroxisomes move in unison over large time windows and show
correlations with microtubule tip positions, indicating rapid microtubule
fluctuations in the longitudinal direction. We report the first high-resolution
measurement of longitudinal microtubule fluctuations performed by tracing such
pairs of co-moving peroxisomes. The resulting picture shows that
motor-dependent longitudinal microtubule oscillations contribute significantly
to cargo movement along microtubules. Thus, contrary to the conventional view,
organelle transport cannot be described solely in terms of cargo movement along
stationary microtubule tracks, but instead includes a strong contribution from
the movement of the tracks.


  This article is interested in the origin of the genetic code, it puts forward
a scenario of a simultaneous selection of the bases and amino acids and setting
up of a correlation between them. Each amino acid is associated with a pair of
its own kind, called the binding pair and each binding pair is associated with
the codon(s) corresponding to the same amino acid. An explanation is also
proposed about the origin of the start and stop codons.


  We discuss an amazing prey--predator model with variable coefficients,
analyze its predictions and the accuracy of the variational iteration method
used to solve the nonlinear equations.


  We develop equilibrium and kinetic theories that describe the assembly of
viral capsid proteins on a charged central core, as seen in recent experiments
in which brome mosaic virus (BMV) capsids assemble around nanoparticles
functionalized with polyelectrolyte. We model interactions between capsid
proteins and nanoparticle surfaces as the interaction of polyelectrolyte
brushes with opposite charge, using the nonlinear Poisson Boltzmann equation.
The models predict that there is a threshold density of functionalized charge,
above which capsids efficiently assemble around nanoparticles, and that light
scatter intensity increases rapidly at early times, without the lag phase
characteristic of empty capsid assembly. These predictions are consistent with,
and enable interpretation of, preliminary experimental data. However, the
models predict a stronger dependence of nanoparticle incorporation efficiency
on functionalized charge density than measured in experiments, and do not
completely capture a logarithmic growth phase seen in experimental light
scatter. These discrepancies may suggest the presence of metastable disordered
states in the experimental system. In addition to discussing future experiments
for nanoparticle-capsid systems, we discuss broader implications for
understanding assembly around charged cores such as nucleic acids.


  Models of protein energetics which neglect interactions between amino acids
that are not adjacent in the native state, such as the Go model, encode or
underlie many influential ideas on protein folding. Implicit in this
simplification is a crucial assumption that has never been critically evaluated
in a broad context: Detailed mechanisms of protein folding are not biased by
non-native contacts, typically imagined as a consequence of sequence design
and/or topology. Here we present, using computer simulations of a well-studied
lattice heteropolymer model, the first systematic test of this oft-assumed
correspondence over the statistically significant range of hundreds of
thousands of amino acid sequences, and a concomitantly diverse set of folding
pathways. Enabled by a novel means of fingerprinting folding trajectories, our
study reveals a profound insensitivity of the order in which native contacts
accumulate to the omission of non-native interactions. Contrary to conventional
thinking, this robustness does not arise from topological restrictions and does
not depend on folding rate. We find instead that the crucial factor in
discriminating among topological pathways is the heterogeneity of native
contact energies. Our results challenge conventional thinking on the
relationship between sequence design and free energy landscapes for protein
folding, and help justify the widespread use of Go-like models to scrutinize
detailed folding mechanisms of real proteins.


  Several of the next generation sequencers are limited in their sample
preparation process by the need to make an absolute measurement of the number
of template molecules in the library to be sequenced. As currently practiced,
the practical effects of this requirement compromise sequencing performance,
both by requiring large amounts of sample DNA and by requiring extra sequencing
runs to be performed. We used digital PCR to quantitate sequencing libraries,
and demonstrated its sensitivity and robustness by preparing and sequencing
libraries from subnanogram amounts of bacterial and human DNA on the 454 and
Solexa sequencing platforms. This assay allows absolute quantitation and
eliminates uncertainties associated with the construction and application of
standard curves. The digital PCR platform consumes subfemptogram amounts of the
sequencing library and gives highly accurate results, allowing the optimal DNA
concentration to be used in setting up sequencing runs without costly and
time-consuming titration techniques. This approach also reduces the input
sample requirement more than 1000-fold: from micrograms of DNA to less than a
nanogram.


  We investigate the evolution of populations of random Boolean networks under
selection for robustness of the dynamics with respect to the perturbation of
the state of a node. The fitness landscape contains a huge plateau of maximum
fitness that spans the entire network space. When selection is so strong that
it dominates over drift, the evolutionary process is accompanied by a slow
increase in the mean connectivity and a slow decrease in the mean fitness.
Populations evolved with higher mutation rates show a higher robustness under
mutations. This means that even though all the evolved populations exist close
to the plateau of maximum fitness, they end up in different regions of network
space.


  The densities of states are evaluated for very short chain molecules made up
of overlapping monomers, using a model which has previously been shown to
produce helical structure. The results of numerical calculations are presented
for tetramers and pentamers. We show that these models demonstrate behaviors
relevant to the behaviors seen in longer, helix forming chains, particularly,
"magic numbers" of the overlap parameter where the derivatives of the densities
of states change discontinuously, and a region of bimodal energy probability
distributions, reminiscent of a first order phase transition in a bulk system.


  This article presents a physical biology approach to understanding
organization and segregation of bacterial chromosomes. The author uses a
"piston" analogy for bacterial chromosomes in a cell, which leads to a phase
diagram for the organization of two athermal chains confined in a closed
geometry characterized by two length scales (length and width). When applied to
rod-shaped bacteria such as Escherichia coli, this phase diagram predicts that,
despite strong confinement, duplicated chromosomes will demix, i.e., there
exists a primordial physical driving force for chromosome segregation. The
author discusses segregation of duplicating chromosomes using the
concentric-shell model, which predicts that newly synthesized DNA will be found
in the periphery of the chromosome during replication. In contrast to
chromosomes, these results suggest that most plasmids will be randomly
distributed inside the cell because of their small sizes. An active
partitioning system is therefore required for accurate segregation of low-copy
number plasmids. Implications of these results are also sketched, e.g., on the
role of proteins, segregation mechanisms for bacteria of diverse shapes, cell
cycle of an artificial cell, and evolution.


  Distance-based approaches in phylogenetics such as Neighbor-Joining are a
fast and popular approach for building trees. These methods take pairs of
sequences from them construct a value that, in expectation, is additive under a
stochastic model of site substitution. Most models assume a distribution of
rates across sites, often based on a gamma distribution. Provided the (shape)
parameter of this distribution is known, the method can correctly reconstruct
the tree. However, if the shape parameter is not known then we show that
topologically different trees, with different shape parameters and associated
positive branch lengths, can lead to exactly matching distributions on pairwise
site patterns between all pairs of taxa. Thus, one could not distinguish
between the two trees using pairs of sequences without some prior knowledge of
the shape parameter. More surprisingly, this can happen for {\em any} choice of
distinct shape parameters on the two trees, and thus the result is not peculiar
to a particular or contrived selection of the shape parameters. On a positive
note, we point out known conditions where identifiability can be restored
(namely, when the branch lengths are clocklike, or if methods such as maximum
likelihood are used).


  We present results contrasting food webs constructed using the same model
where the source of species was either evolution or immigration from a
previously evolved species pool. The overall structure of the webs are
remarkably similar, although we find some important differences which mainly
relate to the percentage of basal and top species. Food webs assembled from
evolved webs also show distinct plateaux in the number of tropic levels as the
resources available to system increase, in contrast to evolved webs. By
equating the resources available to basal species to area, we are able to
examine the species-area curve created by each process separately. They are
found to correspond to different regimes of the tri-phasic species-area curve.


  We present a panoramic view of the utility of coarse-grained (CG) models to
study folding and functions of proteins and RNA. Drawing largely on the methods
developed in our group over the last twenty years, we describe a number of key
applications ranging from folding of proteins with disulfide bonds to functions
of molecular machines. After presenting the theoretical basis that justifies
the use of CG models, we explore the biophysical basis for the emergence of a
finite number of folds from lattice models. The lattice model simulations of
approach to the folded state show that non-native interactions are relevant
only early in the folding process - a finding that rationalizes the success of
structure-based models that emphasize native interactions. Applications of
off-lattice $C_{\alpha}$ and models that explicitly consider side chains
($C_{\alpha}$-SCM) to folding of $\beta$-hairpin and effects of macromolecular
crowding are briefly discussed. Successful application of a new class of
off-lattice model, referred to as the Self-Organized Polymer (SOP), is shown by
describing the response of Green Fluorescent Protein (GFP) to mechanical force.
The utility of the SOP model is further illustrated by applications that
clarify the functions of the chaperonin GroEL and motion of the molecular motor
kinesin. We also present two distinct models for RNA, namely, the Three Site
Interaction (TIS) model and the SOP model, that probe forced unfolding and
force quench refolding of a simple hairpin and {\it Azoarcus} ribozyme. The
predictions based on the SOP model show that force-induced unfolding pathways
of the ribozyme can be dramatically changed by varying the loading rate. We
conclude with a discussion of future prospects for the use of coarse-grained
models in addressing problems of outstanding interest in biology.


  Sex is considered as an evolutionary paradox, since its evolutionary
advantage does not necessarily overcome the two fold cost of sharing half of
one's offspring's genome with another member of the population. Here we
demonstrate that sexual reproduction can be evolutionary stable even when its
Darwinian fitness is twice as low when compared to the fitness of asexual
mutants. We also show that more than two sexes are always evolutionary
unstable. Our approach generalizes the evolutionary game theory to analyze
species whose members are able to sense the sexual state of their conspecifics
and to switch sexes consequently. The widespread emergence and maintenance of
sex follows therefore from its co-evolution with even more widespread
environmental sensing abilities.


  Cyclic dominance of three species is a commonly occurring interaction
dynamics, often denoted the rock-paper-scissors (RPS) game. Such type of
interactions is known to promote species coexistence. Here, we generalize
recent results of Reichenbach et al. (e.g. Nature 448, 1046 (2007)) of a
four-state variant of RPS. We show that spiral formation takes place only
without a conservation law for the total density. Nevertheless, in general fast
diffusion can destroy species coexistence. We also generalize the four-state
model to slightly varying reaction rates. This is shown both analytically and
numerically not to change pattern formation, or the effective wave length of
the spirals, and therefore does not alter the qualitative properties of the
cross-over to extinction.


  The problem of the directionality of genome evolution is studied. Based on
the analysis of C-value paradox and the evolution of genome size we propose
that the function-coding information quantity of a genome always grows in the
course of evolution through sequence duplication, expansion of code, and gene
transfer from outside. The function-coding information quantity of a genome
consists of two parts, p-coding information quantity which encodes functional
protein and n-coding information quantity which encodes other functional
elements except amino acid sequence. The evidences on the evolutionary law
about the function-coding information quantity are listed. The needs of
function is the motive force for the expansion of coding information quantity
and the information quantity expansion is the way to make functional innovation
and extension for a species. So, the increase of coding information quantity of
a genome is a measure of the acquired new function and it determines the
directionality of genome evolution.


  Within large taxonomic assemblages, the number of species with adult body
mass M is characterized by a broad but asymmetric distribution, with the
largest mass being orders of magnitude larger than the typical mass. This
canonical shape can be explained by cladogenetic diffusion that is bounded
below by a hard limit on viable species mass and above by extinction risks that
increase weakly with mass. Here we introduce and analytically solve a
simplified cladogenetic diffusion model. When appropriately parameterized, the
diffusion-reaction equation predicts mass distributions that are in good
agreement with data on 4002 terrestrial mammal from the late Quaternary and
8617 extant bird species. Under this model, we show that a specific tradeoff
between the strength of within-lineage drift toward larger masses (Cope's rule)
and the increased risk of extinction from increased mass is necessary to
produce realistic mass distributions for both taxa. We then make several
predictions about the evolution of avian species masses.


  Tuffley and Steel (1997) proved that Maximum Likelihood and Maximum Parsimony
methods in phylogenetics are equivalent for sequences of characters under a
simple symmetric model of substitution with no common mechanism. This result
has been widely cited ever since. We show that small changes to the model
assumptions suffice to make the two methods inequivalent. In particular, we
analyze the case of bounded substitution probabilities as well as the molecular
clock assumption. We show that in these cases, even under no common mechanism,
Maximum Parsimony and Maximum Likelihood might make conflicting choices. We
also show that if there is an upper bound on the substitution probabilities
which is `sufficiently small', every Maximum Likelihood tree is also a Maximum
Parsimony tree (but not vice versa).


  New models for evolutionary processes of mutation accumulation allow
hypotheses about the age-specificity of mutational effects to be translated
into predictions of heterogeneous population hazard functions. We apply these
models to questions in the biodemography of longevity, including proposed
explanations of Gompertz hazards and mortality plateaus, and use them to
explore the possibility of melding evolutionary and functional models of aging.


  Demographic noise causes unlimited population growth in a broad class of
models which, without noise, would predict a stable finite population. We study
this effect on the example of a stochastic birth-death model which includes
immigration, binary reproduction and death. The unlimited population growth
proceeds as an exponentially slow decay of a metastable probability
distribution (MPD) of the population. We develop a systematic WKB theory,
complemented by the van Kampen system size expansion, for the MPD and for the
decay time. Important signatures of the MPD is a power-law tail (such that all
the distribution moments, except the zeroth one, diverge) and the presence in
the solution of two different WKB modes.


  We develop a stochastic, individual-based model for food web simulation which
in the large-population limit reduces to the well-studied Webworld model, which
has been used to successfully construct model food webs with several realistic
features. We demonstrate that an almost exact match is found between the
population dynamics in fixed food webs, and that the demographic fluctuations
have systematic effects when the new model is used to construct food webs due
to the presence of species with small populations.


  The inhabitants of the bush lupine in coastal California have been the
subject of scientific scrutiny in recent years. Observations of a
host-parasitoid interaction in the shrub's foliage, in which victims are
significantly less motile than their exploiters, record stable spatial patterns
in a fairly homogeneous environment. Though such pattern formation has been
found in reaction-diffusion models, the correspondence of these models to
continuous-time predator-prey interactions does not reflect the reality of the
system being studied. Near the root of the lupine, another host-parasitoid
interaction is also of considerable interest. In some cases this interaction,
which promotes the health of the lupine, has been observed to be much more
persistent than suggested by mathematical models.
  In this work a discrete-time spatial model of the first host-parasitoid
system is introduced. We analyze the model, describing its transient behavior
and finding the conditions under which spatial patterns occur, as well as an
estimate of outbreak size under those conditions. We consider the feasibility
of the necessary conditions in the natural system by modeling the mechanisms
responsible for them, and discuss the effects of variable habitat on pattern
formation. We also explore one possible explanation for the persistence of the
second host-parasitoid system -- the existence of an alternate host. Under
certain surprising conditions and by means quite different from previous models
of similar situations, an alternate host can greatly enhance persistence of the
nematode parasitoid.


  We present a quantitative model for the biological evolution of species body
masses within large groups of related species, e.g., terrestrial mammals, in
which body mass M evolves according to branching (speciation), multiplicative
diffusion, and an extinction probability that increases logarithmically with
mass. We describe this evolution in terms of a convection-diffusion-reaction
equation for ln M. The steady-state behavior is in good agreement with
empirical data on recent terrestrial mammals, and the time-dependent behavior
also agrees with data on extinct mammal species between 95 - 50 Myr ago.


  We propose a universal approach for analysis and fast simulations of stiff
stochastic biochemical kinetics networks, which rests on elimination of fast
chemical species without a loss of information about mesoscopic, non-Poissonian
fluctuations of the slow ones. Our approach, which is similar to the
Born-Oppenheimer approximation in quantum mechanics, follows from the
stochastic path integral representation of the full counting statistics of
reaction events (also known as the cumulant generating function). In
applications with a small number of chemical reactions, this approach produces
analytical expressions for moments of chemical fluxes between slow variables.
This allows for a low-dimensional, interpretable representation of the
biochemical system, that can be used for coarse-grained numerical simulation
schemes with a small computational complexity and yet high accuracy. As an
example, we consider a chain of biochemical reactions, derive its
coarse-grained description, and show that the Gillespie simulations of the
original stiff system, the coarse-grained simulations, and the full analytical
treatment are in an agreement, but the coarse-grained simulations are three
orders of magnitude faster than the Gillespie analogue.


  We consider the application of Kramers theory to the microscopic calculation
of rates of conformational transitions of macromolecules. The main difficulty
in such an approach is to locate the transition state in a huge configuration
space. We present a method which identifies the transition state along the most
probable reaction pathway. It is then possible to microscopically compute the
activation energy, the damping coefficient, the eigenfrequencies at the
transition state and obtain the rate, without any a-priori choice of a reaction
coordinate. Our theoretical results are tested against the results of Molecular
Dynamics simulations for transitions in a 2-dimensional double well and for the
cis-trans isomerization of a linear molecule.


  The two classic theories for the existence of sexual replication are that sex
purges deleterious mutations from a population, and that sex allows a
population to adapt more rapidly to changing environments. These two theories
have often been presented as opposing explanations for the existence of sex.
Here, we develop and analyze evolutionary models based on the asexual and
sexual replication pathways in Saccharomyces cerevisiae (Baker's yeast), and
show that sexual replication can both purge deleterious mutations in a static
environment, as well as lead to faster adaptation in a dynamic environment.
This implies that sex can serve a dual role, which is in sharp contrast to
previous theories.


  We prove that Nakhleh's latest dissimilarity measure for phylogenetic
networks is a metric on the classes of tree-child phylogenetic networks, of
semi-binary time consistent tree-sibling phylogenetic networks, and of
multi-labeled phylogenetic trees. We also prove that it distinguishes
phylogenetic networks with different reduced versions. In this way, it becomes
the dissimilarity measure for phylogenetic networks with the strongest
separation power available so far.


  Market-based instruments such as payments, auctions or tradable permits have
been proposed as flexible and cost-effective instruments for biodiversity
conservation on private lands. Trading the service of conservation requires one
to define a metric that determines the extent to which a conserved site adds to
the regional conservation objective. Yet, while markets for conservation are
widely discussed and increasingly applied, little research has been conducted
on explicitly accounting for spatial ecological processes in the trading. In
this paper, we use a coupled ecological economic simulation model to examine
how spatial connectivity may be considered in the financial incentives created
by a market-based conservation scheme. Land use decisions, driven by changing
conservation costs and the conservation market, are simulated by an agent-based
model of land users. On top of that, a metapopulation model evaluates the
conservational success of the market. We find that optimal spatial incentives
for agents correlate with species characteristics such as the dispersal
distance, but they also depend on the spatio-temporal distribution of
conservation costs. We conclude that a combined analysis of ecological and
socio-economic conditions should be applied when designing market instruments
to protect biodiversity.


  Several aspects of mathematical astrobiology are discussed. It is argued that
around the time of the origin of life the handedness of biomolecules must have
established itself through an instability. Possible pathways of producing a
certain handedness include mechanisms involving either autocatalysis or,
alternatively, epimerization as governing effects. Concepts for establishing
hereditary information are discussed in terms of the theory of hypercycles.
Instabilities toward parasites and possible remedies by invoking spatial extent
are reviewed. Finally, some effects of early life are discussed that
contributed to modifying and regulating atmosphere and climate of the Earth,
and that could have contributed to the highly oxidized state of its crust.


  An agent-based computer simulation of death by inheritable mutations in a
changing environment shows a maximal population, or avoids extinction, at so
intermediate mutation rate of the individuals. Thus death seems needed to al
for evolution of the fittest, as required by a changing environment.


  Arguably the most influential force in human history is the formation of
social coalitions and alliances (i.e., long-lasting coalitions) and their
impact on individual power. In most great ape species, coalitions occur at
individual and group levels and among both kin and non-kin. Nonetheless, ape
societies remain essentially hierarchical, and coalitions rarely weaken social
inequality. In contrast, human hunter-gatherers show a remarkable tendency to
egalitarianism, and human coalitions and alliances occur not only among
individuals and groups, but also among groups of groups. Here, we develop a
stochastic model describing the emergence of networks of allies resulting from
within-group competition for status or mates between individuals utilizing
dyadic information. The model shows that alliances often emerge in a phase
transition-like fashion if the group size, awareness, aggressiveness, and
persuasiveness of individuals are large and the decay rate of individual
affinities is small. With cultural inheritance of social networks, a single
leveling alliance including all group members can emerge in several
generations. Our results suggest that a rapid transition from a hierarchical
society of great apes to an egalitarian society of hunter-gatherers (often
referred to as "egalitarian revolution") could indeed follow an increase in
human cognitive abilities. The establishment of stable group-wide egalitarian
alliances creates conditions promoting the origin of cultural norms favoring
the group interests over those of individuals.


  Mimicry is a resemblance between species that benefits at least one of the
species. It is a ubiquitous evolutionary phenomenon particularly common among
prey species, in which case the advantage involves better protection from
predation. We formulate a mathematical description of mimicry among prey
species, to investigate benefits and disadvantages of mimicry. The basic setup
involves differential equations for quantities representing predator behavior,
namely, the probabilities for attacking prey at the next encounter. Using this
framework, we present new quantitative results, and also provide a unified
description of a significant fraction of the quantitative mimicry literature.
The new results include `temporary' mutualism between prey species, and an
optimal density at which the survival benefit is greatest for the mimic. The
formalism leads naturally to extensions in several directions, such as the
evolution of mimicry, the interplay of mimicry with population dynamics, etc.
We demonstrate this extensibility by presenting some explorations on
spatiotemporal pattern dynamics.


  The $\sigma$ profile is presented as a tool to analyze the organization of
systems at different scales, and how this organization changes in time.
Describing structures at different scales as goal-oriented agents, one can
define $\sigma \in [0,1]$ ("satisfaction") as the degree to which the goals of
each agent at each scale have been met. $\sigma$ reflects the organization
degree at that scale. The $\sigma$ profile of a system shows the satisfaction
at different scales, with the possibility to study their dependencies and
evolution. It can also be used to extend game theoretic models. A general
tendency on the evolution of complexity and cooperation naturally follows from
the $\sigma$ profile. Experiments on a virtual ecosystem are used as
illustration.