  We demonstrate that Protein-Protein Interaction (PPI) networks in several
eucaryotic organisms contain significantly more self-interacting proteins than
expected if such homodimers randomly appeared in the course of the evolution.
We also show that on average homodimers have twice as many interaction partners
than non-self-interacting proteins. More specifically the likelihood of a
protein to physically interact with itself was found to be proportional to the
total number of its binding partners. These properties of dimers are are in
agreement with a phenomenological model in which individual proteins differ
from each other by the degree of their ``stickiness'' or general propensity
towards interaction with other proteins including oneself. A duplication of
self-interacting proteins creates a pair of paralogous proteins interacting
with each other. We show that such pairs occur more frequently than could be
explained by pure chance alone. Similar to homodimers, proteins involved in
heterodimers with their paralogs on average have twice as many interacting
partners than the rest of the network. The likelihood of a pair of paralogous
proteins to interact with each other was also shown to decrease with their
sequence similarity. This all points to the conclusion that most of
interactions between paralogs are inherited from ancestral homodimeric
proteins, rather than established de novo after the duplication. We finally
discuss possible implications of our empirical observations from functional and
evolutionary standpoints.


  Protein one-dimensional (1D) structures such as secondary structure and
contact number provide intuitive pictures to understand how the native
three-dimensional (3D) structure of a protein is encoded in the amino acid
sequence. However, it has not been clear whether a given set of 1D structures
contains sufficient information for recovering the underlying 3D structure.
Here we show that the 3D structure of a protein can be recovered from a set of
three types of 1D structures, namely, secondary structure, contact number and
residue-wise contact order which is introduced here for the first time. Using
simulated annealing molecular dynamics simulations, the structures satisfying
the given native 1D structural restraints were sought for 16 proteins of
various structural classes and of sizes ranging from 56 to 146 residues. By
selecting the structures best satisfying the restraints, all the proteins
showed a coordinate RMS deviation of less than 4\AA{} from the native
structure, and for most of them, the deviation was even less than 2\AA{}. The
present result opens a new possibility to protein structure prediction and our
understanding of the sequence-structure relationship.


  Background: Hidden Markov models (HMM) are powerful machine learning tools
successfully applied to problems of computational Molecular Biology. In a
predictive task, the HMM is endowed with a decoding algorithm in order to
assign the most probable state path, and in turn the class labeling, to an
unknown sequence. The Viterbi and the posterior decoding algorithms are the
most common. The former is very efficient when one path dominates, while the
latter, even though does not guarantee to preserve the automaton grammar, is
more effective when several concurring paths have similar probabilities. A
third good alternative is 1-best, which was shown to perform equal or better
than Viterbi. Results: In this paper we introduce the posterior-Viterbi (PV) a
new decoding which combines the posterior and Viterbi algorithms. PV is a two
step process: first the posterior probability of each state is computed and
then the best posterior allowed path through the model is evaluated by a
Viterbi algorithm.
  Conclusions: We show that PV decoding performs better than other algorithms
first on toy models and then on the computational biological problem of the
prediction of the topology of beta-barrel membrane proteins.


  The lack of specificity in microarray experiments due to non-specific
hybridization raises a serious problem for the analysis of microarray data
because the residual chemical background intensity is not related to the
expression degree of the gene of interest. We analyzed the concentration
dependence of the signal intensity of perfect match (PM) and mismatch (MM)
probes in terms using a microscopic binding model using a combination of mean
hybridization isotherms and single base related affinity terms. The signal
intensities of the PM and MM probes and their difference are assessed with
regard to their sensitivity, specificity and resolution for gene expression
measures. The presented theory implies the refinement of existing algorithms of
probe level analysis to correct microarray data for non-specific background
intensities.


  A simple deterministic model is proposed to represent the basic aspects
concerning the effects of different antiretroviral treatment schedules on HIV
incidence and prevalence of affected populations. The model mimics current
treatment guidelines applied in Brazil. However, the model does not intend to
fit the data with any acceptable degree of accuracy since uncertainties on the
values of the parameters and on the precise effect of the treatment put some
limits on the practical implications of our model from which only orders of
magnitude and some qualitative effects can be deduced. So, this paper intends
to provide a conceptual and mechanistic understanding of the possible long term
effects of treatment on the dynamics of HIV transmission.
  According to the model, the effect of the treatment depends on the level of
sexual activity of the subpopulations considered, being more pronounced on the
subpopulations with the highest sexual activity levels. Also, inefficient
treatment can be prejudicial depending on the level of sexual activity and on
the capacity to provide adequate treatment coverages to the population
affected.


  We study a minimal model for genome evolution whose elementary processes are
single site mutation, duplication and deletion of sequence regions and
insertion of random segments. These processes are found to generate long-range
correlations in the composition of letters as long as the sequence length is
growing, i.e., the combined rates of duplications and insertions are higher
than the deletion rate. For constant sequence length, on the other hand, all
initial correlations decay exponentially. These results are obtained
analytically and by simulations. They are compared with the long-range
correlations observed in genomic DNA, and the implications for genome evolution
are discussed.


  Information capacity of a symbol sequence is a measure of the unexpectedness
of a continuation of given string of symbols. Continuation of a string is
determined through the maximum entropy of the reconstructed frequency
dictionary; the capacity, in turn, is determined through the calculation of
mutual entropy of a real frequency dictionary of a sequence with respect to the
reconstructed one. The capacity does not depend on the length of strings in a
dictionary. The capacity calculated for various genomes exhibits a multi-minima
pattern reflecting an order observed within a sequence.


  We proposed a probabilistic algorithm to solve the Multiple Sequence
Alignment problem. The algorithm is a Simulated Annealing (SA) that exploits
the representation of the Multiple Alignment between $D$ sequences as a
directed polymer in $D$ dimensions. Within this representation we can easily
track the evolution in the configuration space of the alignment through local
moves of low computational cost. At variance with other probabilistic
algorithms proposed to solve this problem, our approach allows for the creation
and deletion of gaps without extra computational cost. The algorithm was tested
aligning proteins from the kinases family. When D=3 the results are consistent
with those obtained using a complete algorithm. For $D>3$ where the complete
algorithm fails, we show that our algorithm still converges to reasonable
alignments. Moreover, we study the space of solutions obtained and show that
depending on the number of sequences aligned the solutions are organized in
different ways, suggesting a possible source of errors for progressive
algorithms.


  The thermosynthesis concept, biological free energy gain from thermal
cycling, is combined with the concept of the RNA World. The resulting overall
origin of life model gives new explanations for the emergence of the genetic
code and the ribosome. The first protein named pF1 obtains the energy to
support the RNA world by a thermal variation of F1 ATP synthase's binding
change mechanism. This pF1 is the single translation product during the
emergence of the genetic machinery. During thermal cycling pF1 condenses many
substrates with broad specificity, yielding NTPs and randomly constituted
protein and RNA libraries that contain (self)-replicating RNA. The smallness of
pF1 permits the emergence of the genetic machinery by selection of RNA that
increases the fraction of pF1s in the protein library: (1) a progenitor of rRNA
that concatenates amino acids bound to (2) a chain of 'positional tRNAs' linked
by mutual recognition, yielding a pF1 (or its main motif); this positional tRNA
set gradually evolves to a set of regular tRNAs functioning according to the
genetic code, with concomitant emergence of (3) an mRNA coding for pF1.


  We investigate the translocation of a single stranded DNA (ssDNA) through a
pore, which fluctuates between two conformations, by using coupled master
equations (ME). The probability density function (PDF) of the first passage
times (FPT) of the translocation process is calculated, displaying a triple,
double or mono-peaked behavior, depending on the system parameters. An
analytical expression for the mean first passage time (MFPT) of the
translocation process is derived, and provides an extensive characterization of
the translocation process.


  Residue-wise contact order (RWCO) is a new kind of one-dimensional protein
structures which represents the extent of long-range contacts. We have recently
shown that a set of three types of one-dimensional structures (secondary
structure, contact number, and RWCO) contains sufficient information for
reconstructing the three-dimensional structure of proteins. Currently, there
exist prediction methods for secondary structure and contact number from amino
acid sequence, but none exists for RWCO. Also, the properties of amino acids
that affect RWCO is not clearly understood. Here, we present a linear
regression-based method to predict RWCO from amino acid sequence, and analyze
the regression parameters to identify the properties that correlates with the
RWCO. The present method achieves the significant correlation of 0.59 between
the native and predicted RWCOs on average. An unusual feature of the RWCO
prediction is the remarkably large optimal half window size of 26 residues. The
regression parameters for the central and near-central residues of the local
sequence segment highly correlate with those of the contact number prediction,
and hence with hydrophobicity.


  Complex dynamical networks consisting of many components that interact and
produce each other are difficult to understand, especially, when new components
may appear. In this paper we outline a theory to deal with such systems. The
theory consists of two parts. The first part introduces the concept of a
chemical organization as a closed and mass-maintaining set of components. This
concept allows to map a complex (reaction) network to the set of organizations,
providing a new view on the system's structure. The second part connects
dynamics with the set of organizations, which allows to map a movement of the
system in state space to a movement in the set of organizations.


  The presence of neighbor dependencies generated a specific pattern of
dinucleotide frequencies in all organisms. Especially, the
CpG-methylation-deamination process is the predominant substitution process in
vertebrates and needs to be incorporated into a more realistic model for
nucleotide substitutions. Based on a general framework of nucleotide
substitutions we develop a method that is able to identify the most relevant
neighbor dependent substitution processes, measure their strength, and judge
their importance to be included into the modeling. Starting from a model for
neighbor independent nucleotide substitution we successively add neighbor
dependent substitution processes in the order of their ability to increase the
likelihood of the model describing given data. The analysis of neighbor
dependent nucleotide substitutions in human, zebrafish and fruit fly is
presented. A web server to perform the presented analysis is publicly
available.


  What are proteins made from, as the working parts of the living cells protein
