  The population dynamics and stability of ecosystems of interacting species is
studied from the perspective of non-equilibrium thermodynamics by assuming that
species, through their biotic and abiotic interactions, are units of entropy
production and exchange in an open thermodynamic system with constant external
constraints. Within the context of the linear theory of irreversible
thermodynamics, such a system will naturally evolve towards a stable stationary
state in which the production of entropy within the ecosystem is at a local
minimum value. It is shown that this extremal condition leads to equations for
the stationary (steady) state population dynamics of interacting species, more
general than those of Lotka-Volterra, and to conditions on the parameters of
the community interaction matrix guaranteeing ecosystem stability. The
paradoxical stability of real complex ecosystems thus has a simple explanation
within the proposed framework. Furthermore, it is shown that the second law of
thermodynamics constrains the inter- and intra-species interaction coefficients
in the sense of maintaining stability during evolution from one stationary
state to another. A firm connection is thus established between the second law
of thermodynamics and natural selection.


  A simple model of macroevolution is proposed exhibiting both the property of
punctuated equilibrium and the dynamics of potentialities for different species
to evolve towards increasingly higher complexity. It is based on the phenomenon
of "fractal evolution" which has been shown to constitute a fundamental
property of nonlinear discretized systems with one memory- or random-based
feedback loop. The latter involves a basic "cognitive" function of each species
given by the power of distinction of states within some predefined resolution.
The introduction of a realistic background noise limiting the range of the
feedback operation yields a pattern signature in fitness space with a
distribution of temporal boost/mutation distances according to a randomized
devil's staircase function. Introducing a further level in the hierarchy of the
system's rules, the possibility of an adaptive evolutionary change of the
resolution itself is implemented, thereby providing a time-dependent measure of
the species' cognitive abilities: an additional feedback loop makes use of the
inevitable intrinsic fluctuations within the fitness landscape to direct the
temporal change of the resolution. Feeding back the small adaptive changes of
resolution into the essentially directionless variations of the patterns'
lifetimes in fitness space effectively leads to a clear tendency towards
increasing evolution potentials for each species ("hierarchically emergent
fractal evolution").


  We offer simple solutions to three kinematic problems that occur in the
folding of proteins. We show how to construct suitably local elementary Monte
Carlo moves, how to close a loop, and how to fold a loop without breaking the
bond that closes it.


  Textual analysis of typical microbial genomes reveals that they have the
statistical characteristics of a DNA sequence of a much shorter length. This
peculiar property supports an evolutionary model in which a genome evolves by
random mutation but primarily grows by random segmental self-copying. That
genomes grew mostly by self-copying is consistent with the observation that
repeat sequences in all genomes are widespread and intragenomic and
intergenomic homologous genes are preponderance across all life forms. The
model predicates the coexistence of the two competing modes of evolution: the
gradual changes of classical Darwinism and the stochastic spurts envisioned in
``punctuated equilibrium''.


  Motivated by recent experimental developments in functional genomics, we
construct and test a numerical technique for inferring it process pathways, in
which one process calls another process, from time series data. We validate
using a case in which data are readily available and formulate an extension,
appropriate for genetic regulatory networks, which exploits Bayesian inference
and in which the present--day undersampling is compensated for by prior
understanding of genetic regulation.


  Genome wide comparisons between enteric bacteria yield large sets of
conserved putative regulatory sites on a gene by gene basis that need to be
clustered into regulons. Using the assumption that regulatory sites can be
represented as samples from weight matrices we derive a unique probability
distribution for assignments of sites into clusters. Our algorithm, 'PROCSE'
(probabilistic clustering of sequences), uses Monte-Carlo sampling of this
distribution to partition and align thousands of short DNA sequences into
clusters. The algorithm internally determines the number of clusters from the
data, and assigns significance to the resulting clusters. We place theoretical
limits on the ability of any algorithm to correctly cluster sequences drawn
from weight matrices (WMs) when these WMs are unknown. Our analysis suggests
that the set of all putative sites for a single genome (e.g. E. coli) is
largely inadequate for clustering. When sites from different genomes are
combined and all the homologous sites from the various species are used as a
block, clustering becomes feasible. We predict 50-100 new regulons as well as
many new members of existing regulons, potentially doubling the number of known
regulatory sites in E. coli.


  Gene function prediction from microarray data is a first step toward better
understanding the machinery of the cell from relatively cheap and
easy-to-produce data. In this paper we investigate whether the knowledge of
many metabolic pathways and their catalyzing enzymes accumulated over the years
can help improve the performance of classifiers for this problem.
  The complex network of known biochemical reactions in the cell results in a
representation where genes are nodes of a graph. Formulating the problem as a
graph-driven features extraction problem, based on the simple idea that
relevant features are likely to exhibit correlation with respect to the
topology of the graph, we end up with an algorithm which involves encoding the
network and the set of expression profiles into kernel functions, and
performing a regularized form of canonical correlation analysis in the
corresponding reproducible kernel Hilbert spaces.
  Function prediction experiments for the genes of the yeast S. Cerevisiae
validate this approach by showing a consistent increase in performance when a
state-of-the-art classifier uses the vector of features instead of the original
expression profile to predict the functional class of a gene.


  We present and review Coupled Two Way Clustering, a method designed to mine
gene expression data. The method identifies submatrices of the total expression
matrix, whose clustering analysis reveals partitions of samples (and genes)
into biologically relevant classes. We demonstrate, on data from colon and
breast cancer, that we are able to identify partitions that elude standard
clustering analysis.


  Single molecule experiments on single- and double stranded DNA have sparked a
renewed interest in the force-extension of polymers. The extensible Freely
Jointed Chain (FJC) model is frequently invoked to explain the observed
behavior of single-stranded DNA. We demonstrate that this model does not
satisfactorily describe recent high-force stretching data. We instead propose a
model (the Discrete Persistent Chain, or ``DPC'') that borrows features from
both the FJC and the Wormlike Chain, and show that it resembles the data more
closely. We find that most of the high-force behavior previously attributed to
stretch elasticity is really a feature of the corrected entropic elasticity;
the true stretch compliance of single-stranded DNA is several times smaller
than that found by previous authors. Next we elaborate our model to allow
coexistence of two conformational states of DNA, each with its own stretch and
bend elastic constants. Our model is computationally simple, and gives an
excellent fit through the entire overstretching transition of nicked,
double-stranded DNA. The fit gives the first values for the elastic constants
of the stretched state. In particular we find the effective bend stiffness for
DNA in this state to be about 10 nm*kbt, a value quite different from either
B-form or single-stranded DNA


  We study the space of all compact structures on a two-dimensional square
lattice of size $N=6\times6$. Each structure is mapped onto a vector in
$N$-dimensions according to a hydrophobic model. Previous work has shown that
the designabilities of structures are closely related to the distribution of
the structure vectors in the $N$-dimensional space, with highly designable
structures predominantly found in low density regions. We use principal
component analysis to probe and characterize the distribution of structure
vectors, and find a non-uniform density with a single peak. Interestingly, the
principal axes of this peak are almost aligned with Fourier eigenvectors, and
the corresponding Fourier eigenvalues go to zero continuously at the
wave-number for alternating patterns ($q=\pi$). These observations provide a
stepping stone for an analytic description of the distribution of structural
points, and open the possibility of estimating designabilities of realistic
structures by simply Fourier transforming the hydrophobicities of the
corresponding sequences.


  By using the Jensen-Shannon divergence, genomic DNA can be divided into
compositionally distinct domains through a standard recursive segmentation
procedure. Each domain, while significantly different from its neighbours, may
however share compositional similarity with one or more distant
(non--neighbouring) domains. We thus obtain a coarse--grained description of
the given DNA string in terms of a smaller set of distinct domain labels. This
yields a minimal domain description of a given DNA sequence, significantly
reducing its organizational complexity. This procedure gives a new means of
evaluating genomic complexity as one examines organisms ranging from bacteria
to human. The mosaic organization of DNA sequences could have originated from
the insertion of fragments of one genome (the parasite) inside another (the
host), and we present numerical experiments that are suggestive of this
scenario.


  Plankton is the productive base of aquatic ecosystems and plays a major role
in the global control of atmospheric carbon dioxide. Nevertheless, after
intensive study, the factors that drive its spatial distribution are still far
from being clear. The models proposed so far show very limited agreement with
actual data as many of their results are not consistent with field
observations. Here we show that fluctuations and turbulent diffusion in
standard prey-predator models are able to accurately and consistently explain
plankton field observations at mesoscales (1-100 km). This includes not only
the spatial pattern but also its temporal evolution. We explicitly elucidate
the interplay between physical and biological factors, suggesting that the form
in which small scale biotic fluctuations are transferred to larger scales may
constitute one of the key elements in determining the spatial distribution of
plankton in the sea.


  A wide range of organisms use circadian clocks to keep internal sense of
daily time and regulate their behavior accordingly. Most of these clocks use
intracellular genetic networks based on positive and negative regulatory
elements. The integration of these "circuits" at the cellular level imposes
strong constraints on their functioning and design. Here we study a recently
proposed model [N. Barkai and S. Leibler, Nature, 403:267--268, 2000] that
incorporates just the essential elements found experimentally. We show that
this type of oscillator is driven mainly by two elements: the concentration of
a repressor protein and the dynamics of an activator protein forming an
inactive complex with the repressor. Thus the clock does not need to rely on
mRNA dynamics to oscillate, which makes it especially resistant to
fluctuations. Oscillations can be present even when the time average of the
number of mRNA molecules goes below one. Under some conditions, this oscillator
is not only resistant to but paradoxically also enhanced by the intrinsic
biochemical noise.


  The isochore concept in human genome sequence was challenged in an analysis
by the International Human Genome Sequencing Consortium (IHGSC). We argue here
that a statement in IGHSC analysis concerning the existence of isochore is
incorrect, because it had applied an inappropriate statistical test. To test
the existence of isochores should be equivalent to a test of homogeneity of
windowed GC%. The statistical test applied in the IHGSC's analysis, the
binomial test, is however a test of a sequence being random on the base level.
For testing the existence of isochore, or homogeneity in GC%, we propose to use
another statistical test: the analysis of variance (ANOVA). It can be shown
that DNA sequences that are rejected by binomial test may not be rejected by
the ANOVA test.


  Artificial Life provides the opportunity to study the emergence and evolution
of simple ecosystems in real time. We give an overview of the advantages and
limitations of such an approach, as well as its relation to individual-based
modeling techniques. The Digital Life system Avida is introduced and prospects
for experiments with ab initio evolution (evolution "from scratch"),
maintenance, as well as stability of ecosystems are discussed.


  The probability that an advantageous mutant rises to fixation in a viral
quasispecies is investigated in the framework of multi-type branching
processes. Whether fixation is possible depends on the overall growth rate of
the quasispecies that will form if invasion is successful, rather than on the
individual fitness of the invading mutant. The exact fixation probability can
only be calculated if the fitnesses of all potential members of the invading
quasispecies are known. Quasispecies fixation has two important
characteristics: First, a sequence with negative selection coefficient has a
positive fixation probability as long as it has the potential to grow into a
quasispecies with an overall growth rate that exceeds the one of the
established quasispecies. Second, the fixation probabilities of sequences with
identical fitnesses can nevertheless vary over many orders of magnitudes. Two
approximations for the probability of fixation are introduced. Both
approximations require only partial knowledge about the potential members of
the invading quasispecies. The performance of these two approximations is
compared to the exact fixation probability on a network of RNA sequences with
identical secondary structure.


  We present a new approach for the analysis of genome-wide expression data.
Our method is designed to overcome the limitations of traditional techniques,
when applied to large-scale data. Rather than alloting each gene to a single
cluster, we assign both genes and conditions to context-dependent and
potentially overlapping transcription modules. We provide a rigorous definition
of a transcription module as the object to be retrieved from the expression
data. An efficient algorithm, that searches for the modules encoded in the data
by iteratively refining sets of genes and conditions until they match this
definition, is established. Each iteration involves a linear map, induced by
the normalized expression matrix, followed by the application of a threshold
function. We argue that our method is in fact a generalization of Singular
Value Decomposition, which corresponds to the special case where no threshold
is applied. We show analytically that for noisy expression data our approach
leads to better classification due to the implementation of the threshold. This
result is confirmed by numerical analyses based on in-silico expression data.
We discuss briefly results obtained by applying our algorithm to expression
data from the yeast S. cerevisiae.


  Ownership and trade emerge from anarchy as evolutionary stable strategies. In
these evolutionary game models, ownership status provides an endogenous
asymmetrizing criterion enabling cheaper resolution of property conflicts.


  Poisson representation techniques provide a powerful method for mapping
master equations for birth/death processes -- found in many fields of physics,
chemistry and biology -- into more tractable stochastic differential equations.
However, the usual expansion is not exact in the presence of boundary terms,
which commonly occur when the differential equations are nonlinear. In this
paper, a gauge Poisson technique is introduced that eliminates boundary terms,
to give an exact representation as a weighted rate equation with stochastic
terms. These methods provide novel techniques for calculating and understanding
the effects of number correlations in systems that have a master equation
description. As examples, correlations induced by strong mutations in genetics,
and the astrophysical problem of molecule formation on microscopic grain
surfaces are analyzed. Exact analytic results are obtained that can be compared
with numerical simulations, demonstrating that stochastic gauge techniques can
give exact results where standard Poisson expansions are not able to.


  We formulate and solve a two-state model for the elasticity of nicked,
double-stranded DNA that borrows features from both the Worm Like Chain and the
Bragg--Zimm model. Our model is computationally simple, and gives an excellent
fit to recent experimental data through the entire overstretching transition.
The fit gives the first value for the bending stiffness of the overstretched
state as about 10 nm*kbt, a value quite different from either B-form or
single-stranded DNA.


  In linear elastic theory, membrane deformation energy caused by inclusions
(e.g., membrane-spanning peptides in lipid bilayers) can be rigorously
presented in terms of coupled harmonic oscillators. The general expressions for
the corresponding "effective spring constants" are derived here for quite
general boundary conditions allowing, e.g., for the azimuthal variation of the
contact slope. The case of two insertions is considered in more details.


  Helicase opening of double-stranded nucleic acids may be "active" (the
helicase directly destabilizes the dsNA to promote opening) or "passive" (the
helicase binds ssNA available due to a thermal fluctuation which opens part of
the dsNA). We describe helicase opening of dsNA, based on helicases which bind
single NA strands and move towards the double-stranded region, using a discrete
``hopping'' model. The interaction between the helicase and the junction where
the double strand opens is characterized by an interaction potential. The form
of the potential determines whether the opening is active or passive. We
calculate the rate of passive opening for the helicase PcrA, and show that the
rate increases when the opening is active. Finally, we examine how to choose
the interaction potential to optimize the rate of strand separation. One
important result is our finding that active opening can increase the unwinding
rate by 7 fold compared to passive opening.


  We present a new method to extract distance and orientation dependent
potentials between amino acid side chains using a database of protein
structures and the standard Boltzmann device. The importance of orientation
dependent interactions is first established by computing orientational order
parameters for proteins with alpha-helical and beta-sheet architecture.
Extraction of the anisotropic interactions requires defining local reference
frames for each amino acid that uniquely determine the coordinates of the
neighboring residues. Using the local reference frames and histograms of the
radial and angular correlation functions for a standard set of non-homologue
protein structures, we construct the anisotropic pair potentials. The
performance of the orientation dependent potentials was studied using a large
database of decoy proteins. The results demonstrate that the new distance and
orientation dependent residue-residue potentials present a significantly
improved ability to recognize native folds from a set of native and decoy
protein structures.


  Clarification of the detailed mechanisms involved in the DNA polymorphism is
an important challenge for computational molecular biophysics. This paper
reports about reversible A/B transitions in DNA observed in silico in a
simulated titration experiment by smooth variation of the size of a water drop
around a double helical solute. The estimated range of hydration numbers
corresponding to the B/A transition roughly agrees with experimental data. The
chain length dependence was studied and it appeared that the transition to
A-form is hindered when the fragment becomes shorter than one helical turn.
Dynamics of the A/B transition at low hydration is cooperative and is driven
mainly by medium range electrostatic interactions of counterions sandwiched
between phosphate strands in the major groove. The correspondence of these
computational observations to common experimental conditions of A/B transitions
is discussed.


  We show that textual analysis of microbial genomes reveal telling footprints
of the early evolution of the genomes. The frequencies of word occurrence of
random DNA sequences considered as texts in their four nucleotides are expected
to obey Poisson distributions. It is noticed that for words less than nine
letters the average width of the distributions for complete microbial genomes
is many times that of a Poisson distribution. We interpret this phenomenon as
follows: the genome is a large system that possesses the statistical
characteristics of a much smaller ``random'' system, and certain textual
statistical properties of genomes we now see are remnants of those of their
ancestral genomes, which were much shorter than the genomes are now. This
interpretation suggests a simple biologically plausible model for the growth of
genomes: the genome first grows randomly to an initial length of approximately
one thousand nucleotides (1k nt), or about one thousandth of its final length,
thereafter mainly grows by random segmental duplication. We show that using
duplicated segments averaging around 25 nt, the model sequences generated
possess statistical properties characteristic of present day genomes. Both the
initial length and the duplicated segment length support an RNA world at the
time duplication began. Random segmental duplication would greatly enhance the
ability of a genome to use its hard-to-acquire codes repeatedly, and a genome
that practiced it would have evolved enormously faster than those that did not.


  We address the question of finding the community structure of a complex
network. In an earlier effort [H. Zhou, {\em Phys. Rev. E} (2003)], the concept
of network random walking is introduced and a distance measure defined. Here we
calculate, based on this distance measure, the dissimilarity index between
nearest-neighboring vertices of a network and design an algorithm to partition
these vertices into communities that are hierarchically organized. Each
community is characterized by an upper and a lower dissimilarity threshold. The
algorithm is applied to several artificial and real-world networks, and
excellent results are obtained. In the case of artificially generated random
modular networks, this method outperforms the algorithm based on the concept of
edge betweenness centrality. For yeast's protein-protein interaction network,
we are able to identify many clusters that have well defined biological
functions.


  In competition experiments between two RNA viruses of equal or almost equal
fitness, often both strains gain in fitness before one eventually excludes the
other. This observation has been linked to the Red Queen effect, which
describes a situation in which organisms have to constantly adapt just to keep
their status quo. I carried out experiments with digital organisms
(self-replicating computer programs) in order to clarify how the competing
strains' location in fitness space influences the Red-Queen effect. I found
that gains in fitness during competition were prevalent for organisms that were
taken from the base of a fitness peak, but absent or rare for organisms that
were taken from the top of a peak or from a considerable distance away from the
nearest peak. In the latter two cases, either neutral drift and loss of the
fittest mutants or the waiting time to the first beneficial mutation were more
important factors. Moreover, I found that the Red-Queen dynamic in general led
to faster exclusion than the other two mechanisms.


  Background: The rate at which fitness declines as an organism's genome
accumulates random mutations is an important variable in several evolutionary
theories. At an intuitive level, it might seem natural that random mutations
should tend to interact synergistically, such that the rate of mean fitness
decline accelerates as the number of random mutations is increased. However, in
a number of recent studies, a prevalence of antagonistic epistasis (the
tendency of multiple mutations to have a mitigating rather than reinforcing
effect) has been observed.
  Results: We studied in silico the net amount and form of epistatic
interactions in RNA secondary structure folding by measuring the fraction of
neutral mutants as a function of mutational distance d. We found a clear
prevalence of antagonistic epistasis in RNA secondary structure folding. By
relating the fraction of neutral mutants at distance d to the average
neutrality at distance d, we showed that this prevalence derives from the
existence of many compensatory mutations at larger mutational distances.
  Conclusions: Our findings imply that the average direction of epistasis in
simple fitness landscapes is directly related to the density with which fitness
peaks are distributed in these landscapes.


  Understanding functional associations among genes discovered in sequencing
projects is a key issue in post-genomic biology. However, reliable
interpretation of the protein interaction data has been difficult. In this
work, we show that if two proteins share significantly larger number of common
interaction partners than random, they have close functional associations.
Analysis of publicly available data from Saccharomyces cerevisiae reveals more
than 2800 reliable functional associations, 29% of which involve at least one
unannotated protein. By further analyzing these associations, we derive
tentative functions for 81 unannotated proteins with high certainty.


  We study the equilibrium behaviour of a deterministic four-state
mutation-selection model as a model for the evolution of a population of
nucleotide sequences in sequence space. The mutation model is the Kimura 3ST
mutation scheme, and the selection scheme is assumed to be invariant under
permutation of sites. Considering the evolution process both forward and
backward in time, we use the ancestral distribution as the stationary state of
the backward process to derive an expression for the mutational loss (as the
difference between ancestral and population mean fitness), and we prove a
maximum principle that determines the population mean fitness in
mutation-selection balance.


  We investigate a multi-locus evolutionary model which is based on the DNA
shuffling protocol widely applied in \textit{in vitro} directed evolution. This
model incorporates selection, recombination and point mutations. The simplicity
of the model allows us to obtain a full analytical treatment of both its
dynamical and equilibrium properties, for the case of an infinite population.
We also briefly discuss finite population size corrections.


  We study the competition between several advantageous mutants in an asexual
population (clonal interference) as a function of the time between the
appearance of the mutants, their selective advantages, and the rate of
deleterious mutations. We find that the overall probability of fixation (the
probability that at least one of the mutants becomes the ancestor of the entire
population) does not depend on the time interval between the appearance of
these mutants, and equals the probability that a genotype bearing all of these
mutations reaches fixation. This result holds also in the presence of
deleterious mutations, and for an arbitrary number of competing mutants. We
also show that if mutations interfere, an increase in the mean number of
fixation events is associated with a decrease in the expected fitness gain of
the population.


  Recent theoretical approaches to community structure and dynamics reveal that
many large-scale features of community structure (such as species-rank
distributions and species-area relations) can be explained by a so-called
neutral model. Using this approach, species are taken to be equivalent and
trophic relations are not taken into account explicitly. Here we provide a
general analytic solution to the local community model of Hubbell's neutral
theory of biodiversity by recasting it as an urn model i.e.a Markovian
description of states and their transitions. Both stationary and time-dependent
distributions are analysed. The stationary distribution -- also called the
zero-sum multinomial -- is given in closed form. An approximate form for the
time-dependence is obtained by using an expansion of the master equation. The
temporal evolution of the approximate distribution is shown to be a good
representation for the true temporal evolution for a large range of parameter
values.


  The ATP-driven Plasma Membrane Calcium pump or Ca(2+)-ATPase (PMCA) is
characterized by a high affinity to calcium and a low transport rate compared
to other transmembrane calcium transport proteins. It plays a crucial role for
calcium extrusion from cells. Calmodulin is an intracellular calcium buffering
protein which is capable in its Ca(2+) liganded form of stimulating the PMCA by
increasing both the affinity to calcium and the maximum calcium transport rate.
We introduce a new model of this stimulation process and derive analytical
expressions for experimental observables in order to determine the model
parameters on the basis of specific experiments. We furthermore develop a model
for the pumping activity. The pumping description resolves the seeming
contradiction of the Ca(2+):ATP stoichiometry of 1:1 during a translocation
step and the observation that the pump binds two calcium ions at the
intracellular site. The combination of the calcium pumping and the stimulation
model correctly describes PMCA function. We find that the processes of
calmodulin-calcium complex attachment to the pump and of stimulation have to be
separated. Other PMCA properties are discussed in the framework of the model.
The presented model can serve as a tool for calcium dynamics simulations and
provides the possibility to characterize different pump isoforms by different
type-specific parameter sets.


  Loop formation in long molecules occurs many places in nature, from solutions
of carbon nanotubes to polymers inside a cell. In this article, we review
theoretical studies of the static and dynamic properties of polymer loops. For
example, long polymers must search many configurations to find a "target"
binding site, while short polymers are stiff and resist bending. In between,
there is an optimal loop size, which balances the entropy of long loops against
the energetic cost of short loops. We show that such simple pictures of loop
formation can explain several long-standing observations in DNA replication,
quantitatively.


  Background. In a number of recent experiments with food-and-mouth disease
virus, a deleterious mutant, was found to avoid extinction and remain in the
population for long periods of time. This observation was called quasispecies
memory. The origin of quasispecies memory is not fully understood.
  Results. We propose and analyze a simple model of complementation between the
wild type virus and a mutant that has an impaired ability of cell entry. The
mutant will go extinct unless it is recreated from the wild type through
mutations. However, under phenotypic mixing-and-hiding as a mechanism of
complementation, the time to extinction in the absence of mutations increases
with increasing multiplicity of infection (m.o.i.). The mutant's frequency at
equilibrium under selection-mutation balance also increases with increasing
m.o.i. At high m.o.i., a large fraction of mutant genomes are encapsidated with
wild-type protein, which enables them to infect cells as efficiently as the
wild type virions, and thus increases their fitness to the wild-type level.
Moreover, even at low m.o.i. the equilibrium frequency of the mutant is higher
than predicted by the standard quasispecies model, because a fraction of mutant
virions generated from wild-type parents will also be encapsidated by wild-type
protein.
  Conclusions. Our model predicts that phenotypic hiding will strongly
influence the population dynamics of viruses, particularly at high m.o.i., and
will also have important effects on the mutation--selection balance at low
m.o.i. The delay in mutant extinction and increase in mutant frequencies at
equilibrium may, at least in part, explain memory in quasispecies populations.


  With the number of sequenced genomes now over one hundred, and the
availability of rough functional annotations for a substantial proportion of
their genes, it has become possible to study the statistics of gene content
across genomes. Here I show that, for many high-level functional categories,
the number of genes in the category scales as a power-law in the total number
of genes in the genome. The occurrence of such scaling laws can be explained
with a simple theoretical model, and this model suggests that the exponents of
the observed scaling laws correspond to universal constants of the evolutionary
process. I discuss some consequences of these scaling laws for our
understanding of organism design.


  We extend our previous stochastic cellular automata based model for areal
aggregation of prion proteins on neuronal surfaces. The new anisotropic model
allow us to simulate both strong beta-sheet and weaker attachment bonds between
proteins. Constraining binding directions allows us to generate aggregate
structures with the hexagonal lattice symmetry found in recently observed in
vitro experiments. We argue that these constraints on rules may correspond to
underlying steric constraints on the aggregation process. We find that monomer
dominated growth of the areal aggregate is too slow to account for some
observed doubling time-to-incubation time ratios inferred from data, and so
consider aggregation dominated by relatively stable but non-infectious
oligomeric intermediates. We compare a kinetic theory analysis of oligomeric
aggregation to spatially explicit simulations of the process. We find that with
suitable rules for misfolding of oligomers, possibly due to water exclusion by
the surrounding aggregate, the resulting oligomeric aggregation model maps onto
our previous monomer aggregation model. Therefore it can produce some of the
same attractive features for the description of prion incubation time data. We
propose experiments to test the oligomeric aggregation model.


  Simple hidden Markov models are proposed for predicting secondary structure
of a protein from its amino acid sequence. Since the length of protein
conformation segments varies in a narrow range, we ignore the duration effect
of length distribution, and focus on inclusion of short range correlations of
residues and of conformation states in the models. Conformation-independent and
-dependent amino acid coarse-graining schemes are designed for the models by
means of proper mutual information. We compare models of different level of
complexity, and establish a practical model with a high prediction accuracy.


  Inspired by recent suggestions that the Alzheimer's amyloid beta peptide (A
beta) can insert into cell membranes and form harmful ion channels, we model
insertion of the 40 and 42 residue forms of the peptide into cell membranes
using a Monte Carlo code which is specific at the amino acid level. We examine
insertion of the regular A-beta peptide as well as mutants causing familial
Alzheimer's disease, and find that all but one of the mutants change the
insertion behavior by causing the peptide to spend more simulation steps in
only one leaflet of the bilayer. We also find that A-beta 42, because of the
extra hydrophobic residues relative to A-beta 40, is more likely to adopt this
conformation than A-beta 40 in both wild-type and mutant forms. We argue
qualitatively why these effects happen. Here, we present our results and
develop the hypothesis that this partial insertion increases the probability of
harmful channel formation. This hypothesis can partly explain why these
mutations are neurotoxic simply due to peptide insertion behavior. We further
apply this model to various artificial A-beta mutants which have been examined
experimentally, and offer testable experimental predictions contrasting the
roles of aggregation and insertion with regard to toxicity of A-beta mutants.
These can be used through further experiments to test our hypothesis.


  Using a cellular automata model, we simulate the British Government Policy
