  Proteins are known to locate their specific targets on DNA up to two orders
of magnitude faster than predicted by the Smoluchowski three-dimensional
diffusion rate. One of the mechanisms proposed to resolve this discrepancy is
termed "intersegment transfer". Many proteins have two DNA binding sites and
can transfer from one DNA segment to another without dissociation to water. We
calculate the target search rate for such proteins in a dense globular DNA,
taking into account intersegment transfer working in conjunction with DNA
motion and protein sliding along DNA. We show that intersegment transfer plays
a very important role in cases where the protein spends most of its time
adsorbed on DNA.


  We calculate analytically the critical connectivity $K_c$ of Random Threshold
Networks (RTN) for homogeneous and inhomogeneous thresholds, and confirm the
results by numerical simulations. We find a super-linear increase of $K_c$ with
the (average) absolute threshold $|h|$, which approaches $K_c(|h|) \sim
h^2/(2\ln{|h|})$ for large $|h|$, and show that this asymptotic scaling is
universal for RTN with Poissonian distributed connectivity and threshold
distributions with a variance that grows slower than $h^2$. Interestingly, we
find that inhomogeneous distribution of thresholds leads to increased
propagation of perturbations for sparsely connected networks, while for densely
connected networks damage is reduced; the cross-over point yields a novel,
characteristic connectivity $K_d$, that has no counterpart in Boolean networks.
Last, local correlations between node thresholds and in-degree are introduced.
Here, numerical simulations show that even weak (anti-)correlations can lead to
a transition from ordered to chaotic dynamics, and vice versa. It is shown that
the naive mean-field assumption typical for the annealed approximation leads to
false predictions in this case, since correlations between thresholds and
out-degree that emerge as a side-effect strongly modify damage propagation
behavior.


  The maintainance of a stable periodicity during the yeast metabolic cycle
involving approximately half of the genome requires a very strict and efficient
control of gene expression. For this reason, the metabolic cycle is a very good
candidate for testing the role of a class of post-transcriptional regulators,
the so called PUF-family, whose genome-wide mRNA binding specificity was
recently experimentally assessed. Here we show that an integrated computational
analysis of gene expression time series during the metabolic cycle and the mRNA
binding specificity of PUF-family proteins allow for a clear demonstration of
the very specific role exerted by selective post-transcriptional mRNA
degradation in yeast metabolic cycle global regulation.


  Local minima and the saddle points separating them in the energy landscape
are known to dominate the dynamics of biopolymer folding. Here we introduce a
notion of a "folding funnel" that is concisely defined in terms of energy
minima and saddle points, while at the same time conforming to a notion of a
"folding funnel" as it is discussed in the protein folding literature.


  Cancer results from genetic alterations that disturb the normal cooperative
behavior of cells. Recent high-throughput genomic studies of cancer cells have
shown that the mutational landscape of cancer is complex and that individual
cancers may evolve through mutations in as many as 20 different
cancer-associated genes. We use data published by Sjoblom et al. (2006) to
develop a new mathematical model for the somatic evolution of colorectal
cancers. We employ the Wright-Fisher process for exploring the basic parameters
of this evolutionary process and derive an analytical approximation for the
expected waiting time to the cancer phenotype. Our results highlight the
relative importance of selection over both the size of the cell population at
risk and the mutation rate. The model predicts that the observed genetic
diversity of cancer genomes can arise under a normal mutation rate if the
average selective advantage per mutation is on the order of 1%. Increased
mutation rates due to genetic instability would allow even smaller selective
advantages during tumorigenesis. The complexity of cancer progression thus can
be understood as the result of multiple sequential mutations, each of which has
a relatively small but positive effect on net cell growth.


  Site-specific recombination on supercoiled circular DNA yields a variety of
knotted or catenated products. We develop a model of this process, and give
extensive experimental evidence that the assumptions of our model are
reasonable. We then characterize all possible knot or catenane products that
arise from the most common substrates. We apply our model to tightly prescribe
the knot or catenane type of previously uncharacterized data.


  In eukaryotes, DNA replication is initiated along each chromosome at multiple
sites called replication origins. Locally, each replication origin is
"licensed", or specified, at the end of the M and the beginning of G1 phases of
the cell cycle. During S phase when DNA synthesis takes place, origins are
activated in stages corresponding to early and late replicating domains. The
staged and progressive activation of replication origins reflects the need to
maintain a strict balance between the number of active replication forks and
the rate at which DNA synthesis procedes. This suggests that origin densities
(frequency of intiation) and replication fork movement (rates of elongation)
must be co-regulated in order to guarantee the efficient and complete
duplication of each subchromosomal domain. Emerging evidence supports this
proposal and suggests that the ATM/ATR intra-S phase checkpoint plays an
important role in the co-regulation of initiation frequencies and rates of
elongation. In the following, we review recent results concerning the
mechanisms governing the global regulation of DNA replication and discuss the
roles these mechanisms play in maintaining genome stability during both a
normal and perturbed S phase.


  The multiple worlds of genetically manipulated laboratory organisms such as
transgenic mice or worms with certain gene mutations are somewhat reminiscent
of parallel worlds in quantum mechanics. So are various models of aging tested
in such organisms. In this context, the tumor suppressor p53 has been found to
either accelerate or delay aging, the latter, for instance, in conjunction with
ARF, another tumor suppressor, as shown very recently. To more easily determine
which of these artificial settings comes closest to real life, I discuss here
their features in the light of my protein structure-based insights that have
led me to propose a physiological anti-aging role for the retinoblastoma tumor
suppressor protein (RB) over the past four years.


  We propose a simple model of pathologic microsatellite expansion, and
describe an inherent self-repairing mechanism working against expansion. We
prove that if the probabilities of elementary expansions and contractions are
equal, microsatellite expansions are always self-repairing. If these
probabilities are different, self-reparation does not work. Mosaicism,
anticipation and reverse mutation cases are discussed in the framework of the
model. We explain these phenomena and provide some theoretical evidence for
their properties, for example the rarity of reverse mutations.


  Recently, structural analysis of the human transferrin and growth hormone
(GH) amino acid sequences has unravelled that they harbor a motif identical to
a pattern found in viral oncoproteins known to bind the primarily nuclear tumor
suppressor retinoblastoma protein (RB). Since related signatures had previously
been identified also in insulin and the two insulin-like growth factors (IGFs),
the aim of the current study has been to investigate whether further hints
substantiating these reported homologies can be found in silico. Here,
additional similarities are presented supporting the notion of an insulin
superfamily of growth-promoting proteins with dual localization in the
extracellular environment and the intracellular space, particularly in the
nucleus, as well as characterized by a tropism for RB.


  We simulate an individual-based model that represents both the phenotype and
genome of digital organisms with predator-prey interactions. We show how
open-ended growth of complexity arises from the invariance of genetic evolution
operators with respect to changes in the complexity, and that the dynamics
which emerges is controlled by a non-equilibrium critical point. The mechanism
is analogous to the development of the cascade in fluid turbulence.


  Twist stiffness and an asymmetric bending stiffness of a polymer or a polymer
bundle is captured by the elastic ribbon model. We investigate the effects a
ring geometry induces to a thermally fluctuating ribbon, finding bend-bend
coupling in addition to twist-bend coupling. Furthermore, due to the geometric
constraint the polymer's effective bending stiffness increases. A new parameter
for experimental investigations of polymer bundles is proposed: the mean square
diameter of a ribbonlike ring, which is determined analytically in the
semiflexible limit. Monte Carlo simulations are performed which affirm the
model's prediction up to high flexibility.


  Using magnetic tweezers to investigate the mechanical response of single
chromatin fibers, we show that fibers submitted to large positive torsion
transiently trap positive turns, at a rate of one turn per nucleosome. A
comparison with the response of fibers of tetrasomes (the (H3-H4)2 tetramer
bound with ~50 bp of DNA) obtained by depletion of H2A-H2B dimers, suggests
that the trapping reflects a nucleosome chiral transition to a metastable form
built on the previously documented righthanded tetrasome. In view of its low
energy, <8 kT, we propose this transition is physiologically relevant and
serves to break the docking of the dimers on the tetramer which in the absence
of other factors exerts a strong block against elongation of transcription by
the main RNA polymerase.


  Recently, several studies have investigated the transcription process
associated to specific genetic regulatory networks. In this work, we present a
stochastic approach for analyzing the dynamics and effect of negative feedback
loops (FBL) on the transcriptional noise. First, our analysis allows us to
identify a bimodal activity depending of the strength of self-repression
coupling D. In the strong coupling region D>>1, the variance of the
transcriptional noise is found to be reduced a 28 % more than described
earlier. Secondly, the contribution of the noise effect to the abundance of
regulating protein becomes manifest when the coefficient of variation is
computed. In the strong coupling region, this coefficient is found to be
independent of all parameters and in fair agreement with the experimentally
observed values. Finally, our analysis reveals that the regulating protein is
significantly induced by the intrinsic and external noise in the strong
coupling region. In short, it indicates that the existence of inherent noise in
FBL makes it possible to produce a basal amount of proteins even though the
repression level D is very strong.


  Chemical reaction systems with a low to moderate number of molecules are
typically modeled as discrete jump Markov processes. These systems are
oftentimes simulated with methods that produce statistically exact sample paths
such as the Gillespie Algorithm or the Next Reaction Method. In this paper we
make explicit use of the fact that the initiation times of the reactions can be
represented as the firing times of independent, unit rate Poisson processes
with internal times given by integrated propensity functions. Using this
representation we derive a modified Next Reaction Method and, in a way that
achieves efficiency over existing approaches for exact simulation, extend it to
systems with time dependent propensities as well as to systems with delays.


  By explicitly representing the reaction times of discrete chemical systems as
