  Protein structures are a very special class among all possible structures. It
was suggested that a ``designability principle'' plays a crucial role in
nature's selection of protein sequences and structures. Here we provide a
theoretical base for such a selection principle, using a novel formulation of
the protein folding problem based on hydrophobic interactions. A structure is
reduced to a string of 0's and 1's which represent the surface and core sites,
respectively, as the backbone is traced. Each structure is therefore associated
with one point in a high dimensional space. Sequences are represented by
strings of their hydrophobicities and thus can be mapped into the same space. A
sequence which lies closer to a particular structure in this space than to any
other structures will have that structure as its ground state. Atypical
structures, namely those far away from other structures in the high dimensional
space, have more sequences which fold into them, and are thermodynamically more
stable. We argue that the most common folds of proteins are the most atypical
in the space of possible structures.


  A new method for sequence optimization in protein models is presented. The
approach, which has inherited its basic philosophy from recent work by Deutsch
and Kurosky [Phys. Rev. Lett. 76, 323 (1996)] by maximizing conditional
probabilities rather than minimizing energy functions, is based upon a novel
and very efficient multisequence Monte Carlo scheme. By construction, the
method ensures that the designed sequences represent good folders
thermodynamically. A bootstrap procedure for the sequence space search is
devised making very large chains feasible. The algorithm is successfully
explored on the two-dimensional HP model with chain lengths N=16, 18 and 32.


  A general strategy is described for finding which amino acid sequences have
native states in a desired conformation (inverse design). The approach is used
to design sequences of 48 hydrophobic and polar aminoacids on three-dimensional
lattice structures. Previous studies employing a sequence-space Monte-Carlo
technique resulted in the successful design of one sequence in ten attempts.
The present work also entails the exploration of conformations that compete
significantly with the target structure for being its ground state. The design
procedure is successful in all the ten cases.


  DNA stretching experiments are usually interpreted using the worm-like chain
model; the persistence length A appearing in the model is then interpreted as
the elastic stiffness of the double helix. In fact the persistence length
obtained by this method is a combination of bend stiffness and intrinsic bend
effects reflecting sequence information, just as at zero stretching force. This
observation resolves the discrepancy between the value of A measured in these
experiments and the larger ``dynamic persistence length'' measured by other
means. On the other hand, the twist persistence length deduced from
torsionally-constrained stretching experiments suffers no such correction. Our
calculation is very simple and analytic; it applies to DNA and other polymers
with weak intrinsic disorder.


  These lectures contain a brief description of evolutionary models inspired by
the statistical mechanics of disordered systems. After an introduction
describing the Darwinian paradigm of evolving populations, the deterministic
quasispecies equation is described, and the simplest fitness landscapes are
discussed. The effect of finite population size is then considered, from the
opposing points of view leading to stochastic escape and to adaptive walks. A
synthesis is attempted. Finally the effects of coevolution are considered, and
the promising models of large-scale inspired by the Bak-Sneppen models are
described.


  A two amino acid (hydrophobic and polar) scheme is used to perform the design
on target conformations corresponding to the native states of twenty single
chain proteins. Strikingly, the percentage of successful identification of the
nature of the residues benchmarked against naturally occurring proteins and
their homologues is around 75 % independent of the complexity of the design
procedure. Typically, the lowest success rate occurs for residues such as
alanine that have a high secondary structure functionality. Using a simple
lattice model, we argue that one possible shortcoming of the model studied may
involve the coarse-graining of the twenty kinds of amino acids into just two
effective types.


  A simple lattice model for proteins that allows for distinct sizes of the
amino acids is presented. The model is found to lead to a significant number of
conformations that are the unique ground state of one or more sequences or
encodable. Furthermore, several of the encodable structures are highly
designable and are the non-degenerate ground state of several sequences. Even
though the native state conformations are typically compact, not all compact
conformations are encodable. The incorporation of the hydrophobic and polar
nature of amino acids further enhances the attractive features of the model.


  We present and discuss a novel approach to the direct and inverse protein
folding problem. The proposed strategy is based on a variational approach that
allows the simultaneous extraction of amino acid interactions and the
low-temperature free energy of sequences of amino acids. The knowledge-based
technique is simple and straightforward to implement even for realistic
off-lattice proteins because it does not entail threading-like procedures. Its
validity is assessed in the context of a lattice model by means of a variety of
stringent checks.


  We consider two- and three-dimensional lattice models of proteins which were
characterized previously. We coarse grain their folding dynamics by reducing it
to transitions between effective states. We consider two methods of selection
of the effective states. The first method is based on the steepest descent
mapping of states to underlying local energy minima and the other involves an
additional projection to maximally compact conformations. Both methods generate
connectivity patterns that allow to distinguish between the good and bad
folders. Connectivity graphs corresponding to the folding funnel have few loops
and are thus tree-like. The Arrhenius law for the median folding time of a
16-monomer sequence is established and the corresponding barrier is related to
easily identifiable kinetic trap states.


  We study folding in 16-monomer heteropolymers on the square lattice. For a
given sequence, thermodynamic properties and stability of the native state are
unique. However, the kinetics of folding depends on the model of dynamics
adopted for the time evolution of the system. We consider three such models:
Rouse-like dynamics with either single monomer moves or with single and double
monomer moves, and the 'slithering snake' dynamics. Usually, the snake dynamics
has poorer folding properties compared to the Rouse-like dynamics, but examples
of opposite behavior can also be found. This behavior relates to which
conformations act as local energy minima when their stability is checked
against the moves of a particular dynamics. A characteristic temperature
related to the combined probability, $P_L$, to stay in the non-native minima
during folding coincides with the temperature of the fastest folding. Studies
of $P_L$ yield an easy numerical way to determine conditions of the optimal
folding.


  A model for the evolution of a finite population in a rugged fitness
landscape is introduced and solved. The population is trapped in an
evolutionary loop, alternating periods of stasis to periods in which it
performs adaptive walks. The dependence of the average rarity of the population
(a quantity related to the fitness of the most adapted individual) and of the
duration of stases on population size and mutation rate is calculated.


  Various bacterial strains exhibit colonial branching patterns during growth
on poor substrates. These patterns reflect bacterial cooperative
self-organization and cybernetic processes of communication, regulation and
control employed during colonial development. One method of modeling is the
continuous, or coupled reaction-diffusion approach, in which continuous time
evolution equations describe the bacterial density and the concentration of the
relevant chemical fields. In the context of branching growth, this idea has
been pursued by a number of groups. We present an additional model which
includes a lubrication fluid excreted by the bacteria. We also add fields of
chemotactic agents to the other models. We then present a critique of this
whole enterprise with focus on the models' potential for revealing new
biological features.


  Cooperativity plays an important role in the action of proteins bound to DNA.