(BGP) in the 2001 foot and mouth epidemic in Great Britain. When clinical
symptoms of the disease appeared on a farm, there is mandatory slaughter
(culling) of all livestock on an infected premise (IP). Those farms that
neighbor an IP (contiguous premise, CP), are also culled, aka nearest neighbor
interaction. Farms where the disease may be prevalent from animal, human,
vehicle or airborne transmission (dangerous contact, DC), are additionally
culled, aka next-to-nearest neighbor iteractions and lightning factor. The
resulting mathematical model possesses a phase transition, whereupon if the
physical disease transmission kernel exceeds a critical value, catastrophic
loss of animals ensues. The non-local disease transport probability can be as
low as .01% per day and the disease can still be in the high mortality phase.
We show that the fundamental equation for sustainable disease transport is the
criticality equation for neutron fission cascade. Finally, we calculate that
the percentage of culled animals that are actually healthy is ~ 30%.


  A new suggestion for pattern formation in biological systems is proposed. The
interaction of two signaling pathways provides complex patterning. The
interaction of two pathways involves few of the biochemical complexities of the
pathways. The 'Wnt' pathway is proposed to be one of these, in interaction with
others, one at a time, perhaps sequentially, such as Hunchback (Hh), Notch,
BMP, etc., those involved in development of metazoa. It is further shown how
such interactions can easily lead to such epithelial movements as gastrulation.
Numerical simulations of the model are given.


  We studied the adsorption of a charged protein onto an oppositely charged