the firing times of independent, unit rate Poisson processes, we develop a new
adaptive tau-leaping procedure. The procedure developed is novel in that
accuracy is guaranteed by performing postleap checks. Because the
representation we use separates the randomness of the model from the state of
the system, we are able to perform the postleap checks in such a way that the
statistics of the sample paths generated will not be biased by the rejections
of leaps. Further, since any leap condition is ensured with a probability of
one, the simulation method naturally avoids negative population values


  Phylogenetic diversity is a measure for describing how much of an
evolutionary tree is spanned by a subset of species. If one applies this to the
(unknown) subset of current species that will still be present at some future
time, then this `future phylogenetic diversity' provides a measure of the
impact of various extinction scenarios in biodiversity conservation. In this
paper we study the distribution of future phylogenetic diversity under a simple
model of extinction (a generalized `field of bullets' model). We show that the
distribution of future phylogenetic diversity converges to a normal
distribution as the number of species grows (under mild conditions, which are
necessary). We also describe an algorithm to compute the distribution
efficiently, provided the edge lengths are integral, and briefly outline the
significance of our findings for biodiversity conservation.


  We perform extensive Monte Carlo simulations of a lattice model and the Go
potential to investigate the existence of folding pathways at the level of
contact cluster formation for two native structures with markedly different
geometries. Our analysis of folding pathways revealed a common underlying
folding mechanism, based on nucleation phenomena, for both protein models.
However, folding to the more complex geometry (i.e. that with more non-local
contacts) is driven by a folding nucleus whose geometric traits more closely
resemble those of the native fold. For this geometry folding is clearly a more
cooperative process.


  In this study we evaluate, at full atomic detail, the folding processes of
two small helical proteins, the B domain of protein A and the Villin headpiece.
Folding kinetics are studied by performing a large number of ab initio Monte
Carlo folding simulations using a single transferable all-atom potential. Using
these trajectories, we examine the relaxation behavior, secondary structure
formation, and transition-state ensembles (TSEs) of the two proteins and
compare our results with experimental data and previous computational studies.
To obtain a detailed structural information on the folding dynamics viewed as
an ensemble process, we perform a clustering analysis procedure based on graph
theory. Moreover, rigorous pfold analysis is used to obtain representative
samples of the TSEs and a good quantitative agreement between experimental and
simulated Fi-values is obtained for protein A. Fi-values for Villin are also
obtained and left as predictions to be tested by future experiments. Our
analysis shows that two-helix hairpin is a common partially stable structural
motif that gets formed prior to entering the TSE in the studied proteins. These
results together with our earlier study of Engrailed Homeodomain and recent
experimental studies provide a comprehensive, atomic-level picture of folding
mechanics of three-helix bundle proteins.


  Stability of the branching structure of an RNA molecule is an important
condition for its function. In this letter we show that the melting
thermodynamics of RNA molecules is very sensitive to their branching geometry
for the case of a molecule whose groundstate has the branching geometry of a
Cayley Tree and whose pairing interactions are described by the Go model.
Whereas RNA molecules with a linear geometry melt via a conventional continuous
phase transition with classical exponents, molecules with a Cayley Tree
geometry are found to have a free energy that seems smooth, at least within our
precision. Yet, we show analytically that this free energy in fact has a
mathematical singularity at the stability limit of the ordered structure. The
correlation length appears to diverge on the high-temperature side of this
singularity.


  Within individual bacteria, we combine force-dependent polymerization
dynamics of individual MreB protofilaments with an elastic model of
protofilament bundles buckled into helical configurations. We use variational
techniques and stochastic simulations to relate the pitch of the MreB helix,
the total abundance of MreB, and the number of protofilaments. By comparing our
simulations with mean-field calculations, we find that stochastic fluctuations
are significant. We examine the quasi-static evolution of the helical pitch
with cell growth, as well as timescales of helix turnover and denovo
establishment. We find that while the body of a polarized MreB helix treadmills
towards its slow-growing end, the fast-growing tips of laterally associated
protofilaments move towards the opposite fast-growing end of the MreB helix.
This offers a possible mechanism for targeted polar localization without
cytoplasmic motor proteins.


  In the present study, we show that ascorbic acid dose-dependently inhibited
interleukin-1beta (IL-1beta)-mediated PGE2 synthesis in the human neuronal cell
line, SK-N-SH. Furthermore, in combination with aspirin, ascorbic acid
augmented the inhibitory effect of aspirin on PGE2 synthesis. However, ascorbic
acid had no synergistic effect along with other COX inhibitors (SC-58125 and
indomethacin). The inhibition of IL-1beta-mediated PGE2 synthesis by ascorbic
acid was not due to the inhibition of the expression of COX-2 or microsomal
prostaglandin E synthase (mPGES-1). Rather, ascorbic acid dose-dependently
(0.1-100 microM) produced a significant reduction in IL-1beta-mediated
production of 8-iso-prostaglandin F2alpha (8-iso-PGF2alpha), a reliable
indicator of free radical formation, suggesting that the effects of ascorbic
acid on COX-2-mediated PGE2 biosynthesis may be the result of the maintenance
of the neuronal redox status since COX activity is known to be enhanced by
oxidative stress. Our results provide in vitro evidence that the
neuroprotective effects of ascorbic acid may depend, at least in part, on its
ability to reduce neuronal COX-2 activity and PGE2 synthesis, owing to its
antioxidant properties. Further, these experiments suggest that a combination
of aspirin with ascorbic acid constitutes a novel approach to render COX-2 more
sensitive to inhibition by aspirin, allowing an anti-inflammatory therapy with
lower doses of aspirin, thereby avoiding the side effects of the usually high
dose aspirin treatment.


  Neurons are subject to various kinds of noise. In addition to synaptic noise,
the stochastic opening and closing of ion channels represents an intrinsic
source of noise that affects the signal processing properties of the neuron. In
this paper, we studied the response of a stochastic Hodgkin-Huxley neuron to
transient input subthreshold pulses. It was found that the average response
time decreases but variance increases as the amplitude of channel noise
increases. In the case of single pulse detection, we show that channel noise
enables one neuron to detect the subthreshold signals and an optimal membrane
area (or channel noise intensity) exists for a single neuron to achieve optimal
performance. However, the detection ability of a single neuron is limited by
large errors. Here, we test a simple neuronal network that can enhance the
pulse detecting abilities of neurons and find dozens of neurons can perfectly
detect subthreshold pulses. The phenomenon of intrinsic stochastic resonance is
also found both at the level of single neurons and at the level of networks. At
the network level, the detection ability of networks can be optimized for the
number of neurons comprising the network.


  A tendency in biological theorizing is to formulate principles above or equal
to Evolution by Variation and Selection of Darwin and Wallace. In this letter I
analyze one such recent proposal which did so for the developmental ascendency.
I show that though the idea of developmental ascendency is brilliant, this is
in wrong order in the hierarchical structure of biological theories and can
easily generate confusing. Several other examples are also briefly discussed in
the note added.


  RNA polymerase (RNAP) is an enzyme that synthesizes a messenger RNA (mRNA)
strand which is complementary to a single-stranded DNA template. From the
perspective of physicists, an RNAP is a molecular motor that utilizes chemical
energy input to move along the track formed by a DNA. In many circumstances,
which are described in this paper, a large number of RNAPs move simultaneously
along the same track; we refer to such collective movements of the RNAPs as
RNAP traffic. Here we develop a theoretical model for RNAP traffic by
incorporating the steric interactions between RNAPs as well as the
mechano-chemical cycle of individual RNAPs during the elongation of the mRNA.
By a combination of analytical and numerical techniques, we calculate the rates
of mRNA synthesis and the average density profile of the RNAPs on the DNA
track. We also introduce, and compute, two new measures of fluctuations in the
synthesis of RNA. Analyzing these fluctuations, we show how the level of {\it
intrinsic noise} in mRNA synthesis depends on the concentrations of the RNAPs
as well as on those of some of the reactants and the products of the enzymatic
reactions catalyzed by RNAP. We suggest appropriate experimental systems and
techniques for testing our theoretical predictions.


  Strong experimental and theoretical evidence shows that transcription factors
and other specific DNA-binding proteins find their sites using a two-mode
search: alternating between 3D diffusion through the cell and 1D sliding along
the DNA. We consider the role spatial effects in the mechanism on two different
scales. First, we reconcile recent experimental findings by showing that the 3D
diffusion of the transcription factor is often local, i.e. the transcription
factor lands quite near its dissociation site. Second, we discriminate between
two types of searches: global searches and local searches. We show that these
searches differ significantly in average search time and the variability of
search time. Using experimentally measured parameter values, we also show that
1D and 3D search is not optimally balanced, leading to much larger estimates of
search time. Together, these results lead to a number of biological
implications including suggestions of how prokaryotes and eukaryotes achieve
rapid gene regulation and the relationship between the search mechanism and
noise in gene expression.


  We propose an order index, phi, which quantifies the notion of ``life at the
edge of chaos'' when applied to genome sequences. It maps genomes to a number
from 0 (random and of infinite length) to 1 (fully ordered) and applies
regardless of sequence length. The 786 complete genomic sequences in GenBank
were found to have phi values in a very narrow range, 0.037+/-0.027. We show
this implies that genomes are halfway towards being completely random, namely,
at the edge of chaos. We argue that this narrow range represents the
neighborhood of a fixed-point in the space of sequences, and genomes are driven
there by the dynamics of a robust, predominantly neutral evolution process.


  We study an evolutionary algorithm that locally adapts thresholds and wiring
in Random Threshold Networks, based on measurements of a dynamical order
parameter. A control parameter $p$ determines the probability of threshold
adaptations vs. link rewiring. For any $p < 1$, we find spontaneous symmetry
breaking into a new class of self-organized networks, characterized by a much
higher average connectivity $\bar{K}_{evo}$ than networks without threshold
adaptation ($p =1$). While $\bar{K}_{evo}$ and evolved out-degree distributions
are independent from $p$ for $p <1$, in-degree distributions become broader
when $p \to 1$, approaching a power-law. In this limit, time scale separation
between threshold adaptions and rewiring also leads to strong correlations
between thresholds and in-degree. Finally, evidence is presented that networks
converge to self-organized criticality for large $N$.


  The diversity in wealth and social status is present not only among humans,