A simple, mechanical mechanism for cooperativity, in the form of a
tension-mediated interaction between proteins bound to DNA at two different
locations is proposed. These proteins are not in direct physical contact. DNA
segments intercalating bound proteins are modeled as a Worm-Like Chain, which
is free to deform in two dimensions. The tension-controlled protein-protein
interaction is the consequence of two effects produced by the protein binding.
The first is the introduction of a bend in the host DNA and the second is the
modification of the bending modulus of the DNA in the immediate vicinity of the
bound protein. The interaction between two bound proteins may be either
attractive or repulsive, depending on their relative orientation on the DNA.
Applied tension controls both the strength and the range of protein-protein
interactions in this model. Properties of the cooperative interaction are
discussed, along with experimental implications.


  An evolving population, in which individual members (`agents') adapt their
behaviour according to past experience, is of central importance to many
disciplines. Because of their limited knowledge and capabilities, agents are
forced to make decisions based on inductive, rather than deductive, thinking.
We show that a population of competing agents with similar capabilities and
knowledge will tend to self-segregate into opposing groups characterized by
extreme behavior. Cautious agents perform poorly and tend to become rare.


  Segregation of populations is a key question in evolution theory. One
important aspect is the relation between spatial organization and the
population's composition. Here we study a specific example -- sectors in
expanding bacterial colonies. Such sectors are spatially segregated
sub-populations of mutants. The sectors can be seen both in disk-shaped
colonies and in branching colonies. We study the sectors using two models we
have used in the past to study bacterial colonies -- a continuous
reaction-diffusion model with non-linear diffusion and a discrete
``Communicating Walkers'' model. We find that in expanding colonies, and
especially in branching colonies, segregation processes are more likely than in
a spatially static population. One such process is the establishment of stable
sub- population having neutral mutation. Another example is the maintenance of
wild-type population along side with sub-population of advantageous mutants.
Understanding such processes in bacterial colonies is an important subject by
itself, as well as a model system for similar processes in other spreading
populations.


  Novel numerical techniques, validated by an analysis of barnase and
chymotrypsin inhibitor, are used to elucidate the paramount role played by the
geometry of the protein backbone in steering the folding to the correct native
state. It is found that, irrespective of the sequence, the native state of a
protein has exceedingly large number of conformations with a given amount of
structural overlap compared to other compact artificial backbones; moreover the
conformational entropies of unrelated proteins of the same length are nearly
equal at any given stage of folding. These results are suggestive of an
extremality principle underlying protein evolution, which, in turn, is shown to
be associated with the emergence of secondary structures.


  We study the role of phylogenetic trees on correlations in mutation
processes. Generally, correlations decay exponentially with the generation
number. We find that two distinct regimes of behavior exist. For mutation rates
smaller than a critical rate, the underlying tree morphology is almost
irrelevant, while mutation rates higher than this critical rate lead to strong
tree-dependent correlations. We show analytically that identical critical
behavior underlies all multiple point correlations. This behavior generally
characterizes branching processes undergoing mutation.


  We study the dynamics of an age-structured population in which the life
expectancy of an offspring may be mutated with respect to that of its parent.
When advantageous mutation is favored, the average fitness of the population
grows linearly with time $t$, while in the opposite case the average fitness is
constant. For no mutational bias, the average fitness grows as t^{2/3}. The
average age of the population remains finite in all cases and paradoxically is
a decreasing function of the overall population fitness.


  The genetic algorithm is an optimization procedure motivated by biological
evolution and is successfully applied to optimization problems in different
areas. A statistical mechanics model for its dynamics is proposed based on the
parent-child fitness correlation of the genetic operators, making it applicable
to general fitness landscapes. It is compared to a recent model based on a
maximum entropy ansatz. Finally it is applied to modeling the dynamics of a
genetic algorithm on the rugged fitness landscape of the NK model.


  We introduce a dynamical annealing schedule for population-based optimization
algorithms with mutation. On the basis of a statistical mechanics formulation
of the population dynamics, the mutation rate adapts to a value maximizing
expected rewards at each time step. Thereby, the mutation rate is eliminated as
a free parameter from the algorithm.


  Scale-free dynamics in physical and biological systems can arise from a
variety of causes. Here, we explore a branching process which leads to such
dynamics. We find conditions for the appearance of power laws and study
quantitatively what happens to these power laws when such conditions are
violated. From a branching process model, we predict the behavior of two
systems which seem to exhibit near scale-free behavior--rank-frequency
distributions of number of subtaxa in biology, and abundance distributions of
genotypes in an artificial life system. In the light of these, we discuss
distributions of avalanche sizes in the Bak-Tang-Wiesenfeld sandpile model.


  The statistical mechanics of heteropolymer structure formation is studied in
the context of RNA secondary structures. A designed RNA sequence biased
energetically towards a particular native structure (a hairpin) is used to
study the transition between the native and molten phase of the RNA as a
function of temperature. The transition is driven by a competition between the
energy gained from the polymer's overlap with the native structure and the
entropic gain of forming random contacts. A simplified Go-like model is
proposed and solved exactly. The predicted critical behavior is verified via
exact numerical enumeration of a large ensemble of similarly designed
sequences.


  We examine the conformations of a model for a short segment of closed DNA.
The molecule is represented as a cylindrically symmetric elastic rod with a
constraint corresponding to a specification of the linking number. We obtain
analytic expressions leading to the spatial configuration of a family of
solutions representing distortions that interpolate between the circular form
of DNA and a figure-eight form that represents the onset of interwinding. We
are also able to generate knotted loops. We suggest ways to use our approach to
produce other configurations relevant to studies of DNA structure. The
stability of the distorted configurations is assessed, along with the effects
of fluctuations on the free energy of the various configurations.


  We study the elastic behaviour of a supercoiled DNA molecule. The simplest
model is that of a rod like chain, involving two elastic constants, the bending
and the twist rigidities. We show that this model is singular and needs a small
distance cut-off, which is a natural length scale giving the limit of validity
of the model, of the order of the double helix pitch. The rod like chain in
presence of the cutoff is able to reproduce quantitatively the experimentally
observed effects of supercoiling on the elongation-force characteristics, in
the small supercoiling regime. An exact solution of the model, using both
transfer matrix techniques and its mapping to a quantum mechanics problem,
allows to extract, from the experimental data,the value of the twist rigidity.
We also analyse the variation of the torque and the writhe to twist ratio
versus supercoiling, showing analytically the existence of a rather sharp
crossover regime which can be related to the excitation of plectonemic-like
structures. Finally we study the extension fluctuations of a stretched and
supercoiled DNA molecule, both at fixed torque and at fixed supercoiling angle,
and we compare the theoretical predictions to some preliminary experimental
data.


  Dynamical connectivity graphs, which describe dynamical transition rates
between local energy minima of a system, can be displayed against the
background of a disconnectivity graph which represents the energy landscape of
the system. The resulting supergraph describes both dynamics and statics of the
system in a unified coarse-grained sense. We give examples of the supergraphs
for several two dimensional spin and protein-related systems. We demonstrate
that disordered ferromagnets have supergraphs akin to those of model proteins
whereas spin glasses behave like random sequences of aminoacids which fold
badly.


  When two binary solutions are separated by a permeable barrier, the
individual species typically diffuse and mix, dissipating their chemical
potential gradients. However, we use model lattice simulations to show that
single-file molecular-sized channels (such biomembrane channels and zeolites)
can exhibit diffusional pumping, where one type of particle uses its entropy of
mixing to drive another up its chemical potential gradient. Quantitative
analyses of rates and efficiencies of transport are plotted as functions of
transmembrane potential, pore length, and particle-pore interactions. Our
results qualitatively explain recent measurements of ``negative'' osmosis and
suggest new, more systematic experiments, particularly in zeolite transport
systems.


  We show that evolution in a population of adaptive agents, repeatedly
competing for a limited resource, can come to an abrupt halt. This transition
from evolutionary to non-evolutionary behavior arises as the global resource
level is changed, and is reminiscent of a phase transition to a frozen state.
Its origin lies in the inductive decision-making of the agents, the limited
global information that they possess and the dynamical feedback inherent in the
system.


  We study the dynamics of the Bak-Sneppen model on small-world networks. For
each site in the network, we define a ``connectance,'' which measures the
distance to all other sites. We find radically different patterns of activity
for different sites, depending on their connectance and also on the topology of
the network. For a given network, the site with the minimal connectance shows
long periods of stasis interrupted by much smaller periods of activity. In
contrast, the activity pattern for the maximally connected site appears uniform
on the same time scale. We discuss the significance of these results for
speciation events.


  We present results relevant to the equilibrium organization of DNA strands of
arbitrary length interacting with a spherical organizing center, suggestive of
DNA-histone complexation in nucleosomes. We obtain a rich phase diagram in
which a wrapping state is transformed into a complex multi-leafed, rosette
structure as the adhesion energy is reduced. The statistical mechanics of the
"melting" of a rosette can be mapped into an exactly soluble one-dimensional
many-body problem.


  It is well known that transcription can induce torsional stress in DNA,
affecting the activity of nearby genes or even inducing structural transitions
in the DNA duplex. It has long been assumed that the generation of significant
torsional stress requires the DNA to be anchored, forming a limited topological
domain, since otherwise it would spin almost freely about its axis. Previous
estimates of the rotational drag have, however, neglected the role of small
natural bends in the helix backbone. We show how these bends can increase the
drag several thousandfold relative to prior estimates, allowing significant
torsional stress even in linear, unanchored DNA. The model helps explain
several puzzling experimental results on structural transitions induced by
transcription of DNA.


  We numerically study a disordered model for the RNA secondary structure and
we find that it undergoes a phase transition, with a breaking of the replica
symmetry in the low temperature region (like in spin glasses). Our results are
based on the exact evaluation of the partition function.


  Scaling of folding properties of proteins is studied in a toy system -- the
lattice Go model with various two- and three- dimensional geometries of the
maximally compact native states. Characteristic folding times grow as power
laws with the system size. The corresponding exponents are not universal.
Scaling of the thermodynamic stability also indicates size-related
deterioration of the folding properties.


  Motivated by present activities in (statistical) physics directed towards
biological evolution, we review the interplay of three evolutionary forces:
mutation, selection, and genetic drift. The review addresses itself to
physicists and intends to bridge the gap between the biological and the
physical literature. We first clarify the terminology and recapitulate the
basic models of population genetics, which describe the evolution of the
composition of a population under the joint action of the various evolutionary
forces. Building on these foundations, we specify the ingredients explicitly,
namely, the various mutation models and fitness landscapes. We then review
recent developments concerning models of mutational degradation. These predict
upper limits for the mutation rate above which mutation can no longer be
controlled by selection, the most important phenomena being error thresholds,
Muller's ratchet, and mutational meltdowns. Error thresholds are deterministic
phenomena, whereas Muller's ratchet requires the stochastic component brought
about by finite population size. Mutational meltdowns additionally rely on an
explicit model of population dynamics, and describe the extinction of
populations. Special emphasis is put on the mutual relationship between these
phenomena. Finally, a few connections with the process of molecular evolution
are established.


  A coarse grained model of a random polypeptide chain, with only discrete
torsional degrees of freedom and Hookean springs connecting pairs of
hydrophobic residues is shown to display stretched exponential relaxation under
Metropolis dynamics at low temperatures with the exponent $\beta\simeq 1/4$, in
agreement with the best experimental results. The time dependent correlation
functions for fluctuations about the native state, computed in the Gaussian
approximation for real proteins, have also been found to have the same
functional form. Our results indicate that the energy landscape exhibits
universal features over a very large range of energies and is relatively
independent of the specific dynamics.


  Intracellular signaling often arises from ligand-induced oligomerization of
cell surface receptors. This oligomerization or clustering process is
fundamentally a cooperative behavior between near-neighbor receptor molecules;
the properties of this cooperative process clearly affects the signal
transduction. Recent investigations have revealed the molecular basis of
receptor-receptor interactions, but a simple theoretical framework for using
this data to predict cluster formation has been lacking. Here, we propose a
simple, coarse-grained, phenomenological model for ligand-modulated receptor
interactions and discuss its equilibrium properties via mean-field theory. The
existence of a first-order transition for this model has immediate implications
regarding the robustness of the cellular signaling response.


  Understanding the mechanism of protein secondary structure formation is an
essential part of protein-folding puzzle. Here we describe a simple model for
the formation of the $\beta$-hairpin, motivated by the fact that folding of a
$\beta$-hairpin captures much of the basic physics of protein folding. We argue
that the coupling of ``primary'' backbone stiffness and ``secondary'' contact
formation (similar to the coupling between the ``secondary'' and ``tertiary''
structure in globular proteins), caused for example by side-chain packing
regularities, is responsible for producing an all-or-none 2-state
$\beta$-hairpin formation. We also develop a recursive relation to compute the
phase diagram and single exponential folding/unfolding rate arising via a
dominant transition state.


  The extension elasticity of rod-coil mutliblock copolymers is analyzed for
two experimentally accessible situations. In the quenched case, when the
architecture is fixed by the synthesis, the force law is distinguished by a
sharp change in the slope. In the annealed case, where interconversion between
rod and coil states is possible, the resulting force law is sigmoid with a
pronounced plateau. This last case is realized, for example, when
homopolypeptides capable of undergoing a helix-coil transition are extended
from a coil state. Both scenarios are relevant to the analysis and design of
experiments involving single molecule mechanical measurements of biopolymers
and synthetic macromolecules.


  A generic model of a random polypeptide chain, with discrete torsional
degrees of freedom and Hookean springs connecting pairs of hydrophobic
residues, reproduces the energy probability distribution of real proteins over
a very large range of energies. We show that this system with harmonic
interactions, under dissipative dynamics driven by random noise, leads to a
distribution of energy states obeying a modified one-dimensional
Ornstein-Uhlenbeck process and giving rise to the so called Wigner
distribution. A tunably fine- or coarse-grained sampling of the energy
landscape yields a family of distributions for the energies and energy
spacings.


  We demonstrate a variant of the Bond Fluctuation lattice Monte Carlo model in
which moves through cis conformations are forbidden. Ring polymers in this
model have a conserved quantity that amounts to a topological linking number.
Increased linking number reduces the radius of gyration mildly. A linking
number of order 0.2 per bond leads to an eight-percent reduction of the radius
for 128-bond chains. This percentage appears to rise with increasing chain
length, contrary to expectation. For ring chains evolving without the
conservation of linking number, we demonstrate a substantial anti-correlation
between the twist and writhe variables whose sum yields the linking number. We
raise the possibility that our observed anti-correlations may have counterparts
in the most important practical polymer that conserves linking number, DNA.


  A theory is presented for the binding of small molecules such as surfactants
to semiflexible polymers. The persistence length is assumed to be large
compared to the monomer size but much smaller than the total chain length. Such
polymers (e.g. DNA) represent an intermediate case between flexible polymers
and stiff, rod-like ones, whose association with small molecules was previously
studied. The chains are not flexible enough to actively participate in the
self-assembly, yet their fluctuations induce long-range attractive interactions
between bound molecules. In cases where the binding significantly affects the
local chain stiffness, those interactions lead to a very sharp, cooperative
association. This scenario is of relevance to the association of DNA with
surfactants and compact proteins such as RecA. External tension exerted on the
chain is found to significantly modify the binding by suppressing the
fluctuation-induced interaction.


  The helix-coil transition is modified by grafting to a surface. This
modification is studied for short peptides capable of forming $\alpha$-helices.
Three factors are involved: (i) the grafting can induced change of the boundary
free energy of the helical domain (ii) the van der Waals attraction between the
helices and (iii) the crowding induced stretching of the coils. As a result the
helix-coil transition acquires ``all or nothing'' characteristics. In addition
the transition temperature is elevated and the transition itself sharpens as
the grafting density increases.


  Reaching a ground state of a spin system is analogous to a protein evolving
into its native state. We study the ``folding'' times for various random Ising
spin systems and determine characteristic temperatures that relate to the
``folding''. Under optimal kinetic conditions, the ``folding'' times scale with
the system size as a power law with a non-universal exponent. This is similar
to what happens in model proteins. On the other hand, the scaling behavior of
the characteristic temperatures is different than in model proteins. Both in
the spin systems and in proteins, the folding properties deteriorate with the
system size.


  We consider six different secondary structures of proteins and construct two
types of Go-type off-lattice models: with the steric constraints and without.
The basic aminoacid-aminoacid potential is Lennard Jones for the native
contacts and a soft repulsion for the non-native contacts. The interactions are
chosen to make the target secondary structure be the native state of the
system. We provide a thorough equilibrium and kinetic characterization of the
sequences through the molecular dynamics simulations with the Langevin noise.
Models with the steric constraints are found to be better folders and to be
more stable, especially in the case of the $\beta$-structures. Phononic spectra
for vibrations around the native states have low frequency gaps that correlate
with the thermodynamic stability. Folding of the secondary structures proceeds
through a well defined sequence of events. For instance, $\alpha$-helices fold
from the ends first. The closer to the native state, the faster establishment
of the contacts. Increasing the system size deteriorates the folding
characteristics. We study the folding times as a function of viscous friction
and find a regime of moderate friction with the linear dependence. We also
consider folding when one end of a structure is pinned which imitates
instantaneous conditions when a protein is being synthesized. We find that,
under such circumstances, folding of helices is faster and of the
$\beta$-sequences slower.


  The protein folding problem has attracted an increasing attention from
physicists. The problem has a flavor of statistical mechanics, but possesses
the most common feature of most biological problems -- the profound effects of
evolution. I will give an introduction to the problem, and then focus on some
recent work concerning the so-called ``designability principle''. The
designability of a structure is measured by the number of sequences that have
that structure as their unique ground state. Structures differ drastically in
terms of their designability; highly designable structures emerge with a number
of associated sequences much larger than the average. These highly designable
structures 1) possess ``proteinlike'' secondary structures and motifs, 2) are
thermodynamically more stable, and 3) fold faster than other structures. These
results suggest that protein structures are selected in nature because they are
readily designed and stable against mutations, and that such selection
simultaneously leads to thermodynamic stability and foldability. According to
this picture, a key to the protein folding problem is to understand the
emergence and the properties of the highly designable structures.


  In modern biology, one of the most important research problems is to
understand how protein sequences fold into their native 3D structures. To
investigate this problem at a high level, one wishes to analyze the protein
landscapes, i.e., the structures of the space of all protein sequences and
their native 3D structures. Perhaps the most basic computational problem at
this level is to take a target 3D structure as input and design a fittest
protein sequence with respect to one or more fitness functions of the target 3D
structure. We develop a toolbox of combinatorial techniques for protein
landscape analysis in the Grand Canonical model of Sun, Brem, Chan, and Dill.
The toolbox is based on linear programming, network flow, and a linear-size
representation of all minimum cuts of a network. It not only substantially
expands the network flow technique for protein sequence design in Kleinberg's
seminal work but also is applicable to a considerably broader collection of
computational problems than those considered by Kleinberg. We have used this
toolbox to obtain a number of efficient algorithms and hardness results. We
have further used the algorithms to analyze 3D structures drawn from the
Protein Data Bank and have discovered some novel relationships between such
native 3D structures and the Grand Canonical model.


  We have recently built a database system called BODHI, intended to store
plant bio-diversity information. It is based on an object-oriented modeling
approach and is developed completely around public-domain software. The unique
feature of BODHI is that it seamlessly integrates diverse types of data,
including taxonomic characteristics, spatial distributions, and genetic
sequences, thereby spanning the entire range from molecular to organism-level
information. A variety of sophisticated indexing strategies are incorporated to
efficiently access the various types of data, and a rule-based query processor
is employed for optimizing query execution. In this paper, we report on our
experiences in building BODHI and on its performance characteristics for a
representative set of queries.


  Conception, design, and implementation of cDNA microarray experiments present
a variety of bioinformatics challenges for biologists and computational
scientists. The multiple stages of data acquisition and analysis have motivated
the design of Expresso, a system for microarray experiment management. Salient
aspects of Expresso include support for clone replication and randomized
placement; automatic gridding, extraction of expression data from each spot,
and quality monitoring; flexible methods of combining data from individual
spots into information about clones and functional categories; and the use of
inductive logic programming for higher-level data analysis and mining. The
development of Expresso is occurring in parallel with several generations of
microarray experiments aimed at elucidating genomic responses to drought stress
in loblolly pine seedlings. The current experimental design incorporates 384
pine cDNAs replicated and randomly placed in two specific microarray layouts.
We describe the design of Expresso as well as results of analysis with Expresso
that suggest the importance of molecular chaperones and membrane transport
proteins in mechanisms conferring successful adaptation to long-term drought
stress.


  A new class of distances appropriate for measuring similarity relations
between sequences, say one type of similarity per distance, is studied. We
propose a new ``normalized information distance'', based on the noncomputable
notion of Kolmogorov complexity, and show that it is in this class and it
minorizes every computable distance in the class (that is, it is universal in
that it discovers all computable similarities). We demonstrate that it is a
metric and call it the {\em similarity metric}. This theory forms the
foundation for a new practical tool. To evidence generality and robustness we
give two distinctive applications in widely divergent areas using standard
compression programs like gzip and GenCompress. First, we compare whole
mitochondrial genomes and infer their evolutionary history. This results in a
first completely automatic computed whole mitochondrial phylogeny tree.
Secondly, we fully automatically compute the language tree of 52 different
languages.


  It is widely accepted that (1) the natural or folded state of proteins is a
global energy minimum, and (2) in most cases proteins fold to a unique state
determined by their amino acid sequence. The H-P (hydrophobic-hydrophilic)
model is a simple combinatorial model designed to answer qualitative questions
about the protein folding process. In this paper we consider a problem
suggested by Brian Hayes in 1998: what proteins in the two-dimensional H-P
model have unique optimal (minimum energy) foldings? In particular, we prove
that there are closed chains of monomers (amino acids) with this property for
all (even) lengths; and that there are open monomer chains with this property
for all lengths divisible by four.


  The idea that there are any large-scale trends in the evolution of biological
organisms is highly controversial. It is commonly believed, for example, that
there is a large-scale trend in evolution towards increasing complexity, but
empirical and theoretical arguments undermine this belief. Natural selection
results in organisms that are well adapted to their local environments, but it
is not clear how local adaptation can produce a global trend. In this paper, I
present a simple computational model, in which local adaptation to a randomly
changing environment results in a global trend towards increasing evolutionary
versatility. In this model, for evolutionary versatility to increase without
bound, the environment must be highly dynamic. The model also shows that
unbounded evolutionary versatility implies an accelerating evolutionary pace. I
believe that unbounded increase in evolutionary versatility is a large-scale
trend in evolution. I discuss some of the testable predictions about organismal
evolution that are suggested by the model.


  Evolvability is the capacity to evolve. This paper introduces a simple
computational model of evolvability and demonstrates that, under certain
conditions, evolvability can increase indefinitely, even when there is no
direct selection for evolvability. The model shows that increasing evolvability
implies an accelerating evolutionary pace. It is suggested that the conditions
for indefinitely increasing evolvability are satisfied in biological and
cultural evolution. We claim that increasing evolvability is a large-scale
trend in evolution. This hypothesis leads to testable predictions about
biological and cultural evolution.


  JohnnyVon is an implementation of self-replicating machines in continuous
two-dimensional space. Two types of particles drift about in a virtual liquid.
The particles are automata with discrete internal states but continuous
external relationships. Their internal states are governed by finite state
machines but their external relationships are governed by a simulated physics
that includes Brownian motion, viscosity, and spring-like attractive and
repulsive forces. The particles can be assembled into patterns that can encode
arbitrary strings of bits. We demonstrate that, if an arbitrary "seed" pattern
is put in a "soup" of separate individual particles, the pattern will replicate
by assembling the individual particles into copies of itself. We also show
that, given sufficient time, a soup of separate individual particles will
eventually spontaneously form self-replicating patterns. We discuss the
implications of JohnnyVon for research in nanotechnology, theoretical biology,
and artificial life.


  We propose a family of very efficient hierarchical indexing schemes for
ungapped, score matrix-based similarity search in large datasets of short (4-12
amino acid) protein fragments. This type of similarity search has importance in
both providing a building block to more complex algorithms and for possible use
in direct biological investigations where datasets are of the order of 60
million objects. Our scheme is based on the internal geometry of the amino acid
alphabet and performs exceptionally well, for example outputting 100 nearest
neighbours to any possible fragment of length 10 after scanning on average less
than one per cent of the entire dataset.


  This paper presents our computational methodology using Genetic Algorithms
(GA) for exploring the nature of RNA editing. These models are constructed
using several genetic editing characteristics that are gleaned from the RNA
editing system as observed in several organisms. We have expanded the
traditional Genetic Algorithm with artificial editing mechanisms as proposed by
(Rocha, 1997). The incorporation of editing mechanisms provides a means for
artificial agents with genetic descriptions to gain greater phenotypic
plasticity, which may be environmentally regulated. Our first implementations
of these ideas have shed some light into the evolutionary implications of RNA
editing. Based on these understandings, we demonstrate how to select proper RNA
editors for designing more robust GAs, and the results will show promising
applications to real-world problems. We expect that the framework proposed will
both facilitate determining the evolutionary role of RNA editing in biology,
and advance the current state of research in Genetic Algorithms.


  Within the literature on non-cooperative game theory, there have been a
number of attempts to propose logorithms which will compute Nash equilibria.
Rather than derive a new algorithm, this paper shows that the family of
algorithms known as Markov chain Monte Carlo (MCMC) can be used to calculate
Nash equilibria. MCMC is a type of Monte Carlo simulation that relies on Markov
chains to ensure its regularity conditions. MCMC has been widely used
throughout the statistics and optimization literature, where variants of this
algorithm are known as simulated annealing. This paper shows that there is
interesting connection between the trembles that underlie the functioning of
this algorithm and the type of Nash refinement known as trembling hand
perfection.


  We present a new method for clustering based on compression. The method
doesn't use subject-specific features or background knowledge, and works as
follows: First, we determine a universal similarity distance, the normalized
compression distance or NCD, computed from the lengths of compressed data files
(singly and in pairwise concatenation). Second, we apply a hierarchical
clustering method. The NCD is universal in that it is not restricted to a
specific application area, and works across application area boundaries. A
theoretical precursor, the normalized information distance, co-developed by one
of the authors, is provably optimal but uses the non-computable notion of
Kolmogorov complexity. We propose precise notions of similarity metric, normal
compressor, and show that the NCD based on a normal compressor is a similarity
metric that approximates universality. To extract a hierarchy of clusters from
the distance matrix, we determine a dendrogram (binary tree) by a new quartet
method and a fast heuristic to implement it. The method is implemented and
available as public software, and is robust under choice of different
compressors. To substantiate our claims of universality and robustness, we
report evidence of successful application in areas as diverse as genomics,
virology, languages, literature, music, handwritten digits, astronomy, and
combinations of objects from completely different domains, using statistical,
dictionary, and block sorting compressors. In genomics we presented new
evidence for major questions in Mammalian evolution, based on
whole-mitochondrial genomic analysis: the Eutherian orders and the Marsupionta
hypothesis against the Theria hypothesis.


  After reviewing unnormalized and normalized information distances based on
incomputable notions of Kolmogorov complexity, we discuss how Kolmogorov
complexity can be approximated by data compression algorithms. We argue that
optimal algorithms for data compression with side information can be
successfully used to approximate the normalized distance. Next, we discuss an
alternative information distance, which is based on relative entropy rate (also
known as Kullback-Leibler divergence), and compression-based algorithms for its
estimation. Based on available biological and linguistic data, we arrive to
unexpected conclusion that in Bioinformatics and Computational Linguistics this
alternative distance is more relevant and important than the ones based on
Kolmogorov complexity.


  Evolutionary role of the separation into two sexes from a cyberneticist's
point of view. [I translated this 1965 article from Russian "Nauka i Zhizn"
(Science and Life) in 1988. In a popular form, the article puts forward several
useful ideas not all of which even today are necessarily well known or widely
accepted. Boris Lubachevsky, bdl@bell-labs.com ]


  In the computational-mechanics structural analysis of one-dimensional
cellular automata the following automata-theoretic analogue of the
\emph{change-point problem} from time series analysis arises: \emph{Given a
string $\sigma$ and a collection $\{\mc{D}_i\}$ of finite automata, identify
the regions of $\sigma$ that belong to each $\mc{D}_i$ and, in particular, the
boundaries separating them.} We present two methods for solving this
\emph{multi-regular language filtering problem}. The first, although providing
the ideal solution, requires a stack, has a worst-case compute time that grows
quadratically in $\sigma$'s length and conditions its output at any point on
arbitrarily long windows of future input. The second method is to
algorithmically construct a transducer that approximates the first algorithm.
In contrast to the stack-based algorithm, however, the transducer requires only
a finite amount of memory, runs in linear time, and gives immediate output for
each letter read; it is, moreover, the best possible finite-state approximation
with these three features.


  In this paper we study the problem of learning phylogenies and hidden Markov
models. We call a Markov model nonsingular if all transition matrices have
determinants bounded away from 0 (and 1). We highlight the role of the
nonsingularity condition for the learning problem. Learning hidden Markov
models without the nonsingularity condition is at least as hard as learning
parity with noise, a well-known learning problem conjectured to be
computationally hard. On the other hand, we give a polynomial-time algorithm
for learning nonsingular phylogenies and hidden Markov models.


  In this paper we propose a new symbolic-numeric algorithm to find positive
equilibria of a n-dimensional dynamical system. This algorithm implies a
symbolic manipulation of ODE in order to give a local approximation of
differential equations with power-law dynamics (S-systems). A numerical
calculus is then needed to converge towards an equilibrium, giving at the same
time a S-system approximating the initial system around this equilibrium. This
algorithm is applied to a real biological example in 14 dimensions which is a
subsystem of a metabolic pathway in Arabidopsis Thaliana.


  Despite much progress over the past decade, current Single Nucleotide
Polymorphism (SNP) genotyping technologies still offer an insufficient degree
of multiplexing when required to handle user-selected sets of SNPs. In this
paper we propose a new genotyping assay architecture combining multiplexed
solution-phase single-base extension (SBE) reactions with sequencing by
hybridization (SBH) using universal DNA arrays such as all $k$-mer arrays. In
addition to PCR amplification of genomic DNA, SNP genotyping using SBE/SBH
assays involves the following steps: (1) Synthesizing primers complementing the
genomic sequence immediately preceding SNPs of interest; (2) Hybridizing these
primers with the genomic DNA; (3) Extending each primer by a single base using
polymerase enzyme and dideoxynucleotides labeled with 4 different fluorescent
dyes; and finally (4) Hybridizing extended primers to a universal DNA array and
determining the identity of the bases that extend each primer by hybridization
pattern analysis. Our contributions include a study of multiplexing algorithms
for SBE/SBH genotyping assays and preliminary experimental results showing the
achievable tradeoffs between the number of array probes and primer length on
one hand and the number of SNPs that can be assayed simultaneously on the
other. Simulation results on datasets both randomly generated and extracted
from the NCBI dbSNP database suggest that the SBE/SBH architecture provides a
flexible and cost-effective alternative to genotyping assays currently used in
the industry, enabling genotyping of up to hundreds of thousands of
user-specified SNPs per assay.


  Gene assembly in ciliates is one of the most involved DNA processings going
on in any organism. This process transforms one nucleus (the micronucleus) into
another functionally different nucleus (the macronucleus). We continue the
development of the theoretical models of gene assembly, and in particular we
demonstrate the use of the concept of the breakpoint graph, known from another
branch of DNA transformation research. More specifically: (1) we characterize
the intermediate gene patterns that can occur during the transformation of a
given micronuclear gene pattern to its macronuclear form; (2) we determine the
number of applications of the loop recombination operation (the most basic of
the three molecular operations that accomplish gene assembly) needed in this
transformation; (3) we generalize previous results (and give elegant
alternatives for some proofs) concerning characterizations of the micronuclear
gene patterns that can be assembled using a specific subset of the three
molecular operations.


  Gene assembly in ciliates is an extremely involved DNA transformation
process, which transforms a nucleus, the micronucleus, to another functionally
different nucleus, the macronucleus. In this paper we characterize which loop
recombination operations (one of the three types of molecular operations that
accomplish gene assembly) can possibly be applied in the transformation of a
given gene from its micronuclear form to its macronuclear form. We also
characterize in which order these loop recombination operations are applicable.
This is done in the abstract and more general setting of so-called legal
strings.


  The greatest challenge in maximizing the use of gene expression data is to
develop new computational tools capable of interconnecting and interpreting the
results from different organisms and experimental settings. We propose an
integrative and comprehensive approach including a super-chip containing data
from microarray experiments collected on different species subjected to hypoxic
and anoxic stress. A data mining technology called Peano count tree (P-trees)
is used to represent genomic data in multidimensions. Each microarray spot is
presented as a pixel with its corresponding red/green intensity feature bands.
Each bad is stored separately in a reorganized 8-separate (bSQ) file format.
Each bSQ is converted to a quadrant base tree structure (P-tree) from which a
superchip is represented as expression P-trees (EP-trees) and repression
P-trees (RP-trees). The use of association rule mining is proposed to derived
to meanigingfully organize signal transduction pathways taking in consideration
evolutionary considerations. We argue that the genetic constitution of an
organism (K) can be represented by the total number of genes belonging to two
groups. The group X constitutes genes (X1,Xn) and they can be represented as 1
or 0 depending on whether the gene was expressed or not. The second group of Y
genes (Y1,Yn) is expressed at different levels. These genes have a very high
repression, high expression, very repressed or highly repressed. However, many
genes of the group Y are specie specific and modulated by the products and
combinations of genes of the group X. In this paper, we introduce the dSQ and
P-tree technology; the biological implications of association rule mining using
X and Y gene groups and some advances in the integration of this information
using the BRAIN architecture.


  A common method of making a theory more understandable, is by comparing it to
another theory which has been better developed. Radical interpretation is a
theory which attempts to explain how communication has meaning. Radical
interpretation is treated as another time-dependent theory and compared to the
time dependent theory of biological evolution. The main reason for doing this
is to find the nature of the time dependence; producing analogs between the two
theories is a necessary prerequisite to this and brings up many problems. Once
the nature of the time dependence is better known it might allow the underlying
mechanism to be uncovered. Several similarities and differences are uncovered,
there appear to be more differences than similarities.


  We apply a recently developed method, multicanonical algorithm, to the
problem of tertiary structure prediction of peptides and proteins. As a simple
example to test the effectiveness of the algorithm, Met-enkephalin is studied
and the ergodicity problem, or multiple-minima problem, is shown to be overcome
by this algorithm. The lowest-energy conformation obtained agrees with that
determined by other efficient methods such as Monte Carlo simulated annealing.
The superiority of the present method to simulated annealing lies in the fact
that the relationship to the canonical ensemble remains exactly controlled.
Once the multicanonical parameters are determined, only one simulation run is
necessary to obtain the lowest-energy conformation and furthermore the results
of this one run can be used to calculate various thermodynamic quantities at
any temperature. The latter point is demonstrated by the calculation of the
average potential energy and specific heat as functions of temperature.


  We have developed a simple, phenomenological, Monte-Carlo code that predicts
the three-dimensional structure of globular proteins from the DNA sequences
that define them. We have applied this code to two small proteins, the villin
headpiece (1VII) and cole1 rop (1ROP). Our code folds both proteins to within 5
A rms of their native structures.


  We propose a mechanism for fixing the velocity of relativistic soliton based
on the breaking of the Lorentz symmetry of the sine-Gordon (SG) model. The
proposal is first elaborated for a molecular chain model, as the simple
pendulum limit of a double pendulums chain. It is then generalized to a full
class of two-dimensional field theories of the sine-Gordon type. From a
phenomenological point of view, the mechanism allows one to select the speed of
a SG soliton just by tuning elastic couplings constants and kinematical
parameters. From a fundamental, field-theoretical point of view we show that
the characterizing features of relativistic SG solitons (existence of conserved
topological charges and stability) may be still preserved even if the Lorentz
symmetry is broken and a soliton of a given speed is selected.


  In the framework of the crystal basis model of the genetic code, where each
codon is assigned to an irreducible representation of $U_{q \to 0}(sl(2) \oplus
sl(2))$, single base mutation matrices are introduced. The strength of the
mutation is assumed to depend on the "distance" between the codons. Preliminary
general predictions of the model are compared with experimental data, with a
satisfactory agreement.


  Traditional chemical kinetics may be inappropriate to describe chemical
reactions in micro-domains involving only a small number of substrate and
reactant molecules. Starting with the stochastic dynamics of the molecules, we
derive a master-diffusion equation for the joint probability density of a
mobile reactant and the number of bound substrate in a confined domain. We use
the equation to calculate the fluctuations in the number of bound substrate
molecules as a function of initial reactant distribution. A second model is
presented based on a Markov description of the binding and unbinding and on the
mean first passage time of a molecule to a small portion of the boundary. These
models can be used for the description of noise due to gating of ionic channels
by random binding and unbinding of ligands in biological sensor cells, such as
olfactory cilia, photo-receptors, hair cells in the cochlea.


  In this work we intend to report on some results obtained by an analytical
modelling of biomacromolecular structures as driven by the study of Steiner
points and Steiner trees with an Euclidean definition of distance.


  New periods can emerge from data as a byproduct of incorrect processing or
even the method applied. In one such recent instance, a new life cycle with a
62+-3 Myr period was reportedly found (about trend) in genus variations from
the Sepkoski compendium, the world most complete fossil record. The approach
that led to reporting this period was based on Fourier method of spectral
analysis. I show here that no such period is found when the original data set
is considered rigorously and processed in the Gauss-Vanicek spectral analysis.
I also demonstrate that data altering can boost spectral power up to a nearly
100 percent increase in the signal range, thus introducing artificial, "99
percent significant" periods as seen in the corresponding variance-spectra of
noise. Besides geology and paleontology, virtually all science and engineering
disciplines could benefit from the approach described here. The main general
advantages of the Gauss-Vanicek spectral analysis lay in period detection from
gapped records and in straightforward testing of statistical null hypothesis.
The main advantage of the method for physical sciences is its use as a field
descriptor for accurate simultaneous detection of eigenfrequencies and relative
dynamics. Besides analyzing incomplete records, researchers might also want to
remove less-trustworthy data from any time series before analyzing it with the
Gauss-Vanicek method. This could increase both the accuracy and reliability of
spectral analyses in general.


  Karlin and Altschul in their statistical analysis for multiple high-scoring
segments in molecular sequences introduced a distribution function which gives
the probability there are at least r distinct and consistently ordered segment
pairs all with score at least x. For long sequences this distribution can be
expressed in terms of the distribution of the length of the longest increasing
subsequence in a random permutation. Within the past few years, this last
quantity has been extensively studied in the mathematics literature. The
purpose of these notes is to summarize these new mathematical developments in a
form suitable for use in computational biology.


  For supercritical multitype branching processes in continuous time, we
investigate the evolution of types along those lineages that survive up to some
time t. We establish almost-sure convergence theorems for both time and
population averages of ancestral types (conditioned on non-extinction), and
identify the mutation process describing the type evolution along typical
lineages. An important tool is a representation of the family tree in terms of
a suitable size-biased tree with trunk. As a by-product, this representation
allows a `conceptual proof' (in the sense of Kurtz, Lyons, Pemantle, Peres
1997) of the continuous-time version of the Kesten-Stigum theorem.


  We derive theoretical upper and lower bounds on the maximum size of DNA codes