membrane, composed of mobile phospholipids of differing valence, using a
statistical-thermodynamical approach. A two-block model was employed, one block
corresponding to the protein-affected region on the membrane, referred to as
the adsorption domain, and the other to the unaffected remainder of the
membrane. We calculated the protein-induced lipid rearrangement in the
adsorption domain as arising from the interplay between the electrostatic
interactions in the system and the mixing entropy of the lipids. Equating the
electrochemical potentials of the lipids in the two blocks yields an expression
for the relations among the various lipid fractions in the adsorption domain,
indicating a sensitive dependence of lipid fraction on valence. This expression
is a result of the two-block picture but does not depend on further details of
the protein-membrane interaction. We subsequently calculated the lipid
fractions themselves using the Poisson-Boltzmann theory. We examined the
dependence of lipid enrichment, i.e., the ratio between the lipid fractions
inside and outside the adsorption domain, on various parameters such as ionic
strength and lipid valence. Maximum enrichment was found for lipid valence of
about (-3) to (-4) in physiological conditions. Our results are in qualitative
agreement with recent experimental studies on the interactions between peptides
having a domain of basic residues and membranes containing a small fraction of
the polyvalent phosphatidylinositol 4,5-bisphosphate (PIP2). This study
provides theoretical support for the suggestion that proteins adsorbed onto
membranes through a cluster of basic residues may sequester PIP2 and other
polyvalent lipids.


  For a substance diffusing on a curved surface, we obtain an explicit relation