but throughout the animal world. We account for this observation by generating
random variables that determ ine the social diversity of players engaging in
the prisoner's dilemma game. Here the term social diversity is used to address
extrinsic factors that determine the mapping of game pay offs to individual
fitness. These factors may increase or decrease the fitness of a player
depending on its location on the spatial grid. We consider different
distributions of extrin sic factors that determine the social diversity of
players, and find that the power-law distribution enables the best promotion of
cooperation. The facilitation of the cooperative str ategy relies mostly on the
inhomogeneous social state of players, resulting in the formation of
cooperative clusters which are ruled by socially high-ranking players that are
able to prevail against the defectors even when there is a large temptation to
defect. To confirm this, we also study the impact of spatially correlated
social diversity and find that coopera tion deteriorates as the spatial
correlation length increases. Our results suggest that the distribution of
wealth and social status might have played a crucial role by the evolution of
cooperation amongst egoistic individuals.


  To establish the relationship between locomotory behavior and dynamics of
neural circuits in the nematode C. elegans we combined molecular and
theoretical approaches. In particular, we quantitatively analyzed the motion of
C. elegans with defective synaptic GABA and acetylcholine transmission,
defective muscle calcium signaling, and defective muscles and cuticle
structures, and compared the data with our systems level circuit model. The
major experimental findings are: (i) anterior-to-posterior gradients of body
bending flex for almost all strains both for forward and backward motion, and
for neuronal mutants, also analogous weak gradients of undulatory frequency,
(ii) existence of some form of neuromuscular (stretch receptor) feedback, (iii)
invariance of neuromuscular wavelength, (iv) biphasic dependence of frequency
on synaptic signaling, and (v) decrease of frequency with increase of the
muscle time constant. Based on (i) we hypothesize that the Central Pattern
Generator (CPG) is located in the head both for forward and backward motion.
Points (i) and (ii) are the starting assumptions for our theoretical model,
whose dynamical patterns are qualitatively insensitive to the details of the
CPG design if stretch receptor feedback is sufficiently strong and slow. The
model reveals that stretch receptor coupling in the body wall is critical for
generation of the neuromuscular wave. Our model agrees with our behavioral
data(iii), (iv), and (v), and with other pertinent published data, e.g., that
frequency is an increasing function of muscle gap-junction coupling.


  We present a Brownian inchworm model of a self-propelled elastic dimer in the
absence of an external potential. Nonequilibrium noise together with a
stretch-dependent damping form the propulsion mechanism. Our model connects
three key nonequilibrium features -- position-velocity correlations, a nonzero
mean internal force, and a drift velocity. Our analytical results, including
striking current reversals, compare very well with numerical simulations. The
model unifies the propulsion mechanisms of DNA helicases, polar rods on a
vibrated surface, crawling keratocytes and Myosin VI. We suggest experimental
realizations and tests of the model.


  This paper describes a design of a molecular propagation system in molecular
communication. Molecular communication is a new communication paradigm where
biological and artificially-created nanomachines communicate over a short
distance using molecules. A molecular propagation system in molecular
communication directionally transports molecules from a sender to a receiver.
In the design described in this paper, protein filaments glide over immobilized
motor proteins along preconfigured microlithographic tracks, and the gliding
protein filaments carry and transport molecules from a sender to a receiver. In
the design, DNA hybridization is used to load and unload the molecules onto and
from the carriers at a sender and a receiver. In the design,
loading/transporting/unloading processes are autonomous and require no external
control.


  Complex biological systems are very robust to genetic and environmental
changes at all levels of organization. Many biological functions of Escherichia
coli metabolism can be sustained against single-gene or even multiple-gene
mutations by using redundant or alternative pathways. Thus, only a limited
number of genes have been identified to be lethal to the cell. In this regard,
the reaction-centric gene deletion study has a limitation in understanding the
metabolic robustness. Here, we report the use of flux-sum, which is the
summation of all incoming or outgoing fluxes around a particular metabolite
under pseudo-steady state conditions, as a good conserved property for
elucidating such robustness of E. coli from the metabolite point of view. The
functional behavior, as well as the structural and evolutionary properties of
metabolites essential to the cell survival, was investigated by means of a
constraints-based flux analysis under perturbed conditions. The essential
metabolites are capable of maintaining a steady flux-sum even against severe
perturbation by actively redistributing the relevant fluxes. Disrupting the
flux-sum maintenance was found to suppress cell growth. This approach of
analyzing metabolite essentiality provides insight into cellular robustness and
concomitant fragility, which can be used for several applications, including
the development of new drugs for treating pathogens.


  Thermal shape fluctuations of grafted microtubules were studied using high
resolution particle tracking of attached fluorescent beads. First mode
relaxation times were extracted from the mean square displacement in the
transverse coordinate. For microtubules shorter than 10 um, the relaxation
times were found to follow an L^2 dependence instead of L^4 as expected from
the standard wormlike chain model. This length dependence is shown to result
from a complex length dependence of the bending stiffness which can be
understood as a result of the molecular architecture of microtubules. For
microtubules shorter than 5 um, high drag coefficients indicate contributions
from internal friction to the fluctuation dynamics.


  For most of the important processes in DNA metabolism, a protein has to reach
a specific binding site on the DNA. The specific binding site may consist of
just a few base pairs while the DNA is usually several millions of base pairs
long. How does the protein search for the target site? What is the most
efficient mechanism for a successful search? Motivated by these fundamental
questions on intracellular biological processes, we have developed a model for
searching a specific site on a model DNA by a single protein. We have made a
comparative quantitative study of the efficiencies of sliding, inter-segmental
hoppings and detachment/re-attachments of the particle during its search for
the specific site on the DNA. We also introduce some new quantitative measures
of {\it efficiency} of a search process by defining a relevant quantity, which
can be measured in {\it in-vitro} experiments.


  The very insightful Trivers-Willard hypothesis, proposed in the early 1970s,
states that females in good physiological conditions are more likely to produce
male offspring, when the variance of reproductive success amongst males is
high. A number of studies, aimed at its experimental verification, have found
adequate supportive evidence in its favour. Theoretical investigations,
however, have been few, perhaps because formulating a population-genetic model
for describing the Trivers-Willard hypothesis turns out to be surprisingly
complex. The present study describes a minimal population genetic model to
explore one specific scenario, viz. how is the preference for a male offspring
by females in good condition altered when 'g', the proportion of such females
in the population changes from a low to a high value. As expected, when the
proportion of such females is low, i.e., for low values of 'g', the
Trivers-Willard (TW) strategy goes to fixation against the equal investment
strategy. This holds true up to gmax, a critical value of 'g', above which the
two strategies coexist, but the proportion of the TW strategy steadily
decreases as 'g' increases to unity. Similarly, when the effect of well-endowed
males attaining disproportionately high number of matings is more pronounced,
the TW strategy is more likely to go to fixation. Interestingly, the success of
the TW strategy has a complex dependence on the variance in the physiological
condition of females. If the difference in the two types of conditions is not
large, TW strategy is favoured, and its success is more likely as the
difference increases. However, beyond a critical value of the difference, the
TW strategy is found to be less and less likely to succeed as the difference
becomes larger. Possible reasons for these effects are discussed.


  Across all kingdoms of biological life, protein-coding genes exhibit unequal
usage of synonmous codons. Although alternative theories abound, translational
selection has been accepted as an important mechanism that shapes the patterns
of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns
of codon usage across 74 diverse bacteriophages that infect E. coli, P.
aeruginosa and L. lactis as their primary host. We introduce the concept of a
`genome landscape,' which helps reveal non-trivial, long-range patterns in
codon usage across a genome. We develop a series of randomization tests that
allow us to interrogate the significance of one aspect of codon usage, such a
GC content, while controlling for another aspect, such as adaptation to
host-preferred codons. We find that 33 phage genomes exhibit highly non-random
patterns in their GC3-content, use of host-preferred codons, or both. We show
that the head and tail proteins of these phages exhibit significant bias
towards host-preferred codons, relative to the non-structural phage proteins.
Our results support the hypothesis of translational selection on viral genes
for host-preferred codons, over a broad range of bacteriophages.


  E. Coli. dihydrofolate reductase (DHFR) undergoes conformational transitions
between the closed (CS) and occluded (OS) states which, respectively, describe
whether the active site is closed or occluded by the Met20 loop. A
sequence-based approach is used to identify a network of residues that
represents the allostery wiring diagram. We also use a self-organized polymer
model to monitor the kinetics of the CS->OS and the reverse transitions. a
sliding motion of Met20 loop is observed. The residues that facilitate the
Met20 loop motion are part of the network of residues that transmit allosteric
signals during the CS->OS transition.


  This work is concerned with the development of a well-founded, theoretically
justified, and least complicated metric for the classification of proteins with
reference to enzymes. As the signature of an enzyme family, a catalytic domain
is easily fingerprinted. Given that the classification problem has so far
seemed intractable, a classification schema derived from the catalytic domain
would be satisfying. Here I show that there exists a natural ab initio if
nonobvious basis to theorize that the catalytic domain of an enzyme is uniquely
informative about its regulation. This annotates its function. Based on this
hypothesis, a method that correctly classifies potassium ion channels into
their respective subfamilies is described. To put the principle on firmer
ground, extra validation was sought and obtained through co-evolutionary
analyses. The co-evolutionary analyses reveal a departure from the notion that
potassium ion channel proteins are functionally modular. This finding is
discussed in light of the prevailing notion of domain. These studies establish
that significant co-evolution of the catalytic domain of a gene with its
conjoint domain is a specialized, necessary process following fusion and
swapping events in evolution. Instances of this discovery are likely to be
found pervasive in protein science.


  We analyse a maximum-likelihood approach for combining phylogenetic trees
into a larger `supertree'. This is based on a simple exponential model of
phylogenetic error, which ensures that ML supertrees have a simple
combinatorial description (as a median tree, minimising a weighted sum of
distances to the input trees). We show that this approach to ML supertree
reconstruction is statistically consistent (it converges on the true species
supertree as more input trees are combined), in contrast to the widely-used MRP
method, which we show can be statistically inconsistent under the exponential
error model. We also show that this statistical consistency extends to an ML
approach for constructing species supertrees from gene trees. In this setting,
incomplete lineage sorting (due to coalescence rates of homologous genes being
lower than speciation rates) has been shown to lead to gene trees that are
frequently different from species trees, and this can confound efforts to
reconstruct the species phylogeny correctly.


  Species coexistence is one of the central themes in modern ecology.