of length n with constant GC-content w and minimum Hamming distance d, both
with and without the additional constraint that the minimum Hamming distance
between any codeword and the reverse-complement of any codeword be at least d.
We also explicitly construct codes that are larger than the best
previously-published codes for many choices of the parameters n, d and w.


  This paper explains the uniqueness of positive steady state of general
Lotka-Volterra competition model of two species of animals in the same
environment.


  In this paper, we study the uniqueness of coexistence state for multiple
species of animals in the same environment.


  We consider two versions of a simple evolutionary algorithm model for protein
folding at temperature zero: the (1+1)-EA on the LeadingOnes problem. In this
schematic model, the structure of the protein, which is encoded as a bit-string
of length n, is evolved to its native conformation through a stochastic pathway
of sequential contact bindings.
  We study the asymptotic behavior of the hitting time, in the mean case
scenario, under two different mutations: the one flip which flips a unique bit
chosen uniformly at random in the bit-string, and the Bernoulli flip which
flips each bit in the bit-string independently with probability c/n. For each
algorithm we prove a law of large numbers, a central limit theorem and compare
the performance of the two models.


  The purpose of this note is to provide proofs for some facts about the NK
model of evolution proposed by Kauffman and Levin. In the case of normally
distributed fitness summands, some of these facts have been previously
conjectured and heuristics given. In particular, we provide rigorous asymptotic
estimates for the number of local fitness maxima in the case when K is
unbounded. We also examine the role of the individual fitness distribution and
find the model to be quite robust with respect to this.


  We consider the optimal spacing between marine reserves for maximising the