machines? To answer this question, we need a technology to disassemble proteins
onto elementary func-tional details and to prepare lumped description of such
details. This lumped description might have a multiple material realization (in
amino acids). Our hypothesis is that informational approach to this problem is
possible. We propose a way of hierarchical classification that makes the
primary structure of protein maximally non-random. The first steps of the
suggested research program are realized: the method and the analysis of optimal
informational protein binary alphabet. The general method is used to answer
several specific questions, for example: (i) Is there a syntactic difference
between Globular and Membrane proteins? (ii) Are proteins random sequences of
amino acids (a long discussion)? For these questions, the answers are as
follows: (i) There exists significant syntactic difference between Globular and
Membrane proteins, and this difference is described; (ii) Amino acid sequences
in proteins are definitely not random.


  This study presents the first global, 1 Mbp level analysis of patterns of
nucleotide substitutions along the human lineage. The study is based on the
analysis of a large amount of repetitive elements deposited into the human
genome since the mammalian radiation, yielding a number of results that would
have been difficult to obtain using the more conventional comparative method of
analysis. This analysis revealed substantial and consistent variability of
rates of substitution, with the variability ranging up to 2-fold among
different regions. The rates of substitutions of C or G nucleotides with A or T
nucleotides vary much more sharply than the reverse rates suggesting that much
of that variation is due to differences in mutation rates rather than in the
probabilities of fixation of C/G vs. A/T nucleotides across the genome. For all
types of substitution we observe substantially more hotspots than coldspots,
with hotspots showing substantial clustering over tens of Mbp's. Our analysis
revealed that GC-content of surrounding sequences is the best predictor of the
rates of substitution. The pattern of substitution appears very different near
telomeres compared to the rest of the genome and cannot be explained by the
genome-wide correlations of the substitution rates with GC content or exon
density. The telomere pattern of substitution is consistent with natural
selection or biased gene conversion acting to increase the GC-content of the
sequences that are within 10-15 Mbp away from the telomere.


  In this paper we introduce a class of stochastic population models based on
"patch dynamics". The size of the patch may be varied, and this allows one to
quantify the departures of these stochastic models from various mean field
theories, which are generally valid as the patch size becomes very large. These
models may be used to formulate a broad range of biological processes in both
spatial and non-spatial contexts. Here, we concentrate on two-species
competition. We present both a mathematical analysis of the patch model, in
which we derive the precise form of the competition mean field equations (and
their first order corrections in the non-spatial case), and simulation results.
These mean field equations differ, in some important ways, from those which are
normally written down on phenomenological grounds. Our general conclusion is
that mean field theory is more robust for spatial models than for a single
isolated patch. This is due to the dilution of stochastic effects in a spatial
setting resulting from repeated rescue events mediated by inter-patch
diffusion. However, discrete effects due to modest patch sizes lead to striking
deviations from mean field theory even in a spatial setting.


  In this paper we present the simplest individual level model of predator-prey
dynamics and show, via direct calculation, that it exhibits cycling behavior.
The deterministic analogue of our model, recovered when the number of
individuals is infinitely large, is the Volterra system (with density-dependent
prey reproduction) which is well-known to fail to predict cycles. This
difference in behavior can be traced to a resonant amplification of demographic
fluctuations which disappears only when the number of individuals is strictly
infinite. Our results indicate that additional biological mechanisms, such as
predator satiation, may not be necessary to explain observed predator-prey
cycles in real (finite) populations.


  Based on experimental observations it is known that various biological cells