Coexistence is a prerequisite of biological diversity. However, the question
arises how biodiversity can be reconciled with the statement of competition
theory, which asserts that competing species cannot coexist. To solve this
problem natural selection theory is rejected because it contradicts kinetic
models of interacting populations. Biological evolution is presented as a
process equivalent to a chemical reaction. The main point is that interactions
occur between self-replicating units. Under these assumptions biodiversity is
possible if and only if species are identical with respect to the patterns of
energy flow in which individuals are involved.


This paper has been withdrawn.


  In this paper we made a review of some papers about probabilistic regulatory
networks (PRN), in particular we introduce our concept of homomorphisms of PRN
with an example of projection of a regulatory network to a smaller one. We
apply the model PRN (or Probabilistic Boolean Network) to the immune system,
the PRN works with two functions. The model called ""The B/T-cells
interaction"" is Boolean, so we are really working with a Probabilistic Boolean
Network. Using Markov Chains we determine the state of equilibrium of the
immune response.


  We study how the dynamic equilibrium of the reversible protein-protein
binding network in yeast Saccharomyces cerevisiae responds to large changes in
abundances of individual proteins. The magnitude of shifts between free and
bound concentrations of their immediate and more distant neighbors in the
network is influenced by such factors as the network topology, the distribution
of protein concentrations among its nodes, and the average binding strength.
Our primary conclusion is that, on average, the effects of a perturbation are
strongly localized and exponentially decay with the network distance away from
the perturbed node, which explains why, despite globally connected topology,
individual functional modules in such networks are able to operate fairly
independently. We also found that under specific favorable conditions, realized
in a significant number of paths in the yeast network, concentration
perturbations can selectively propagate over considerable network distances (up
to four steps). Such "action-at-a-distance" requires high concentrations of
heterodimers along the path as well as low free (unbound) concentration of
intermediate proteins.


  A model is presented to describe the nucleotide and repeat addition
processivity by the telomerase. In the model, the processive nucleotide
addition is implemented on the basis of two requirements: One is that stem IV
loop stimulates the chemical reaction of nucleotide incorporation, and the
other one is the existence of an ssRNA-binding site adjacent to the polymerase
site that has a high affinity for the unpaired base of the template. The
unpairing of DNA:RNA hybrid after the incorporation of the nucleotide paired
with the last base on the template, which is the prerequisite for repeat
addition processivity, is caused by a force acting on the primer. The force is
resulted from the unfolding of stem III pseudoknot that is induced by the
swinging of stem IV loop towards the nucleotide-bound polymerase site. Based on
the model, the dynamics of processive nucleotide and repeat additions by
Tetrahymena telomerase are quantitatively studied, which give good explanations
to the previous experimental results. Moreover, some predictions are presented.
In particular, it is predicted that the repeat addition processivity is mainly
determined by the difference between the free energy required to disrupt the
DNA:RNA hybrid and that required to unfold the stem III pseudoknot, with the
large difference corresponding to a low repeat addition processivity while the
small one corresponding to a high repeat addition processivity.


  The critical properties of the Frank model of spontaneous chiral synthesis
are discussed by applying results from the field theoretic renormalization
group (RG). The long time and long wavelength features of this microscopic
reaction scheme belong to the same universality class as multi-colored directed
percolation processes. Thus, the following RG fixed points (FP) govern the
critical dynamics of the Frank model for d<4: one unstable FP that corresponds
to complete decoupling between the two enantiomers, a saddle-point that
corresponds to symmetric interspecies coupling, and two stable FPs that
individually correspond to unidirectional couplings between the two chiral
molecules. These latter two FPs are associated with the breakdown of mirror or
chiral symmetry. In this simplified model of molecular synthesis, homochirality
is a natural consequence of the intrinsic reaction noise in the critical
regime, which corresponds to extremely dilute chemical systems.


  Cooperation plays a key role in the evolution of complex systems. However,
the level of cooperation extensively varies with the topology of agent networks
in the widely used models of repeated games. Here we show that cooperation
remains rather stable by applying the reinforcement learning strategy adoption
rule, Q-learning on a variety of random, regular, small-word, scale-free and
modular network models in repeated, multi-agent Prisoners Dilemma and Hawk-Dove
games. Furthermore, we found that using the above model systems other long-term
learning strategy adoption rules also promote cooperation, while introducing a
low level of noise (as a model of innovation) to the strategy adoption rules
makes the level of cooperation less dependent on the actual network topology.
Our results demonstrate that long-term learning and random elements in the
strategy adoption rules, when acting together, extend the range of network
topologies enabling the development of cooperation at a wider range of costs
and temptations. These results suggest that a balanced duo of learning and
innovation may help to preserve cooperation during the re-organization of
real-world networks, and may play a prominent role in the evolution of
self-organizing, complex systems.


  With the continued improvement of sequencing technologies, the prospect of
genome-based medicine is now at the forefront of scientific research. To
realize this potential, however, we need a revolutionary sequencing method for
the cost-effective and rapid interrogation of individual genomes. This
capability is likely to be provided by a physical approach to probing DNA at
the single nucleotide level. This is in sharp contrast to current techniques
and instruments which probe, through chemical elongation, electrophoresis, and
optical detection, length differences and terminating bases of strands of DNA.
In this Colloquium we review several physical approaches to DNA detection that
have the potential to deliver fast and low-cost sequencing. Center-fold to
these approaches is the concept of nanochannels or nanopores which allow for
the spatial confinement of DNA molecules. In addition to their possible impact
in medicine and biology, the methods offer ideal test beds to study open
scientific issues and challenges in the relatively unexplored area at the
interface between solids, liquids, and biomolecules at the nanometer length
scale. We emphasize the physics behind these methods and ideas, critically
describe their advantages and drawbacks, and discuss future research
opportunities in this field.


  Tetracyclines are an old class of molecules that constitute a broad-spectrum
antibiotics. Since the first member of tetracycline family were isolated, the
clinical importance of these compounds as therapeutic and prophylactic agents
against a wide range of infections has stimulated efforts to define their mode
of action as inhibitors of bacterial reproduction. We used three SAR
methodologies for the analysis of biological activity of a set of 104
tetracycline compounds. Our calculation were carried out using the
semi-empirical Austin Method One (AM1) and Parametric Method 3 (PM3).
Electronic Indices Methodology (EIM), Principal Component Analysis (PCA) and
Artificial Neural Networks (ANN) were applied to the classification of 14 old
and 90 new proposed derivatives of tetracyclines. Our results make evident the
importance of EIM descriptors in pattern recognition and also show that the EIM
can be effectively used to predict the biological activity of Tetracyclines.


  In an equilibrium chemical reaction mixture, the number of molecules present
obeys a Poisson distribution. We ask when the same is true of the steady state
of a nonequilibrium reaction network and obtain an essentially complete answer.
In particular, we show that networks with certain topological features must
have a Poisson distribution, whatever the reaction rates. Such driven systems
also obey an analog of the fluctuation-dissipation theorem. Our results may be
relevant to biological systems and to the larger question of how equilibrium
concepts might apply to nonequilibrium systems.


  The evolution of languages closely resembles the evolution of haploid
organisms. This similarity has been recently exploited \cite{GA,GJ} to
construct language trees. The key point is the definition of a distance among
all pairs of languages which is the analogous of a genetic distance. Many
methods have been proposed to define these distances, one of this, used by
glottochronology, compute distance from the percentage of shared ``cognates''.
Cognates are words inferred to have a common historical origin, and subjective
judgment plays a relevant role in the identification process. Here we push
closer the analogy with evolutionary biology and we introduce a genetic
distance among language pairs by considering a renormalized Levenshtein
distance among words with same meaning and averaging on all the words contained
in a Swadesh list \cite{Sw}. The subjectivity of process is consistently
reduced and the reproducibility is highly facilitated. We test our method
against the Indo-European group considering fifty different languages and the
two hundred words of the Swadesh list for any of them. We find out a tree which
closely resembles the one published in \cite{GA} with some significant
differences.


  A major hindrance to studies of microbial diversity has been that the vast
majority of microbes cannot be cultured in the laboratory and thus are not
amenable to traditional methods of characterization. Environmental shotgun
sequencing (ESS) overcomes this hurdle by sequencing the DNA from the organisms
present in a microbial community. The interpretation of this metagenomic data
can be greatly facilitated by associating every sequence read with its source
organism. We report the development of CompostBin, a DNA composition-based
algorithm for analyzing metagenomic sequence reads and distributing them into
taxon-specific bins. Unlike previous methods that seek to bin assembled contigs
and often require training on known reference genomes, CompostBin has the
ability to accurately bin raw sequence reads without need for assembly or
training. It applies principal component analysis to project the data into an
informative lower-dimensional space, and then uses the normalized cut
clustering algorithm on this filtered data set to classify sequences into
taxon-specific bins. We demonstrate the algorithm's accuracy on a variety of
simulated data sets and on one metagenomic data set with known species
assignments. CompostBin is a work in progress, with several refinements of the
algorithm planned for the future.


  Cichlid fishes are one of the best model system for the study of evolution of