viability of a species occupying a reserve network. The closer the networks are
placed together, the higher the probability of colonisation of an empty reserve
by an occupied reserve, thus increasing population viability. However, the
closer the networks are placed together, the higher the probability that a
catastrophe will cause extinction of the species in both reserves, thus
decreasing population viability. Using a simple discrete-time Markov chain
model for the presence or absence of the species in each reserve we determine
the distance between the two reserves which provides the optimal trade-off
between these processes, resulting in maximum viability of the species.


  The objective of this paper is the study of the equilibrium behavior of a
population on the hierarchical group $\Omega_N$ consisting of families of
individuals undergoing critical branching random walk and in addition these
families also develop according to a critical branching process. Strong
transience of the random walk guarantees existence of an equilibrium for this
two-level branching system. In the limit $N\to\infty$ (called the hierarchical
mean field limit), the equilibrium aggregated populations in a nested sequence
of balls $B^{(N)}_\ell$ of hierarchical radius $\ell$ converge to a backward
Markov chain on $\mathbb{R_+}$. This limiting Markov chain can be explicitly
represented in terms of a cascade of subordinators which in turn makes possible
a description of the genealogy of the population.


  The tree-metric theorem provides a necessary and sufficient condition for a
dissimilarity matrix to be a tree metric, and has served as the foundation for
numerous distance-based reconstruction methods in phylogenetics. Our main
result is an extension of the tree-metric theorem to more general dissimilarity
maps. In particular, we show that a tree with n leaves is reconstructible from
the weights of the m-leaf subtrees provided that n \geq 2m-1.


  Weakly activated signaling cascades can be modeled as linear systems. The