valid for very small values of the time, between the local concentration, the
diffusion coefficient, the intrinsic spatial curvature and the time. We recover
the known solution of Fick's law of diffusion in the flat space limit. In the
biological context, this result would be useful in understanding the variations
in the diffusion rates of integral proteins and other molecules on membranes.


  Group theoretical concepts are invoked in a specific model to explain how
only twenty amino acids occur in nature out of a possible sixty four. The
methods we use enable us to justify the occurrence of the recently discovered
twenty first amino acid selenocysteine, and also enables us to predict the
possible existence of two more, as yet undiscovered amino acids.


  Statistical analysis of distributions of occurrence frequencies of short
words in 108 microbial complete genomes reveals the existence of a set of
universal "root-sequence lengths" shared by all microbial genomes. These
lengths and their universality give powerful clues to the way microbial genomes
are grown. We show that the observed genomic properties are explained by a
model for genome growth in which primitive genomes grew mainly by maximally
stochastic duplications of short segments from an initial length of about 200
nucleotides (nt) to a length of about one million nt typical of microbial
genomes. The relevance of the result of this study to the nature of
simultaneous random growth and information acquisition by genomes, to the
so-called RNA world in which life evolved before the rise of proteins and
enzymes and to several other topics are discussed.


  Although the auger-like 'swimming' motility of the African trypanosome was
described upon its discovery over one hundred years ago, the precise
biomechanical and biophysical properties of trypanosome flagellar motion has
not been elucidated. In this study, we describe five different modes of
flagellar beat/wave patterns in African trypanosomes by microscopically
examining the flagellar movements of chemically tethered cells. The dynamic
nature of the different beat/wave patterns suggests that flagellar motion in
Trypanosoma brucei is a complex mixture of oscillating waves, rigid bends,
helical twists and non-linear waves. Interestingly, we have observed
soliton-like depression waves along the flagellar membrane, suggesting a
nonlinear mechanism for the dynamics of this system. The physical model is
inspired by the 2-dimensional elastic dynamics of a beam, and by taking into
account uniform distribution of molecular motors torque and nonlinear terms in
the curvature.


  RNA secondary structures of increasing complexity are probed combining single
molecule stretching experiments and stochastic unfolding/refolding simulations.
We find that force-induced unfolding pathways cannot usually be interpretated
by solely invoking successive openings of native helices. Indeed, typical
force-extension responses of complex RNA molecules are largely shaped by
stretching-induced, long-lived intermediates including non-native helices. This
is first shown for a set of generic structural motifs found in larger RNA
structures, and then for Escherichia coli's 1540-base long 16S ribosomal RNA,
which exhibits a surprisingly well-structured and reproducible unfolding
pathway under mechanical stretching. Using out-of-equilibrium stochastic
simulations, we demonstrate that these experimental results reflect the slow
relaxation of RNA structural rearrangements. Hence, micromanipulations of
single RNA molecules probe both their native structures and long-lived
intermediates, so-called "kinetic traps", thereby capturing -at the single
molecular level- the hallmark of RNA folding/unfolding dynamics.


  Ab initio RNA secondary structure predictions have long dismissed helices
interior to loops, so-called pseudoknots, despite their structural importance.
Here, we report that many pseudoknots can be predicted through long time scales
RNA folding simulations, which follow the stochastic closing and opening of
individual RNA helices. The numerical efficacy of these stochastic simulations
relies on an O(n^2) clustering algorithm which computes time averages over a
continously updated set of n reference structures. Applying this exact
stochastic clustering approach, we typically obtain a 5- to 100-fold simulation
speed-up for RNA sequences up to 400 bases, while the effective acceleration
can be as high as 100,000-fold for short multistable molecules (<150 bases). We
performed extensive folding statistics on random and natural RNA sequences, and
found that pseudoknots are unevenly distributed amongst RNAstructures and
account for up to 30% of base pairs in G+C rich RNA sequences (Online RNA
folding kinetics server including pseudoknots : http://kinefold.u-strasbg.fr/
).


  We discuss the possibility of controlling biological systems, by exciting in
the near infrared region {\it hybrid} metallic nanotube ropes, dressed with
proteins and embedded in the biosystems. If one nanotube, in a double-tube
rope, is filled with metallofullerenes and the other is empty, the two tubes
change their opposite equilibrium charging during the irradiation. The
resulting change of the local electric field can deform proteins attached to
the tubes, and change their catalytic properties.


  We investigate the effects of spatial discreteness of molecules in
reaction-diffusion systems. It is found that discreteness within the so called
Kuramoto length can lead to a localization of molecules, resulting in novel
steady states that do not exist in the continuous case. These novel states are
analyzed theoretically as the fixed points of accelerated localized reactions,
an approach that was verified to be in good agreement with stochastic particle
simulations. The relevance of this discreteness-induced state to biological
intracellular processes is discussed.


  A nonlinear autocatalysis of a chiral substance is shown to achieve
homochirality in a closed system, if the back-reaction is included. Asymmetry
in the concentration of two enantiomers or the enantiometric excess increases
due to the nonlinear autocatalysis. Furthermore, when the back-reaction is
taken into account, the reactant supplied by the decomposition of the
enantiomers is recycled to produce more and more the dominant one, and
eventually the homochirality is established.


  We report the reproducible first-principles folding of the 40 amino acid,
three-helix headpiece of the HIV accessory protein in a recently developed
all-atom free-energy forcefield. Six of twenty simulations using an adapted
basin-hopping method converged to better than 3 \AA backbone RMS deviation to
the experimental structure. Using over 60,000 low-energy conformations of this
protein, we constructed a decoy tree that completely characterizes its folding
funnel.


  By defining the phenotype of a biopolymer by its active three-dimensional
shape, and its genotype by its primary sequence, we propose a model that
predicts and characterizes the statistical distribution of a population of
biopolymers with a specific phenotype, that originated from a given genotypic
sequence by a single mutational event. Depending on the ratio g0 that
characterizes the spread of potential energies of the mutated population with
respect to temperature, three different statistical regimes have been
identified. We suggest that biopolymers found in nature are in a critical
regime with g0 in the range 1-6, corresponding to a broad, but not too broad,
phenotypic distribution resembling a truncated Levy flight. Thus the biopolymer
phenotype can be considerably modified in just a few mutations. The proposed
model is in good agreement with the experimental distribution of activities
determined for a population of single mutants of a group I ribozyme.


  We explore the consequences of very high dimensionality in the dynamical