the species. Inspired by them, in this paper we simulated the splitting of a
single species into two separate ones via random mutations, with both
populations living together in sympatry, sharing the same habitat. We study the
ecological, mating and genetic conditions needed to reproduce the
polychromatism and polymorphism of three species of the Midas Cichlid species
complex. Our results show two scenarios for the A. Citrinellus speciation
process, one with and the other without disruptive natural selection. In the
first scenario, the ecological and genetic conditions are sufficient to create
two new species, while in the second the mating and genetic conditions must be
synchronized in order to control the velocity of genetic drift.


  In this paper we study $k$-noncrossing RNA structures with arc-length $\ge
3$, i.e. RNA molecules in which for any $i$, the nucleotides labeled $i$ and
$i+j$ ($j=1,2$) cannot form a bond and in which there are at most $k-1$
mutually crossing arcs. Let ${\sf S}_{k,3}(n)$ denote their number. Based on a
novel functional equation for the generating function $\sum_{n\ge 0}{\sf
S}_{k,3}(n)z^n$, we derive for arbitrary $k\ge 3$ exponential growth factors
and for $k=3$ the subexponential factor. Our main result is the derivation of
the formula ${\sf S}_{3,3}(n) \sim \frac{6.11170\cdot 4!}{n(n-1)...(n-4)}
4.54920^n$.


  The p53 protein is well-known for its tumour suppressor function. The
p53-MDM2 negative feedback loop constitutes the core module of a network of
regulatory interactions activated under cellular stress. In normal cells, the
level of p53 proteins is kept low by MDM2, i.e. MDM2 negatively regulates the
activity of p53. In the case of DNA damage,the p53-mediated pathways are
activated leading to cell cycle arrest and repair of the DNA. If repair is not
possible due to excessive damage, the p53-mediated apoptotic pathway is
activated bringing about cell death. In this paper, we give an overview of our
studies on the p53-MDM2 module and the associated pathways from a systems
biology perspective. We discuss a number of key predictions, related to some
specific aspects of cell cycle arrest and cell death, which could be tested in
experiments.


  We report on a theoretical study of point mutations effects on charge
transfer properties in the DNA sequence of the tumor-suppressor p53 gene. On
the basis of effective single-strand or double-strand tight-binding models
which simulate hole propagation along the DNA, a statistical analysis of charge
transmission modulations associated with all possible point mutations is
performed. We find that in contrast to non-cancerous mutations, mutation
hotspots tend to result in significantly weaker {\em changes of transmission
properties}. This suggests that charge transport could play a significant role
for DNA-repairing deficiency yielding carcinogenesis.


  Phylogenetic networks are a generalization of phylogenetic trees that allow
for the representation of non-treelike evolutionary events, like recombination,
hybridization, or lateral gene transfer. In this paper, we present and study a
new class of phylogenetic networks, called tree-child phylogenetic networks,
where every non-extant species has some descendant through mutation. We provide
an injective representation of these networks as multisets of vectors of
natural numbers, their path multiplicity vectors, and we use this
representation to define a distance on this class and to give an alignment
method for pairs of these networks. To the best of our knowledge, they are
respectively the first true distance and the first alignment method defined on
a meaningful class of phylogenetic networks strictly extending the class of
phylogenetic trees. Simple, polynomial algorithms for reconstructing a
tree-child phylogenetic network from its path multiplicity vectors, for
computing the distance between two tree-child phylogenetic networks, and for
aligning a pair of tree-child phylogenetic networks, are provided, and they
have been implemented as a Perl package and a Java applet, and they are
available at http://bioinfo.uib.es/~recerca/phylonetworks/mudistance


  Assuming the deleterious mutations in the Penna ageing model to affect mainly
the young ages, we get an enhanced mortality at very young age, followed by a
minimum of the mortality, and then the usual exponential increase of mortality
with age.


  In this paper a quantum mechanical description of the assembly/disassembly
process for microtubules is proposed. We introduce creation and annihilation
operators that raise or lower the microtubule length by a tubulin layer.
Following that, the Hamiltonian and corresponding equations of motion for the
quantum fields are derived that describe the dynamics of microtubules. These
Heisenberg-type equations are then transformed to semi-classical equations
using the method of coherent structures. We find that the dynamics of a
microtubule can be mathematically expressed via a cubic-quintic nonlinear
Schr\"{o}dinger (NLS) equation. We show that a vortex filament, a generic
solution of the NLS equation, exhibits linear growth/shrinkage in time as well
as temporal fluctuations about some mean value which is qualitatively similar
to the dynamic instability of microtubules.


  Push-pull networks are ubiquitous in signal transduction pathways in both
prokaryotic and eukaryotic cells. They allow cells to strongly amplify signals
via the mechanism of zero-order ultrasensitivity. In a push-pull network, two
antagonistic enzymes control the activity of a protein by covalent
modification. These enzymes are often uniformly distributed in the cytoplasm.
They can, however, also be colocalized in space, for instance, near the pole of
the cell. Moreover, it is increasingly recognized that these enzymes can also
be spatially separated, leading to gradients of the active form of the
messenger protein. Here, we investigate the consequences of the spatial
distributions of the enzymes for the amplification properties of push-pull
networks. Our calculations reveal that enzyme localization by itself can have a
dramatic effect on the gain. The gain is maximized when the two enzymes are
either uniformly distributed or colocalized in one region in the cell.
Depending on the diffusion constants, however, the sharpness of the response
can be strongly reduced when the enzymes are spatially separated. We discuss
how our predictions could be tested experimentally.


  In this short note we sketch the statistical physics framework of the replica
exchange technique when applied to molecular dynamics simulations. In
particular, we draw attention to generalized move sets that allow a variety of
optimizations as well as new applications of the method.


  We consider a general incompressible finite model protein of size M in its
environment, which we represent by a semiflexible copolymer consisting of amino
acid residues classified into only two species (H and P, see text) following
Lau and Dill. We allow various interactions between chemically unbonded
residues in a given sequence and the solvent (water), and exactly enumerate the
number of conformations W(E) as a function of the energy E on an infinite
lattice under two different conditions: (i) we allow conformations that are
restricted to be compact (known as Hamilton walk conformations), and (ii) we
allow unrestricted conformations that can also be non-compact. It is easily
demonstrated using plausible arguments that our model does not possess any
energy gap even though it is supposed to exhibit a sharp folding transition in
the thermodynamic limit. The enumeration allows us to investigate exactly the
effects of energetics on the native state(s), and the effect of small size on
protein thermodynamics and, in particular, on the differences between the
microcanonical and canonical ensembles. We find that the canonical entropy is
much larger than the microcanonical entropy for finite systems. We investigate
the property of self-averaging and conclude that small proteins do not
self-average. We also present results that (i) provide some understanding of
the energy landscape, and (ii) shed light on the free energy landscape at
different temperatures.


  The effects of demographic stochasticity in the long term behaviour of
endemic infectious diseases have been considered for long as a necessary
addition to an underlying deterministic theory. The latter would explain the
regular behaviour of recurrent epidemics, and the former the superimposed noise
of observed incidence patterns. Recently, a stochastic theory based on a
mechanism of resonance with internal noise has shifted the role of
stochasticity closer to the center stage, by showing that the major dynamic
patterns found in the incidence data can be explained as resonant fluctuations,
whose behaviour is largely independent of the amplitude of seasonal forcing,
and by contrast very sensitive to the basic epidemiological parameters. Here we
elaborate on that approach, by adding an ingredient which is missing in
standard epidemic models, the 'mixing network' through which infection may
propagate. We find that spatial correlations have a major effect in the
enhancement of the amplitude and the coherence of the resonant stochastic
fluctuations, providing the ordered patterns of recurrent epidemics, whose
period may differ significantly from that of the small oscillations around the
deterministic equilibrium. We also show that the inclusion of a more realistic,
time correlated, recovery profile instead of exponentially distributed
infectious periods may, even in the random-mixing limit, contribute to the same
effect.


  Considering a basic enzyme-catalysed reaction, in which the rate of input of
the substrate varies periodically in time, we give a necessary and sufficient
condition for the existence of a periodic solution of the reaction equations.
The proof employs the Leray-Schauder degree, applied to an appropriately
constructed homotopy.


  The determination of a patient's DNA sequence can, in principle, reveal an
increased risk to fall ill with particular diseases [1,2] and help to design
"personalized medicine" [3]. Moreover, statistical studies and comparison of
genomes [4] of a large number of individuals are crucial for the analysis of
mutations [5] and hereditary diseases, paving the way to preventive medicine
[6]. DNA sequencing is, however, currently still a vastly time-consuming and
very expensive task [4], consisting of pre-processing steps, the actual
sequencing using the Sanger method, and post-processing in the form of data
analysis [7]. Here we propose a new approach that relies on functionalized
nanopore-embedded electrodes to achieve an unambiguous distinction of the four
nucleic acid bases in the DNA sequencing process. This represents a significant
improvement over previously studied designs [8,9] which cannot reliably
distinguish all four bases of DNA. The transport properties of the setup
investigated by us, employing state-of-the-art density functional theory
together with the non-equilibrium Green's Function method, leads to current
responses that differ by at least one order of magnitude for different bases
and can thus provide a much more robust read-out of the base sequence. The
implementation of our proposed setup could thus lead to a viable protocol for
rapid DNA sequencing with significant consequences for the future of genome
related research in particular and health care in general.


  In this article, we study the kinetics of reversible ligand binding to