input-to-output transfer function and the internal gain of a linear system,
provide natural measures for the propagation of the input signal down the
cascade and for the characterization of the final outcome. The most efficient
design of a cascade for generating sharp signals, is obtained by choosing all
the off rates equal, and a ``universal'' finite optimal length.


  We consider the following interacting particle system: There is a ``gas'' of
particles, each of which performs a continuous-time simple random walk on
$\mathbb{Z}^d$, with jump rate $D_A$. These particles are called $A$-particles
and move independently of each other. They are regarded as individuals who are
ignorant of a rumor or are healthy. We assume that we start the system with
$N_A(x,0-)$ $A$-particles at $x$, and that the $N_A(x,0-),x\in\mathbb{Z}^d$,
are i.i.d., mean-$\mu_A$ Poisson random variables. In addition, there are
$B$-particles which perform continuous-time simple random walks with jump rate
$D_B$. We start with a finite number of $B$-particles in the system at time 0.
$B$-particles are interpreted as individuals who have heard a certain rumor or
who are infected. The $B$-particles move independently of each other. The only
interaction is that when a $B$-particle and an $A$-particle coincide, the
latter instantaneously turns into a $B$-particle. We investigate how fast the
rumor, or infection, spreads. Specifically, if
$\widetilde{B}(t):=\{x\in\mathbb{Z}^d:$ a $B$-particle visits $x$ during
$[0,t]\}$ and $B(t)=\widetilde{B}(t)+[-1/2,1/2]^d$, then we investigate the
asymptotic behavior of $B(t)$. Our principal result states that if $D_A=D_B$
(so that the $A$- and $B$-particles perform the same random walk), then there
exist constants $0<C_i<\infty$ such that almost surely
$\mathcal{C}(C_2t)\subset B(t)\subset \mathcal{C}(C_1t)$ for all large $t$,
where $\mathcal{C}(r)=[-r,r]^d$. In a further paper we shall use the results
presented here to prove a full ``shape theorem,'' saying that $t^{-1}B(t)$
converges almost surely to a nonrandom set $B_0$, with the origin as an
interior point, so that the true growth rate for $B(t)$ is linear in $t$. If
$D_A\ne D_B$, then we can only prove the upper bound $B(t)\subset
\mathcal{C}(C_1t)$ eventually.


  We analyze a class of models for unequal crossover (UC) of sequences