exhibit a persistent random walk during migration on flat substrates. The
persistent random walk is characterized by `stop-and-go' movements :
unidirectional motions over distances of the order of several cell diameter are
separated by localized short time erratic movements. Using computer simulations
the reasons for this phenomena had been unveiled and shown to be attributed to
two antagonistic nucleation processes during the polymerization of the cell's
actin cytoskeleton : the (ordinary) spontaneous nucleation and the dendritic
nucleation processes. Whereas spontaneous nucleations generate actin filaments
growing in different directions and hence create motions in random directions,
dendritic nucleations provide a unidirectional growth. Since dendritic growth
exhibits stochastic fluctuations, spontaneous nucleation may eventually compete
or even dominate, which results in a reorientation of filament growth and hence
a new direction of cell motion. The event of reorientation takes place at
instants of vanishing polarity of the actin skeleton.


  We present a simple model of adaptive radiations in evolution based on
species competition. Competition is found to promote species divergence and
branching, and to dampen the net species production. In the model simulations,
high taxonomic diversification and branching take place during the beginning of
the radiation. The results show striking similarities with empirical data and
highlight the mechanism of competition as an important driving factor for
accelerated evolutionary transformation.


  Network of packages with regulatory interactions (dependences and conflicts)
from Debian GNU/Linux operating system is compiled and used as analogy of a
gene regulatory network. Using a trace-back algorithm we assembly networks from
the potential pool of packages for both scale-free and exponential topology
from real and a null model data, respectively. We calculate the maximum number
of packages that can be functionally installed in the system (i.e., the active
network size). We show that scale-free regulatory networks allow a larger
active network size than random ones. Small genomes with scale-free regulatory
topology could allow much more functionality than large genomes with an
exponential one, with implications on its dynamics, robustness and evolution.


  We consider the evolutionary trajectories traced out by an infinite
population undergoing mutation-selection dynamics in static, uncorrelated
random fitness landscapes. Starting from the population that consists of a
single genotype, the most populated genotype \textit{jumps} from a local
fitness maximum to another and eventually reaches the global maximum. We use a
strong selection limit, which reduces the dynamics beyond the first time step
to the competition between independent mutant subpopulations, to study the
dynamics of this model and of a simpler one-dimensional model which ignores the
geometry of the sequence space. We find that the fit genotypes that appear
along a trajectory are a subset of suitably defined fitness \textit{records},
and exploit several results from the record theory for non-identically
distributed random variables. The genotypes that contribute to the trajectory
are those records that are not \textit{bypassed} by superior records arising
further away from the initial population. Several conjectures concerning the
statistics of bypassing are extracted from numerical simulations. In
particular, for the one-dimensional model, we propose a simple relation between
the bypassing probability and the dynamic exponent which describes the scaling
of the typical evolution time with genome size. The latter can be determined
exactly in terms of the extremal properties of the fitness distribution.


  Trajectories of on-off events are the output of many single molecule
experiments. Usually, one describes the underlying mechanism that generates the
trajectory using a kinetic scheme, and by analyzing the trajectory aims at
deducing this scheme. In a previous work [O. Flomenbom, J. Klafter, and A.
Szabo, submitted (2004)], we showed that when successive events along a
trajectory are uncorrelated, all the information in the trajectory is contained
in two basic functions, which are the waiting time probability functions (PDFs)
of the on state and of the off state. The kinetic schemes that lead to such
uncorrelated trajectories were termed reducible. Here we discuss the reasons
that lead to reducible schemes. In particular, the topology of reducible
schemes is characterized and proven.


  In molecular phylogeny, relationships among organisms are reconstructed using
DNA or protein sequences and are displayed as trees. A linear increase in the
number of sequences results in an exponential increase of possible trees. Thus,
inferring trees from molecular data was shown to be NP-hard. This causes
problems, if large data sets are used. This review gives an introduction to
molecular phylogenetic methods and to the problems biologists are facing in
molecular phylogenetic analyses.


  In this paper the heat signaling in microtubules (MT) is investigated. It is
argued that for the description of the heat signaling phenomena in MT, the
hyperbolic heat transport (HHT) equation must be used. It is shown that HHT is
the Klein-Gordon (K-G) equation. The general solution for the K-G equation for
MT is obtained. For the undistorted signal propagation in MT the Heisenberg
uncertainty principle is formulated and discussed.
  Key words: Microtubules; Heat signaling; Klein-Gordon equation; Heisenberg
principle.


  Motor proteins that specifically interact with the ends of cytoskeletal
filaments can induce filament depolymerization. A phenomenological description
of this process is presented. We show that under certain conditions motors
dynamically accumulate at the filament ends. We compare simulations of two
microscopic models to the phenomenological description. The depolymerization
rate can exhibit maxima and dynamic instabilities as a function of the bulk
motor density for processive depolymerization. We discuss our results in
relation to experimental studies of Kin-13 family motor proteins.


  We model the competition between recombination and point mutation in
microbial genomes, and present evidence for two distinct phases, one uniform,
the other genetically diverse. Depending on the specifics of homologous
recombination, we find that global sequence divergence can be mediated by
fronts propagating along the genome, whose characteristic signature on genome
structure is elucidated, and apparently observed in closely-related {\it
Bacillus} strains. Front propagation provides an emergent, generic mechanism
for microbial "speciation", and suggests a classification of microorganisms on
the basis of their propensity to support propagating fronts.


  It was shown by Gillespie (1974) that if two genotypes produce the same
average number of offspring on but have a different variance associated within
each generation, the genotype with a lower variance will have a higher
effective fitness. Specifically, the effective fitness is W(effective)=w-var/N,
where w is the mean fitness, var is the variance in offspring number, and N is
the total population size. The model also predicts that if a strategy has a
higher arithmetic mean fitness and a higher variance than the competitor, the
outcome of selection will depend on the population size (with larger population
sizes favoring the high variance, high mean genotype). This suggests that for
metapopulations with large numbers of (relatively) small demes, a strategy with
lower variance and lower mean may be favored if the migration rate is low while
higher migration rates (consistent with a larger effective population size)
favor the opposite strategy. Individual based simulation confirms that this is
indeed the case for an island model of migration, though the effect of
migration differs greatly depending on whether migration precedes or follows
selection. It is noted in the appendix that while Gillespie 1974 does seem to
be heuristically accurate, it is not clear that the definition of effective
fitness follows from his derivation.


  Motivated by the finding that there is some biological universality in the
relationship between school geometry and school biomass of various pelagic
fishes in various conditions, I here establish a scaling law for school
dimensions: the school diameter increases as a power-law function of school
biomass. The power-law exponent is extracted through the data collapse, and is
close to 3/5. This value of the exponent implies that the mean packing density
decreases as the school biomass increases, and the packing structure displays a
mass-fractal dimension of 5/3. By exploiting an analogy between school geometry
and polymer chain statistics, I examine the behavioral algorithm governing the
swollen conformation of large-sized schools of pelagics, and I explain the
value of the exponent.


  By starting from the four DNA bases order in the Boolean lattice, a novel Lie
Algebra of the genetic code is proposed. Here, the principal partitions of the
genetic code table were obtained as equivalent classes of quotient subspaces of
the genetic code vector space over the Galois field of the four DNA bases. The
new algebraic structure shows strong connections among algebraic relationships,
codon assignment and physicochemical properties of amino acids. Moreover, a
distance function defined between codons in the Lie algebra was demonstrated to
have a linear behavior respect to physical variables such as the mean of amino
acids energy in proteins. It was also noticed that the distance between wild
type and mutant codons approach smaller values in mutational variants of four
genes, i.e, human phenylalanine hydroxylase, human beta-globin, HIV-1 protease
and HIV-1 reverse transcriptase. These results strongly suggest that
deterministic rules in genetic code origin must be involved.


  Interactions between genes and gene products give rise to complex circuits
that enable cells to process information and respond to external signals.
Theoretical studies often describe these interactions using continuous,
stochastic, or logical approaches. We propose a new modeling framework for gene
regulatory networks, that combines the intuitive appeal of a qualitative
description of gene states with a high flexibility in incorporating
stochasticity in the duration of cellular processes. We apply our methods to
the regulatory network of the segment polarity genes, thus gaining novel
insights into the development of gene expression patterns. For example, we show
that very short synthesis and decay times can perturb the wild type pattern. On
the other hand, separation of timescales between pre- and posttranslational
processes and a minimal prepattern ensure convergence to the wild type
expression pattern regardless of fluctuations.


  The promise of discovering a functional blueprint of a cellular system from
large-scale and high-throughput sequence and experimental data is predicated on
the belief that the same top-down investigative approach that proved successful
in other biological problems (e.g. DNA sequencing) will be as effective when it
comes to inferring more complex intracellular processes. The results in this
paper address this fundamental issue in the specific context of transcription
regulatory networks. Although simple recurring regulatory motifs have been
identified in the past, due to the size and complexity of the connectivity
structure, the subdivision of such networks into larger, and possibly
inter-connected, regulatory modules is still under investigation. Specifically,
it is unclear whether functionally well-characterized transcriptional
sub-networks can be identified by solely analyzing the connectivity structure
of the overall network topology. In this paper, we show that transcriptional
regulatory networks can be systematically partitioned into communities whose
members are consistently functionally related. We applied the partitioning
method to the transcriptional regulatory networks of the yeast Saccharomyces
cerevisiae; the resulting communities of gene and transcriptional regulators
can be associated to distinct functional units, such as amino acid metabolism,
cell cycle regulation, protein biosynthesis and localization, DNA replication
and maintenance, lipid catabolism, stress response and so on. Moreover, the
observation of inter-community connectivity patterns provides a valuable tool
for elucidating the inter-dependency between the discovered regulatory modules.


  We studied the effects of epidermal growth factor receptor (EGFR) density on
tumor growth dynamics, both on the sub- and the multi-cellular level using our
previously developed model. This algorithm simulates the growth of a brain
tumor using a multi-scale two-dimensional agent-based approach with an
integrated transforming growth factor alpha (TGFalpha) induced
EGFR-gene-protein interaction network. The results confirm that increasing cell
receptor density correlates with an acceleration of the tumor system's
spatio-temporal expansion dynamics. This multicellular behavior cannot be
explained solely on the basis of spatial sub-cellular dynamics, which remain
qualitatively similar amongst the three glioma cell lines investigated here in
silico. Rather, we find that cells with higher EGFR density show an early
increase in the phenotypic switching activity between proliferative and
migratory traits, linked to a higher level of initial auto-stimulation by the
PLCgamma-mediated TGFalpha-EGFR autocrine network. This indicates a more active
protein level interaction in these chemotactically acting tumor systems and
supports the role of post-translational regulation for the implemented EGFR
pathway. Implications of these results for experimental cancer research are
discussed.


  A population of complete subgraphs or cliques in a network evolving via
duplication-divergence is considered. We find that a number of cliques of each
size scales linearly with the size of the network. We also derive a clique
population distribution that is in perfect agreement with both the simulation
results and the clique statistic of the protein-protein binding network of the
fruit fly. In addition, we show that such features as fat-tail degree
distribution, various rates of average degree growth and non-averaging,
revealed recently for only the particular case of a completely asymmetric
divergence, are present in a general case of arbitrary divergence.


  A time trajectory of an observable that fluctuates between two values (say,
on and off), stemming from some unknown multi-substate kinetic scheme, is the
output of many single molecule experiments. Here we show that when all
successive waiting times along the trajectory are uncorrelated the on and the
off waiting time probability density functions (PDFs) contain all the
information. By relating the lack of correlation in the trajectory to the
topology of kinetic schemes, we can immediately specify those kinetic schemes
that are equally consistent with experiment, which means that it is impossible
to differentiate between them by any sophisticated analyses of the trajectory.
Correlated trajectories, however, contain additional information about the
underlying kinetic scheme, and we consider the strategy that one should use to
extract it. An example is given on correlations in the activity of individual
lipase molecules.


  We present a mathematical analysis of the speciation model for food-web
structure, which had in previous work been shown to yield a good description of
empirical data of food-web topology. The degree distributions of the network
are derived. Properties of the speciation model are compared to those of other
models that successfully describe empirical data. It is argued that the
speciation model unifies the underlying ideas of previous theories. In
particular, it offers a mechanistic explanation for the success of the niche
model of Williams and Martinez and the frequent observation of intervality in
empirical food webs.


  A fully self-contained model of homochirality is presented that contains the
effects of both polymerization and dissociation. The dissociation fragments are
assumed to replenish the substrate from which new monomers can grow and undergo
new polymerization. The mean length of isotactic polymers is found to grow
slowly with the normalized total number of corresponding building blocks.
Alternatively, if one assumes that the dissociation fragments themselves can
polymerize further, then this corresponds to a strong source of short polymers,
and an unrealistically short average length of only 3. By contrast, without
dissociation, isotactic polymers becomes infinitely long.


  The human genome contains repetitive DNA at different level of sequence
length, number and dispersion. Highly repetitive DNA is particularly rich in
homo-- and di--nucleotide repeats, while middle repetitive DNA is rich of
families of interspersed, mobile elements hundreds of base pairs (bp) long,
among which the Alu families. A link between homo- and di-polymeric tracts and
mobile elements has been recently highlighted. In particular, the mobility of
Alu repeats, which form 10% of the human genome, has been correlated with the
length of poly(A) tracts located at one end of the Alu. These tracts have a
rigid and non-bendable structure and have an inhibitory effect on nucleosomes,
which normally compact the DNA. We performed a statistical analysis of the
genome-wide distribution of lengths and inter--tract separations of poly(X) and
poly(XY) tracts in the human genome. Our study shows that in humans the length
distributions of these sequences reflect the dynamics of their expansion and
DNA replication. By means of general tools from linguistics, we show that the
latter play the role of highly-significant content-bearing terms in the DNA
text. Furthermore, we find that such tracts are positioned in a non-random
fashion, with an apparent periodicity of 150 bases. This allows us to extend
the link between repetitive, highly mobile elements such as Alus and
low-complexity words in human DNA. More precisely, we show that Alus are
sources of poly(X) tracts, which in turn affect in a subtle way the combination
and diversification of gene expression and the fixation of multigene families.


  Observations of the evolution of species groups in nature, such as well
recognized Galapagos finches, have motivated much theoretical research aimed at
understanding the processes associated with such radiations. The Penna model is
one such model and has been widely used to study aging. In this paper we use
the basic Penna model to investigate the process of sympatric speciation in a
simple food web model. Initially our web consists of a primary food source and
a single herbivore species that feeds on this resource. Subsequently we
introduce a predator that feeds on the herbivore. In both instances we directly
manipulate the food source, that is, its size distribution, and monitore the
changes in the populations structures. Sympatric speciation is obtained for the
consumer species in both webs, and our results confirm that the speciation
velocity depends on how far up,in the food chain, the focus population is
feeding. Simulations are done with three different sexual imprinting-like
mechanisms, in order to discuss adaptation by natural selection.


  We determine both barrier heights and prefactors for protein folding by
applying constraints determined from experimental rate measurements to a
Kramers theory for folding rate. The theoretical values are required to match
the experimental values at two conditions of temperature and denaturant that
induce the same stability. Several expressions for the prefactor in the Kramers
rate equation are examined: a random energy approximation, a correlated energy
approximation, and an approximation using a single Arrhenius activation energy.
Barriers and prefactors are generally found to be large as a result of
implementing this recipe, i.e. the folding landscape is cooperative and smooth.
Interestingly, a prefactor with a single Arrhenius activation energy admits no
formal solution.


  We demonstrate a novel all-optical switch consisting of two molecules: a
primary fluorophore that can be switched between a fluorescent and a dark state
by light of different wavelengths, and a secondary chromophore that facilitates
switching. The interaction between the two molecules exhibits a distance
dependence much steeper than that of Forster resonance energy transfer. This
enables the switch to act as a ruler with the capability to probe distances
difficult to access by other spectroscopic methods, thus presenting a new tool
for the study of biomolecules at the single-molecule level.


  Many growing networks possess accelerating statistics where the number of
links added with each new node is an increasing function of network size so the
total number of links increases faster than linearly with network size. In
particular, biological networks can display a quadratic growth in regulator
number with genome size even while remaining sparsely connected. These features
are mutually incompatible in standard treatments of network theory which
typically require that every new network node possesses at least one
connection. To model sparsely connected networks, we generalize existing
approaches and add each new node with a probabilistic number of links to
generate either accelerating, hyper-accelerating, or even decelerating network
statistics in different regimes. Under preferential attachment for example,
slowly accelerating networks display stationary scale-free statistics
relatively independent of network size while more rapidly accelerating networks
display a transition from scale-free to exponential statistics with network
growth. Such transitions explain, for instance, the evolutionary record of
single-celled organisms which display strict size and complexity limits.


  We extend replica exchange simulation in two ways, and apply our approaches
to biomolecules. The first generalization permits exchange simulation between
models of differing resolution -- i.e., between detailed and coarse-grained
models. Such ``resolution exchange'' can be applied to molecular systems or
spin systems. The second extension is to ``pseudo-exchange'' simulations, which
require little CPU usage for most levels of the exchange ladder and also
substantially reduces the need for overlap between levels. Pseudo exchanges can
be used in either replica or resolution exchange simulations. We perform
efficient, converged simulations of a 50-atom peptide to illustrate the new
approaches.


  In the last years, tens of thousands gene expression profiles for cells of
several organisms have been monitored. Gene expression is a complex
transcriptional process where mRNA molecules are translated into proteins,
which control most of the cell functions. In this process, the correlation
among genes is crucial to determine the specific functions of genes. Here, we
propose a novel multi-dimensional stochastic approach to deal with the gene
correlation phenomena. Interestingly, our stochastic framework suggests that
the study of the gene correlation requires only one theoretical assumption
-Markov property- and the experimental transition probability, which
characterizes the gene correlation system. Finally, a gene expression
experiment is proposed for future applications of the model.


  In the past years, the folding kinetics of many small single-domain proteins
has been characterized by mutational Phi-value analysis. In this article, a
simple, essentially parameter-free model is introduced which derives folding
routes from native structures by minimizing the entropic loop-closure cost
during folding. The model predicts characteristic folding sequences of
structural elements such as helices and beta-strand pairings. Based on few
simple rules, the kinetic impact of these structural elements is estimated from
the routes and compared to average experimental Phi-values for the helices and
strands of 15 small, well-characterized proteins. The comparison leads on
average to a correlation coefficient of 0.62 for all proteins with polarized
Phi-value distributions, and 0.74 if distributions with negative average
Phi-values are excluded. The diffuse Phi-value distributions of the remaining
proteins are reproduced correctly. The model shows that Phi-value
distributions, averaged over secondary structural elements, can often be traced
back to entropic loop-closure events, but also indicates energetic preferences
in the case of a few proteins governed by parallel folding processes.


  Identity, response to external stimuli, and spatial architecture of a living
system are central topics of molecular biology. Presently, they are largely
seen as a result of the interplay between a gene repertoire and the regulatory
machinery of the cell. At the transcriptional level, the cis-regulatory regions
establish sets of interdependencies between transcription factors and genes,
including other transcription factors. These ``transcription networks'' are too
large to be approached globally with a detailed dynamical model. In this paper,
we describe an approach to this problem that focuses solely on the
compatibility between gene expression patterns and signal integration
functions, discussing calculations carried on the simplest, Boolean,
realization of the model, and a first application to experimental data sets.


  Thanks to recent progress in high-throughput experimental techniques, the
datasets of large-scale protein interactions of prototypical multicellular
species, the nematode worm Caenorhabditis elegans and the fruit fly Drosophila
melanogaster, have been assayed. The datasets are obtained mainly by using the
yeast hybrid method, which contains false-positive and false-negative
simultaneously. Accordingly, while it is desirable to test such datasets
through further wet experiments, here we invoke recent developed network theory
to test such high throughput datasets in a simple way. Based on the fact that
the key biological processes indispensable to maintaining life are universal
across eukaryotic species, and the comparison of structural properties of the
protein interaction networks (PINs) of the two species with those of the yeast
PIN, we find that while the worm and the yeast PIN datasets exhibit similar
structural properties, the current fly dataset, though most comprehensively
screened ever, does not reflect generic structural properties correctly as it
is. The modularity is suppressed and the connectivity correlation is lacking.
Addition of interlogs to the current fly dataset increases the modularity and
enhances the occurrence of triangular motifs as well. The connectivity
correlation function of the fly, however, remains distinct under such interlogs
addition, for which we present a possible scenario through an in silico
modeling.


  The importance of understanding the mechanism of protein aggregation into
insoluble amyloid fibrils relies not only on its medical consequences, but also
on its more basic properties of self--organization. The discovery that a large
number of uncorrelated proteins can form, under proper conditions, structurally
similar fibrils has suggested that the underlying mechanism is a general
feature of polypeptide chains. In the present work, we address the early events
preceeding amyloid fibril formation in solutions of zinc--free human insulin
incubated at low pH and high temperature. Aside from being a easy--to--handle
model for protein fibrillation, subcutaneous aggregation of insulin after
injection is a nuisance which affects patients with diabetes. Here, we show by
time--lapse atomic force microscopy (AFM) that a steady-state distribution of
protein oligomers with an exponential tail is reached within few minutes after
heating. This metastable phase lasts for few hours until aggregation into
fibrils suddenly occurs. A theoretical explanation of the oligomer
pre--fibrillar distribution is given in terms of a simple
coagulation--evaporation kinetic model, in which concentration plays the role
of a critical parameter. Due to high resolution and sensitivity of AFM
technique, the observation of a long-lasting latency time should be considered
an actual feature of the aggregation process, and not simply ascribed to
instrumental inefficency. These experimental facts, along with the kinetic
model used, claim for a critical role of thermal concentration fluctuations in
the process of fibril nucleation.


  Motivation: Standard algorithms for pairwise protein sequence alignment make
the simplifying assumption that amino acid substitutions at neighboring sites
are uncorrelated. This assumption allows implementation of fast algorithms for
pairwise sequence alignment, but it ignores information that could conceivably
increase the power of remote homolog detection. We examine the validity of this
assumption by constructing extended substitution matrixes that encapsulate the
observed correlations between neighboring sites, by developing an efficient and
rigorous algorithm for pairwise protein sequence alignment that incorporates
these local substitution correlations, and by assessing the ability of this
algorithm to detect remote homologies. Results: Our analysis indicates that
local correlations between substitutions are not strong on the average.
Furthermore, incorporating local substitution correlations into pairwise
alignment did not lead to a statistically significant improvement in remote
homology detection. Therefore, the standard assumption that individual residues
within protein sequences evolve independently of neighboring positions appears
to be an efficient and appropriate approximation.


  This is the first of two papers where we discuss the limits imposed by
competition to the biodiversity of species communities. In this first paper we
study the coexistence of competing species at the fixed point of population
dynamic equations. For many simple models, this imposes a limit on the width of
the productivity distribution, which is more severe the more diverse the
ecosystem is (Chesson, 1994). Here we review and generalize this analysis,
beyond the ``mean-field''-like approximation of the competition matrix used in
previous works, and extend it to structured food webs. In all cases analysed,
we obtain qualitatively similar relations between biodiversity and competition:
the narrower the productivity distribution is, the more species can stably
coexist. We discuss how this result, considered together with environmental
fluctuations, limits the maximal biodiversity that a trophic level can host.


  This is the second of two papers dedicated to the relationship between
population models of competition and biodiversity. Here we consider species
assembly models where the population dynamics is kept far from fixed points
through the continuous introduction of new species, and generalize to such
models thecoexistence condition derived for systems at the fixed point. The
ecological overlap between species with shared preys, that we define here,
provides a quantitative measure of the effective interspecies competition and
of the trophic network topology. We obtain distributions of the overlap from
simulations of a new model based both on immigration and speciation, and show
that they are in good agreement with those measured for three large natural
food webs. As discussed in the first paper, rapid environmental fluctuations,
interacting with the condition for coexistence of competing species, limit the
maximal biodiversity that a trophic level can host. This horizontal limitation
to biodiversity is here combined with either dissipation of energy or growth of
fluctuations, which in our model limit the length of food webs in the vertical
direction. These ingredients yield an effective model of food webs that produce
a biodiversity profile with a maximum at an intermediate trophic level, in
agreement with field studies.


  The stationary states of systems with many molecular motors are studied
theoretically for uniaxial and centered (aster-like) arrangements of
cytoskeletal filaments using Monte Carlo simulations and a two-state model.
Mutual exclusion of motors from binding sites of the filaments is taken into
account. For small overall motor concentration, the density profiles are
exponential and algebraic in uniaxial and centered filament systems,
respectively. For uniaxial systems, exclusion leads to the coexistence of
regions of high and low densities of bound motors corresponding to motor
traffic jams, which grow upon increasing the overall motor concentration. These
jams are insensitive to the motor behavior at the end of the filament. In
centered systems, traffic jams remain small and an increase in the motor
concentration leads to a flattening of the profile, if the motors move inwards,
and to the build-up of a concentration maximum in the center of the aster if
motors move outwards. In addition to motors density patterns, we also determine
the corresponding patterns of the motor current.


  We study a stochastic model proposed recently in the genetic literature to
explain the heterogeneity of cell populations or of gene products. Cells are
located in two colonies, whose sizes fluctuate as birth and migration processes
in switching environments. We prove that there is a range of parameters where
heterogeneity induces a larger mean fitness


  The Penna model is a strategy to simulate the genetic dynamics of
age-structured populations, in which the individuals genomes are represented by
bit-strings. It provides a simple metaphor for the evolutionary process in
terms of the mutation accumulation theory. In its original version, an
individual dies due to inherited diseases when its current number of
accumulated mutations, n, reaches a threshold value, T. Since the number of
accumulated diseases increases with age, the probability to die is zero for
very young ages (n < T) and equals 1 for the old ones (n >= T). Here, instead
of using a step function to determine the genetic death age, we test several
other functions that may or may not slightly increase the death probability at
young ages (n < T), but that decreases this probability at old ones. Our
purpose is to study the oldest old effect, that is, a plateau in the mortality
curves at advanced ages. Imposing certain conditions, it has been possible to
obtain a clear plateau using the Penna model. However, a more realistic one
appears when a modified version, that keeps the population size fixed without
fluctuations, is used. We also find a relation between the birth rate, the
age-structure of the population and the death probability.


  Preparations of arachidonic acid binding and non-binding proteins from bovine
brain, four acidic proteins (alpha-casein, phosvitin, beta-lactoglobulin A and
B), the peptide polyglutamate, and two polyanions (heparin, dextran sulfate)
enhanced both basal and glucose 6-phosphate induced solubilization of rat brain
mitochondrial hexokinase (ATP:D-hexose 6-phosphotransferase, EC 2.7.1.1). In
contrast, three other acidic proteins, had little (alpha-lactalbumin) or no
effect (bovine serum albumin, ovalbumin) and five basic proteins inhibited
release of the enzyme. Solubilizing activity therefore appears to require a net
negative charge and one or more of the following structural features: extended
conformation, random coil, and unordered or beta-structure, in the latter case,
as the beta-barrel in the fatty acid binding proteins and beta-lactoglobulins.
It is of interest that a difference of a single negative charge between
beta-lactoglobulin A and B, resulted in a statistically significant difference
in the stimulation of hexokinase release. Possible physiological and
pathological roles of this hexokinase solubilizing effect are discussed
briefly.


  Studies on distribution, abundance and diversity of species revealed
fascinating universalities in macroecology. Many of these patterns, like the
species-area and range-abundance relationship or the year-to-year fluctuations
in population sizes are expressed as power-law distributions, and indicate thus
scale-invariance. The species abundance distribution (SAD) apparently shows
this scale-free nature only for rare species, and its mathematical form is much
debated. In the present work we propose a new mathematical expression for SAD
which describes reasonable well most of the presently available large-scale
experimental data and the results of the neutral models. This distribution
function leads to an interesting relation between the total number of
individuals, total number of species and the size of the most abundant species
of the meta-community. This novel scaling relation is confirmed by computer
simulations on neutral models.


  We use the oscillating Min proteins of Escherichia coli as a prototype system
to illustrate the current state and potential of modeling protein dynamics in
space and time. We demonstrate how a theoretical approach has led to striking
new insights into the mechanisms of self-organization in bacterial cells and
indicate how these ideas may be applicable to more complex structure formation
in eukaryotic cells.


  The structure of a genetic network is uncovered by studying its response to
external stimuli (input signals). We present a theory of propagation of an
input signal through a linear stochastic genetic network. It is found that
there are important advantages in using oscillatory signals over step or
impulse signals, and that the system may enter into a pure fluctuation
resonance for a specific input frequency.


  It is widely recognized that the cleaving rate of a restriction enzyme on
target DNA sequences is several orders of magnitude faster than the maximal one
calculated from the diffusion--limited theory. It was therefore commonly
assumed that the target site interaction of a restriction enzyme with DNA has
to occur via two steps: one--dimensional diffusion along a DNA segment, and
long--range jumps coming from association/dissociation events. We propose here
a stochastic model for this reaction which comprises a series of 1D diffusions
of a restriction enzyme on non-specific DNA sequences interrupted by 3D
excursions in the solution until the target sequence is reached. This model
provides an optimal finding strategy which explains the fast association rate.
Modeling the excursions by uncorrelated random jumps, we recover the expression
of the mean time required for target site association to occur given by Berg &
al. \cite{berg81}, and we explicitly give several physical quantities
describing the stochastic pathway of the enzyme. For competitive target sites
we calculate two quantities: processivity and preference. By comparing these
theoretical expressions to recent experimental data obtained for
\textit{Eco}RV--DNA interaction, we quantify: i) the mean residence time per
binding event of \textit{Eco}RV on DNA for a representative 1D diffusion
coefficient, ii) the average lengths of DNA scanned during the 1D diffusion
(during one binding event and during the overall process), iii) the mean time
and the mean number of visits needed to go from one target site to the other.
Further, we evaluate the dynamics of DNA cleavage with regard to the
probability for the restriction enzyme to perform another 1D diffusion on the
same DNA substrate following a 3D excursion.


  High-throughput techniques are leading to an explosive growth in the size of
biological databases and creating the opportunity to revolutionize our
understanding of life and disease. Interpretation of these data remains,
however, a major scientific challenge. Here, we propose a methodology that
enables us to extract and display information contained in complex networks.
Specifically, we demonstrate that one can (i) find functional modules in
complex networks, and (ii) classify nodes into universal roles according to
their pattern of intra- and inter-module connections. The method thus yields a
``cartographic representation'' of complex networks. Metabolic networks are
among the most challenging biological networks and, arguably, the ones with
more potential for immediate applicability. We use our method to analyze the
metabolic networks of twelve organisms from three different super-kingdoms. We
find that, typically, 80% of the nodes are only connected to other nodes within
their respective modules, and that nodes with different roles are affected by
different evolutionary constraints and pressures. Remarkably, we find that
low-degree metabolites that connect different modules are more conserved than
hubs whose links are mostly within a single module.


  We investigate the effect of time-dependent noise on the shape of a morphogen
gradient in a developing embryo. Perturbation theory is used to calculate the
deviations from deterministic behavior in a simple reaction-diffusion model of
robust gradient formation, and the results are confirmed by numerical
simulation. It is shown that such deviations can disrupt robustness for
sufficiently high noise levels, and the implications of these findings for more
complex models of gradient-shaping pathways are discussed.


  This paper deals with a new model for clonal network dynamics. We describe in
detail this model and derive special equations governing immune system dynamics
based on the general gradient type principles that can be inherent to a wide
class of real living objects. A special clonal network is modeled by two
symmetric projector matrix variables simultaneously taking into account both
asymmetry of the interaction to each other and adaptation states that can be
realized owing to possible idiotypic clonal suppresions. We perform computer
simulations of the model dynamics for some simple cases of relatively low
dimension, paying special attention to the dynamics of amounts of activated
receptor strings within clonal network.


  We analyse the effect of harvesting in a resource dependent age structured
population model, deriving the conditions for the existence of a stable steady
state as a function of fertility coefficients, harvesting mortality and
carrying capacity of the resources. Under the effect of proportional harvest,
we give a sufficient condition for a population to extinguish, and we show that
the magnitude of proportional harvest depends on the resources available to the
population. We show that the harvesting yield can be periodic, quasi-periodic
or chaotic, depending on the dynamics of the harvested population. For
populations with large fertility numbers, small harvesting mortality leads to
abrupt extinction, but larger harvesting mortality leads to controlled
population numbers by avoiding over consumption of resources. Harvesting can be
a strategy in order to stabilise periodic or quasi-periodic oscillations in the
number of individuals of a population.


  Identical objects, regularly assembled, form a helix, which is the principal
motif of nucleic acids, proteins, and viral capsids.


  The formation of DNA loops by proteins and protein complexes that bind at
distal DNA sites plays a central role in many cellular processes, such as
transcription, recombination, and replication. Here we review the basic
thermodynamic concepts underlying the assembly of macromolecular complexes on
looped DNA and the effects that this process has in the properties of gene
regulation. Beyond the traditional view of DNA looping as a mechanism to
increase the affinity of regulatory molecules for their cognate sites, recent
developments indicate that DNA looping can also lead to the suppression of
cell-to-cell variability, the control of transcriptional noise, and the
activation of cooperative interactions on demand.


  A-to-I RNA editing by ADARs is a post-transcriptional mechanism for expanding
the proteomic repertoire. Genetic recoding by editing was so far observed for
only a few mammalian RNAs that are predominantly expressed in nervous tissues.
However, as these editing targets fail to explain the broad and severe
phenotypes of ADAR1 knockout mice, additional targets for editing by ADARs were
always expected. Using comparative genomics and expressed sequence analysis, we
identified and experimentally verified four additional candidate human
substrates for ADAR-mediated editing: FLNA, BLCAP, CYFIP2 and IGFBP7.
Additionally, editing of three of these substrates was verified in the mouse
while two of them were validated in chicken. Interestingly, none of these
substrates encodes a receptor protein but two of them are strongly expressed in
the CNS and seem important for proper nervous system function. The editing
pattern observed suggests that some of the affected proteins might have altered
physiological properties leaving the possibility that they can be related to
the phenotypes of ADAR1 knockout mice.


  Fluctuations in diversity and extinction sizes are discussed and compared for
two different, individual-based models of biological coevolution. Both models
display power-law distributions for various quantities of evolutionary
interest, such as the lifetimes of individual species, the quiet periods
between evolutionary upheavals larger than a given cutoff, and the sizes of
extinction events. Time series of the diversity and measures of the size of
extinctions give rise to flicker noise. Surprisingly, the power-law behaviors
of the probability densities of quiet periods in the two models differ, while
the distributions of the lifetimes of individual species are the same.


  Axons of retinal ganglion cells establish orderly projections to the superior
colliculus of the midbrain. Axons of neighboring cells terminate proximally in
the superior colliculus thus forming a topographically precise representation
of the visual world. Coordinate axes are encoded in retina and in the target
through graded expression of chemical labels. Mapping based on chemical labels
alone does not yield required specificity of connections. Additional sharpening
is provided by electrical activity, which is correlated between neighboring
axons. Here we propose a quantitative model which allows combining the effects
of chemical labels and correlated activity in a single approach. Using this
model we study a complete structure of two-dimensional topographic maps in
mutant mice, in which one of the labels, encoding horizontal retinal coordinate
(ephrin-A), is reduced/eliminated. We show that topographic maps in ephrin-A
deficient mice display granular structure, with the regions of smooth mapping
separated by linear discontinuities reminiscent of fractures observed in the
maps of preferred orientation. We also demonstrate the presence of point
singularities in topographic maps.


  The ambitious and ultimate research purpose in Systems Biology is the
understanding and modelling of the cell's system. Although a vast number of
models have been developed in order to extract biological knowledge from
complex systems composed of basic elements as proteins, genes and chemical
compounds, a need remains for improving our understanding of dynamical features
of the systems (i.e., temporal-dependence).
  In this article, we analyze the gene expression dynamics (i.e., how the genes
expression fluctuates in time) by using a new constructive approach. This
approach is based on only two fundamental ingredients: symmetry and the Markov
property of dynamics. First, by using experimental data of human and yeast gene
expression time series, we found a symmetry in short-time transition
probability from time $t$ to time $t+1$. We call it self-similarity symmetry
(i.e., surprisingly, the gene expression short-time fluctuations contain a
repeating pattern of smaller and smaller parts that are like the whole, but
different in size). Secondly, the Markov property of dynamics reflects that the
short-time fluctuation governs the full-time behaviour of the system. Here, we
succeed in reconstructing naturally the global behavior of the observed
distribution of gene expression (i.e., scaling-law) and the local behaviour of
the power-law tail of this distribution, by using only these two ingredients:
symmetry and the Markov property of dynamics. This approach may represent a
step forward toward an integrated image of the basic elements of the whole
cell.


  In this work, the dynamics of fluctuations in gene expression time series is
investigated. By using collected data of gene expression from yeast and human
organisms, we found that the fluctuations of gene expression level and its
average value over time are strongly correlated and obey a scaling law. As this
feature is found in yeast and human organisms, it suggests that probably this
coupling is a common dynamical organizing property of all living systems. To
understand these observations, we propose a stochastic model which can explain
these collective fluctuations, and predict the scaling exponent. Interestingly,
our results indicate that the observed scaling law emerges from the
self-similarity symmetry embedded in gene expression fluctuations.


  In contrast to engineering applications, in which the structure of control
laws are designed to satisfy prescribed function requirements, in biology it is
often necessary to infer gene-circuit function from incomplete data on
gene-circuit structure. By using the feed-forward loop as a model system, this
paper introduces a technique for classifying gene-circuit function given
gene-circuit structure. In simulations performed on a comprehensive set of
models that span a broad range of parameter space, some designs are robust,
producing one unique type of functional response regardless of parameter
selection. Other designs may exhibit a variety of functional responses,
depending upon parameter values. We conclude that, although some feed-forward
loop models have designs that lend themselves to unique function inference,
others have designs for which the function type may be uncertain.


  We study the mechanism underlying the attraction between nucleosomes, the
fundamental packaging units of DNA inside the chromatin complex. We introduce a
simple model of the nucleosome, the eight-tail colloid, consisting of a charged
sphere with eight oppositely charged, flexible, grafted chains that represent
the terminal histone tails. We demonstrate that our complexes are attracted via
the formation of chain bridges and that this attraction can be tuned by
changing the fraction of charged monomers on the tails. This suggests a
physical mechanism of chromatin compaction where the degree of DNA condensation
can be controlled via biochemical means, namely the acetylation and
deacetylation of lysines in the histone tails.


  Protein interaction networks aim to summarize the complex interplay of
proteins in an organism. Early studies suggested that the position of a protein
in the network determines its evolutionary rate but there has been considerable
disagreement as to what extent other factors, such as protein abundance, modify
this reported dependence.
  We compare the genomes of Saccharomyces cerevisiae and Caenorhabditis elegans
with those of closely related species to elucidate the recent evolutionary
history of their respective protein interaction networks. Interaction and
expression data are studied in the light of a detailed phylogenetic analysis.
The underlying network structure is incorporated explicitly into the
statistical analysis.
  The increased phylogenetic resolution, paired with high-quality interaction
data, allows us to resolve the way in which protein interaction network
structure and abundance of proteins affect the evolutionary rate. We find that
expression levels are better predictors of the evolutionary rate than a
protein's connectivity. Detailed analysis of the two organisms also shows that
the evolutionary rates of interacting proteins are not sufficiently similar to
be mutually predictive.
  It appears that meaningful inferences about the evolution of protein
interaction networks require comparative analysis of reasonably closely related
species. The signature of protein evolution is shaped by a protein's abundance
in the organism and its function and the biological process it is involved in.
Its position in the interaction networks and its connectivity may modulate this
but they appear to have only minor influence on a protein's evolutionary rate.


  Information valuable words are the strings with the significant deviation of
real frequency from the expected one. The expected frequency is determined
through the maximum entropy principle of the reconstructed (extended) frequency
dictionary of strings composed from the shorter words. The information valuable
words are found to be the complementary palindromes: they are read equally in
opposite directions, if nucleotides are changed for the complementary ones (A
<--> T; C <--> G) in one of them. Some properties of such symmetric words are
discussed.


  We propose a general method to predict functions of vertices where: 1. The
wiring of the network is somehow related to the vertex functionality. 2. A
fraction of the vertices are functionally classified. The method is influenced
by role-similarity measures of social network analysis. The two versions of our
prediction scheme is tested on model networks were the functions of the
vertices are designed to match their network surroundings. We also apply these
methods to the proteome of the yeast Saccharomyces cerevisiae and find the
results compatible with more specialized methods.


  In this paper we derive a statistical law of Life. It governs the probability
of death, or complementary of survival, of the living organisms. We have
deduced such a law coupling the widely used Weibull statistics, developed for
describing the distribution of the strength of solids, with the universal model
for ontogenetic growth only recently proposed by West and co-authors. The main
idea presented in this paper is that cracks can propagate in solids and cause
their failure as sick cells in living organisms can cause their death. Making a
rough analogy, living organisms are found to behave as growing mechanical
components under cyclic, i.e., fatigue, loadings and composed by a dynamic
evolutionary material that, as an ineluctable fate, deteriorates. The
implications on biological scaling laws are discussed. As an example of
application, we apply such a statistical law to large data collections on human
deaths due to cancer of various types recorded in Italy: a relevant agreement
is observed.


  The effects of monovalent (Na+, K+) and divalent (Mg2+, Ca2+, Mn2+) ions on
the interaction between DNA and histone are studied using the molecular combing
technique. Lamda-DNA molecules and DNA-histone complexes incubated with metal
cations (Na+, K+, Mg2+, Ca2+, Mn2+) are stretched on hydrophobic surfaces, and
directly observed by fluorescence microscopy. The results indicate that when
these cations are added into the DNA solution, the fluorescence intensities of
the stained DNA are reduced differently. The monovalent cations (Na+, K+)
inhibit binding of histone to DNA. The divalent cations (Mg2+, Ca2+, Mn2+)
enhance significantly the binding of histone to DNA and the binding of the
DNA-histone complex to the hydrophobic surface. Mn2+ also induces condensation
and aggregation of the DNA-histone complex.


  This manuscript uses a statistical mechanical approach to study the effect of
the adhesion, through MOCA protein, on cell locomotion. The MOCA protein
regulates cell-cell adhesion, and we explore its potential role in the cell
movement. We present a series of statistical descriptions of the motion in
order to characterize the cell movement, and found that MOCA affects the
statistical scenario of cell locomotion. In particular, we observe that MOCA
enhances the tendency of joint motion, inhibits super-diffusion, and decreases
overall cell motion. These facts are compatible with the hypothesis that the
cells move faster in a less cohesive environment. Furthermore, we observe that
velocity distribution tails are longer than those predicted by
Maxwell-Boltzmann in all cases studied here, indicating that cell movement is
more complex than that of a liquid.


  The assumption that similar structures have similar folding probabilities
($p_{fold}$) leads naturally to a procedure to evaluate $p_{fold}$ for every
snapshot saved along an equilibrium folding-unfolding trajectory of a
structured peptide or protein. The procedure utilizes a structurally
homogeneous clustering and does not require any additional simulation. It can
be used to detect multiple folding pathways as shown for a three-stranded
antiparallel $\beta$-sheet peptide investigated by implicit solvent molecular
dynamics simulations.


  This third part of the lecture series deals with the question: Who will pay
for your retirement? For Western Europe the answer may be ``nobody'', but for
Algeria the demography looks more promising.


  PDZ (Post-synaptic density-95/discs large/zonula occludens-1) domains are
relatively small (80 to 120 residues) protein binding modules central in the
organization of receptor clusters and in the association of cellular proteins.
Their main function is to bind C-terminals of selected proteins that are
recognized through specific amino-acids in their carboxyl end. Binding is
associated with a deformation of the PDZ native structure and is responsible
for dynamical changes in regions not in direct contact with the target. We
investigate how this deformation is related to the harmonic dynamics of the PDZ
structure and show that one low-frequency collective normal mode, characterized
by the concerted movements of different secondary structures, is involved in
the binding process. Our results suggest that even minimal structural changes
are responsible of communication between distant regions of the protein, in
agreement with recent Nuclear Magnetic Resonance (NMR) experiments. Thus PDZ
domains are a very clear example of how collective normal modes are able to
characterize the relation between function and dynamics of proteins, and to
provide indications on the precursors of binding/unbonding events.


  We discuss probabilistic methods for predicting protein functions from
protein-protein interaction networks. Previous work based on Markov Randon
Fields is extended and compared to a general machine-learning theoretic
approach. Using actual protein interaction networks for yeast from the MIPS
database and GO-SLIM function assignments, we compare the predictions of the
different probabilistic methods and of a standard support vector machine. It
turns out that, with the currently available networks, the simple methods based
on counting frequencies perform as well as the more sophisticated approaches.


  Twenty years ago, after analysing palaeontological data, Raup and Sepkoski
suggested that mass extinctions on Earth appear cyclically in time with a
period of approximately 26 million years (My). To explain the 26My period, a
number of proposals were made involving, e.g., astronomical effects, increased
volcanic activity, or the Earth's magnetic field reversal, none of which,
however, has been confirmed. Here we study a spatially extended discrete model
of an ecosystem and show that the periodicity of mass extinctions might be a
natural feature of the ecosystem's dynamics and not the result of a periodic
external perturbation. In our model, periodic changes of the diversity of an
ecosystem and some of its other characteristics are induced by the coevolution
of species. In agreement with some palaeontological data, our results show that
the longevity of a species depends on the evolutionary stage at which the
species is created. Possible further tests of our model are also discussed.


  We study theoretically the effects of confinement on active polar gels such
as the actin network of eukaryotic cells. Using generalized hydrodynamics
equations derived for active gels, we predict, in the case of quasi
one-dimensional geometry, a spontaneous flow transition from a homogeneously
polarized immobile state for small thicknesses, to a perturbed flowing state
for larger thicknesses. The transition is not driven by an external field but
by the activity of the system. We suggest several possible experimental
realizations.


  Human behavior is one of the main problems for evolution, as it is often the
case that human actions are disadvantageous for the self and advantageous for
other people. Behind this puzzle are our beliefs about rational behavior, based
on game theory. Here we show that by going beyond the standard game-theoretical
conventions, apparently altruistic behavior can be understood as
self-interested. We discuss in detail an example related to the so called
Ultimatum game and illustrate the appearance of altruistic behavior induced by
fluctuations. In addition, we claim that in general settings, fluctuations play
a very relevant role, and we support this claim by considering a completely
different example, namely the Stag-Hunt game.


  An analytical approximation is derived for the Zero Sum Multinomial
distribution which gives the Species Abundance Distribution in Neutral
Community Models. The obtained distribution function describes well computer
simulation results on the model, and leads to an interesting relation between
the total number of individuals, total number of species and the size of the
most abundant species of the considered metacommunity. Computer simulations on
neutral community models, proves also the validity of this scaling relation.


  We present a verified computational model of the SH3 domain transition state
(TS) ensemble. This model was built for three separate SH3 domains using
experimental s in all-atom protein folding simulations. While averaging over
all conformations incorrectly considers non-TS conformations as transition
states, quantifying structures as pre-TS, TS, and post-TS by measurement of
their transmission coefficient (pfold, or probability to fold) allows for
rigorous conclusions regarding the structure of the folding nucleus and a full
mechanistic analysis of the folding process. Through analysis of the TS, we
observe a highly polarized nucleus in which many residues are solvent-exposed.
Mechanistic analysis suggests the hydrophobic core forms largely after an early
nucleation step. SH3 presents an ideal system for studying the
nucleation-condensation mechanism and highlights the synergistic relationship
between experiment and simulation in the study of protein folding.


  This series presents an approach to mathematical biology which makes precise
the function of biological molecules. Because biological systems compute, the
theory is a general purpose computer language. I build a language for
efficiently representing the function of protein-like molecules in a cell. The
first paper only presents the kinematic part of the formalism, but this is
already useful for representing large-scale protein networks. The full
formalism allows us to investigate the properties of protein interaction
models, ultimately yielding an estimate of the random-access memory of the
proteins, a measure of their capacity for computation.


  Mechanical unfolding trajectories, generated by applying constant force in
optical tweezer experiments, show that RNA hairpins and the P5abc subdomain of
the group I intron unfold reversibly. We use coarse-grained Go-like models for
RNA hairpins to explore forced-unfolding over a broad range of temperatures. A
number of predictions that are amenable to experimental tests are made. At the
critical force the hairpin jumps between folded and unfolded conformations
without populating any discernible intermediates. The phase diagram in the
force-temperature (f,T) plane shows that the hairpin unfolds by an all-or-none
process. The cooperativity of the unfolding transition increases dramatically
at low temperatures. Free energy of stability, obtained from time averages of
mechanical unfolding trajectories, coincide with ensemble averages which
establishes ergodicity. The hopping time between the the native basin of
attraction (NBA) and the unfolded basin increases dramatically along the phase
boundary. Thermal unfolding is stochastic whereas mechanical unfolding occurs
in "quantized steps" with great variations in the step lengths. Refolding
times, upon force quench, from stretched states to the NBA is "at least an
order of magnitude" greater than folding times by temperature quench. Upon
force quench from stretched states the NBA is reached in at least three stages.
In the initial stages the mean end-to-end distance decreases nearly
continuously and only in the last stage there is a sudden transition to the
NBA. Because of the generality of the results we propose that similar behavior
should be observed in force quench refolding of proteins.


  We introduce a new measure of antigenic distance between influenza A vaccine
and circulating strains. The measure correlates well with efficacies of the
H3N2 influenza A component of the annual vaccine between 1971 and 2004, as do
results of a theory of the immune response to influenza following vaccination.
This new measure of antigenic distance is correlated with vaccine efficacy to a
greater degree than are current state-of-the-art phylogenetic sequence analyzes
or ferret antisera inhibition assays. We suggest that this new measure of
antigenic distance be used in the design of the annual influenza vaccine and in
the interpretation of vaccine efficacy monitoring.


  We look at the problem of signal transduction by extracellular agonist
binding to a receptor protein at the membrane (sensor) via binding of
G-proteins (effectors) to a highly integrative target molecule, such as the
second messenger cAMP (target). We explore the effects of binding times,
effector assignment and effector pool size on the shape of the output signal
under different input scenarios. We conclude that low rates of information
transfer may sometimes coincide with a high probability or efficiency of
plasticity induction.


  Prediction of one-dimensional protein structures such as secondary structures
and contact numbers is useful for the three-dimensional structure prediction
and important for the understanding of sequence-structure relationship. Here we
present a new machine-learning method, critical random networks (CRNs), for
predicting one-dimensional structures, and apply it, with position-specific
scoring matrices, to the prediction of secondary structures (SS), contact
numbers (CN), and residue-wise contact orders (RWCO). The present method
achieves, on average, $Q_3$ accuracy of 77.8% for SS, correlation coefficients
of 0.726 and 0.601 for CN and RWCO, respectively. The accuracy of the SS
prediction is comparable to other state-of-the-art methods, and that of the CN
prediction is a significant improvement over previous methods. We give a
detailed formulation of critical random networks-based prediction scheme, and
examine the context-dependence of prediction accuracies. In order to study the
nonlinear and multi-body effects, we compare the CRNs-based method with a
purely linear method based on position-specific scoring matrices. Although not
superior to the CRNs-based method, the surprisingly good accuracy achieved by
the linear method highlights the difficulty in extracting structural features
of higher order from amino acid sequence beyond that provided by the
position-specific scoring matrices.


  We describe the results obtained from an improved model for protein folding.
We find that a good agreement with the native structure of a 46 residue long,
five-letter protein segment is obtained by carefully tuning the parameters of
the self-avoiding energy. In particular we find an improved free-energy
profile. We also compare the efficiency of the multidimensional replica
exchange method with the widely used parallel tempering.


  We show that spectral fluctuation of interaction matrices of yeast a core
protein interaction network and a metabolic network follows the description of
the Gaussian orthogonal ensemble (GOE) of random matrix theory (RMT).
Furthermore, we demonstrate that while the global biological networks evaluated
belong to GOE, removal of interactions between constituents transitions the
networks to systems of isolated modules described by the Poisson statistics of
RMT. Our results indicate that although biological networks are very different
from other complex systems at the molecular level, they display the same
statistical properties at large scale. The transition point provides a new
objective approach for the identification of functional modules.


  We consider a neutral haploid population whose generations are not
overlapping and whose size is large and constantly of $N$ individuals. Any
generation is replaced by a new one and any individual has a single parent. We
do not choose the stochastic rule which assigns the number of offsprings to any
individual since results do not depend on the details of the dynamics, and, as
a consequence, the model is parameter free. The genealogical tree is very
complex, and distances between individuals (number of generations from the
common ancestor) are distributed according to probability density which remains
random in the thermodynamic limit (large population). We give a theoretical and
numerical description of this distribution and we also consider the dynamical
aspects of the problem describing the time evolution of the maximum and mean
distances in a single population.


  We use single-particle tracking to study the elastic properties of single
microtubules grafted to a substrate. Thermal fluctuations of the free
microtubule's end are recorded, in order to measure position distribution
functions from which we calculate the persistence length of microtubules with
contour lengths between 2.6 and 48 micrometers. We find the persistence length
to vary by more than a factor of 20 over the total range of contour lengths.
Our results support the hypothesis that shearing between protofilaments
contributes significantly to the mechanics of microtubules.


  Understanding design principles of complex cellular organization is one of
the major challenges in biology. Recent analysis of the large-scale cellular
organization has revealed the scale-free nature and robustness of metabolic and
protein networks. However, the underlying evolutional process that creates such
a cellular organization is not fully elucidated. To approach this problem, we
analyzed the metabolic networks of 126 organisms, whose draft or complete
genome sequences have been published. This analysis has revealed that the
evolutional process of metabolic networks follows the same and surprisingly
simple principles in Archaea, Bacteria and Eukaryotes; where highly linked
metabolites change their chemical links more dynamically than less linked
metabolites. Here we demonstrate that this rich-travel-more mechanism rather
than the previously proposed rich-get-richer mechanism can generate the
observed scale-free organization of metabolic networks. These findings
illustrate universal principles in evolution of metabolic networks and suggest
marked flexibility of metabolic network throughout evolution.


  Optimal structure of proteins is described by linear stochastic differential
equation with mean decrease of free energy and volatility. Structure
determining strategy is given by a twin of stochastic variables for which
empirical conditions are not postulated. Optimal structure determination will
be deformed to be adoptive to trading strategy employing martingale property
where stochastic integral w.r.t. analytical solution of stochastic differential
equation will be employed.


  The discovery of two fundamental laws concerning cellular dynamics with
recursive growth is reported. First, the chemical abundances measured over many
cells are found to obey a log-normal distribution and second, the relationship
between the average and standard deviation of the abundances is found to be
linear. The ubiquity of the laws is explored both theoretically and
experimentally. First by means of a model with a catalytic reaction network,
the laws are shown to appear near the critical state with efficient
self-reproduction. Second by measuring distributions of fluorescent proteins in
bacteria cells the ubiquity of log-normal distribution of protein abundances is
confirmed. Relevance of these findings to cellular function and biological
plasticity is briefly discussed.


  This article discusses how the individual morphological properties of basic
objects (e.g. neurons, molecules and aggregates), jointly with their particular
spatial distribution, can determine the connectivity and dynamics of systems
composed by those objects. This problem is characterized as a particular case
of the more general shape and function paradigm, which emphasizes the interplay
between shape and function in nature and evolution. Five key issues are
addressed: (a) how to measure shapes; (b) how to obtain stochastic models of
classes of shapes; (c) how to simulate morphologically realistic systems of
multiple objects; (d) how to characterize the connectivity and topology of such
systems in terms of complex network concepts and measurements; and (e) how the
dynamics of such systems can be ultimately affected, and even determined, by
the individual morphological features of the basic objects. Although emphasis
is placed on neuromorphic systems, the presented concepts and methods are
useful also for several other multiple object systems, such as protein-protein
interaction, tissues, aggregates and polymers.


  We construct a continuum model for biological aggregations in which
individuals experience long-range social attraction and short range dispersal.
For the case of one spatial dimension, we study the steady states analytically
and numerically. There exist strongly nonlinear states with compact support and
steep edges that correspond to localized biological aggregations, or clumps.
These steady state clumps are approached through a dynamic coarsening process.
In the limit of large population size, the clumps approach a constant density
swarm with abrupt edges. We use energy arguments to understand the nonlinear
selection of clump solutions, and to predict the internal density in the large
population limit. The energy result holds in higher dimensions as well, and is
demonstrated via numerical simulations in two dimensions.


  In biology experiments, oligonucleotide microarrays are contacted with a
solution of long nucleic acid (NA) targets. The hybridized probes thus carry
long tails. When the surface density of the oligonucleotide probes is high
enough, the progress of hybridization leads to the formation of a
polyelectrolyte brush due to mutual crowding of the NA tails. The free energy
penalty associated with the brush modifies both the hybridization isotherms and
the rate equations: the attainable hybridization is lowered significantly as is
the hybridization rate. While the equilibrium hybridization fraction, $x_{eq}$,
is low, the hybridization follows a Langmuir type isotherm, $x_{eq}/(1-x_{eq})
= c_t K$ where $c_t$ is the target concentration and $K$ is the equilibrium
constant smaller than its bulk value by a factor $(n/N)^{2/5}$ due to wall
effects where $n$ and $N$ denote the number of bases in the probe and the
target. At higher $x_{eq}$, when the brush is formed, the leading correction is
$x_{eq}/(1-x_{eq}) = c_t K \exp [ - const' (x_{eq}^{2/3} - x_B^{2/3})]$ where
$x_B$ corresponds to the onset of the brush regime. The denaturation rate
constant in the two regimes are identical. However, the hybridization rate
constant in the brush regime is lower, the leading correction being $\exp [-
const' (x^{2/3} - x_B^{2/3})]$.


  We propose to model the dynamics of metabolic networks from a systems biology
point of view by four dynamical structure elements: potential function,
transverse matrix, degradation matrix, and stochastic force. These four
elements are balanced to determine the network dynamics, which gives arise to a
special stochastic differential equation supplemented by a relationship between
the stochastic force and the degradation matrix. Important network behaviors
can be obtained from the potential function without explicitly solving for the
time-dependent solution. The existence of such a potential function suggests a
global optimization principle, and the existence stochastic force corresponds
natural to the hierarchical structure in metabolic networks. We provide
theoretical evidences to justify our proposal by discussing its connections to
others large-scale biochemical systems approaches, such as the network
thermodynamics theory, biochemical systems theory, metabolic control analysis,
and flux balance analysis. Experimental data displaying stochasticity are also
pointed out.


  We study the electronic properties of DNA by way of a tight-binding model
applied to four particular DNA sequences. The charge transfer properties are
presented in terms of localisation lengths, crudely speaking the length over
which electrons travel. Various types of disorder, including random potentials,
are employed to account for different real environments. We have performed
calculations on poly(dG)-poly(dC), telomeric-DNA, random-ATGC DNA and
lambda-DNA. We find that random and lambda-DNA have localisation lengths
allowing for electron motion among a few dozen base pairs only. A novel
enhancement of localisation lengths is observed at particular energies for an
increasing binary backbone disorder. We comment on the possible biological
relevance of sequence dependent charge transfer in DNA.


  A fundamental task in developmental biology is to identify the mechanisms
which drive morphogenesis. In many cases, pattern formation is driven by the
positional information determined by both the gradient of maternal factors and
hard-wired mechanisms embedded in the genome. Alternative mechanisms of
positional information that contribute to patterning are the influence of
signals derived from surrounding tissues. In this paper, we show that the
interplay of geometrical constraints imposed by tissue shapes and hard-wired
mechanisms, computationally implemented by a gene network model, can able to
induce stable complex patterns. The rise of these patterns depends strongly on
the geometrical constraints such as the shape of the surrounding tissues.


  We perform simulations based on the Penna model for biological ageing, now
with the purpose of studying sympatric speciation, that is, the division of a
single species into two or more populations, reproductively isolated, but
without any physical barrier separating them. For that we introduce a new kind
of competition among the individuals, using a modified Verhulst factor. The new
competition depends on some specific phenotypic characteristic of each
individual, which is represented by a pair of bitstrings. These strings are
read in parallel and have no age structure. In this way, each individual genome
consists of two parts. The first one has an age-structure and is related to the
appearance of inherited diseases; the second part is not structured and takes
into account the competition for the available resources. We also introduce
sexual selection into the model, making use of another non-structured and
independent pair of bitstrings.
  In this thesis we present three different models; two of them use, besides
the competition, a sudden change in the ecology to obtain speciation. They were
motivated by the speciation process observed in the Darwin finches, a family of
birds that inhabits the Galapagos Islands, and also by that observed in the
cichlids, a family of fish that lives in the Nicaragua Lakes and in the Vitoria
Lake, in Africa. The third model does not use any ecological change: sympatric
speciation is obtained depending only on the strength of competition among
individuals with similar phenotypic characteristics.


  L-selectin mediated tethers result in leukocyte rolling only above a
threshold in shear. Here we present biophysical modeling based on recently
published data from flow chamber experiments (Dwir et al., J. Cell Biol. 163:
649-659, 2003) which supports the interpretation that L-selectin mediated
tethers below the shear threshold correspond to single L-selectin carbohydrate
bonds dissociating on the time scale of milliseconds, whereas L-selectin
mediated tethers above the shear threshold are stabilized by multiple bonds and
fast rebinding of broken bonds, resulting in tether lifetimes on the timescale
of $10^{-1}$ seconds. Our calculations for cluster dissociation suggest that
the single molecule rebinding rate is of the order of $10^4$ Hz. A similar
estimate results if increased tether dissociation for tail-truncated L-selectin
mutants above the shear threshold is modeled as diffusive escape of single
receptors from the rebinding region due to increased mobility. Using computer
simulations, we show that our model yields first order dissociation kinetics
and exponential dependence of tether dissociation rates on shear stress. Our
results suggest that multiple contacts, cytoskeletal anchorage of L-selectin
and local rebinding of ligand play important roles in L-selectin tether
stabilization and progression of tethers into persistent rolling on endothelial
surfaces.