receptors on a spherical cell surface using a self-consistent stochastic
theory. Binding, dissociation, diffusion and rebinding of ligands are
incorporated into the theory in a systematic manner. We derive explicitly the
time evolution of the ligand-bound receptor fraction p(t) in various regimes .
Contrary to the commonly accepted view, we find that the well-known
Berg-Purcell scaling for the association rate is modified as a function of
time. Specifically, the effective on-rate changes non-monotonically as a
function of time and equals the intrinsic rate at very early as well as late
times, while being approximately equal to the Berg-Purcell value at
intermediate times. The effective dissociation rate, as it appears in the
binding curve or measured in a dissociation experiment, is strongly modified by
rebinding events and assumes the Berg-Purcell value except at very late times,
where the decay is algebraic and not exponential. In equilibrium, the ligand
concentration everywhere in the solution is the same and equals its spatial
mean, thus ensuring that there is no depletion in the vicinity of the cell.
Implications of our results for binding experiments and numerical simulations
of ligand-receptor systems are also discussed.


  We analyze a model of interacting agents (e.g. prebiotic chemical species)
which are represended by nodes of a network, whereas their interactions are
mapped onto directed links between these nodes. On a fast time scale, each
agent follows an eigendynamics based on catalytic support from other nodes,
whereas on a much slower time scale the network evolves through selection and
mutation of its nodes-agent. In the first part of the paper, we explain the
dynamics of the model by means of characteristic snapshots of the network
evolution and confirm earlier findings on crashes an recoveries in the network
structure. In the second part, we focus on the aggregate behavior of the
network dynamics. We show that the disruptions in the network structure are
smoothed out, so that the average evolution can be described by a growth regime
followed by a saturation regime, without an initial random regime. For the
saturation regime, we obtain a logarithmic scaling between the average
connectivity per node $\mean{l}_{s}$ and a parameter $m$, describing the
average incoming connectivity, which is independent of the system size $N$.


  In this work we suggest a simple mathematical model for the dynamics of the
population of children and adolescents without problematic behavior (criminal
activities etc.). This model represents a typical population growth equation
but with time dependent (linearly decreasing) population growth coefficient.
Given equation admits definition of the half-life time of the non-problematic
children behavior as well as a criterion for estimation of the social
regulation of the children behavior.


  The use of mathematical models has helped to shed light on countless
phenomena in chemistry and biology. Often, though, one finds that systems of
interest in these fields are dauntingly complex. In this paper, we attempt to
synthesize and expand upon the body of mathematical results pertaining to the
theory of multiple equilibria in chemical reaction networks (CRNs), which has
yielded surprising insights with minimal computational effort. Our central
focus is a recent, cycle-based theorem by Gheorghe Craciun and Martin Feinberg,
which is of significant interest in its own right and also serves, in a
somewhat restated form, as the basis for a number of fruitful connections among
related results.


  Biodiversity is essential to the viability of ecological systems. Species