containing sections with repeated units that may differ in length. In these,
the probability of an `imperfect' alignment, in which the shorter sequence has
d units without a partner in the longer one, scales like q^d as compared to
`perfect' alignments where all these copies are paired. The class is
parameterized by this penalty factor q. An effectively infinite population size
and thus deterministic dynamics is assumed. For the extreme cases q=0 and q=1,
and any initial distribution whose moments satisfy certain conditions, we prove
the convergence to one of the known fixed points, uniquely determined by the
mean copy number, in both discrete and continuous time. For the intermediate
parameter values, the existence of fixed points is shown.


  Qian, Luscombe and Gerstein [J. Molecular Biol. 313 (2001) 673--681]
introduced a model of the diversification of protein folds in a genome that we
may formulate as follows. Consider a multitype Yule process starting with one
individual in which there are no deaths and each individual gives birth to a
new individual at rate 1. When a new individual is born, it has the same type
as its parent with probability $1-r$ and is a new type, different from all
previously observed types, with probability $r$. We refer to individuals with
the same type as families and provide an approximation to the joint
distribution of family sizes when the population size reaches $N$. We also show
that if $1\ll S\ll N^{1-r}$, then the number of families of size at least $S$
is approximately $CNS^{-1/(1-r)}$, while if $N^{1-r}\ll S$ the distribution
decays more rapidly than any power.


  Phylogenetic algebraic geometry is concerned with certain complex projective
algebraic varieties derived from finite trees. Real positive points on these
varieties represent probabilistic models of evolution. For small trees, we
recover classical geometric objects, such as toric and determinantal varieties
and their secant varieties, but larger trees lead to new and largely unexplored
territory. This paper gives a self-contained introduction to this subject and
offers numerous open problems for algebraic geometers.


  Motivated as a null model for comparison with data, we study the following
model for a phylogenetic tree on $n$ extant species. The origin of the clade is
a random time in the past, whose (improper) distribution is uniform on
$(0,\infty)$. After that origin, the process of extinctions and speciations is
a continuous-time critical branching process of constant rate, conditioned on
having the prescribed number $n$ of species at the present time. We study
various mathematical properties of this model as $n \to \infty$ limits: time of
origin and of most recent common ancestor; pattern of divergence times within
lineage trees; time series of numbers of species; number of extinct species in
total, or ancestral to extant species; and "local" structure of the tree
itself. We emphasize several mathematical techniques: associating walks with
trees, a point process representation of lineage trees, and Brownian limits.


  The general Markov model of the evolution of biological sequences along a
tree leads to a parameterization of an algebraic variety. Understanding this
variety and the polynomials, called phylogenetic invariants, which vanish on
it, is a problem within the broader area of Algebraic Statistics.
  For an arbitrary trivalent tree, we determine the full ideal of invariants
for the 2-state model, establishing a conjecture of Pachter-Sturmfels. For the
$\kappa$-state model, we reduce the problem of determining a defining set of
polynomials to that of determining a defining set for a 3-leaved tree.
  Along the way, we prove several new cases of a conjecture of
Garcia-Stillman-Sturmfels on certain statistical models on star trees, and
reduce their conjecture to a family of subcases.


  When a beneficial mutation occurs in a population, the new, favored allele
may spread to the entire population. This process is known as a selective
sweep. Suppose we sample $n$ individuals at the end of a selective sweep. If we
focus on a site on the chromosome that is close to the location of the
beneficial mutation, then many of the lineages will likely be descended from
the individual that had the beneficial mutation, while others will be descended
from a different individual because of recombination between the two sites. We
introduce two approximations for the effect of a selective sweep. The first one
is simple but not very accurate: flip $n$ independent coins with probability
$p$ of heads and say that the lineages whose coins come up heads are those that
are descended from the individual with the beneficial mutation. A second
approximation, which is related to Kingman's paintbox construction, replaces
the coin flips by integer-valued random variables and leads to very accurate
results.


  When an advantageous mutation occurs in a population, the favorable allele
may spread to the entire population in a short time, an event known as a
selective sweep. As a result, when we sample $n$ individuals from a population
and trace their ancestral lines backwards in time, many lineages may coalesce
almost instantaneously at the time of a selective sweep. We show that as the
population size goes to infinity, this process converges to a coalescent
process called a coalescent with multiple collisions. A better approximation
for finite populations can be obtained using a coalescent with simultaneous
multiple collisions. We also show how these coalescent approximations can be
used to get insight into how beneficial mutations affect the behavior of
statistics that have been used to detect departures from the usual Kingman's
coalescent.


  This article is not a research paper, but a little note on the history of
combinatorics: We present here a tentative short biography of Henri Delannoy,
and a survey of his most notable works. This answers to the question raised in
the title, as these works are related to lattice paths enumeration, to the
so-called Delannoy numbers, and were the first general way to solve Ballot-like
problems. These numbers appear in probabilistic game theory, alignments of DNA
sequences, tiling problems, temporal representation models, analysis of
algorithms and combinatorial structures.


  In the evolution of a genome, the gene sequence is sometimes rearranged, for
example by transposition of two adjacent gene blocks. In biocombinatorics, one
tries to reconstruct these rearrangement incidents from the resulting
permutation. It seems that the algorithms used are too effective and find a
shorter path than the real one. For the simplified case of adjacent
transpositions, we give expressions for the expected number of inversions after
t random moves. This average can be much smaller than t, a fact that has
largely been neglected so far.


  We study two systems of tangle equations that arise when modeling the action
of the Integrase family of proteins on DNA. These two systems--direct and
inverted repeats--correspond to two different possibilities for the initial DNA
sequence. We present one new class of solutions to the tangle equations. In the
case of inverted repeats we prove that any solution not in these, or 2
previously known classes, would have to belong to one specific class. In the
case of direct repeats we prove that the three classes are are the only
solutions possible.


  Integrase proteins acting on circular double-stranded DNA often change its
topology by transforming unknotted circles into torus knots and links. Two
systems of tangle equations--corresponding to the two initial DNA
sequences--arise when modelling this transformation: direct and inverted.
  With no a priori assumptions on the constituent tangles, we utilize Dehn
surgery arguments to completely classify the tangle solutions for each of the
two systems. A key step is to combine work of our previous paper with recent
results of Kronheimer, Mrowka, Ozsvath and Szabo, and Ernst to show a certain
prime tangle must in fact be a generalized Montesinos tangle.
  These tangle solutions are divided into three classes, common to both