landscape of protein folding. Consideration of both typical range of
stabilising interactions, and folding rates themselves, leads to a model of the
energy hypersurface that is characterised by the structure of diffusive
"hypergutters" as well as the familiar "funnels". Several general predictions
result: (1) intermediate subspaces of configurations will always be visited;
(2) specific but non-native interactions are important in stabilising these
low-dimensional diffusive searches on the folding pathway; (3) sequential
barriers will commonly be found, even in "two-state"proteins; (4) very early
times will show charactreristic departures from single-exponential kinetics;
(5) contributions of non-native interactions to phi-values are calculable, and
may be significant. The example of a three-helix bundle is treated in more
detail as an illustration. The model also shows that high-dimensional
structures provide conceptual relations between the "folding funnel",
"diffusion-collision", "nucleation-condensation" and "topomer search" models of
protein folding. It suggests that kinetic strategies for fast folding may be
encoded rather generally in non-native, rather than native interactions. The
predictions are related to very recent findings in experiment and simulation.


  We present a general expression for the association rate for partially
diffusion-controlled reactions between spherical molecules with an asymmetric
reactive patch on each surface. Reaction can occur only if the two patches are
in contact and properly aligned to within specified angular tolerances. This
extends and generalizes previous approaches that considered only axially
symmetric patches; the earlier solutions are shown to be limiting cases of our
general expression. Previous numerical results on the rate of protein-protein
association with high steric specificity are in very good agreement with the
value computed from our analytic expression. Using the new expression, we
investigate the influence of orientational constraints on the rate constant. We
find that for angular constraints of ~5-15 degrees, a typical range for example
in the case of protein-protein interactions, the reaction rate is about 2 to 3
orders of magnitude higher than expected from a simple geometric model.


  A common mechanism for intracellular transport is the use of controlled
deformations of the membrane to create spherical or tubular buds. While the
basic physical properties of homogeneous membranes are relatively well-known,
the effects of inhomogeneities within membranes are very much an active field
of study. Membrane domains enriched in certain lipids in particular are
attracting much attention, and in this Letter we investigate the effect of such
domains on the shape and fate of membrane tubes. Recent experiments have
demonstrated that forced lipid phase separation can trigger tube fission, and
we demonstrate how this can be understood purely from the difference in elastic
constants between the domains. Moreover, the proposed model predicts timescales
for fission that agree well with experimental findings.


  DNA migration in topologically structured microchannels with periodic
cavities is investigated experimentally and with Brownian dynamics simulations
of a simple bead-spring model. The results are in very good agreement with one
another. In particular, the experimentally observed migration order of Lambda-
and T2-DNA molecules is reproduced by the simulations. The simulation data
indicate that the mobility may depend on the chain length in a nonmonotonic way
at high electric fields. This is found to be the signature of a nonequilibrium
phase transition between two different migration states, a slow one and a fast
one, which can also be observed experimentally under appropriate conditions.


  We develop a general theory for active viscoelastic materials made of polar
filaments. This theory is motivated by the dynamics of the cytoskeleton. The
continuous consumption of a fuel generates a non equilibrium state
characterized by the generation of flows and stresses. Our theory can be
applied to experiments in which cytoskeletal patterns are set in motion by
active processes such as those which are at work in cells.


  An anomalous diffusion model for ion channel gating is put forward. This
scheme is able to describe non-exponential, power-law like distributions of
residence time intervals in several types of ion channels. Our method presents
a generalization of the discrete diffusion model by Millhauser, Salpeter and
Oswald [Proc. Natl. Acad. Sci. USA 85, 1503 (1988)] to the case of a
continuous, anomalous slow conformational diffusion. The corresponding
generalization is derived from a continuous time random walk composed of
nearest neighbor jumps which in the scaling limit results in a fractional
diffusion equation. The studied model contains three parameters only: the mean
residence time, a characteristic time of conformational diffusion, and the
index of subdiffusion. A tractable analytical expression for the characteristic
function of the residence time distribution is obtained. In the limiting case
of normal diffusion, our prior findings [Proc. Natl. Acad. Sci. USA 99, 3552
(2002)] are reproduced. Depending on the chosen parameters, the fractional
diffusion model exhibits a very rich behavior of the residence time
distribution with different characteristic time-regimes. Moreover, the
corresponding autocorrelation function of conductance fluctuations displays
nontrivial features. Our theoretical model is in good agreement with
experimental data for large conductance potassium ion channels.


  Recently the first chemical systems which show the amplification of
enantiometric excess (ee) was found. Inspired by these experiments, we propose
a few chemical reaction models in a closed system. The reactions consist of
autocatalytic production of chiral molecules and their decomposition processes.
By means of flow diagrams in the phase space of enatiometric concentrations,
chiral symmetry breaking is studied.
  Without any destruction processes, the quadratic autocatalysis leads to the
amplification of ee, but the flow is attracted to a fixed line, which is only
marginally stable. The destruction processes alter the attractor structure
completely. With a cross inhibition the linear autocatalysis is sufficient for
the chiral symmetry breaking with a complete homochirality. With a back
reaction, the quadratic autocatalysis is necessary to bring about the chiral
symmetry breaking. The final ee is not complete if there is a racemization
process.


  A self-consistent equation to derive a discreteness-induced stochastic steady
state is presented for reaction-diffusion systems. For this formalism, we use
the so-called Kuramoto length, a typical distance over which a molecule
diffuses in its lifetime, as was originally introduced to determine if local
fluctuations influence globally the whole system. We show that this Kuramoto
length is also relevant to determine whether the discreteness of molecules is
significant or not. If the number of molecules of a certain species within the
Kuramoto length is small and discrete, localization of some other chemicals is
brought about, which can accelerate certain reactions. When this acceleration
influences the concentration of the original molecule species, it is shown that
a novel, stochastic steady state is induced that does not appear in the
continuum limit. A theory to obtain and characterize this state is introduced,
based on the self-consistent equation for chemical concentrations. This
stochastic steady state is confirmed by numerical simulations on a certain
reaction model, which agrees well with the theoretical estimation. Formation
and coexistence of domains with different stochastic states are also reported,
which is maintained by the discreteness. Relevance of our result to
intracellular reactions is briefly discussed.


  We have used kinetic Monte Carlo simulations to study the kinetics of
unfolding of cross-linked polymer chains under mechanical loading. As the ends
of a chain are pulled apart, the force transmitted by each crosslink increases
until it ruptures. The stochastic crosslink rupture process is assumed to be
governed by first order kinetics with a rate that depends exponentially on the
transmitted force. We have performed random searches to identify optimal
crosslink configurations whose unfolding requires a large applied force
(measure of strength) and/or large dissipated energy (measure of toughness). We
found that such optimal chains always involve cross-links arranged to form
parallel strands. The location of those optimal strands generally depends on
the loading rate. Optimal chains with a small number of cross-links were found
to be almost as strong and tough as optimal chains with a large number of
cross-links. Furthermore, optimality of chains with a small number of
cross-links can be easily destroyed by adding cross-links at random. The
present findings are relevant for the interpretation of single molecule force
probe spectroscopy studies of the mechanical unfolding of load-bearing
proteins, whose native topology often involves parallel strand arrangements
similar to the optimal configurations identified in the study.


  Allometric scaling is one of the most pervasive laws in biology. Its origin,
however, is still a matter of dispute. Recent studies have established that
maximum metabolic rate scales with an exponent larger than that found for basal
metabolism. This unpredicted result sets a challenge that can decide which of
the concurrent hypotheses is the correct theory. Here we show that both scaling
laws can be deduced from a single network model. Besides the 3/4-law for basal
metabolism, the model predicts that maximum metabolic rate scales as $M^{6/7}$,
maximum heart rate as $M^{-1/7}$, and muscular capillary density as $M^{-1/7}$,
in agreement with data.


  Nonlinear control techniques by means of a software sensor that are commonly
used in chemical engineering could be also applied to genetic regulation
processes. We provide here a realistic formulation of this procedure by
introducing an additive white Gaussian noise, which is usually found in
experimental data. Besides, we include model errors, meaning that we assume we
do not know the nonlinear regulation function of the process. In order to
illustrate this procedure, we employ the Goodwin dynamics of the concentrations
[B.C. Goodwin, Temporal Oscillations in Cells, (Academic Press, New York,
1963)] in the simple form recently applied to single gene systems and some
operon cases [H. De Jong, J. Comp. Biol. 9, 67 (2002)], which involves the
dynamics of the mRNA, given protein, and metabolite concentrations. Further, we
present results for a three gene case in co-regulated sets of transcription
units as they occur in prokaryotes. However, instead of considering their full
dynamics, we use only the data of the metabolites and a designed software
sensor. We also show, more generally, that it is possible to rebuild the
complete set of nonmeasured concentrations despite the uncertainties in the
regulation function or, even more, in the case of not knowing the mRNA
dynamics. In addition, the rebuilding of concentrations is not affected by the
perturbation due to the additive white Gaussian noise and also we managed to
filter the noisy output of the biological system


  Channel proteins, that selectively conduct molecules across cell membranes,
often exhibit an asymmetric structure. By means of a stochastic model, we argue
that channel asymmetry in the presence of non-equilibrium fluctuations, fueled
by the cell's metabolism as observed recently, can dramatically influence the
transport through such channels by a ratchet-like mechanism. For an
aquaglyceroporin that conducts water and glycerol we show that a previously
determined asymmetric glycerol potential leads to enhanced inward transport of
glycerol, but for unfavorably high glycerol concentrations also to enhanced
outward transport that protects a cell against poisoning.


  A scaling analysis within a model of hierarchically constrained dynamics is
shown to reproduce the main features of non-exponential relaxation observed in
kinetic studies of carbonmonoxymyoglobin.


  We present a method to search low energy configurations of polar molecules in
the complex potential energy surfaces associated with dense fluids. The search
is done in the configurational space of the translational and rotational
degrees of freedom of the molecule, combining steepest-descent and
Newton-Raphson steps which embed information on the average sizes of the
potential energy wells obtained from prior inspection of the liquid structure.
We perform a molecular dynamics simulation of a liquid water shell which
demonstrates that the method enables fast and energy-controlled water molecule
insertion in aqueous environments. The algorithm finds low energy
configurations of incoming water molecules around three orders of magnitude
faster than direct random insertion.
  This method is an important step towards dynamic simulations of open systems