diversity in ecosystems is promoted by cyclic, non-hierarchical interactions
among competing populations. Such non-transitive relations lead to an evolution
with central features represented by the `rock-paper-scissors' game, where rock
crushes scissors, scissors cut paper, and paper wraps rock. In combination with
spatial dispersal of static populations, this type of competition results in
the stable coexistence of all species and the long-term maintenance of
biodiversity. However, population mobility is a central feature of real
ecosystems: animals migrate, bacteria run and tumble. Here, we observe a
critical influence of mobility on species diversity. When mobility exceeds a
certain value, biodiversity is jeopardized and lost. In contrast, below this
critical threshold all subpopulations coexist and an entanglement of travelling
spiral waves forms in the course of temporal evolution. We establish that this
phenomenon is robust, it does not depend on the details of cyclic competition
or spatial environment. These findings have important implications for
maintenance and evolution of ecological systems and are relevant for the
formation and propagation of patterns in excitable media, such as chemical
kinetics or epidemic outbreaks.


  The Battle of the Sexes describes asymmetric conflicts in mating behavior of
males and females. Males can be philanderer or faithful, while females are
either fast or coy, leading to a cyclic dynamics. The adjusted replicator
equation predicts stable coexistence of all four strategies. In this situation,
we consider the effects of fluctuations stemming from a finite population size.
We show that they unavoidably lead to extinction of two strategies in the
population. However, the typical time until extinction occurs strongly prolongs
with increasing system size. In the meantime, a quasi-stationary probability
distribution forms that is anomalously flat in the vicinity of the coexistence
state. This behavior originates in a vanishing linear deterministic drift near
the fixed point. We provide numerical data as well as an analytical approach to
the mean extinction time and the quasi-stationary probability distribution.


  We analytically derive the lower bound of the total conformational energy of
a protein structure by assuming that the total conformational energy is well
approximated by the sum of sequence-dependent pairwise contact energies. The
condition for the native structure achieving the lower bound leads to the
contact energy matrix that is a scalar multiple of the native contact matrix,
i.e., the so-called Go potential. We also derive spectral relations between
contact matrix and energy matrix, and approximations related to one-dimensional
protein structures. Implications for protein structure prediction are
discussed.


  Occupancy of new habitats through dispersion is a central process in nature.
In particular, long range dispersal is involved in the spread of species and
epidemics, although it has not been previously related with cancer invasion, a
process that involves spread to new tissues. We show that the early spread of
cancer cells is similar to the species individuals spread and that both
processes are represented by a common spatio-temporal signature, characterized
by a particular fractal geometry of the boundaries of patches generated, and a
power law-scaled, disrupted patch size distribution. We show that both
properties are a direct result of long-distance dispersal, and that they
reflect homologous ecological processes of population self-organization. Our
results are significant for processes involving long-range dispersal like
biological invasions, epidemics and cancer metastasis.


  Inference of evolutionary trees and rates from biological sequences is
commonly performed using continuous-time Markov models of character change. The
Markov process evolves along an unknown tree while observations arise only from
the tips of the tree. Rate heterogeneity is present in most real data sets and
is accounted for by the use of flexible mixture models where each site is
allowed its own rate. Very little has been rigorously established concerning
the identifiability of the models currently in common use in data analysis,
although non-identifiability was proven for a semi-parametric model and an
incorrect proof of identifiability was published for a general parametric model
(GTR+Gamma+I). Here we prove that one of the most widely used models
(GTR+Gamma) is identifiable for generic parameters, and for all parameter
choices in the case of 4-state (DNA) models. This is the first proof of
identifiability of a phylogenetic model with a continuous distribution of
rates.


  An endogenous molecular-cellular network for both normal and abnormal
functions is assumed to exist. This endogenous network forms a nonlinear
stochastic dynamical system, with many stable attractors in its functional
landscape. Normal or abnormal robust states can be decided by this network in a
manner similar to the neural network. In this context cancer is hypothesized as
one of its robust intrinsic states. This hypothesis implies that a nonlinear
stochastic mathematical cancer model is constructible based on available
experimental data and its quantitative prediction is directly testable. Within
such model the genesis and progression of cancer may be viewed as stochastic
transitions between different attractors. Thus it further suggests that
progressions are not arbitrary. Other important issues on cancer, such as
genetic vs epigenetics, double-edge effect, dormancy, are discussed in the
light of present hypothesis. A different set of strategies for cancer
prevention, cure, and care, is therefore suggested.


  "Leaping" methods show great promise for significantly accelerating
stochastic simulations of complex biochemical reaction networks. However, few
practical applications of leaping have appeared in the literature to date.
Here, we address this issue using the "partitioned leaping algorithm" (PLA)
[L.A. Harris and P. Clancy, J. Chem. Phys. 125, 144107 (2006)], a
recently-introduced multiscale leaping approach. We use the PLA to investigate
stochastic effects in two model biochemical reaction networks. The networks
that we consider are simple enough so as to be accessible to our intuition but
sufficiently complex so as to be generally representative of real biological
systems. We demonstrate how the PLA allows us to quantify subtle effects of
stochasticity in these systems that would be difficult to ascertain otherwise
as well as not-so-subtle behaviors that would strain commonly-used "exact"
stochastic methods. We also illustrate bottlenecks that can hinder the approach
and exemplify and discuss possible strategies for overcoming them. Overall, our
aim is to aid and motivate future applications of leaping by providing stark
illustrations of the benefits of the method while at the same time elucidating
obstacles that are often encountered in practice.


  Intracellular pathogens such as Listeria monocytogenes and Rickettsia
rickettsii move within a host cell by polymerizing a comet-tail of actin fibers
that ultimately pushes the cell forward. This dense network of cross-linked
actin polymers typically exhibits a striking curvature that causes bacteria to
move in gently looping paths. Theoretically, tail curvature has been linked to
details of motility by considering force and torque balances from a finite
number of polymerizing filaments. Here we track beads coated with a prokaryotic
activator of actin polymerization in three dimensions to directly quantify the
curvature and torsion of bead motility paths. We find that bead paths are more
likely to have low rather than high curvature at any given time. Furthermore,
path curvature changes very slowly in time, with an autocorrelation decay time
of 200 seconds. Paths with a small radius of curvature, therefore, remain so
for an extended period resulting in loops when confined to two dimensions. When
allowed to explore a 3D space, path loops are less evident. Finally, we
quantify the torsion in the bead paths and show that beads do not exhibit a
significant left- or right-handed bias to their motion in 3D. These results
suggest that paths of actin-propelled objects may be attributed to slow changes
in curvature rather than a fixed torque.


  The architecture of biological networks has been reported to exhibit high
level of modularity, and to some extent, topological modules of networks
overlap with known functional modules. However, how the modular topology of the
molecular network affects the evolution of its member proteins remains unclear.
In this work, the functional and evolutionary modularity of Homo sapiens (H.
sapiens) metabolic network were investigated from a topological point of view.
Network decomposition shows that the metabolic network is organized in a highly
modular core-periphery way, in which the core modules are tightly linked
together and perform basic metabolism functions, whereas the periphery modules
only interact with few modules and accomplish relatively independent and
specialized functions. Moreover, over half of the modules exhibit
co-evolutionary feature and belong to specific evolutionary ages. Peripheral
modules tend to evolve more cohesively and faster than core modules do. The
correlation between functional, evolutionary and topological modularity
suggests that the evolutionary history and functional requirements of metabolic
systems have been imprinted in the architecture of metabolic networks. Such
systems level analysis could demonstrate how the evolution of genes may be
placed in a genome-scale network context, giving a novel perspective on
molecular evolution.


  In a recent paper [C. Marr, M. Mueller-Linow, and M.-T. Huett, Phys. Rev. E
75, 041917 (2007)] we discuss the pronounced potential of real metabolic
network topologies, compared to randomized counterparts, to regularize complex
binary dynamics. In their comment [P. Holme and M. Huss, arXiv:0705.4084v1],
Holme and Huss criticize our approach and repeat our study with more realistic
dynamics, where stylized reaction kinetics are implemented on sets of pairwise
reactions. The authors find no dynamic difference between the reaction sets
recreated from the metabolic networks and randomized counterparts. We reproduce
the author's observation and find that their algorithm leads to a dynamical
fragmentation and thus eliminates the topological information contained in the
graphs. Hence, their approach cannot rule out a connection between the topology
of metabolic networks and the ubiquity of steady states.


  In this work we propose a novel algorithmic strategy that allows for an
efficient characterization of the whole set of stable fluxes compatible with
the metabolic constraints. The algorithm, based on the well-known Bethe
approximation, allows the computation in polynomial time of the volume of a non
full-dimensional convex polytope in high dimensions. The result of our
algorithm match closely the prediction of Monte Carlo based estimations of the
flux distributions of the Red Blood Cell metabolic network but in incomparably
shorter time. We also analyze the statistical properties of the average fluxes
of the reactions in the E-Coli metabolic network and finally to test the effect
of gene knock-outs on the size of the solution space of the E-Coli central
metabolism.


  . Genechip oligonucleotide microarrays have been used widely for
transcriptional profiling of a large number of genes in a given paradigm. Gene
expression estimation precedes biological inference and is given as a complex
combination of atomic entities on the array called probes. These probe
intensities are further classified into perfect-match (PM) and mis-match (MM)
probes. While former is a measure of specific binding, the lat-ter is a measure
of non-specific binding. The behavior of the MM probes has especially proven to
be elusive. The present study investigates qualita-tive similarities in the
distributional signatures and local correlation struc-tures/patchiness between
the PM and MM probe intensities. These qualita-tive similarities are
established on publicly available microarrays generated across laboratories
investigating the same paradigm. Persistence of these similarities across raw
as well as background subtracted probe intensities is also investigated. The
results presented raise fundamental concerns in inter-preting Genechip
oligonucleotide microarray data.


  This paper develops a simplified set of models describing asexual and sexual
replication in unicel- lular diploid organisms. The models assume organisms
whose genomes consist of two chromosomes, where each chromosome is assumed to
be functional if it is equal to some master sequence $ \sigma_0 $, and
non-functional otherwise. The first-order growth rate constant, or fitness, of
an organism, is determined by whether it has zero, one, or two functional
chromosomes in its genome. For a population replicating asexually, a given cell
replicates both of its chromosomes, and splits its genetic material evenly
between the two cells. For a population replicating sexually, a given cell
first divides into two haploids, which enter a haploid pool, fuse into
diploids, and then divide via the normal mitotic process. Haploid fusion is
modeled as a second-order rate process. When the cost for sex is small, as
measured by the ratio of the characteristic haploid fusion time to the
characteristic growth time, we find that sexual replication with random haploid
fusion leads to a greater mean fitness for the population than a purely asexual
strategy. However, independently of the cost for sex, we find that sexual
replication with a selective mating strategy leads to a higher mean fitness
than the random mating strategy. This result is based on the assumption that a
selective mating strategy does not have any additional time or energy costs
over the random mating strategy, an assumption that is discussed in the paper.
The results of this paper are consistent with previous studies suggesting that
sex is favored at intermediate mutation rates, for slowly replicating
organisms, and at high population densities.


  Conformational fluctuations are believed to play an important role in the
process by which transcription factor proteins locate and bind their target
site on the genome of a bacterium. Using a simple model, we show that the
binding time can be minimized, under selective pressure, by adjusting the
spectrum of conformational states so that the fraction of time spent in more
mobile conformations is matched with the target recognition rate. The
associated optimal binding time is then within an order of magnitude of the
limiting binding time imposed by thermodynamics, corresponding to an idealized
protein with instant target recognition. Numerical estimates suggest that
typical bacteria operate in this regime of optimized conformational
fluctuations.


  Based on our experience in kinetic modeling of coupled multiple metabolic
pathways we propose a generic rate equation for the dynamical modeling of
metabolic kinetics. Its symmetric form makes the kinetic parameters (or
functions) easy to relate to values in database and to use in computation. In
addition, such form is workable to arbitrary number of substrates and products
with different stoichiometry. We explicitly show how to obtain such rate
equation exactly for various binding mechanisms. Hence the proposed rate
equation is formally rigorous. Various features of such a generic rate equation
are discussed. For irreversible reactions, the product inhibition which
directly arise from enzymatic reaction is eliminated in a natural way. We also
discuss how to include the effects of modifiers and cooperativity.


  We study a continuous-time dynamical system that models the evolving
distribution of genotypes in an infinite population where genomes may have
infinitely many or even a continuum of loci, mutations accumulate along
lineages without back-mutation, added mutations reduce fitness, and
recombination occurs on a faster time scale than mutation and selection. Some
features of the model, such as existence and uniqueness of solutions and
convergence to the dynamical system of an approximating sequence of discrete
time models, were presented in earlier work by Evans, Steinsaltz, and Wachter
for quite general selective costs. Here we study a special case where the
selective cost of a genotype with a given accumulation of ancestral mutations
from a wild type ancestor is a sum of costs attributable to each individual
mutation plus successive interaction contributions from each $k$-tuple of
mutations for $k$ up to some finite ``degree''. Using ideas from complex
chemical reaction networks and a novel Lyapunov function, we establish that the
phenomenon of mutation-selection balance occurs for such selection costs under
mild conditions. That is, we show that the dynamical system has a unique
equilibrium and that it converges to this equilibrium from all initial
conditions.


  Biological functions in each animal cell depend on coordinated operations of
a wide variety of molecular motors. Some of the these motors transport cargo to
their respective destinations whereas some others are mobile workshops which
synthesize macromolecules while moving on their tracks. Some other motors are
designed to function as packers and movers. All these motors require input
energy for performing their mechanical works and operate under conditions far
from thermodynamic equilibrium. The typical size of these motors and the forces
they generate are of the order of nano-meters and pico-Newtons, respectively.
They are subjected to random bombardments by the molecules of the surrounding
aqueous medium and, therefore, follow noisy trajectories. Because of their
small inertia, their movements in the viscous intracellular space exhibits
features that are characteristics of hydrodynamics at low Reynold's number. In
this article we discuss how theoretical modeling and computer simulations of
these machines by physicists are providing insight into their mechanisms which
engineers can exploit to design and control artificial nano-motors.


  Adaptive networks appear in many biological applications. They combine
topological evolution of the network with dynamics in the network nodes.
Recently, the dynamics of adaptive networks has been investigated in a number
of parallel studies from different fields, ranging from genomics to game
theory. Here we review these recent developments and show that they can be
viewed from a unique angle. We demonstrate that all these studies are
characterized by common themes, most prominently: complex dynamics and robust
topological self-organization based on simple local rules.


  Genetic recombination can produce heterogeneous phylogenetic histories within
a set of homologous genes. Delineating recombination events is important in the
study of molecular evolution, as inference of such events provides a clearer
picture of the phylogenetic relationships among different gene sequences or
genomes. Nevertheless, detecting recombination events can be a daunting task,
as the performance of different recombinationdetecting approaches can vary,
depending on evolutionary events that take place after recombination. We
recently evaluated the effects of postrecombination events on the prediction
accuracy of recombination-detecting approaches using simulated nucleotide
sequence data. The main conclusion, supported by other studies, is that one
should not depend on a single method when searching for recombination events.
In this paper, we introduce a two-phase strategy, applying three statistical
measures to detect the occurrence of recombination events, and a Bayesian
phylogenetic approach in delineating breakpoints of such events in nucleotide
sequences. We evaluate the performance of these approaches using simulated
data, and demonstrate the applicability of this strategy to empirical data. The
two-phase strategy proves to be time-efficient when applied to large datasets,
and yields high-confidence results.


  Numerical models indicate that collective animal behaviour may emerge from
simple local rules of interaction among the individuals. However, very little
is known about the nature of such interaction, so that models and theories
mostly rely on aprioristic assumptions. By reconstructing the three-dimensional
position of individual birds in airborne flocks of few thousands members, we
prove that the interaction does not depend on the metric distance, as most
current models and theories assume, but rather on the topological distance. In
fact, we discover that each bird interacts on average with a fixed number of
neighbours (six-seven), rather than with all neighbours within a fixed metric
distance. We argue that a topological interaction is indispensable to maintain
flock's cohesion against the large density changes caused by external
perturbations, typically predation. We support this hypothesis by numerical
simulations, showing that a topological interaction grants significantly higher
cohesion of the aggregation compared to a standard metric one.


  We discuss the relevance of studying ecology within the framework of
Complexity Science from a statistical mechanics approach. Ecology is concerned
with understanding how systems level properties emerge out of the multitude of
interactions amongst large numbers of components, leading to ecosystems that
possess the prototypical characteristics of complex systems. We argue that
statistical mechanics is at present the best methodology available to obtain a
quantitative description of complex systems, and that ecology is in urgent need
of ``integrative'' approaches that are quantitative and non-stationary. We
describe examples where combining statistical mechanics and ecology has led to
improved ecological modelling and, at the same time, broadened the scope of
statistical mechanics.


  The transfer of genetic materials across species (lateral genetic transfer,
LGT) contributes to genomic and physiological innovation in prokaryotes. The
extent of LGT in prokaryotes has been examined in a number of studies, but the
unit of transfer has not been studied in a rigorous manner. Using a rigorous
phylogenetic approach, we analysed the units of LGT within families of
single-copy genes obtained from 144 fully sequenced prokaryote genomes. A total
of 30.3% of these gene families show evidence of LGT. We found that the
transfer of gene fragments has been more frequent than the transfer of entire
genes, suggesting the extent of LGT has been underestimated. We found little
functional bias between within-gene (fragmentary) and whole-gene
(non-fragmentary) genetic transfer, but non-fragmentary transfer has been more
frequent into pathogens than into non-pathogens. As gene families that contain
probable paralogs were excluded from the current study, our results may still
underestimate the extent of LGT; nonetheless this is the most-comprehensive
study to date of the unit of LGT among prokaryote genomes.


  Genomes evolve as modules. In prokaryotes (and some eukaryotes), genetic
material can be transferred between species and integrated into the genome via
homologous or illegitimate recombination. There is little reason to imagine
that the units of transfer correspond to entire genes; however, such units have
not been rigorously characterized. We examined fragmentary genetic transfers in
single-copy gene families from 144 prokaryotic genomes and found that
breakpoints are located significantly closer to the boundaries of genomic
regions that encode annotated structural domains of proteins than expected by
chance, particularly when recombining sequences are more divergent. This
correlation results from recombination events themselves and not from
differential nucleotide substitution. We report the first systematic study
relating genetic recombination to structural features at the protein level.


  The widespread presence of antibiotic resistance and virulence among
Staphylococcus isolates has been attributed to lateral genetic transfer (LGT)
between different strains or species. However, there has been very little study
of the extent of LGT in Staphylococcus species using a phylogenetic approach,
particularly of the units of such genetic transfer. Here we report the first
systematic study of the units of genetic transfer in 13 Staphylococcus genomes,
using a rigorous phylogenetic approach. We found clear evidence of LGT in 26.1%
of the 1354 homologous gene families examined, and possibly more in another
17.9% of the total families. Within-gene and whole-gene transfer contribute
almost equally to the discordance of these gene families against a reference
phylogeny. Comparing genetic transfer in single-copy and in multi-copy gene
families, we found little functional bias in cases of within-gene (fragmentary)
genetic transfer but substantial functional bias in cases of whole-gene
(non-fragmentary) genetic transfer, and we observed a higher frequency of LGT
in multi-copy gene families. Our results demonstrate that LGT and gene
duplication play an important part among the factors that contribute to
functional innovation in staphylococcal genomes.


  Previous research has shown a strong correlation of protein folding rates to
the native state geometry, yet a complete explanation for this dependence is
still lacking. Here we study the rate-geometry relationship with a simple
statistical physics model, and focus on two classes of model geometries,
representing ideal parallel and antiparallel structures. We find that the
logarithm of the rate shows an almost perfect linear correlation with the
"absolute contact order", but the slope depends on the particular class
considered. We discuss these findings in the light of experimental results.


  The authors address the problem of downhill protein folding in the framework
of a simple statistical mechanical model, which allows an exact solution for
the equilibrium and a semianalytical treatment of the kinetics. Focusing on
protein 1BBL, a candidate for downhill folding behavior, and comparing it to
the WW domain of protein PIN1, a two-state folder of comparable size, the
authors show that there are qualitative differences in both the equilibrium and
kinetic properties of the two molecules. However, the barrierless scenario
which would be expected if 1BBL were a true downhill folder is observed only at
low enough temperature.


  One of the mainstays of the controversial "rare Earth" hypothesis is the
"Goldilocks problem" regarding various parameters describing a habitable
planet, partially involving the role of mass extinctions and other catastrophic
processes in biological evolution. Usually, this is construed as support for
the uniqueness of the Earth's biosphere and intelligent human life. Here I
argue that this is a misconstrual and that, on the contrary,
observation-selection effects, when applied to catastrophic processes, make it
very difficult for us to discern whether the terrestrial biosphere and
evolutionary processes which created it are exceptional in the Milky Way or
not. In particular, an anthropic overconfidence bias related to the temporal
asymmetry of evolutionary processes appears when we try to straightforwardly
estimate catastrophic risks from the past records on Earth. This agnosticism,
in turn, supports the validity and significance of practical astrobiological
and SETI research.


  In simulations of sexual reproduction with diploid individuals, we introduce
that female haploid gametes recognize one specific allele of the genomes as a
marker of the male haploid gametes. They fuse to zygotes preferrably with male
gametes having a different marker than their own. This gamete recognition
enhances the advantage of complementary bit-strings in the simulated diploid
individuals, at low recombination rates. Thus with rare recombinations the
bit-string evolve to be complementary; with recombination rate above about 0.1
instead they evolve under Darwinian purification selection, with few bits
mutated.


  We introduce and analyze a waiting time model for the accumulation of genetic
changes. The continuous time conjunctive Bayesian network is defined by a
partially ordered set of mutations and by the rate of fixation of each
mutation. The partial order encodes constraints on the order in which mutations
can fixate in the population, shedding light on the mutational pathways
underlying the evolutionary process. We study a censored version of the model
and derive equations for an EM algorithm to perform maximum likelihood
estimation of the model parameters. We also show how to select the maximum
likelihood poset. The model is applied to genetic data from different cancers
and from drug resistant HIV samples, indicating implications for diagnosis and
treatment.


  The evolutionary force of recombination is lacking in asexually reproducing
populations. As a consequence, the population can suffer an irreversible
accumulation of deleterious mutations, a phenomenon known as Muller's ratchet.
We formulate discrete and continuous time versions of Muller's ratchet.
Inspired by Haigh's (1978) analysis of a dynamical system which arises in the
limit of large populations, we identify the parameter gamma =
N*lambda/(Ns*log(N*lambda)) as most important for the speed of accumulation of
deleterious mutations. Here N is population size, s is the selection
coefficient and lambda is the deleterious mutation rate. For large parts of the
parameter range, measuring time in units of size N, deleterious mutations
accumulate according to a power law in N*lambda with exponent gamma if
gamma>0.5. For gamma<0.5 mutations cannot accumulate. We obtain diffusion
approximations for three different parameter regimes, depending on the speed of
the ratchet. Our approximations shed new light on analyses of Stephan et al.
(1993) and Gordo & Charlesworth (2000). The heuristics leading to the
approximations are supported by simulations.


  The goal of this paper is to outline a scenario of emerging stochasticity in
high-dimensional highly nonlinear systems, such as genetic regulatory networks
(GRN). We focus attention on the fact that in such systems confluence of all
the factors necessary for gene expression is a comparatively rare event, and
only massive redundancy makes such events sufficiently frequent. An immediate
consequence of this rareness is burstiness in mRNA and protein copy numbers, a
well known experimentally observed effect. We introduce the concept of
stochastic cooperativity and show that this phenomenon is a natural consequence
of high dimensionality coupled with highly nonlinearity of a dynamical system.
In mathematical terms, burstiness is associated with heavy-tailed probability
distributions of stochastic processes describing the dynamics of the system.
The sequence of stochastic cooperativity events allows for transition from
continuous deterministic dynamics expressed in terms of ordinary differential
equations (ODE) to discrete stochastic dynamics expressed in terms of Langevin
and Fokker-Plank equations. We demonstrate also that high-dimensional nonlinear
systems, even in the absence of explicit mechanisms for suppressing inherent
instability, may nevertheless reside in a state of stationary pseudo-random
fluctuations which for all practical purposes may be regarded as stochastic
process. This type of stochastic behavior is an inherent property of such
systems and requires neither an external random force, nor highly specialized
conditions of bistability.


  Phylogenetic invariants are certain polynomials in the joint probability
distribution of a Markov model on a phylogenetic tree. Such polynomials are of
theoretical interest in the field of algebraic statistics and they are also of
practical interest--they can be used to construct phylogenetic trees. This
paper is a self-contained introduction to the algebraic, statistical, and
computational challenges involved in the practical use of phylogenetic
invariants. We survey the relevant literature and provide some partial answers
and many open problems.


  Biological activity gives rise to non-equilibrium fluctuations in the
cytoplasm of cells; however, there are few methods to directly measure these
fluctuations. Using a reconstituted actin cytoskeleton, we show that the
bending dynamics of embedded microtubules can be used to probe local stress
fluctuations. We add myosin motors that drive the network out of equilibrium,
resulting in an increased amplitude and modified time-dependence of microtubule
bending fluctuations. We show that this behavior results from step-like forces
on the order of 10 pN driven by collective motor dynamics.


  We generalize the Poland-Scheraga model to consider DNA denaturation in the
presence of an external stretching force. We demonstrate the existence of a
force-induced DNA denaturation transition and obtain the temperature-force
phase diagram. The transition is determined by the loop exponent $c$ for which
we find the new value $c=4\nu-1/2$ such that the transition is second order
with $c=1.85<2$ in $d=3$. We show that a finite stretching force $F$
destabilizes DNA, corresponding to a lower melting temperature $T(F)$, in
agreement with single-molecule DNA stretching experiments.


  We study the Susceptible-Infected-Susceptible model of the spread of an
endemic infection. We calculate an exact expression for the mean number of
transmissions for all values of the population and the infectivity. We derive
the large-N asymptotic behavior for the infectivitiy below, above, and in the
critical region. We obtain an analytical expression for the probability
distribution of the number of transmissions, n, in the critical region. We show
that this distribution has a $n^3/2$ singularity for small n and decays
exponentially for large n. The exponent decreases with the distance from
threshold, diverging to infinity far below and approaching zero far above.