systems, plus a fourth class for the inverted system that contains the sole
generalized Montesinos tangle. We discuss the possible biological implications
of our classification, and of this novel solution.


  We use Dirichlet form methods to construct and analyze a reversible Markov
process, the stationary distribution of which is the Brownian continuum random
tree. This process is inspired by the subtree prune and regraft (SPR) Markov
chains that appear in phylogenetic analysis. A key technical ingredient in this
work is the use of a novel Gromov--Hausdorff type distance to metrize the space
whose elements are compact real trees equipped with a probability measure.
Also, the investigation of the Dirichlet form hinges on a new path
decomposition of the Brownian excursion.


  For a genetic locus carrying a strongly beneficial allele which has just
fixed in a large population, we study the ancestry at a linked neutral locus.
During this ``selective sweep'' the linkage between the two loci is broken up
by recombination and the ancestry at the neutral locus is modeled by a
structured coalescent in a random background. For large selection coefficients
$\alpha$ and under an appropriate scaling of the recombination rate, we derive
a sampling formula with an order of accuracy of $\mathcal{O}((\log
\alpha)^{-2})$ in probability. In particular we see that, with this order of
accuracy, in a sample of fixed size there are at most two nonsingleton families
of individuals which are identical by descent at the neutral locus from the
beginning of the sweep. This refines a formula going back to the work of
Maynard Smith and Haigh, and complements recent work of Schweinsberg and
Durrett on selective sweeps in the Moran model.


  We extend in two directions our previous results about the sampling and the
empirical measures of immortal branching Markov processes. Direct applications
to molecular biology are rigorous estimates of the mutation rates of polymerase
chain reactions from uniform samples of the population after the reaction.
First, we consider nonhomogeneous processes, which are more adapted to real
reactions. Second, recalling that the first moment estimator is analytically
known only in the infinite population limit, we provide rigorous confidence
intervals for this estimator that are valid for any finite population. Our
bounds are explicit, nonasymptotic and valid for a wide class of nonhomogeneous
branching Markov processes that we describe in detail. In the setting of
polymerase chain reactions, our results imply that enlarging the size of the
sample becomes useless for surprisingly small sizes. Establishing confidence
intervals requires precise estimates of the second moment of random samples.
The proof of these estimates is more involved than the proofs that allowed us,
in a previous paper, to deal with the first moment. On the other hand, our
method uses various, seemingly new, monotonicity properties of the harmonic
moments of sums of exchangeable random variables.


  Maximum likelihood is one of the most widely used techniques to infer
evolutionary histories. Although it is thought to be intractable, a proof of
its hardness has been lacking. Here, we give a short proof that computing the
maximum likelihood tree is NP-hard by exploiting a connection between
likelihood and parsimony observed by Tuffley and Steel.


  We introduce simple models of genetic regulatory networks and we proceed to
the mathematical analysis of their dynamics. The models are discrete time
dynamical systems generated by piecewise affine contracting mappings whose
variables represent gene expression levels. When compared to other models of
regulatory networks, these models have an additional parameter which is
identified as quantifying interaction delays. In spite of their simplicity,
their dynamics presents a rich variety of behaviours. This phenomenology is not
limited to piecewise affine model but extends to smooth nonlinear discrete time
models of regulatory networks. In a first step, our analysis concerns general
properties of networks on arbitrary graphs (characterisation of the attractor,
symbolic dynamics, Lyapunov stability, structural stability, symmetries, etc).
In a second step, focus is made on simple circuits for which the attractor and
its changes with parameters are described. In the negative circuit of 2 genes,
a thorough study is presented which concern stable (quasi-)periodic
oscillations governed by rotations on the unit circle -- with a rotation number
depending continuously and monotonically on threshold parameters. These regular
oscillations exist in negative circuits with arbitrary number of genes where
they are most likely to be observed in genetic systems with non-negligible
delay effects.


  In nature the three-dimensional structure of a protein is encoded in the
corresponding gene. In this paper we describe a new method for encoding the
three-dimensional structure of a protein into a binary sequence. The feature of
the method is the correspondence between protein-folding and ``integration''. A
protein is approximated by a folded tetrahedron sequence. And the binary code
of a protein is obtained as the ``second derivative'' of the shape of the
folded tetrahedron sequence. With this method at hand, we can extract static
structural information of a protein from its gene. And we can describe the
distribution of three-dimensional structures of proteins without any subjective
hierarchical classification.


  The dynamics of recombination in genetics leads to an interesting nonlinear
differential equation, which has a natural generalization to a measure valued
version. The latter can be solved explicitly under rather general
circumstances. It admits a closed formula for the semigroup of nonlinear
positive operators that emerges from the forward flow and is, in general,
embedded in a multi-parameter semigroup.


  A major task of evolutionary biology is the reconstruction of phylogenetic
trees from molecular data. The evolutionary model is given by a Markov chain on
a tree. Given samples from the leaves of the Markov chain, the goal is to
reconstruct the leaf-labelled tree.
  It is well known that in order to reconstruct a tree on $n$ leaves, sample
sequences of length $\Omega(\log n)$ are needed. It was conjectured by M. Steel
that for the CFN/Ising evolutionary model, if the mutation probability on all
edges of the tree is less than $p^{\ast} = (\sqrt{2}-1)/2^{3/2}$, then the tree
can be recovered from sequences of length $O(\log n)$. The value $p^{\ast}$ is
given by the transition point for the extremality of the free Gibbs measure for
the Ising model on the binary tree. Steel's conjecture was proven by the second
author in the special case where the tree is "balanced." The second author also
proved that if all edges have mutation probability larger than $p^{\ast}$ then
the length needed is $n^{\Omega(1)}$. Here we show that Steel's conjecture
holds true for general trees by giving a reconstruction algorithm that recovers
the tree from $O(\log n)$-length sequences when the mutation probabilities are
discretized and less than $p^\ast$. Our proof and results demonstrate that
extremality of the free Gibbs measure on the infinite binary tree, which has
been studied before in probability, statistical physics and computer science,
determines how distinguishable are Gibbs measures on finite binary trees.


  We prove that a wide class of models of Markov neighbor-dependent
substitution processes on the integer line is solvable. This class contains
some models of nucleotide substitutions recently introduced and studied
empirically by molecular biologists. We show that the polynucleotide
frequencies at equilibrium solve explicit finite-size linear systems. Finally,
the dynamics of the process and the distribution at equilibrium exhibit some
stringent, rather unexpected, independence properties. For example, nucleotide
sites at distance at least three evolve independently, and the sites, if
encoded as purines and pyrimidines, evolve independently.


  We study some stochastic models of physical mapping of genomic sequences. Our
starting point is a global construction of the process of the clones and of the
process of the anchors which are used to map the sequence. This yields explicit
formulas for the moments of the proportion occupied by the anchored clones,
even in inhomogeneous models. This also allows to compare, in this respect,
inhomogeneous models to homogeneous ones. Finally, for homogeneous models, we
provide nonasymptotic bounds of the variance and we prove functional invariance
results.


  We study the iterated Galton-Watson process (IGW), possibly with thinning,
introduced by Gawe{\l}and Kimmel to model the number of repeats of DNA triplets
during some genetic disorders. If the process involves some thinning, then
extinction and explosion can have positive probability simultaneously. If the
underlying (simple) Galton-Watson process is nondecreasing with mean m, then,
conditionally on the explosion, the logarithm of the population of the IGW at
time n+1 is equivalent to log(m) times the population at time n, almost surely.
This simplifies arguments of Gawe{\l}and Kimmel, and confirms and extends a
conjecture of Pakes.


  We are interested in modelling Darwinian evolution, resulting from the
interplay of phenotypic variation and natural selection through ecological
interactions. Our models are rooted in the microscopic, stochastic description
of a population of discrete individuals characterized by one or several
adaptive traits. The population is modelled as a stochastic point process whose
generator captures the probabilistic dynamics over continuous time of birth,
mutation, and death, as influenced by each individual's trait values, and
interactions between individuals. An offspring usually inherits the trait
values of her progenitor, except when a mutation causes the offspring to take
an instantaneous mutation step at birth to new trait values. We look for
tractable large population approximations. By combining various scalings on
population size, birth and death rates, mutation rate, mutation step, or time,
a single microscopic model is shown to lead to contrasting macroscopic limits,
of different nature: deterministic, in the form of ordinary, integro-, or
partial differential equations, or probabilistic, like stochastic partial
differential equations or superprocesses. In the limit of rare mutations, we
show that a possible approximation is a jump process, justifying rigorously the
so-called trait substitution sequence. We thus unify different points of view
concerning mutation-selection evolutionary models.


  The alpha model, a parametrized family of probabilities on cladograms (rooted
binary leaf labeled trees), is introduced. This model is Markovian
self-similar, deletion-stable (sampling consistent), and passes through the
Yule, Uniform and Comb models. An explicit formula is given to calculate the
probability of any cladogram or tree shape under the alpha model. Sackin's and
Colless' index are shown to be $O(n^{1+\alpha})$ with asymptotic covariance
equal to 1. Thus the expected depth of a random leaf with $n$ leaves is
$O(n^\alpha)$. The number of cherries on a random alpha tree is shown to be
asymptotically normal with known mean and variance. Finally the shape of
published phylogenies is examined, using trees from Treebase.


  We present the formalism of sequential and asynchronous processes defined in
terms of random or quantum grammars and argue that these processes have
relevance in genomics. To make the article accessible to the
non-mathematicians, we keep the mathematical exposition as elementary as
possible, focusing on some general ideas behind the formalism and stating the
implications of the known mathematical results. We close with a set of open
challenging problems.


  Consider a haploid population which has evolved through an exchangeable