and it may also prove useful for energy-biased ensemble average calculations of
the chemical potential.


  B-DNA solutions of suitable concentration form left-handed chiral nematic
phases (cholesterics). Such phases have also been observed in solutions of
other stiff or semiflexible chiral polymers; magnitude and handedness of the
cholesteric pitch are uniquely related to the molecular features. In this work
we present a theoretical method and a numerical procedure which, starting from
the structure of polyelectrolytes, lead to the prediction of the cholesteric
pitch. Molecular expressions for the free energy of the system are obtained on
the basis of steric and electrostatic interactions between polymers; the former
are described in terms of excluded volume, while a mean field approximation is
used for the latter. Calculations have been performed for 130 bp fragments of
B-DNA. The theoretical predictions provide an explanation for the experimental
behavior, by showing the counteracting role played by shape and charge
chirality of the molecule.


  The rate-limiting step of some enzymatic reactions is a physical step, i.e.
diffusion. The efficiency of such reactions can be improved through an increase
in the arrival rate of the substrate molecules, e.g. by a directed passage of
substrate (ligand) to active site after its random encounter with the enzyme
surface. Herein, we introduce a cellular automata model simulating the ligand
passage over the protein surface to its destined active site. The system is
simulated using the lattice gas automata with probabilistic transition rules.
Different distributions of amino acids over the protein surface are examined.
For each distribution, the hydration pattern is achieved and the mean number of
iteration steps needed for the ligand to arrive at the active site calculated.
Comparison of results indicates that the rate at which ligand arrives at the
active site is clearly affected by the distribution of amino acids outside the
active side. Such a process can facilitate the ligand diffusion towards the
active site thereby enhancing the efficiency of the enzyme action.


  We consider a model of active Brownian agents interacting via a harmonic
attractive potential in a two-dimensional system in the presence of noise. By
numerical simulations, we show that this model possesses a noise-induced
transition characterized by the breakdown of translational motion and the onset
of swarm rotation as the noise intensity is increased. Statistical properties
of swarm dynamics in the weak noise limit are further analytically
investigated.


  One of the most important steps in the developmental process of the bacteria
cell at the cellular level is the determination of the middle of the cell and
the proper placement of the septum, these being essential to the division of
the cell. In E. coli, this step depends on the proteins MinC, MinD, and MinE.
Exposure to a constant electric field may cause the bacteria cell division
mechanism to change, resulting in an abnormal cytokinesis. To see the effects
of an external field e.g., an electric or magnetic field on this process, we
have solved a set of deterministic reaction diffusion equations, which
incorporate the influence of an electric field. We have found some changes in
the dynamics of the oscillations of the min proteins from pole to pole. The
numerical results show some interesting effects, which are qualitatively in
good agreement with some experimental results.


  To study the dynamics of chemical processes, we often adopt rate equations to
observe the change in chemical concentrations. However, when the number of the
molecules is small, the fluctuations cannot be neglected. We often study the
effects of fluctuations with the help of stochastic differential equations.
  Chemicals are composed of molecules on a microscopic level. In principle, the
number of molecules must be an integer, which must only change discretely.
However, in analysis using stochastic differential equations, the fluctuations
are regarded as continuous changes. This approximation can only be valid if
applied to fluctuations that involve a sufficiently large number of molecules.
In the case of extremely rare chemical species, the actual discreteness of the
molecules may critically affect the dynamics of the system.
  To elucidate the effects of the discreteness, we study an autocatalytic
system consisting of several interacting chemical species with a small number
of molecules through stochastic particle simulations. We found novel states,
which were characterized as an extinction of molecule species, due to the
discrete nature of the molecules. We also observed a strong dependence of the
chemical concentrations on the size of the system, which was caused by
transitions to the novel states.


  Spreading of bacteria in a highly advective, disordered environment is
examined. Predictions of super-diffusive spreading for a simplified
reaction-diffusion equation are tested. Concentration profiles display
anomalous growth and super-diffusive spreading. A perturbation analysis yields
a crossover time between diffusive and super-diffusive behavior. The time's
dependence on the convection velocity and disorder is tested. Like the
simplified equation, the full linear reaction-diffusion equation displays
super-diffusive spreading perpendicular to the convection. However, for mean
positive growth rates the full nonlinear reaction-diffusion equation produces
symmetric spreading with a Fisher wavefront, whereas net negative growth rates
cause an asymmetry, with a slower wavefront velocity perpendicular to the
convection.


  As an attempt to understand the homochirality of organic molecules in life, a
chemical reaction model is proposed where the production of chiral monomers
from achiral substrate is catalyzed by the polymers of the same enatiomeric
type. This system has to be open because in a closed system the enhanced
production of chiral monomers by enzymes is compensated by the associated
enhancement in back reaction, and the chiral symmetry is conserved. Open flow
without cross inhibition is shown to lead to the chirality selection in a
general model. In polymerization, the influx of substrate from the ambience and
the efflux of chiral products for purposes other than the catalyst production
make the system necessarily open. The chiral symmetry is found to be broken if
the influx of substrate lies within a finite interval. As the efficiency of the
enzyme increases, the maximum value of the enantiomeric excess approaches unity
so that the chirality selection becomes complete.


  We present a mechanochemical model for myosin V, a two-headed processive
motor protein. We derive the properties of a dimer from those of an individual
head, which we model both with a 4-state cycle (detached, attached with ADP.Pi,
attached with ADP and attached without nucleotide) and alternatively with a
5-state cycle (where the power stroke is not tightly coupled to the phosphate
release). In each state the lever arm leaves the head at a different, but
fixed, angle. The lever arm itself is described as an elastic rod. The chemical
cycles of both heads are coordinated exclusively by the mechanical connection
between the two lever arms. The model explains head coordination by showing
that the lead head only binds to actin after the power stroke in the trail head
and that it only undergoes its power stroke after the trail head unbinds from
actin. Both models (4- and 5-state) reproduce the observed hand-over-hand
motion and fit the measured force-velocity relations. The main difference
between the two models concerns the load dependence of the run length, which is
much weaker in the 5-state model. We show how systematic processivity
measurement under varying conditions could be used to distinguish between both
models and to determine the kinetic parameters.


  We describe a model of cytoskeletal mechanics based on the force-induced
conformational change of protein cross-links in a stressed polymer network.
Slow deformation of simulated networks containing cross-links that undergo
repeated, serial domain unfolding leads to an unusual state--with many
cross-links accumulating near the critical force for further unfolding. Thermal
activation of these links gives rise to power-law rheology resembling the
previously unexplained mechanical response of living cells. Moreover, we
hypothesize that such protein cross-links function as biochemical
mechano-sensors of cytoskeletal deformation.


  We present a model of bi-phasic vesicle in the limit of large surface
tension. In this regime, the vesicle is completely stretched and well described
by two spherical caps with a fold which concentrates the membrane stress. The
conservation laws and geometric constraints restrict the space of possible
shapes to a pair of solutions labeled by a parameter $\tau$ given by {\it line
tension/pressure}. For a given $\tau$ value, the two solutions differ by the
length of the interface between domains. For a critical value $\tau\_c$, the
two vesicle shapes become identical and no solution exists above this critical
value. This model sheds new light on two proposed mechanisms (osmotic shocks
and molecule absorption) to explain the budding and the fission in recent
experiments.


  The type IV pilus retraction motor is found in many important bacterial
pathogens. It is the strongest known linear motor protein and is required for
bacterial infectivity. We characterize the dynamics of type IV pilus retraction
in terms of a stochastic chemical reaction model. We find that a two state
model can describe the experimental force velocity relation and qualitative
dependence of ATP concentration. The results indicate that the dynamics is
limited by an ATP-dependent step at low load and a force-dependent step at high
load, and that at least one step is effectively irreversible in the measured
range of forces. The irreversible nature of the sub-step(s) lead to interesting
predictions for future experiments: We find different parameterizations with
mathematically identical force velocity relations but different fluctuations
(diffusion constant). We also find a longer elementary step compared to an
earlier analysis, which agrees better with known facts about the structure of
the pilus filament and energetic considerations. We conclude that more
experimental data is needed, and that further retraction experiments are likely
to resolve interesting details and give valuable insights into the PilT
machinery. In light of our findings, the fluctuations of the retraction
dynamics emerge as a key property to be studied in future experiments.


  We study the statistical mechanics of RNA secondary structures designed to