reproduction dynamics, and in which all individuals alive at time $t$ have a
most recent common ancestor (MRCA) who lived at time $A_t$, say. As time goes
on, not only the population but also its genealogy evolves: some families will
get lost from the population and eventually a new MRCA will be established. For
a time-stationary situation and in the limit of infinite population size $N$
with time measured in $N$ generations, i.e. in the scaling of population
genetics which leads to Fisher-Wright diffusions and Kingman's coalescent, we
study the process $\mathcal A = (A_t)$ whose jumps form the point process of
time pairs $(E,B)$ when new MRCAs are established and when they lived. By
representing these pairs as the entrance and exit time of particles whose
trajectories are embedded in the look-down graph of Donnelly and Kurtz (1999)
we can show by exchangeability arguments that the times $E$ as well as the
times $B$ from a Poisson process. Furthermore, the particle representation
helps to compute various features of the MRCA process, such as the distribution
of the coalescent at the instant when a new MRCA is established, and the
distribution of the number of MRCAs to come that live in today's past.


  We consider an interacting particle Markov process for Darwinian evolution in
an asexual population with non-constant population size, involving a linear
birth rate, a density-dependent logistic death rate, and a probability $\mu$ of
mutation at each birth event. We introduce a renormalization parameter $K$
scaling the size of the population, which leads, when $K\to+\infty$, to a
deterministic dynamics for the density of individuals holding a given trait. By
combining in a non-standard way the limits of large population ($K\to+\infty$)
and of small mutations ($\mu\to 0$), we prove that a time scales separation
between the birth and death events and the mutation events occurs and that the
interacting particle microscopic process converges for finite dimensional
distributions to the biological model of evolution known as the ``monomorphic
trait substitution sequence'' model of adaptive dynamics, which describes the
Darwinian evolution in an asexual population as a Markov jump process in the
trait space.


  We investigate projective varieties which are geometric models of binary
symmetric phylogenetic 3-valent trees. We prove that these varieties have
Gorenstein terminal singularities (with small resolution) and they are Fano
varieties of index 4. Moreover any two such varieties associated to trees with
the same number of leaves are deformation equivalent, that is, they are in the
same connected component of the Hilbert scheme of the projective space. As an
application we provide a simple formula for computing their Hilbert-Ehrhard
polynomial.


  The biological theory of adaptive dynamics proposes a description of the
long-term evolution of a structured asexual population. It is based on the
assumptions of large population, rare mutations and small mutation steps, that
lead to a deterministic ODE describing the evolution of the dominant type,
called the ``canonical equation of adaptive dynamics.'' Here, in order to
include the effect of stochasticity (genetic drift), we consider self-regulated
randomly fluctuating populations subject to mutation, so that the number of
coexisting types may fluctuate. We apply a limit of rare mutations to these
populations, while keeping the population size finite. This leads to a jump
process, the so-called ``trait substitution sequence,'' where evolution
proceeds by successive invasions and fixations of mutant types. Then we apply a
limit of small mutation steps (weak selection) to this jump process, that leads
to a diffusion process that we call the ``canonical diffusion of adaptive
dynamics,'' in which genetic drift is combined with directional selection
driven by the gradient of the fixation probability, also interpreted as an
invasion fitness. Finally, we study in detail the particular case of multitype
logistic branching populations and seek explicit formulae for the invasion
fitness of a mutant deviating slightly from the resident type. In particular,
second-order terms of the fixation probability are products of functions of the
initial mutant frequency, times functions of the initial total population size,
called the invasibility coefficients of the resident by increased fertility,
defence, aggressiveness, isolation or survival.


  Coalescents with multiple collisions, also known as $\Lambda$-coalescents,
were introduced by Pitman and Sagitov in 1999. These processes describe the
evolution of particles that undergo stochastic coagulation in such a way that
several blocks can merge at the same time to form a single block. In the case
that the measure $\Lambda$ is the $\operatorname {Beta}(2-\alpha,\alpha)$
distribution, they are also known to describe the genealogies of large
populations where a single individual can produce a large number of offspring.
Here, we use a recent result of Birkner et al. to prove that Beta-coalescents
can be embedded in continuous stable random trees, about which much is known
due to the recent progress of Duquesne and Le Gall. Our proof is based on a
construction of the Donnelly--Kurtz lookdown process using continuous random
trees, which is of independent interest. This produces a number of results
concerning the small-time behavior of Beta-coalescents. Most notably, we
recover an almost sure limit theorem of the present authors for the number of
blocks at small times and give the multifractal spectrum corresponding to the
emergence of blocks with atypical size. Also, we are able to find exact
asymptotics for sampling formulae corresponding to the site frequency spectrum
and the allele frequency spectrum associated with mutations in the context of
population genetics.


  Bacteria are able to respond to environmental signals by changing their rules
of movement. When we take into account chemical signals in the environment,
this behaviour is often called chemotaxis. At the individual-level, chemotaxis
consists of several steps. First, the cell detects the extracellular signal
using receptors on its membrane. Then, the cell processes the signal
information through the intracellular signal transduction network, and finally
it responds by altering its motile behaviour accordingly. At the population
level, chemotaxis can lead to aggregation of bacteria, travelling waves or
pattern formation, and the important task is to explain the population-level
behaviour in terms of individual-based models. It has been previously shown
that the transport equation framework is suitable for connecting different
levels of modelling of bacterial chemotaxis. In this paper, we couple the
transport equation for bacteria with the (parabolic/elliptic) equation for the
extracellular signals. We prove global existence of solutions for the general
hyperbolic chemotaxis models of cells which process the information about the
extracellular signal through the intracellular biochemical network and interact
by altering the extracellular signal as well. The conditions for global
existence in terms of the properties of the signal transduction model are
given.


  We consider the so called Moran process with frequency dependent fitness
given by a certain pay-off matrix. For finite populations, we show that the
final state must be homogeneous, and show how to compute the fixation
probabilities. Next, we consider the infinite population limit, and discuss the
appropriate scalings for the drift-diffusion limit. In this case, a degenerated
parabolic PDE is formally obtained that, in the special case of frequency
independent fitness, recovers the celebrated Kimura equation in population
genetics. We then show that the corresponding initial value problem is well
posed and that the discrete model converges to the PDE model as the population
size goes to infinity. We also study some game-theoretic aspects of the
dynamics and characterize the best strategies, in an appropriate sense.


  This paper has been withdrawn by the author due to a crucial error in the
definition of homomorphism.


  In this paper we study finite dynamical systems with $n$ functions acting on
the same set $X$, and probabilities assigned to these functions, that it is
called Probabilistic Regulatory Gene Networks (PRN. his concept is the same or
a natural generalization of the concept Probabilistic Boolean Networks (PBN),
introduced by I. Shmulevich, E. Dougherty, and W. Zhang, particularly the model
PBN has been using to describe genetic networks and has therapeutic
applications. In PRNs the most important question is to describe the steady
states of the systems, so in this paper we pay attention to the idea of
transforming a network to another without lost all the properties, in
particular the probability distribution. Following this objective we develop
the concepts of homomorphism and $\epsilon$-homomorphism of probabilistic
regulatory networks, since these concepts bring the properties from one
networks to another. Projections are special homomorphisms, and they always
induce invariant subnetworks that contain all cycles and steady states in the
network.


  In this paper we study homomorphisms of Probabilistic Regulatory Gene
Networks(PRN) introduced in arXiv:math.DS/0603289 v1 13 Mar 2006. The model PRN
is a natural generalization of the Probabilistic Boolean Networks (PBN),
introduced by I. Shmulevich, E. Dougherty, and W. Zhang in 2001, that has been
using to describe genetic networks and has therapeutic applications. In this
paper, our main objectives are to apply the concept of homomorphism and
$\epsilon$-homomorphism of probabilistic regulatory networks to the dynamic of
the networks. The meaning of $\epsilon$ is that these homomorphic networks have
similar distributions and the distance between the distributions is upper
bounded by $\epsilon$. Additionally, we prove that the class of PRN together
with the homomorphisms form a category with products and coproducts.
Projections are special homomorphisms, and they always induce invariant
subnetworks that contain all the cycles and steady states in the network. Here,
it is proved that the $\epsilon$-homomorphism for $0<\epsilon<1$ produce
simultaneous Markov Chains in both networks, that permit to introduce the
concept of $\epsilon$-isomorphism of Markov Chains, and similar networks.


  The interplay between space and evolution is an important issue in population
dynamics, that is in particular crucial in the emergence of polymorphism and
spatial patterns. Recently, biological studies suggest that invasion and
evolution are closely related. Here we model the interplay between space and
evolution starting with an individual-based approach and show the important
role of parameter scalings on clustering and invasion. We consider a stochastic
discrete model with birth, death, competition, mutation and spatial diffusion,
where all the parameters may depend both on the position and on the trait of
individuals. The spatial motion is driven by a reflected diffusion in a bounded
domain. The interaction is modelled as a trait competition between individuals
within a given spatial interaction range. First, we give an algorithmic
construction of the process. Next, we obtain large population approximations,
as weak solutions of nonlinear reaction-diffusion equations with Neumann's
boundary conditions. As the spatial interaction range is fixed, the
nonlinearity is nonlocal. Then, we make the interaction range decrease to zero
and prove the convergence to spatially localized nonlinear reaction-diffusion
equations, with Neumann's boundary conditions. Finally, simulations based on
the microscopic individual-based model are given, illustrating the strong
effects of the spatial interaction range on the emergence of spatial and
phenotypic diversity (clustering and polymorphism) and on the interplay between
invasion and evolution. The simulations focus on the qualitative differences
between local and nonlocal interactions.


  We establish the exact threshold for the reconstruction problem for a binary
asymmetric channel on the b-ary tree, provided that the asymmetry is
sufficiently small. This is the first exact reconstruction threshold obtained
in roughly a decade. We discuss the implications of our result for Glauber
dynamics, phylogenetic reconstruction, and so-called ``replica symmetry
breaking'' in spin glasses and random satisfiability problems.