have an attraction between two different types of structures as a model system
for heteropolymer aggregation. The competition between the branching entropy of
the secondary structure and the energy gained by pairing drives the RNA to
undergo a `temperature independent' second order phase transition from a molten
to an aggregated phase'. The aggregated phase thus obtained has a
macroscopically large number of contacts between different RNAs. The partition
function scaling exponent for this phase is \theta ~ 1/2 and the crossover
exponent of the phase transition is \nu ~ 5/3. The relevance of these
calculations to the aggregation of biological molecules is discussed.


  Light harvesting complexes 2 (LH2) from Rhodospirillum (Rs.) molischianum and
Rhodopseudomonas (Rps.) acidophila form ring complexes out of eight or nine
identical subunits, respectively. Here, we investigate computationally what
factors govern the different ring sizes. Starting from the crystal structure
geometries, we embed two subunits of each species into their native
lipid-bilayer/water environment. Using molecular dynamics simulations with
umbrella sampling and steered molecular dynamics, we probe the free energy
profiles along two reaction coordinates, the angle and the distance between two
subunits. We find that two subunits prefer to arrange at distinctly different
angles, depending on the species, at about 42.5 deg for Rs. molischianum and at
about 38.5 deg for Rps. acidophila, which is likely to be an important factor
contributing to the assembly into different ring sizes. Our calculations
suggest a key role of surface contacts within the transmembrane domain in
constraining these angles, whereas the strongest interactions stabilizing the
subunit dimers are found in the C-, and to a lesser extent, N-terminal domains.
The presented computational approach provides a promising starting point to
investigate the factors contributing to the assembly of protein complexes, in
particular if combined with modeling of genetic variants.


  Theory and experiment have established that F-actin filaments are strongly
attached to the intracellular parasites (such as Listeria) they propel with
``comet tails''. We consider the implications of these observations for
propulsion. By calculating the motion produced in various models of attachment
and comparing to experiment we demonstrate that the attachment must be sliding
rather than hinged. By modeling experiments on ActA-coated spheres we draw
conclusions regarding the interaction between F-actin and their surfaces that
may also be applicable to living systems.


  Motivated by experiments on single-headed kinesin KIF1A, we develop a model
of intra-cellular transport by interacting molecular motors. It captures
explicitly not only the effects of ATP hydrolysis, but also the ratchet
mechanism which drives individual motors. Our model accounts for the
experimentally observed single molecule properties in the low density limit and
also predicts a phase diagram that shows the influence of hydrolysis and
Langmuir kinetics on the collective spatio-temporal organization of the motors.
Finally, we provide experimental evidence for the existence of domain walls in
our {\it in-vitro} experiment with fluorescently labeled KIF1A.


  Many complex systems in nature and society can be described in terms of
networks capturing the intricate web of connections among the units they are
made of. A key question is how to interpret the global organization of such
networks as the coexistence of their structural subunits (communities)
associated with more highly interconnected parts. Identifying these a priori
unknown building blocks (such as functionally related proteins, industrial
sectors and groups of people) is crucial to the understanding of the structural
and functional properties of networks. The existing deterministic methods used
for large networks find separated communities, whereas most of the actual
networks are made of highly overlapping cohesive groups of nodes. Here we
introduce an approach to analysing the main statistical features of the
interwoven sets of overlapping communities that makes a step towards uncovering
the modular structure of complex systems. After defining a set of new
characteristic quantities for the statistics of communities, we apply an
efficient technique for exploring overlapping communities on a large scale. We
find that overlaps are significant, and the distributions we introduce reveal
universal features of networks. Our studies of collaboration, word-association
and protein interaction graphs show that the web of communities has non-trivial
correlations and specific scaling properties.


  A number of recent works have concentrated on a few statistical properties of
complex networks, such as the clustering, the right-skewed degree distribution
and the community, which are common to many real world networks. In this paper,
we address the hierarchy property sharing among a large amount of networks.
Based upon the eigenvector centrality (EC) measure, a method is proposed to
reconstruct the hierarchical structure of a complex network. It is tested on
the Santa Fe Institute collaboration network, whose structure is well known. We
also apply it to a Mathematicians' collaboration network and the protein
interaction network of Yeast. The method can detect significantly hierarchical
structures in these networks.


  Autonomous oscillations in biological systems may have a biochemical origin
or result from an interplay between force-generating and visco-elastic
elements. In molecular motor assemblies the force-generating elements are
molecular engines and the visco-elastic elements are stiff cytoskeletal
polymers. The physical mechanism leading to oscillations depends on the
particular architecture of the assembly. Existing models can be grouped into
two distinct categories: systems with a delayed force activation and anomalous
force-velocity relations. We discuss these systems within phase plane analysis
known from the theory of dynamic systems and by adopting methods from control
theory, the Nyquist stability criterion.


  We present computer-assisted methods for analyzing stochastic models of gene
regulatory networks. The main idea that underlies this equation-free analysis
is the design and execution of appropriately-initialized short bursts of
stochastic simulations; the results of these are processed to estimate
coarse-grained quantities of interest, such as mesoscopic transport
coefficients. In particular, using a simple model of a genetic toggle switch,
we illustrate the computation of an effective free energy and of a
state-dependent effective diffusion coefficient that characterize an
unavailable effective Fokker-Planck equation. Additionally we illustrate the
linking of equation-free techniques with continuation methods for performing a
form of stochastic "bifurcation analysis"; estimation of mean switching times
in the case of a bistable switch is also implemented in this equation-free
context. The accuracy of our methods is tested by direct comparison with
long-time stochastic simulations. This type of equation-free analysis appears
to be a promising approach to computing features of the long-time,
coarse-grained behavior of certain classes of complex stochastic models of gene
regulatory networks, circumventing the need for long Monte Carlo simulations.


  Helicases are molecular motors which unwind double-stranded nucleic acids
(dsNA) in cells. Many helicases move with directional bias on single-stranded
(ss) nucleic acids, and couple their directional translocation to strand
separation. A model of the coupling between translocation and unwinding uses an
interaction potential to represent passive and active helicase mechanisms. A
passive helicase must wait for thermal fluctuations to open dsNA base pairs
before it can advance and inhibit NA closing. An active helicase directly
destabilizes dsNA base pairs, accelerating the opening rate. Here we extend
this model to include helicase unbinding from the nucleic-acid strand. The
helicase processivity depends on the form of the interaction potential. A
passive helicase has a mean attachment time which does not change between ss
translocation and ds unwinding, while an active helicase in general shows a
decrease in attachment time during unwinding relative to ss translocation. In
addition, we describe how helicase unwinding velocity and processivity vary if
the base-pair binding free energy is changed.


  We study the influence of disorder in the helical actin structure on the
myosin V step size, predicted from the elastic lever arm model. We show that
fluctuations of +-5 degrees per actin subunit, as proposed by Egelman et al.,
significantly alter the distribution of step sizes and improve the agreement
with experimental data.


  Fluctuations of cell state, e.g., abundances of some proteins, have attracted
much attention both theoretically and experimentally. The distribution of such
state over cells, however, is not only a result of intracellular stochastic
process, but is also influenced by the growth in cell numbers that depends on
the state. By incorporating the growth-death process into the standard
Fokker--Planck equation for the probability distribution, a nonlinear temporal
evolution equation of distribution is obtained that includes a self-consistent
growth term. The derived equation is generally solved analytically by means of
eigenfunction expansions. By focusing on the case with linear relaxation, two
examples are considered as applications of the proposed general formalism.
First, by assuming that the growth rate of a cell increases linearly with the
state value $x$, the shift of the average state value $x$ due to the growth
effect is shown to be proportional to the variance of the state $x$ and the
relaxation time, similarly with the biological fluctuation- response
relationship. Second, when there is a gap in the growth rate at some threshold
value for the state $x$, existence of a critical gap value is demonstrated,
beyond which the average growth rate starts to increase. This critical value is
again obtained in terms of the relaxation time and the variance of $x$, all of
which are experimentally measurable quantities. The relevance of the results to
the analysis of biological data on the distribution of cell states, as obtained
for example by flow cytometry, is discussed.


  I introduce the Fermi Paradox and some of its solutions. Then I present my
own solution which includes two proposals called the Subanthropic Principle and
the Undetectability Conjecture. After discussing some consequences of this
solution, I make some comments about brane world scenarios and their potential
to strengthen the Fermi Paradox. Finally, in the appendix I have included some
questions and answers that came up during this Forum. --
  En primer lugar introduzco la Paradoja de Fermi y algunas de sus soluciones.
Luego presento la solucion que yo he propuesto, que incluye dos hipotesis que
denomino el Principio Subantropico y la Conjetura de Indetectabilidad. Despues
de discutir algunas consecuencias de esta solucion, paso a hacer algunos
comentarios sobre las teorias de universos branas y su potencial para reforzar
la Paradoja de Fermi. Finalmente, en el apendice incluyo algunas preguntas y
respuestas que surgieron durante este Forum.


  Topological properties of native folds are obtained from statistical analysis
of 160 low homology proteins covering the four structural classes. This is done
analysing one, two and three-vertex joint distribution of quantities related to
the corresponding network of amino acid residues. Emphasis on the amino acid
residue hydrophobicity leads to the definition of their center of mass as
vertices in this contact network model with interactions represented by edges.
The network analysis helps us to interpret experimental results such as
hydrophobic scales and fraction of buried accessible surface area in terms of
the network connectivity. To explore the vertex type dependent correlations, we
build a network of hydrophobic and polar vertices. This procedure presents the
wiring diagram of the topological structure of globular proteins leading to the
following attachment probabilities between hydrophobic-hydrophobic 0.424(5),
hydrophobic-polar 0.419(2) and polar-polar 0.157(3) residues.


  Single-molecule experiments in which force is applied to DNA or RNA molecules
have enabled important discoveries of nucleic acid properties and nucleic
acid-enzyme interactions. These experiments rely on a model of the polymer
force-extension behavior to calibrate the experiments; typically the
experiments use the worm-like chain (WLC) theory for double-stranded DNA and
RNA. This theory agrees well with experiments for long molecules. Recent
single-molecule experiments have used shorter molecules, with contour lengths
in the range of 1-10 persistence lengths. Most WLC theory calculations to date
have assumed infinite molecule lengths, and do not agree well with experiments
on shorter chains. Key physical effects that become important when shorter
molecules are used include (i) boundary conditions which constrain the allowed
fluctuations at the ends of the molecule and (ii) rotational fluctuations of
the bead to which the polymer is attached, which change the apparent extension
of the molecule. We describe the finite worm-like chain (FWLC) theory, which
takes into account these effects. We show the FWLC predictions diverge from the
classic WLC solution for molecules with contour lengths a few times the
persistence length. Thus the FWLC will allow more accurate experimental
calibration for relatively short molecules, facilitating future discoveries in
single-molecule force microscopy.


  Models of the convergence of opinion in social systems have been the subject
of a considerable amount of recent attention in the physics literature. These
models divide into two classes, those in which individuals form their beliefs
based on the opinions of their neighbors in a social network of personal
acquaintances, and those in which, conversely, network connections form between
individuals of similar beliefs. While both of these processes can give rise to
realistic levels of agreement between acquaintances, practical experience
suggests that opinion formation in the real world is not a result of one
process or the other, but a combination of the two. Here we present a simple
model of this combination, with a single parameter controlling the balance of
the two processes. We find that the model undergoes a continuous phase
transition as this parameter is varied, from a regime in which opinions are
arbitrarily diverse to one in which most individuals hold the same opinion. We
characterize the static and dynamical properties of this transition.


  The TATA-box Binding Protein (TBP) is required by all three eukaryotic RNA
polymerases for the initiation of transcription from most promoters. TBP
recognizes, binds to, and bends promoter sequences called ``TATA-boxes'' in the
DNA. We present results from the study of individual Saccharomyces cerevisia
TBPs interacting with single DNA molecules containing a TATA-box. Using video
microscopy, we observed the Brownian motion of beads tethered by short
surface-bound DNA. When TBP binds to and bends the DNA, the conformation of the
DNA changes and the amplitude of Brownian motion of the tethered bead is
reduced compared to that of unbent DNA. We detected individual binding and
dissociation events and derived kinetic parameters for the process.
Dissociation was induced by increasing the salt concentration or by directly
pulling on the tethered bead using optical tweezers. In addition to the
well-defined free and bound classes of Brownian motion, we observed another two
classes of motion. These extra classes were identified with intermediate states
on a three-step, linear binding pathway. Biological implications of the
intermediate states are discussed.


  Many ribosomes simultaneously move on the same messenger RNA (mRNA), each
separately synthesizing the protein coded by the mRNA. Earlier models of
ribosome traffic represent each ribosome by a ``self-propelled particle'' and
capture the dynamics by an extension of the totally asymmetric simple exclusion
process (TASEP). In contrast, here we develope a theoretical model that not
only incorporates the mutual exclusions of the interacting ribosomes, but also
describes explicitly the mechano-chemistry of each of these individual cyclic
machines during protein synthesis. Using analytical and numerical techniques of
non-equilibrium statistical mechanics, we analyze this model and illustrate its
power by making experimentally testable predictions on the rate of protein
synthesis in real time and the density profile of the ribosomes on some mRNAs
in E-Coli.


  We present a simple model for the formation of pairs in multi-agent
populations of type A and B which move freely on a spatial network. Each agent
of population A (and B) is labeled as Ai (and Bj) with i=1,.. NA (and j=1,..NB)
and carries its own individual list of characteristics or 'phenotype'. When
agents from opposite populations encounter one another on the network, they can
form a relationship if not already engaged in one. The length of time for which
any given pair stays together depends on the compatibility of the two
constituent agents. Possible applications include the human dating scenario,
and the commercial domain where two types of businesses A and B have members of
each type looking for a business partner, i.e. Ai+Bj-->Rij. The pair Rij then
survives for some finite time before dissociating Rij-->Ai+Bj. There are many
possible generalizations of this basic setup. Here we content ourselves with
some initial numerical results for the simplest of network topologies, together
with some accompanying analytic analysis.


  Molecular motors do not work in isolation {\it in-vivo}. We highlight some of
the coordinations, cooperations and competitions that determine the collective
properties of molecular motors in eukaryotic cells. In the context of
traffic-like movement of motors on a track, we emphasize the importance of
single-motor bio-chemical cycle and enzymatic activity on their collective
spatio-temporal organisation. Our modelling strategy is based on a synthesis-
the same model describes the single-motor mechano-chemistry at sufficiently low
densities whereas at higher densities it accounts for the collective flow
properties and the density profiles of the motors. We consider two specific
examples, namely, traffic of single-headed kinesin motors KIF1A on a
microtubule track and ribosome traffic on a messenger RNA track.


  Many {\it ribosomes} simultaneously move on the same messenger RNA (mRNA),
each synthesizing separately a copy of the same protein. In contrast to the
earlier models, here {\it we develop a ``unified'' theoretical model} that not
only incorporates the {\it mutual exclusions} of the interacting ribosomes, but
also describes explicitly the mechano-chemistry of each of these macromolecular
machines during protein synthesis. Using analytical and numerical techniques of
non-equilibrium statistical mechanics, we analyze the rates of protein
synthesis and the spatio-temporal organization of the ribosomes in this model.
We also predict how these properties would change with the changes in the rates
of the various chemo-mechanical processes in each ribosome. Finally, we
illustrate the power of this model by making experimentally testable
predictions on the rates of protein synthesis and the density profiles of the
ribosomes on some mRNAs in {\it E-coli}.


  The time spent by an interacting Brownian molecule inside a bounded
microdomain has many applications in cellular biology, because the number of
bounds is a quantitative signal, which can initiate a cascade of chemical
reactions and thus has physiological consequences. In the present article, we
propose to estimate the mean time spent by a Brownian molecule inside a
microdomain $\Omega$ which contains small holes on the boundary and agonist
molecules located inside. We found that the mean time depends on several
parameters such as the backward binding rate (with the agonist molecules), the
mean escape time from the microdomain and the mean time a molecule reaches the
binding sites (forward binding rate). In addition, we estimate the mean and the
variance of the number of bounds made by a molecule before it exits $\Omega$.
These estimates rely on a boundary layer analysis of a conditional mean first
passage time, solution of a singular partial differential equation. In
particular, we apply the present results to obtain an estimate of the mean time
spent (Dwell time) by a Brownian receptor inside a synaptic domain, when it
moves freely by lateral diffusion on the surface of a neuron and interacts
locally with scaffolding molecules.


  We address the problem of understanding the variable abundance of 3-node and
4-node subgraphs (motifs) in complex networks from a dynamical point of view.
As a criterion in the determination of the functional significance of a n-node
subgraph, we propose an analytic method to measure the stability of the
synchronous state (SSS) the subgraph displays. We show that, for undirected
graphs, the SSS is correlated with the relative abundance, while in directed
graphs the correlation exists only for some specific motifs.


  We consider the inverse-folding problem for RNA secondary structures: for a
given (pseudo-knot-free) secondary structure find a sequence that has that
structure as its ground state. If such a sequence exists, the structure is
called designable. We implemented a branch-and-bound algorithm that is able to
do an exhaustive search within the sequence space, i.e., gives an exact answer
whether such a sequence exists. The bound required by the branch-and-bound
algorithm are calculated by a dynamic programming algorithm. We consider
different alphabet sizes and an ensemble of random structures, which we want to
design. We find that for two letters almost none of these structures are
designable. The designability improves for the three-letter case, but still a
significant fraction of structures is undesignable. This changes when we look
at the natural four-letter case with two pairs of complementary bases:
undesignable structures are the exception, although they still exist. Finally,
we also study the relation between designability and the algorithmic complexity
of the branch-and-bound algorithm. Within the ensemble of structures, a high
average degree of undesignability is correlated to a long time to prove that a
given structure is (un-)designable. In the four-letter case, where the
designability is high everywhere, the algorithmic complexity is highest in the
region of naturally occurring RNA.


  From a nanoscience perspective, cellular processes and their reduced in vitro
imitations provide extraordinary examples for highly robust few or single
molecule reaction pathways. A prime example are biochemical reactions involving
DNA molecules, and the coupling of these reactions to the physical
conformations of DNA. In this review, we summarise recent results on the
following phenomena: We investigate the biophysical properties of DNA-looping
and the equilibrium configurations of DNA-knots, whose relevance to biological
processes are increasingly appreciated. We discuss how random DNA-looping may
be related to the efficiency of the target search process of proteins for their
specific binding site on the DNA molecule. And we dwell on the spontaneous
formation of intermittent DNA nanobubbles and their importance for biological
processes, such as transcription initiation. The physical properties of DNA may
indeed turn out to be particularly suitable for the use of DNA in nanosensing
applications.


  We introduce a minimalistic model based on dynamic node deletion and node
duplication with heterodimerisation. The model is intended to capture the
essential features of the evolution of protein interaction networks. We derive
an exact two-step rate equation to describe the evolution of the degree
distribution. We present results for the case of a fixed-size network. The
results are based on the exact numerical solution to the rate equation which
are consistent with Monte Carlo simulations of the model's dynamics. Power-law
degree distributions with apparent exponents $<1$ were observed for generic
parameter choices. However, a proper finite-size scaling analysis revealed that
the actual critical exponent in such cases is equal to 1. We present a
mean-field argument to determine the asymptotic value of the average degree,
illustrating the existence of an attractive fixed point, and corroborate this
result with numerical simulations of the first moment of the degree
distribution as described by the two-step rate equation. Using the above
results, we show that the apparent exponent is determined by the
heterodimerisation probability. Our preliminary results are consistent with
empirical data for a wide range of organisms, and we believe that through
implementing some of the suggested modifications, the model could be
well-suited to other types of biological and non-biological networks.


  Halting a computer or biological virus outbreak requires a detailed
understanding of the timing of the interactions between susceptible and
infected individuals. While current spreading models assume that users interact
uniformly in time, following a Poisson process, a series of recent measurements
indicate that the inter-contact time distribution is heavy tailed,
corresponding to a temporally inhomogeneous bursty contact process. Here we
show that the non-Poisson nature of the contact dynamics results in prevalence
decay times significantly larger than predicted by the standard Poisson process
based models. Our predictions are in agreement with the detailed time resolved
prevalence data of computer viruses, which, according to virus bulletins, show
a decay time close to a year, in contrast with the one day decay predicted by
the standard Poisson process based models.