  The statistics of steps and dwell times in reversible molecular motors differ
from those of cycle completion in enzyme kinetics. The reason is that a step is
only one of several transitions in the mechanochemical cycle. As a result,
theoretical results for cycle completion in enzyme kinetics do not apply to
stepping data. To allow correct parameter estimation, and to guide data
analysis and experiment design, a theoretical treatment is needed that takes
this observation into account. In this paper, we model the distribution of
dwell times and number of forward and backward steps using first passage
processes, based on the assumption that forward and backward steps correspond
to different directions of the same transition. We extend recent results for
systems with a single cycle and consider the full dwell time distributions as
well as models with multiple pathways, detectable substeps, and detachments.
Our main results are a symmetry relation for the dwell time distributions in
reversible motors, and a relation between certain relative step frequencies and
the free energy per cycle. We demonstrate our results by analyzing recent
stepping data for a bacterial flagellar motor, and discuss the implications for
the efficiency and reversibility of the force-generating subunits. Key words:
motor proteins; single molecule kinetics; enzyme kinetics; flagellar motor;
Markov process; non-equilibrium fluctuations.


  We consider how transfer of genetic information between individuals
influences the phase diagram and mean fitness of both the Eigen and the
parallel, or Crow-Kimura, models of evolution. In the absence of genetic
transfer, these physical models of evolution consider the replication and point
mutation of the genomes of independent individuals in a large population. A
phase transition occurs, such that below a critical mutation rate an
identifiable quasispecies forms. We generalize these models of quasispecies
evolution to include horizontal gene transfer. We show how transfer of genetic
information changes the phase diagram and mean fitness and introduces
metastability in quasispecies theory, via an analytic field theoretic mapping.


  RNA secondary structure is an important computational model to understand how
genetic variation maps into phenotypic (structural) variation. Evolutionary
innovation in RNA structures is facilitated by neutral networks, large
connected sets of RNA sequences that fold into the same structure. Our work
extends and deepens previous studies on neutral networks. First, we show that
even the 1-mutant neighborhood of a given sequence (genotype) G0 with structure
(phenotype) P contains many structural variants that are not close to P. This
holds for biological and generic RNA sequences alike. Second, we analyze the
relation between new structures in the 1-neighborhoods of genotypes Gk that are
only a moderate Hamming distance k away from G0, and the structure of G0
itself, both for biological and for generic RNA structures. Third, we analyze
the relation between mutational robustness of a sequence and the distances of
structural variants near this sequence. Our findings underscore the role of
neutral networks in evolutionary innovation, and the role that high robustness
can play in diminishing the potential for such innovation.


  The scale free structure p(k)~k^{-gamma} of protein-protein interaction
networks can be reproduced by a static physical model in simulation. We inspect
the model theoretically, and find the key reason for the model to generate
apparent scale free degree distributions. This explanation provides a generic
mechanism of "scale free" networks. Moreover, we predict the dependence of
gamma on experimental protein concentrations or other sensitivity factors in
detecting interactions, and find experimental evidence to support the
prediction.


  The folding of the alpha-helice domain hbSBD of the mammalian mitochondrial
branched-chain alpha-ketoacid dehydrogenase (BCKD) complex is studied by the
circular dichroism technique in absence of urea. Thermal denaturation is used
to evaluate various thermodynamic parameters defining the equilibrium
unfolding, which is well described by the two-state model with the folding
temperature T_f = 317.8 K and the enthalpy change Delta H_g = 19.67 kcal/mol.
The folding is also studied numerically using the off-lattice coarse-grained Go
model and the Langevin dynamics. The obtained results, including the population
of the native basin, the free energy landscape as a function of the number of
native contacts and the folding kinetics, also suggest that the hbSBD domain is
a two-state folder. These results are consistent with the biological function
of hbSBD in BCKD.


  This paper (parts I and II) provides an expository introduction to monotone
and near-monotone dynamical systems associated to biochemical networks, those
whose graphs are consistent or near-consistent. Many conclusions can be drawn
from signed network structure, associated to purely stoichiometric information
and ignoring fluxes. In particular, monotone systems respond in a predictable
fashion to perturbations and have robust and ordered dynamical characteristics,
making them reliable components of larger networks. Interconnections of
monotone systems may be fruitfully analyzed using tools from control theory, by
viewing larger systems as interconnections of monotone subsystems. This allows
one to obtain precise bifurcation diagrams without appeal to explicit knowledge
of fluxes or of kinetic constants and other parameters, using merely
"input/output characteristics" (steady-state responses or DC gains). The
procedure may be viewed as a "model reduction" approach in which monotone
subsystems are viewed as essentially one-dimensional objects. The possibility
of performing a decomposition into a small number of monotone components is
closely tied to the question of how "near" a system is to being monotone. We
argue that systems that are "near monotone" may be more biologically more
desirable than systems that are far from being monotone. Indeed, there are
indications that biological networks may be much closer to being monotone than
random networks that have the same numbers of vertices and of positive and
negative edges.


  MicroRNAs (miRNAs) are endogenous 22-nucleotide RNAs, which suppress gene
expression by selectively binding to the 3-noncoding region of specific message
RNAs through base-pairing. Given the diversity and abundance of miRNA targets,
miRNAs appear to functionally interact with various components of many cellular
networks. By analyzing the interactions between miRNAs and a human cellular
signaling network, we found that miRNAs predominantly target positive
regulatory motifs, highly connected scaffolds and most downstream network
components such as signaling transcription factors, but less frequently target
negative regulatory motifs, common components of basic cellular machines and
most upstream network components such as ligands. In addition, when an adaptor
has potential to recruit more downstream components, these components are more
frequently targeted by miRNAs. This work uncovers the principles of miRNA
regulation of signal transduction networks and implies a potential function of
miRNAs for facilitating robust transitions of cellular response to
extracellular signals and maintaining cellular homeostasis.


  We have extended our previously developed 3D multi-scale agent-based brain
tumor model to simulate cancer heterogeneity and to analyze its impact across
the scales of interest. While our algorithm continues to employ an epidermal
growth factor receptor (EGFR) gene-protein interaction network to determine the
cells' phenotype, it now adds an explicit treatment of tumor cell adhesion
related to the model's biochemical microenvironment. We simulate a simplified
tumor progression pathway that leads to the emergence of five distinct glioma
cell clones with different EGFR density and cell 'search precisions'. The in
silico results show that microscopic tumor heterogeneity can impact the tumor
system's multicellular growth patterns. Our findings further confirm that EGFR
density results in the more aggressive clonal populations switching earlier
from proliferation-dominated to a more migratory phenotype. Moreover, analyzing
the dynamic molecular profile that triggers the phenotypic switch between
proliferation and migration, our in silico oncogenomics data display spatial
and temporal diversity in documenting the regional impact of tumorigenesis, and
thus support the added value of multi-site and repeated assessments in vitro
and in vivo. Potential implications from this in silico work for experimental
and computational studies are discussed.


  Relation of genome sizes to organisms complexity is still described rather
equivocally. Neither the number of genes (G-value), nor the total amount of DNA
(C-value) correlates consistently with phenotype complexity. Using information
theory considerations we developed a model that allows a quantative estimate
for the amount of functional information in a genomic sequence. This model
easily answers the long-standing question of why GC content is increased in
functional regions. The model allows consistent estimate of genome
complexities, resolving the major discrepancies of G- and C-values. For related
organisms with similarly complex phenotypes, this estimate provides biological
insights into their niches complexities. This theoretical framework suggests
that biological information can rapidly evolve on demand from environment,
mainly in non-coding genomic sequence and explains the role of duplications in
the evolution of biological information. Knowing the approximate amount of
functionality in a genomic sequence is useful for many applications such as
phylogenetics analyses, in-silico functional elements discovery or prioritising
targets for genotyping and sequencing.


  We analyze the dependence of cooperativity of the thermal denaturation
transition and folding rates of globular proteins on the number of amino acid
residues, $N$, using lattice models with side chains,off-lattice Go models and
the available experimental data. A dimensionless measure of cooperativity,
$\Omega_c$ ($0 < \Omega_c < \infty$), scales as $\Omega_c \sim N^{\zeta}$. The
results of simulations and the analysis of experimental data further confirm
the earlier prediction that $\zeta$ is universal with $\zeta = 1 +\gamma$,
where exponent $\gamma$ characterizes the susceptibility of a self-avoiding
walk. This finding suggests that the structural characteristics in the
denaturated state are manifested in the folding cooperativity at the transition
temperature. The folding rates $k_F$ for the Go models and a dataset of 69
proteins can be fit using $k_F = k_F^0 \exp(-cN^\beta)$. Both $\beta = 1/2$ and
2/3 provide a good fit of the data. We find that $k_F = k_F^0
\exp(-cN^{{1/2}})$, with the average (over the dataset of proteins) $k_F^0
\approx (0.2\mu s)^{-1}$ and $c \approx 1.1$, can be used to estimate folding
rates to within an order of magnitude in most cases. The minimal models give
identical $N$ dependence with $c \approx 1$. The prefactor for off-lattice Go
models is nearly four orders of magnitude larger than the experimental value.


  The didemnins represent a versatile class of depsipeptides of marine origin
and hold a great deal of potential for biomedical application. The biological
and geographical origins of the didemnins are reviewed in addition to the
chemical structures of the major didemnins. The biological mechanisms behind
the antiviral and anticancer effects of selected didemnins are summarized and
the special case of dehydrodidemnin B (Aplidin) is expounded upon including
structural characteristics, synthesis, pharmacological mechanism and a
discussion of its current clinical trials as an anticancer agent.


  We quantify the influence of the topology of a transcriptional regulatory
network on its ability to process environmental signals. By posing the problem
in terms of information theory, we may do this without specifying the function
performed by the network. Specifically, we study the maximum mutual information
between the input (chemical) signal and the output (genetic) response
attainable by the network in the context of an analytic model of particle
number fluctuations. We perform this analysis for all biochemical circuits,
including various feedback loops, that can be built out of 3 chemical species,
each under the control of one regulator. We find that a generic network,
constrained to low molecule numbers and reasonable response times, can
transduce more information than a simple binary switch and, in fact, manages to
achieve close to the optimal information transmission fidelity. These
high-information solutions are robust to tenfold changes in most of the
networks' biochemical parameters; moreover they are easier to achieve in
networks containing cycles with an odd number of negative regulators (overall
negative feedback) due to their decreased molecular noise (a result which we
derive analytically). Finally, we demonstrate that a single circuit can support
multiple high-information solutions. These findings suggest a potential
resolution of the "cross-talk" dilemma as well as the previously unexplained
observation that transcription factors which undergo proteolysis are more
likely to be auto-repressive.


  We are interested in the subtle variations of function among the members of a
protein family. A protein family is usually subdivided into subfamilies based
on functional differences. Existence of this functional diversity is essential
for the successful performance of physiological roles expected of the family.
This presents a unique problem: there must be preservation of the active site;
simultaneously there should be specificity of protein action according to
subfamily. Though the classification into subfamilies is by no means a
formalized one, it is most times based on the character of regulation of the
primary function. The function of a subfamily is a modification of when the
protein performs its function, for example, by changing the protein's
sensitivity to regulatory factors. Rarely, a subfamily possesses a function
completely different for its family. A study of these details is necessary for
understanding the fine-tuning of protein function. I describe a theory for
studying subfamily-based functional specificity and then validate it with an
example application to deciphering the residue-level basis of fine functional
variations in the diverse set of K+-channel subfamilies. I provide specific
results that will be useful to channel physiologists, whereas the strategy
developed will be widely applicable to problems in comparative and functional
genomics.


  Microarray-based genotyping is based on the high discrimination capability of
oligonucleotide probes. For detection of Single Nucleotide Polymorphisms (SNPs)
single-base discrimination is required. We investigate how various
point-mutations, comprising single base mismatches (MMs), insertions and
deletions, affect hybridization of DNA-DNA oligonucleotide duplexes. Employing
light-directed in situ synthesis we fabricate DNA microarrays with
comprehensive sets of cognate point-mutated probes, allowing us to
systematically investigate the influence of defect type, position and nearest
neighbor effects. Defect position has been identified as the dominating
influential factor. This positional effect which is almost identical for the
different point-mutation types, is biased from the local sequence environment.
The impact of the MM type is largely determined by the type of base pair
(either AT or CG) affected by the mismatch. We observe that single base
insertions next to like-bases result in considerably larger hybridization
signals than insertions next to nonidentical bases. The latter as well as the
distinct position dependence could be explained by a kinetic zipper model in
which point defects represent a barrier for the rapid closure of the DNA
duplex.


  Recent advances in high-throughput genomics technologies have resulted in the
sequencing of large numbers of (near) complete genomes. These genome sequences
are being mined for important functional elements, such as genes. They are also
being compared and contrasted in order to identify other functional sequences,
such as those involved in the regulation of genes. In cases where DNA sequences
from different organisms can be determined to have originated from a common
ancestor, it is natural to try to infer the an- cestral sequences. The
reconstruction of ancestral genomes can lead to insights about genome
evolution, and the origins and diversity of function. There are a number of
interesting foundational questions associated with reconstructing ancestral
genomes: Which statistical models for evolution should be used for making
inferences about ancestral sequences? How should extant genomes be compared in
order to facilitate ancestral reconstruction? Which portions of ancestral
genomes can be reconstructed reliably, and what are the limits of ancestral
reconstruction? We discuss recent progress on some of these questions, offer
some of our own opinions, and highlight interesting mathematics, statistics,
and computer science problems.


  The coexistence between different informational molecules has been the
preferred mode to circumvent the limitation posed by imperfect replication on
the amount of information stored by each of these molecules. Here we reexamine
a classic package model in which distinct information carriers or templates are
forced to coexist within vesicles, which in turn can proliferate freely through
binary division. The combined dynamics of vesicles and templates is described
by a multitype branching process which allows us to write equations for the
average number of the different types of vesicles as well as for their
extinction probabilities. The threshold phenomenon associated to the extinction
of the vesicle population is studied quantitatively using finite-size scaling
techniques. We conclude that the resultant coexistence is too frail in the
presence of parasites and so confinement of templates in vesicles without an
explicit mechanism of cooperation does not resolve the information crisis of
prebiotic evolution.


TThe paper had many errors.


  Over the past few years, microRNAs (miRNAs) have emerged as a new prominent
class of gene regulatory factors that negatively regulate expression of
approximately one-third of the genes in animal genomes at post-transcriptional
level. However, it is still unclear why some genes are regulated by miRNAs but
others are not, i.e. what principles govern miRNA regulation in animal genomes.
In this study, we systematically analyzed the relationship between
transcription factors (TFs) and miRNAs in gene regulation. We found that the
genes with more TF-binding sites have a higher probability of being targeted by
miRNAs and have more miRNA-binding sites on average. This observation reveals
that the genes with higher cis-regulation complexity are more coordinately
regulated by TFs at the transcriptional level and by miRNAs at the
post-transcriptional level. This is a potentially novel discovery of mechanism
for coordinated regulation of gene expression. Gene ontology analysis further
demonstrated that such coordinated regulation is more popular in the
developmental genes.


  Even under constant external conditions, the expression levels of genes
fluctuate. Much emphasis has been placed on the components of this noise that
are due to randomness in transcription and translation; here we analyze the
role of noise associated with the inputs to transcriptional regulation, the
random arrival and binding of transcription factors to their target sites along
the genome. This noise sets a fundamental physical limit to the reliability of
genetic control, and has clear signatures, but we show that these are easily
obscured by experimental limitations and even by conventional methods for
plotting the variance vs. mean expression level. We argue that simple, global
models of noise dominated by transcription and translation are inconsistent
with the embedding of gene expression in a network of regulatory interactions.
Analysis of recent experiments on transcriptional control in the early
Drosophila embryo shows that these results are quantitatively consistent with
the predicted signatures of input noise, and we discuss the experiments needed
to test the importance of input noise more generally.


  Thermodynamic constraints on reactions directions are inherent in the
structure of a given biochemical network. However, concrete procedures for
determining feasible reaction directions for large-scale metabolic networks are
not well established. This work introduces a systematic approach to compute
reaction directions, which are constrained by mass balance and thermodynamics,
for genome-scale networks. In addition, it is shown that the nonconvex solution
space constrained by physicochemical constraints can be approximated by a set
of linearized subspaces in which mass and thermodynamic balance are guaranteed.
The developed methodology can be used to {\it ab initio} predict reaction
directions of genome-scale networks based solely on the network stoichoimetry.


  Starting from the chemical master equation, we employ field theoretic
techniques to derive Langevin-type equations that exactly describe the
stochastic dynamics of the Frank chiral amplification model with spatial
diffusion. The intrinsic multiplicative noise properties are completely and
rigorously derived by this procedure. We carry out numerical simulations in two
spatial dimensions. When the inherent spatio-temporal fluctuations are properly
included, then complete chiral amplification results from a purely racemic
initial configuration. Phase separation can also arise in which the enantiomers
coexist in spatially segregated domains separated by a sharp racemic interface
or boundary.


  Here we introduce a model in which individuals differ in the rate at which
they seek new interactions with others, making rational decisions modeled as
general symmetric two-player games. Once a link between two individuals has
formed, the productivity of this link is evaluated. Links can be broken off at
different rates. We provide analytic results for the limiting cases where
linking dynamics is much faster than evolutionary dynamics and vice-versa, and
show how the individual capacity of forming new links or severing inconvenient
ones maps into the problem of strategy evolution in a well-mixed population
under a different game. For intermediate ranges, we investigate numerically the
detailed interplay determined by these two time-scales and show that the scope
of validity of the analytical results extends to a much wider ratio of time
scales than expected.


  Discrete temporal transitions occur in a variety of domains, but this work is
mainly motivated by applications in molecular biology: explaining and analyzing
observed transcriptome and proteome time series by literature and database
knowledge. The starting point of a formal concept analysis model is presented.
The objects of a formal context are states of the interesting entities, and the
attributes are the variable properties defining the current state (e.g.
observed presence or absence of proteins). Temporal transitions assign a
relation to the objects, defined by deterministic or non-deterministic
transition rules between sets of pre- and postconditions. This relation can be
generalized to its transitive closure, i.e. states are related if one results
from the other by a transition sequence of arbitrary length. The focus of the
work is the adaptation of the attribute exploration algorithm to such a
relational context, so that questions concerning temporal dependencies can be
asked during the exploration process and be answered from the computed stem
base. Results are given for the abstract example of a game and a small gene
regulatory network relevant to a biomedical question.


  We develop a simple, chemostat-based model illustrating how a process
analogous to associative learning can occur in a biochemical network.
Associative learning is a form of learning whereby a system "learns" to
associate two stimuli with one another. In our model, two types of replicating
molecules, denoted A and B, are present in some initial concentration in the
chemostat. Molecules A and B are stimulated to replicate by some growth
factors, denoted GA and GB, respectively. It is also assumed that A and B can
covalently link, and that the conjugated molecule can be stimulated by either
the GA or GB growth factors (and can be degraded). We show that, if the
chemostat is stimulated by both growth factors for a certain time, followed by
a time gap during which the chemostat is not stimulated at all, and if the
chemostat is then stimulated again by only one of the growth factors, then
there will be a transient increase in the number of molecules activated by the
other growth factor. Therefore, the chemostat bears the imprint of earlier,
simultaneous stimulation with both growth factors, which is indicative of
associative learning. It is interesting to note that the dynamics of our model
is consistent with various aspects of Pavlov's original series of associative
learning experiments in dogs. We discuss how associative learning can
potentially be performed in vitro within RNA, DNA, or peptide networks. We also
highlight how such a mechanism could potentially be involved in genomic
evolution, and suggest bioinformatics studies that could be used to find
evidence for associative learning processes at work inside living cells.


  The ability to learn and respond to recurrent events depends on the capacity
to remember transient biological signals received in the past. Moreover, it may
be desirable to remember or ignore these transient signals conditioned upon
other signals that are active at specific points in time or in unique
environments. Here, we propose a simple genetic circuit in bacteria that is
capable of conditionally memorizing a signal in the form of a transcription
factor concentration. The circuit behaves similarly to a "data latch" in an
electronic circuit, i.e. it reads and stores an input signal only when
conditioned to do so by a "read command". Our circuit is of the same size as
the well-known genetic toggle switch (an unconditional latch) which consists of
two mutually repressing genes, but is complemented with a "regulatory front
end" involving protein heterodimerization as a simple way to implement
conditional control. Deterministic and stochastic analysis of the circuit
dynamics indicate that an experimental implementation is feasible based on
well-characterized genes and proteins. It is not known, to which extent
molecular networks are able to conditionally store information in natural
contexts for bacteria. However, our results suggest that such sequential logic
elements may be readily implemented by cells through the combination of
existing protein-protein interactions and simple transcriptional regulation.


  In this paper, we investigate the emergence of a ratio-dependent
predator-prey system with Michaelis-Menten-type functional response and
reaction-diffusion. We derive the conditions for Hopf, Turing and Wave
bifurcation on a spatial domain. Furthermore, we present a theoretical analysis
of evolutionary processes that involves organisms distribution and their
interaction of spatially distributed population with local diffusion. The
results of numerical simulations reveal that the typical dynamics of population
density variation is the formation of isolated groups, i.e., stripelike or
spotted or coexistence of both. Our study shows that the spatially extended
model has not only more complex dynamic patterns in the space, but also chaos
and spiral waves. It may help us better understand the dynamics of an aquatic
community in a real marine environment.


  Escherichia coli bacteria respond to DNA damage by a highly orchestrated
series of events known as the SOS response, regulated by transcription factors,
protein-protein binding and active protein degradation. We present a dynamical
model of the UV-induced SOS response, incorporating mutagenesis by the
error-prone polymerase, Pol V. In our model, mutagenesis depends on a
combination of two key processes: damage counting by the replication forks and
a long term memory associated with the accumulation of UmuD'. Together, these
provide a tight regulation of mutagenesis resulting, we show, in a "digital"
turn-on and turn-off of Pol V. Our model provides a compact view of the
topology and design of the SOS network, pinpointing the specific functional
role of each of the regulatory processes. In particular, we suggest that the
recently observed second peak in the activity of promoters in the SOS regulon
(Friedman et al., 2005, PLoS Biol. 3, e238) is the result of a positive
feedback from Pol V to RecA filaments.


  The structure of DNA in the nucleosome core particle is studied using an
elastic model that incorporates anisotropy in the bending energetics and
twist-bend coupling. Using the experimentally determined structure of
nucleosomal DNA [T.J. Richmond and C.A. Davey, Nature {\bf 423}, 145 (2003)],
it is shown that elastic correlations exist between twist, roll, tilt, and
stretching of DNA, as well as the distance between phosphate groups. The
twist-bend coupling term is shown to be able to capture these correlations to a
large extent, and a fit to the experimental data yields a new estimate of G=25
nm for the value of the twist-bend coupling constant.


  Many viruses have the cunning ability to enter a hibernative or off state,
termed latency or lysogeny. When in a latent state, the virus is unable to
replicate, and its gene expression program is largely shut down. This facility
for lying dormant typically ensures lifelong persistence of the virus in the
host; it is also a particularly problematic obstacle in the treatment of HIV.
For most viruses, the molecular regulation of entry into latency is not
completely understood, but it is believed that viral gene expression must be
deactivated in some way. In this study, we introduce a new regulatory motif,
the feedback resistor, that enables a genetic circuit to shut off without the
need for an active repressor molecule. We first show that many animal viruses
might encode feedback resistors in their regulatory circuits. Then, by using a
combination of mathematical theory and single-cell real-time imaging
experiments, we show that a feedback resistor in the HIV Tat transcriptional
circuit likely allows the HIV virus to enter into latency. We postulate that
feedback resistors may give increased stability and control in the complex
noisy signaling environment of the cell.


  This paper argues that self-awareness is a learned behavior that emerges in
organisms whose brains have a sufficiently integrated, complex ability for
associative learning and memory. Continual sensory input of information related
to the organism causes the organism's brain to learn the physical
characteristics of the organism, in the sense that neural pathways are produced
that are reinforced by, and therefore recognize, various features associated
with the organism. This results in the formation of a set of associations that
may be termed an organismal self-image, which provides a mechanistic basis for
the emergence of various behaviors that are associated with self-awareness,
such as self-recognition. In humans, self-awareness includes additional
behaviors such as recognition of self-awareness, the concept of I, and various
existential and religious questions. This paper shows how associative memory
and learning, combined with an organismal self-image and, in the case of
humans, language, leads to the emergence of these various behaviors. This paper
also discusses various tautologies that invariably emerge when discussing
self-awareness, that ultimately prevent an unambiguous resolution to the
various existential issues that arise. We continue with various speculations on
manipulating self-awareness, and discuss how concepts from set and logic may
provide a highly useful set of tools in computational neuroscience for
understanding the emergence of higher cognitive functions in complex organisms.
The existence of other types of awareness, and the role of mirror neurons in
the emergence of self-awareness, are also briefly discussed.


  Oscillatory behavior is a key property of many biological systems. The
Small-Gain Theorem (SGT) for input/output monotone systems provides a
sufficient condition for global asymptotic stability of an equilibrium and
hence its violation is a necessary condition for the existence of periodic
solutions. One advantage of the use of the monotone SGT technique is its
robustness with respect to all perturbations that preserve monotonicity and
stability properties of a very low-dimensional (in many interesting examples,
just one-dimensional) model reduction. This robustness makes the technique
useful in the analysis of molecular biological models in which there is large
uncertainty regarding the values of kinetic and other parameters. However,
verifying the conditions needed in order to apply the SGT is not always easy.
This paper provides an approach to the verification of the needed properties,
and illustrates the approach through an application to a classical model of
circadian oscillations, as a nontrivial ``case study,'' and also provides a
theorem in the converse direction of predicting oscillations when the SGT
conditions fail.


  We propose the concepts of distributed robustness and r-robustness, well
adapted to functional genetics. Then we discuss the robustness of the
relaxation time using a chemical reaction description of genetic and signalling
networks. First, we obtain the following result for linear networks: for large
multiscale systems with hierarchical distribution of time scales the variance
of the inverse relaxation time (as well as the variance of the stationary rate)
is much lower than the variance of the separate constants. Moreover, it can
tend to 0 faster than 1/n, where n is the number of reactions. We argue that
similar phenomena are valid in the nonlinear case as well. As a numerical
illustration we use a model of signalling network that can be applied to
important transcription factors such as NFkB.


  We present MEDUSA, an integrative method for learning motif models of
transcription factor binding sites by incorporating promoter sequence and gene
expression data. We use a modern large-margin machine learning approach, based
on boosting, to enable feature selection from the high-dimensional search space
of candidate binding sequences while avoiding overfitting. At each iteration of
the algorithm, MEDUSA builds a motif model whose presence in the promoter
region of a gene, coupled with activity of a regulator in an experiment, is
predictive of differential expression. In this way, we learn motifs that are
functional and predictive of regulatory response rather than motifs that are
simply overrepresented in promoter sequences. Moreover, MEDUSA produces a model
of the transcriptional control logic that can predict the expression of any
gene in the organism, given the sequence of the promoter region of the target
gene and the expression state of a set of known or putative transcription
factors and signaling molecules. Each motif model is either a $k$-length
sequence, a dimer, or a PSSM that is built by agglomerative probabilistic
clustering of sequences with similar boosting loss. By applying MEDUSA to a set
of environmental stress response expression data in yeast, we learn motifs
whose ability to predict differential expression of target genes outperforms
motifs from the TRANSFAC dataset and from a previously published candidate set
of PSSMs. We also show that MEDUSA retrieves many experimentally confirmed
binding sites associated with environmental stress response from the
literature.


  Evolution of life on earth was governed, primarily, by natural selection,
with major contribution of other evolutionary processes, such as neutral
variation, exaptation, and gene duplication. However, for biological evolution
to take off, a certain minimal degree of complexity is required such that a
replicating genome encodes means for its own replication with sufficient rate
and fidelity. In all existing life forms, this is achieved by dedicated
proteins, polymerases (replicases), that are produced by the elaborate
translation system. However, evolution of the coupled system of replication and
translation does not appear possible without pre-existing efficient
replication; hence a chicken-egg type paradox. I argue that the
many-worlds-in-one version of the cosmological model of eternal inflation
implies that emergence of replication and translation, as well as the major
protein folds, by chance alone, as opposed to biological evolution, is a
realistic possibility and could provide for the onset of biological evolution.
Hence an RNA world, as it is currently conceived, might have never existed
although catalytic activities of RNA were, probably, critical for the onset of
biological evolution and its early stages.


  We study the Susceptible-Infected-Recovered model of epidemics in the
vicinity of the threshold infectivity. We derive the distribution of total
outbreak size in the limit of large population size $N$. This is accomplished
by mapping the problem to the first passage time of a random walker subject to
a drift that increases linearly with time. We recover the scaling results of
Ben-Naim and Krapivsky that the effective maximal size of the outbreak scales
as $N^{2/3}$, with the average scaling as $N^{1/3}$, with an explicit form for
the scaling function.


  Boolean networks with canalizing functions are used to model gene regulatory
networks. In order to learn how such networks may behave under evolutionary
forces, we simulate the evolution of a single Boolean network by means of an
adaptive walk, which allows us to explore the fitness landscape. Mutations
change the connections and the functions of the nodes. Our fitness criterion is
the robustness of the dynamical attractors against small perturbations. We find
that with this fitness criterion the global maximum is always reached and that
there is a huge neutral space of 100% fitness. Furthermore, in spite of having
such a high degree of robustness, the evolved networks still share many
features with "chaotic" networks.


  It has been demonstrated that the observed ratio of phenotypes of marker
enzymes in some sugarbeet plants produced by mitotic agamospermy can be
explained by different degrees of endoreduplication of chromosomes carrying
different alleles of the enzyme loci. In these plants, different patterns of
variability of the enzymes controlled by the linked loci suggest different
degrees of endoreduplication of different chromosomal regions. A concept of
multidimensional encoding of inherited information in eukaryotes has been
proposed.


  We develop a theoretical approach to the protein folding problem based on
out-of-equilibrium stochastic dynamics. Within this framework, the
computational difficulties related to the existence of large time scale gaps in
the protein folding problem are removed and simulating the entire reaction in
atomistic details using existing computers becomes feasible. In addition, this
formalism provides a natural framework to investigate the relationships between
thermodynamical and kinetic aspects of the folding. For example, it is possible
to show that, in order to have a large probability to remain unchanged under
Langevin diffusion, the native state has to be characterized by a small
conformational entropy. We discuss how to determine the most probable folding
pathway, to identify configurations representative of the transition state and
to compute the most probable transition time. We perform an illustrative
application of these ideas, studying the conformational evolution of alanine
di-peptide, within an all-atom model based on the empiric GROMOS96 force field.


  It has been observed that mutualistic bipartite networks have a nested
structure of interactions. In addition, the degree distributions associated
with the two guilds involved in such networks (e.g. plants & pollinators or
plants & seed dispersers) approximately follow a truncated power law. We show
that nestedness and truncated power law distributions are intimately linked,
and that any biological reasons for such truncation are superimposed to finite
size effects . We further explore the internal organization of bipartite
networks by developing a self-organizing network model (SNM) that reproduces
empirical observations of pollination systems of widely different sizes. Since
the only inputs to the SNM are numbers of plant and animal species, and their
interactions (i.e., no data on local abundance of the interacting species are
needed), we suggest that the well-known association between species frequency
of interaction and species degree is a consequence rather than a cause, of the
observed network structure.


  Molecular chaperones play a prominent role in signaling and transcriptional
regulatory networks of the cell. Recent advances uncovered that chaperones act
as genetic buffers stabilizing the phenotype of various cells and organisms and
may serve as potential regulators of evolvability. Chaperones have weak links,
connect hubs, are in the overlaps of network modules and may uncouple these
modules during stress, which gives an additional protection for the cell at the
network-level. Moreover, after stress chaperones are essential to re-build
inter-modular contacts by their low affinity sampling of the potential
interaction partners in different modules. This opens the way to the
chaperone-regulated modular evolution of cellular networks, and helps us to
design novel therapeutic and anti-aging strategies.


  A deterministic model of an age-structured population with genetics analogous
to the discrete time Penna model of genetic evolution is constructed on the
basis of the Lotka-Volterra scheme. It is shown that if, as in the Penna model,
genetic information is represented by the fraction of defective genes in the
population, the population numbers for each specific individual's age are
represented by exactly the same functions of age in both models. This gives us
a new possibility to consider multi-species evolution without using detailed
microscopic Penna model.
  We discuss a particular case of the predator-prey system representing an
ecosystem consisting of a limited amount of energy resources consumed by the
age-structured species living in this ecosystem. Then, the increase in number
of the individuals in the population under consideration depends on the
available energy resources, the shape of the distribution function of defective
genes in the population and the fertility age. We show that these parameters
determine the trend toward equilibrium of the whole ecosystem.


  The apparent stability of population oscillations in ecological systems is a
long-standing puzzle. A generic solution for this problem is suggested here.
The stabilizing mechanism involves the combined effect of spatial migration,
amplitude-dependent frequency, and noise: small differences between spatial
patches, induced by the noise, are amplified by the frequency gradient.
Migration among desynchronized regions then stabilizes the oscillations in the
vicinity of the homogenous manifold. A simple model of diffusively coupled
oscillators allows the derivation of quantitative results, like the dependence
of the desynchronization upon diffusion strength and frequency differences. For
an unstable system, a noise-induced transition is demonstrated, from extinction
for small noise to stability if the noise exceeds some threshold. The coupled
oscillator model is shown to reproduce all previously suggested stabilizing
mechanisms. Accordingly, we suggest using this model as a standard tool for
identification and classification of population oscillations on spatial
domains.


  The E.coli transcription network has an essentially feedforward structure,
with, however, abundant feedback at the level of self-regulations. Here, we
investigate how these properties emerged during evolution. An assessment of the
role of gene duplication based on protein domain architecture shows that (i)
transcriptional autoregulators have mostly arisen through duplication, while
(ii) the expected feedback loops stemming from their initial cross-regulation
are strongly selected against. This requires a divergent coevolution of the
transcription factor DNA-binding sites and their respective DNA cis-regulatory
regions. Moreover, we find that the network tends to grow by expansion of the
existing hierarchical layers of computation, rather than by addition of new
layers. We also argue that rewiring of regulatory links due to
mutation/selection of novel transcription factor/DNA binding interactions
appears not to significantly affect the network global hierarchy, and that
horizontally transferred genes are mainly added at the bottom, as new target
nodes. These findings highlight the important evolutionary roles of both
duplication and selective deletion of crosstalks between autoregulators in the
emergence of the hierarchical transcription network of E.coli.


  We propose a two-level stochastic context-free grammar (SCFG) architecture
for parametrized stochastic modeling of a family of RNA sequences, including
their secondary structure. A stochastic model of this type can be used for
maximum a posteriori estimation of the secondary structure of any new sequence
in the family. The proposed SCFG architecture models RNA subsequences
comprising paired bases as stochastically weighted Dyck-language words, i.e.,
as weighted balanced-parenthesis expressions. The length of each run of
unpaired bases, forming a loop or a bulge, is taken to have a phase-type
distribution: that of the hitting time in a finite-state Markov chain. Without
loss of generality, each such Markov chain can be taken to have a bounded
complexity. The scheme yields an overall family SCFG with a manageable number
of parameters.


  Pathological folding and oligomer formation of the amyloid beta-protein
(Abeta) are widely perceived as central to Alzheimer's disease (AD).
Experimental approaches to study Abeta self-assembly are problematic, because
most relevant aggregates are quasi-stable and inhomogeneous. We apply a
discrete molecular dynamics (DMD) approach combined with a four-bead protein
model to study oligomer formation of the amyloid beta-protein (Abeta). We
address the differences between the two most common Abeta alloforms, Abeta40
and Abeta42, which oligomerize differently in vitro. We study how the presence
of electrostatic interactions (EIs) between pairs of charged amino acids
affects Abeta40 and Abeta42 oligomer formation. Our results indicate that EIs
promote formation of larger oligomers in both Abeta40 and Abeta42. The Abeta40
size distribution remains unimodal, whereas the Abeta42 distribution is
trimodal, as observed experimentally. Abeta42 folded structure is characterized
by a turn in the C-terminus that is not present in Abeta40. We show that the
same C-terminal region is also responsible for the strongest intermolecular
contacts in Abeta42 pentamers and larger oligomers. Our results suggest that
this C-terminal region plays a key role in the formation of Abeta42 oligomers
and the relative importance of this region increases in the presence of EIs.
These results suggest that inhibitors targeting the C-terminal region of
Abeta42 oligomers may be able to prevent oligomer formation or structurally
modify the assemblies to reduce their toxicity.


  Quantitative computational models play an increasingly important role in
modern biology. Such models typically involve many free parameters, and
assigning their values is often a substantial obstacle to model development.
Directly measuring \emph{in vivo} biochemical parameters is difficult, and
collectively fitting them to other data often yields large parameter
uncertainties. Nevertheless, in earlier work we showed in a
growth-factor-signaling model that collective fitting could yield
well-constrained predictions, even when it left individual parameters very
poorly constrained. We also showed that the model had a `sloppy' spectrum of
parameter sensitivities, with eigenvalues roughly evenly distributed over many
decades. Here we use a collection of models from the literature to test whether
such sloppy spectra are common in systems biology. Strikingly, we find that
every model we examine has a sloppy spectrum of sensitivities. We also test
several consequences of this sloppiness for building predictive models. In
particular, sloppiness suggests that collective fits to even large amounts of
ideal time-series data will often leave many parameters poorly constrained.
Tests over our model collection are consistent with this suggestion. This
difficulty with collective fits may seem to argue for direct parameter
measurements, but sloppiness also implies that such measurements must be
formidably precise and complete to usefully constrain many model predictions.
We confirm this implication in our signaling model. Our results suggest that
sloppy sensitivity spectra are universal in systems biology models. The
prevalence of sloppiness highlights the power of collective fits and suggests
that modelers should focus on predictions rather than on parameters.


  Despite decades of conservation efforts on the nesting beaches, the critical
status of leatherback turtles shows that their survival predominantly depends
on our ability to reduce at-sea mortality. Although areas where leatherbacks
meet fisheries have been identified during the long distance movements between
two consecutive nesting seasons, hotspots of lethal interactions are still
poorly defined within the nesting season, when individuals concentrate close to
land. Here we report movements of satellite-tracked gravid leatherback turtles
during the nesting season in Western Central Africa, South America and
Caribbean Sea, accounting for about 70% of the world population. We show that
during, and at the end of, the nesting season leatherback turtles have the
propensity to remain over the continental shelf, yet sometimes perform extended
movements and may even nest in neighbouring countries. Leatherbacks exploit
coastal commercial fishing grounds and face substantial accidental capture by
regional coastal fisheries (e.g. at least 10% in French Guiana). This
emphasises the need for regional conservation strategies to be developed at the
ocean scale, both at sea and on land, to ensure the survival of the last
leatherback turtles.


  Calcium (Ca) plays an important role in regulating various cellular
processes. In a variety of cell types, Ca signaling occurs within microdomains
where channels deliver localized pulses of Ca which activate a nearby
collection of Ca-sensitive receptors. The small number of channels involved
ensures that the signaling process is stochastic. The aggregate response of
several thousand of these microdomains yields a whole-cell response which
dictates the cell behavior. Here, we study analytically the statistical
properties of a population of these microdomains in response to a trigger
signal. We apply these results to understand the relationship between Ca influx
and Ca release in cardiac cells. In this context, we use a first passage time
approach to show analytically how Ca release in the whole cell depends on the
single channel kinetics of Ca channels and the properties of microdomains.
Using these results, we explain the underlying mechanism for the graded
relationship between Ca influx and Ca release in cardiac cells.


  We analyze general two-species stochastic models, of the kind generally used
for the study of population dynamics. We show that the conditions for the
stochastic (microscopic) model to display approximate sustained oscillatory
behavior are governed by the parameters of the corresponding deterministic
(macroscopic) model. We provide a quantitative criterion for the quality of the
stochastic oscillation, using a dimensionless parameter that depends only on
the deterministic model. When this parameter is small, the oscillations are
clear, and the frequencies of the stochastic and deterministic oscillations are
close, for all stochastic models compatible with the same deterministic one. On
the other hand, when it is large, the oscillations cannot be distinguished from
a noise.


  In the bacterium Escherichia coli, selection of the division site involves
pole-to-pole oscillations of the proteins MinD and MinE. Different oscillation
mechanisms based on cooperative effects between Min-proteins and on the
exchange of Min-proteins between the cytoplasm and the cytoplasmic membrane
have been proposed. The parameters characterizing the dynamics of the
Min-proteins in vivo are not known. It has therefore been difficult to compare
the models quantitatively with experiments. Here, we present in vivo
measurements of the mobility of MinD and MinE using fluorescence correlation
spectroscopy. Two distinct time-scales are clearly visible in the correlation
curves. While the faster time-scale can be attributed to cytoplasmic diffusion,
the slower time-scale could result from diffusion of membrane-bound proteins or
from protein exchange between the cytoplasm and the membrane. We determine the
diffusion constant of cytoplasmic MinD to be approximately 16\mu^{2}/s, while
for MinE we find about 10\mu^{2}/s, independently of the processes responsible
for the slower time-scale. Implications of the measured values for the
oscillation mechanism are discussed.


  Tracks of marine animals in the wild, now increasingly acquired by electronic
tagging of individuals, are of prime interest not only to identify habitats and
high-risk areas, but also to gain detailed information about the behaviour of
these animals. Using recent satellite-derived current estimates and leatherback
turtle (Dermochelys coriacea) tracking data, we demonstrate that oceanic
currents, usually neglected when analysing tracking data, can substantially
distort the observed trajectories. Consequently, this will affect several
important results deduced from the analysis of tracking data, such as the
evaluation of the orientation skills and the energy budget of animals or the
identification of foraging areas. We conclude that currents should be
systematically taken into account to ensure the unbiased interpretation of
tracking data, which now play a major role in marine conservation biology.


  The Laessig-Wiese (LW) field theory for the freezing transition of random RNA
secondary structures is generalized to the situation of an external force. We
find a second-order phase transition at a critical applied force f = f_c. For f
< f_c forces are irrelevant. For f > f_c, the extension L as a function of
pulling force f scales as (f-f_c)^(1/gamma-1). The exponent gamma is calculated
in an epsilon-expansion: At 1-loop order gamma = epsilon/2 = 1/2, equivalent to
the disorder-free case. 2-loop results yielding gamma = 0.6 are briefly
mentioned. Using a locking argument, we speculate that this result extends to
the strong-disorder phase.


  The Darwinian paradigm of biological evolution is based on the separability
of the variation and selection processes. As a result, the population thinking
had always been an integral part of the Darwinian approach. I propose an
alternative scheme of biological adaptation. It is based on appreciation of
limits of what we can observe considering an individual biological object. This
leads to a possibility for the adaptation process to occur on the level of a
single object, as a 'selection among virtual states of the organism'. I discuss
the application of this idea to the phenomenon of adaptive mutations in
bacteria.


  Although huddling was shown to be the key by which emperor penguins
(Aptenodytes forsteri) save energy and sustain their breeding fast during the
Antarctic winter, the intricacies of this social behavior have been poorly
studied. We recorded abiotic variables with data loggers glued to the feathers
of eight individually marked emperor penguins to investigate their
thermoregulatory behavior and to estimate their "huddling time budget"
throughout the breeding season (pairing and incubation period). Contrary to the
classic view, huddling episodes were discontinuous and of short and variable
duration, lasting 1.6+/-1.7 (S.D.) h on average. Despite heterogeneous huddling
groups, birds had equal access to the warmth of the huddles. Throughout the
breeding season, males huddled for 38+/-18% (S.D.) of their time, which raised
the ambient temperature that birds were exposed to above 0 degrees C (at
average external temperatures of -17 degrees C). As a consequence of tight
huddles, ambient temperatures were above 20 degrees C during 13+/-12% (S.D.) of
their huddling time. Ambient temperatures increased up to 37.5 degrees C, close
to birds' body temperature. This complex social behavior therefore enables all
breeders to get a regular and equal access to an environment which allows them
to save energy and successfully incubate their eggs during the Antarctic
winter.


  We study the stable attractors of a class of continuous dynamical systems
that may be idealized as networks of Boolean elements, with the goal of
determining which Boolean attractors, if any, are good approximations of the
attractors of generic continuous systems. We investigate the dynamics in simple
rings and rings with one additional self-input. An analysis of switching
characteristics and pulse propagation explains the relation between attractors
of the continuous systems and their Boolean approximations. For simple rings,
"reliable" Boolean attractors correspond to stable continuous attractors. For
networks with more complex logic, the qualitative features of continuous
attractors are influenced by inherently non-Boolean characteristics of
switching events.


  In this paper we develop a theory of general selection systems with discrete
time and explore the evolution of selection systems, in particular,
inhomogeneous populations. We show that the knowledge of the initial
distribution of the selection system allows us to determine explicitly the
system distribution at the entire time interval. All statistical
characteristics of interest, such as mean values of the fitness or any trait
can be predicted effectively for indefinite time and these predictions
dramatically depend on the initial distribution. The Fisher Fundamental theorem
of natural selection (FTNS) and more general the Price equations are the famous
results of the mathematical selection theory. We show that the problem of
dynamic insufficiency for the Price equations and for the FTNS can be resolved
within the framework of selection systems. Effective formulas for solutions of
the Price equations and for the FTNS are derived. Applications of the developed
theory to some other problems of mathematical biology (dynamics of
inhomogeneous logistic and Ricker model, selection in rotifer populations) are
also given. Complex behavior of the total population size, the mean fitness (in
contrast to the plain FTNS) and other traits is possible for inhomogeneous
populations with density-dependent fitness. The temporary dynamics of these
quantities can be investigated with the help of suggested methods.


  We present a method that compares the protein interaction networks of two
species to detect functionally similar (conserved) protein modules between
them. The method is based on an algorithm we developed to identify matching
subgraphs between two graphs. Unlike previous network comparison methods, our
algorithm has provable guarantees on correctness and efficiency. Our algorithm
framework also admits quite general connectivity and local matching criteria
that define when two subgraphs match and constitute a conserved module.
  We apply our method to pairwise comparisons of the yeast protein network with
the human, fruit fly and nematode worm protein networks, using a lenient
criterion based on connectedness and matching edges, coupled with a betweenness
clustering heuristic. We evaluate the detected conserved modules against
reference yeast protein complexes using sensitivity and specificity measures.
In these evaluations, our method performs competitively with and sometimes
better than two previous network comparison methods. Further under some
conditions (proper homolog and species selection), our method performs better
than a popular single-species clustering method. Beyond these evaluations, we
discuss the biology of a couple of conserved modules detected by our method. We
demonstrate the utility of network comparison for transferring annotations from
yeast proteins to human ones, and validate the predicted annotations.


  We consider two versions of stochastic population models with mutation and
selection. The first approach relies on a multitype branching process; here,
individuals reproduce and change type (i.e., mutate) independently of each
other, without restriction on population size. We analyze the equilibrium
behaviour of this model, both in the forward and in the backward direction of
time; the backward point of view emerges if the ancestry of individuals chosen
randomly from the present population is traced back into the past.
  The second approach is the Moran model with selection. Here, the population
has constant size N. Individuals reproduce (at rates depending on their types),
the offspring inherits the parent's type, and replaces a randomly chosen
individual (to keep population size constant). Independently of the
reproduction process, individuals can change type. As in the branching model,
we consider the ancestral lines of single individuals chosen from the
equilibrium population. We use analytical results of Fearnhead (2002) to
determine the explicit properties, and parameter dependence, of the ancestral
distribution of types, and its relationship with the stationary distribution in
forward time.


  Genetic switch systems with mutual repression of two transcription factors
are studied using deterministic methods (rate equations) and stochastic methods
(the master equation and Monte Carlo simulations). These systems exhibit
bistability, namely two stable states such that spontaneous transitions between
them are rare. Induced transitions may take place as a result of an external
stimulus. We study several variants of the genetic switch and examine the
effects of cooperative binding, exclusive binding, protein-protein interactions
and degradation of bound repressors. We identify the range of parameters in
which bistability takes place, enabling the system to function as a switch.
Numerous studies have concluded that cooperative binding is a necessary
condition for the emergence of bistability in these systems. We show that a
suitable combination of network structure and stochastic effects gives rise to
bistability even without cooperative binding. The average time between
spontaneous transitions is evaluated as a function of the biological
parameters.


  Many fundamental questions concerning the emergence and subsequent evolution
of eukaryotic exon-intron organization are still unsettled. Genome-scale
comparative studies, which can shed light on crucial aspects of eukaryotic
evolution, require adequate computational tools.
  We describe novel computational methods for studying spliceosomal intron
evolution. Our goal is to give a reliable characterization of the dynamics of
intron evolution. Our algorithmic innovations address the identification of
orthologous introns, and the likelihood-based analysis of intron data. We
discuss a compression method for the evaluation of the likelihood function,
which is noteworthy for phylogenetic likelihood problems in general. We prove
that after $O(nL)$ preprocessing time, subsequent evaluations take $O(nL/\log
L)$ time almost surely in the Yule-Harding random model of $n$-taxon
phylogenies, where $L$ is the input sequence length.
  We illustrate the practicality of our methods by compiling and analyzing a
data set involving 18 eukaryotes, more than in any other study to date. The
study yields the surprising result that ancestral eukaryotes were fairly
intron-rich. For example, the bilaterian ancestor is estimated to have had more
than 90% as many introns as vertebrates do now.


  The complexity of the cells can be described and understood by a number of
networks such as protein-protein interaction, cytoskeletal, organelle,
signalling, gene transcription and metabolic networks. All these networks are
highly dynamic producing continuous rearrangements in their links, hubs,
network-skeleton and modules. Here we describe the adaptation of cellular
networks after various forms of stress causing perturbations, congestions and
network damage. Chronic stress decreases link-density, decouples or even
quarantines modules, and induces an increased competition between network hubs
and bridges. Extremely long or strong stress may induce a topological phase
transition in the respective cellular networks, which switches the cell to a
completely different mode of cellular function. We summarize our initial
knowledge on network restoration after stress including the role of molecular
chaperones in this process. Finally, we discuss the implications of
stress-induced network rearrangements in diseases and ageing, and propose
therapeutic approaches both to increase the robustness and help the repair of
cellular networks.


  We analytically address disease outbreaks in large, random networks with
heterogeneous infectivity and susceptibility. The transmissibility $T_{uv}$
(the probability that infection of $u$ causes infection of $v$) depends on the
infectivity of $u$ and the susceptibility of $v$. Initially a single node is
infected, following which a large-scale epidemic may or may not occur. We use a
generating function approach to study how heterogeneity affects the probability
that an epidemic occurs and, if one occurs, its attack rate (the fraction
infected). For fixed average transmissibility, we find upper and lower bounds
on these. An epidemic is most likely if infectivity is homogeneous and least
likely if the variance of infectivity is maximized. Similarly, the attack rate
is largest if susceptibility is homogeneous and smallest if the variance is
maximized. We further show that heterogeneity in infectious period is
important, contrary to assumptions of previous studies. We confirm our
theoretical predictions by simulation. Our results have implications for
control strategy design and identification of populations at higher risk from
an epidemic.


  Intestinal hexose absorption and gluconeogenesis have been studied in
relation to refeeding after two different fasting phases: a long period of
protein sparing during which energy expenditure is derived from lipid oxidation
(phase II), and a later phase characterized by a rise in plasma corticosterone
triggering protein catabolism (phase III). Such a switch in body fuel uses,
leading to changes in body reserves and gluconeogenic precursors, could
modulate intestinal gluconeogenesis and glucose transport. The gene and protein
levels, and the cellular localization of the sodium-glucose cotransporter
SGLT1, and of GLUT5 and GLUT2, as well as that of the key gluconeogenic enzymes
phosphoenolpyruvate carboxykinase (PEPCK) and glucose-6-phosphatase (Glc6Pase)
were measured. PEPCK and Glc6Pase activities were also determined. In phase III
fasted rats, SGLT1 was up-regulated and intestinal glucose uptake rates were
higher than in phase II fasted and fed rats. PEPCK and Glc6Pase mRNA, protein
levels and activities also increased in phase III. GLUT5 and GLUT2 were
down-regulated throughout the fast, but increased after refeeding, with GLUT2
recruited to the apical membrane. The increase in SGLT1 expression during phase
III may allow glucose absorption at low concentrations as soon as food is
available. Furthermore, an increased epithelial permeability due to fasting may
induce a paracellular movement of glucose. In the absence of intestinal GLUT2
during fasting, Glc6Pase could be involved in glucose release to the
bloodstream via membrane trafficking. Finally, refeeding triggered GLUT2 and
GLUT5 synthesis and apical recruitment of GLUT2, to absorb larger amounts of
hexoses.


  Biological systems are typically very complex and need to be reduced before
they are amenable to a thorough analysis. Also, they often possess functionally
important dynamic features like bistability. In model reduction, it is
sometimes more desirable to preserve the dynamic features only than to recover
a good quantitative approximation. We present an approach to reduce the order
of a bistable dynamical system significantly while preserving bistability and
the switching threshold. These properties are important for the operation of
the system in the context of a larger network. As an application example, a
bistable model for caspase activation in apoptosis is considered.


  In the presence of ATP, kinesin proceeds along the protofilament of
microtubule by alternated binding of two motor domains on the tubulin binding
sites. Since the processivity of kinesin is much higher than other motor
proteins, it has been speculated that there exists a mechanism for allosteric
regulation between the two monomers. Recent experiments suggest that ATP
binding to the leading head domain in kinesin is regulated by the rearward
strain built on the neck-linker. We test this hypothesis by explicitly modeling
a $C_{\alpha}$-based kinesin structure whose both motor domains are bound on
the tubulin binding sites. The equilibrium structures of kinesin on the
microtubule show disordered and ordered neck-linker configurations for the
leading and the trailing head, respectively. The comparison of the structures
between the two heads shows that several native contacts present at the
nucleotide binding site in the leading head are less intact than those in the
binding site of the rear head. The network of native contacts obtained from
this comparison provides the internal tension propagation pathway, which leads
to the disruption of the nucleotide binding site in the leading head. Also,
using an argument based on polymer theory, we estimate the internal tension
built on the neck-linker to be f~(12-15) pN. Both of these conclusions support
the experimental hypothesis.


  Simple stochastic models for phylogenetic trees on species have been well
studied. But much paleontology data concerns time series or trees on
higher-order taxa, and any broad picture of relationships between extant groups
requires use of higher-order taxa. A coherent model for trees on (say) genera
should involve both a species-level model and a model for the classification
scheme by which species are assigned to genera. We present a general framework
for such models, and describe three alternate classification schemes. Combining
with the species-level model of Aldous-Popovic (2005), one gets models for
higher-order trees, and we initiate analytic study of such models. In
particular we derive formulas for the lifetime of genera, for the distribution
of number of species per genus, and for the offspring structure of the tree on
genera.


  The interpretation of recent environmental genomics data exposes the
far-reaching influence of horizontal gene transfer, and is changing our basic
concepts of organism, species and evolution itself.


  We present an approach for constructing dynamic models for the simulation of
gene regulatory networks from simple computational elements. Each element is
called a ``gene gate'' and defines an input/output-relationship corresponding
to the binding and production of transcription factors. The proposed reaction
kinetics of the gene gates can be mapped onto stochastic processes and the
standard ode-description. While the ode-approach requires fixing the system's
topology before its correct implementation, expressing them in stochastic
pi-calculus leads to a fully compositional scheme: network elements become
autonomous and only the input/output relationships fix their wiring. The
modularity of our approach allows to pass easily from a basic first-level
description to refined models which capture more details of the biological
system. As an illustrative application we present the stochastic repressilator,
an artificial cellular clock, which oscillates readily without any cooperative
effects.


  Immunity is hypothesized to share limited resources with other physiological
functions and this may partly account for the fitness costs of reproduction.
Previous studies have shown that the acquired immunity of female common eider
ducks (Somateria mollissima) is suppressed during their incubation, during
which they entirely fast. Corticosterone was proposed to be an underlying
physiological mechanism for such immunosuppression. Therefore, the current
study aimed to assess the effects of exogenous corticosterone on acquired
immunity in captive eiders. To this end, females were implanted with
corticosterone pellets at different stages of their incubation fast. We
measured total immunoglobulin levels, T-cell-mediated immune response, body
mass and corticosterone levels in these females and compared them with those of
control females prior to and after manipulation (i.e. corticosterone pellet
implantation). To mimic corticosterone effects on body mass, we experimentally
extended fasting duration in a group of females termed ;late fasters'...


  Intestinal apoptosis and expression of apoptosis inducers - the cytokines
TNFalpha, TGFbeta1 - and the intestinal transcription factor Cdx2, were studied
according to two different metabolic and hormonal phases which characterize
long-term fasting: the long period of protein sparing during which energy
expenditure is derived from lipid oxidation (phase II), and the later phase
characterized by a rise in body protein utilization and plasma corticosterone
(phase III). Apoptosis was further studied in 2, 6, and 24 h refed rats.
Morphological apoptotic events were observed by environmental and conventional
scanning electron microscopy and a TUNEL test was used to characterize the
final stages of apoptotic death. The gene and protein expressions of TNFalpha,
TGFbeta1, and Cdx2 were measured. Apoptotic events and TNFalpha, TGFbeta1, and
Cdx2 gene and protein expressions did not vary significantly during phase II as
compared to the normally fed animals. However, a phase III fasting induced a
delay in intestinal epithelial apoptosis, along with a 92, 58, and 25% decrease
in TNFalpha, TGFbeta1, and Cdx2 mRNAs, respectively. The amounts of TNFalpha,
TGFbeta1, and Cdx2 proteins decreased by 70, 36, and 25%, respectively.
Apoptosis was restored rapidly after a 2 h refeeding following the phase III,
accompanied by a significant increase in TNFalpha, TGFbeta1, and Cdx2 mRNA and
the protein levels, compared to the phase III fasting values. The concomitant
decreases in cytokines and Cdx2 and in apoptotic cells during phase III suggest
the preservation of enterocytes during this critical fasting period in order to
optimize nutrient absorption as soon as food is available and thus, to rapidly
restore body mass.


  Most primates live in social groups which survival and stability depend on
individuals' abilities to create strong social relationships with other group
members. The existence of those groups requires to identify individuals and to
assign to each of them a social status. Individual recognition can be achieved
through vocalizations but also through faces. In humans, an efficient system
for the processing of own species faces exists. This specialization is achieved
through experience with faces of conspecifics during development and leads to
the loss of ability to process faces from other primate species. We hypothesize
that a similar mechanism exists in social primates. We investigated face
processing in one Old World species (genus Macaca) and in one New World species
(genus Cebus). Our results show the same advantage for own species face
recognition for all tested subjects. This work suggests in all species tested
the existence of a common trait inherited from the primate ancestor: an
efficient system to identify individual faces of own species only.


  It is shown that the dual to the linear programming problem that arises in
constraint-based models of metabolism can be given a thermodynamic
interpretation in which the shadow prices are chemical potential analogues, and
the objective is to minimise free energy consumption given a free energy drain
corresponding to growth. The interpretation is distinct from conventional
non-equilibrium thermodynamics, although it does satisfy a minimum entropy
production principle. It can be used to motivate extensions of constraint-based
modelling, for example to microbial ecosystems.


  The discovery of motifs underlying gene expression is a challenging one. Some
of these motifs are known transcription factors, but sequence inspection often
provides valuable clues, even discovery of novel motifs with uncharacterized
function in gene expression. Coupled with the complexity underlying
tissue-specific gene expression, there are several motifs that are putatively
responsible for expression in a certain cell type. This has important
implications in understanding fundamental biological processes, such as
development and disease progression. In this work, we present an approach to
the principled selection of motifs (not necessarily transcription factor sites)
and examine its application to several questions in current bioinformatics
research.
  There are two main contributions of this work: Firstly, we introduce a new
metric for variable selection during classification, and secondly, we
investigate a problem of finding specific sequence motifs that underlie tissue
specific gene expression. In conjunction with the SVM classifier we find these
motifs and discover several novel motifs which have not yet been attributed
with any particular functional role (eg: TFBS binding motifs). We hypothesize
that the discovery of these motifs would enable the large-scale investigation
for the tissue specific regulatory potential of any conserved sequence element
identified from genome-wide studies.
  Finally, we propose the utility of this developed framework to not only aid
discovery of discriminatory motifs, but also to examine the role of any motif
of choice in co-regulation or co-expression of gene groups.


  In primate groups, collective movements are typically described as processes
dependent on leadership mechanisms. However, in some species, decision-making
includes negotiations and distributed leadership. These facts suggest that
simple underlying processes may explain certain decision mechanisms during
collective movements. To study such processes, we have designed experiments on
white-faced capuchin monkeys (Cebus capucinus) during which we provoked
collective movements involving a binary choice. These experiments enabled us to
analyse the spatial decisions of individuals in the group.We found that the
underlying process includes anonymous mimetism, which means that each
individual may influence all members of the group. To support this result, we
created a mathematical model issued from our experimental data. A totally
anonymous model does not fit perfectly with our experimental distribution. A
more individualised model, which takes into account the specific behaviour of
social peripheral individuals, revealed the validity of the mimetism
hypothesis. Even though white-faced capuchins have complex cognitive abilities,
a coexistence of anonymous and social mechanisms appears to influence their
choice of direction during collective movements. The present approach may offer
vital insights into the relationships between individual behaviours and their
emergent collective acts.


  The study of synchronization among genetic oscillators is essential for the
understanding of the rhythmic phenomena of living organisms at both molecular
and cellular levels. Genetic networks are intrinsically noisy due to natural
random intra- and inter-cellular fluctuations. Therefore, it is important to
study the effects of noise perturbation on the synchronous dynamics of genetic
oscillators. From the synthetic biology viewpoint, it is also important to
implement biological systems that minimizing the negative influence of the
perturbations. In this paper, based on systems biology approach, we provide a
general theoretical result on the synchronization of genetic oscillators with
stochastic perturbations. By exploiting the specific properties of many genetic
oscillator models, we provide an easy-verified sufficient condition for the
stochastic synchronization of coupled genetic oscillators, based on the Lur'e
system approach in control theory. A design principle for minimizing the
influence of noise is also presented. To demonstrate the effectiveness of our
theoretical results, a population of coupled repressillators is adopted as a
numerical example. In summary, we present an efficient theoretical method for
analyzing the synchronization of genetic oscillator networks, which is helpful
for understanding and testing the synchronization phenomena in biological
organisms. Besides, the results are actually applicable to general oscillator
networks.


  All living cells transport molecules and ions across membranes, often against
concentration gradients. This active transport requires continual energy
expenditure and is clearly a nonequilibrium process for which standard
equilibrium thermodynamics is not rigorously applicable. Here we derive a
nonequilibrium effective potential that evaluates the per particle transport
energy invested by the membrane. A novel method is used whereby a Hamiltonian
function is constructed using particle concentrations as generalized
coordinates. The associated generalized momenta are simply related to the
individual particle energy from which we identify the effective potential.
Examples are given and the formalism is compared with the equilibrium Gibb's
free energy.


  We generalize a sampling algorithm for lattice animals (connected clusters on
a regular lattice) to a Monte Carlo algorithm for `graph animals', i.e.
connected subgraphs in arbitrary networks. As with the algorithm in [N. Kashtan
et al., Bioinformatics 20, 1746 (2004)], it provides a weighted sample, but the
computation of the weights is much faster (linear in the size of subgraphs,
instead of super-exponential). This allows subgraphs with up to ten or more
nodes to be sampled with very high statistics, from arbitrarily large networks.
Using this together with a heuristic algorithm for rapidly classifying
isomorphic graphs, we present results for two protein interaction networks
obtained using the TAP high throughput method: one of Escherichia coli with 230
nodes and 695 links, and one for yeast (Saccharomyces cerevisiae) with roughly
ten times more nodes and links. We find in both cases that most connected
subgraphs are strong motifs (Z-scores >10) or anti-motifs (Z-scores <-10) when
the null model is the ensemble of networks with fixed degree sequence. Strong
differences appear between the two networks, with dominant motifs in E. coli
being (nearly) bipartite graphs and having many pairs of nodes which connect to
the same neighbors, while dominant motifs in yeast tend towards completeness or
contain large cliques. We also explore a number of methods that do not rely on
measurements of Z-scores or comparisons with null models. For instance, we
discuss the influence of specific complexes like the 26S proteasome in yeast,
where a small number of complexes dominate the $k$-cores with large k and have
a decisive effect on the strongest motifs with 6 to 8 nodes. We also present
Zipf plots of counts versus rank. They show broad distributions that are not
power laws, in contrast to the case when disconnected subgraphs are included.


  Mechanosensation is a key part of the sensory repertoire of a vast array of
different cells and organisms. The molecular dissection of the origins of
mechanosensation is rapidly advancing as a result of both structural and
functional studies. One intriguing mode of mechanosensation results from
tension in the membrane of the cell (or vesicle) of interest. The aim of this
review is to catalogue recent work that uses a mix of continuum and statistical
mechanics to explore the role of the lipid bilayer in the function of
mechanosensitive channels that respond to membrane tension. The role of bilayer
deformation will be explored in the context of the well known mechanosensitive
channel MscL. Additionally, we make suggestions for bridging gaps between our
current theoretical understanding and common experimental techniques.


  Biological membranes are elastic media in which the presence of a
transmembrane protein leads to local bilayer deformation. The energetics of
deformation allow two membrane proteins in close proximity to influence each
other's equilibrium conformation via their local deformations, and spatially
organize the proteins based on their geometry. We use the mechanosensitive
channel of large conductance (MscL) as a case study to examine the implications
of bilayer-mediated elastic interactions on protein conformational statistics
and clustering. The deformations around MscL cost energy on the order of 10 kT
and extend ~3nm from the protein edge, as such elastic forces induce
cooperative gating and we propose experiments to measure these effects.
Additionally, since elastic interactions are coupled to protein conformation,
we find that conformational changes can severely alter the average separation
between two proteins. This has important implications for how conformational
changes organize membrane proteins into functional groups within membranes.


  Instructive influence of environment on heredity has been a debated topic for
centuries. Darwin's identification of natural selection coupled to chance
variation as the driving force for evolution, against a formal interpretation
proposed by Lamarck, convinced most scientists that environment does not
specifically instruct evolution in an oriented direction. This is true for
multicellular organisms. In contrast, bacteria were long thought of as prone to
receive oriented influences from their environment, although much was in favour
of the Darwinian route (1). In this context Cairns et al. raised a passionate
debate by suggesting that bacteria generate mutations oriented by the
environmental conditions (2). Several independent pieces of work subsequently
demonstrated that mutations overcoming specific defects arised as a consequence
of cultivation on specific media (3-7). Two diametrically opposed
interpretations were proposed to explain these observations : either induction
of mutations instructed by the environment (e.g. by a process involving a
putative reverse transcription) or selection of variants among a large set of
mutant bacteria generated when stress conditions are present. The experiments
presented below indicate that the Darwinian paradigm is the most plausible.


  We study the properties of Lotka-Volterra competitive models in which the
intensity of the interaction among species depends on their position along an
abstract niche space through a competition kernel. We show analytically and
numerically that the properties of these models change dramatically when the
Fourier transform of this kernel is not positive definite, due to a pattern
forming instability. We estimate properties of the species distributions, such
as the steady number of species and their spacings, for different types of
kernels.


  This work analyzed genome-wide nucleotide distribution patterns in ten insect
genomes. Two internal measures were applied: (i) GC variation and (ii) third
codon nucleotide preference. Although the genome size and overall GC level did
not show any correlation with insect order, the internal measures usually
displayed higher levels of consistency. GC variations in genomes of
hymenopteran insects, honeybee and wasp, ranked highest among all eukaryotic
genomes analyzed by us. Genomes of honeybee and beetle, insects of different
orders with similar overall GC levels, showed significant internal differences.
Honeybee genome stood out as unusual due to its high GC variation and
'left-handed' gene locations.


  Among all insect genomes, honeybee displays one of the most unusual patterns
with interspersed long AT and GC-rich segments. Nearly 75% of the
protein-coding genes are located in the AT-rich segments of the genome, but the
biological significance of the GC-rich regions is not well understood. Based on
an observation that the bee miRNAs, actins and tubulins are located in the
GC-rich segments, this work investigated whether other highly conserved genomic
regions show similar preferences. Sequences ultraconserved between the genomes
of honeybee and Nasonia, another hymenopteran insect, were determined. They
showed strong preferences towards locating in the GC-rich regions of the bee
genome.


  The management of endangered species as metapopulations is becoming
increasingly common. Diverse aspects of metapopulation dynamics and management
have received attention in recent years. In particular, translocation of
individuals between subpopulations of a metapopulation is practiced, or
envisaged, for a variety of reasons and requires careful consideration.
Linklater (2003) proposed that the number of individuals of each sex
translocated into a target population for the purposes of maintaining genetic
diversity could be chosen on the basis of parental investment theory. In this
paper, following basic ideas in the parental investment literature, I propose a
simple model to capture Linklater's proposal and provide a thorough
mathematical analysis of the model. Granted the necessary species-specific
biological information which would determine the model parameters in any
instance of application, the analysis indicates that a practical algorithm can
be constructed to generate the model's predictions for the optimal
translocation.


  The induction of the lac operon follows cooperative kinetics.The first
mechanistic model of these kinetics is the de facto standard in the modeling
literature (Yagil & Yagil, Biophys J, 11, 11-27, 1971). Yet, subsequent studies
have shown that the model is based on incorrect assumptions. Specifically, the
repressor is a tetramerwith four (not two) inducer-binding sites, and the
operon contains two auxiliary operators (in addition to the main operator).
Furthermore, these structural features are crucial for the formation of DNA
loops, the key determinants of lac repression and induction. Indeed, the
repression is determined almost entirely (>95%) by the looped complexes (Oehler
et al, EMBO J, 13, 3348, 1990), and the pronounced cooperativity of the
induction curve hinges upon the existence of the looped complexes (Oehler et
al, Nucleic Acids Res, 34, 606, 2006). Here, we formulate a model of lac
induction taking due account of the tetrameric structure of the repressor and
the existence of looped complexes. We show that: (1) The kinetics are
significantly more cooperative than those predicted by the Yagil & Yagil model.
(2) The model provides good fits to the repression data for cells containing
tetrameric (or mutant dimeric) repressor, as well as the induction curves for 6
different strains of E. coli. It also implies that the ratios of certain looped
and non-looped complexes are independent of inducer and repressor levels, a
conclusion that can be rigorously tested by gel electrophoresis. (3) Repressor
overexpression dramatically increases the cooperativity of the induction curve.
This suggests that repressor overexpression can induce bistability in systems,
such as growth of E. coli on lactose, that are otherwise monostable.


  Time evolution of number of species (genera, families, and others),
population of them, and size distribution of present ones and life times are
studied in terms of a new model, where population of each genetic taxon
increases by a (random) rate and decreases by (random) division
(fragmentation). Each of log-normal, exponential and power law distributions in
various empirical data about mentioned events are considered for gradual,
intermittent and abrupt evolution cases. Results come out in good agreement
with several theoretical literatures and with available empirical data. Sum of
exponential functions with independent random amplitude and exponent are shown
to be useful to express many aspects of the present competition, and they are
shown here to involve power law with power -1.


  The population in the sexual Penna ageing model is first separated into
several reproductively isolated groups. Then, after equilibration, sexual
mixing between the groups is allowed. We study the changes in the population
size due to this mixing and interpret them through a counterplay of purifying
selection and of haplotype complementarity.


  Biological information processing as implemented by regulatory and signaling
networks in living cells requires sufficient specificity of molecular
interaction to distinguish signals from one another, but much of regulation and
signaling involves somewhat fuzzy and promiscuous recognition of molecular
sequences and structures, which can leave systems vulnerable to crosstalk. This
paper examines a simple computational model of protein-protein interactions
which reveals both a sharp onset of crosstalk and a fragmentation of the
neutral network of viable solutions as more proteins compete for regions of
sequence space, revealing intrinsic limits to reliable signaling in the face of
promiscuity. These results suggest connections to both phase transitions in
constraint satisfaction problems and coding theory bounds on the size of
communication codes.


  In this paper, we have studied, within a discrete WLC model, the spatial
constraints in magnetic tweezers used in single molecule experiments. Two
elements are involved: first, the fixed plastic slab on which is stuck the
initial strand, second, the magnetic bead which pulls (or twists) the attached
molecule free end. We have shown that the bead surface can be replaced by its
tangent plane at the anchoring point, when it is close to the bead south pole
relative to the force. We are led to a model with two parallel repulsive
plates: the fixed anchoring plate and a fluctuating plate, simulating the bead,
in thermal equilibrium with the system. The bead effect is a slight upper shift
of the elongation, about four times smaller than the similar effect induced by
the fixed plate. This rather unexpected result, has been qualitatively
confirmed within the soluble Gaussian model. A study of the molecule elongation
versus the countour length exhibits a significant non-extensive behaviour. The
curve for short molecules (with less than 2 kbp) is well fitted by a straight
line, with a slope given by the WLC model, but it does not go through the
origin. The non-extensive offset gives a 15% upward shift to the elongation of
a 2 kbp molecule stretched by a 0.3 pN force.


  Motivation: In microarray analysis, special consideration must be given to
the issues of multiple statistical tests and typically p-values are adjusted to
control family-wise error rate (FWER) or false discovery rate (FDR). FDR
metrics have been suggested for controlling false positives, however, genes
with p-values close to the threshold typically have a higher chance of being
false positives than genes with very low p-values. The local FDR (fdr) metric
gives the number of false positives in the vicinity of a test statistic. We
propose a new fdr estimator that uses windows instead of binsand define
heuristics that use the fluctuations in the estimator to determine robust
thresholds for classifying differential expression.
  Results: Our fdr approach estimates the false discovery rate within a window
of p-values. We present heuristics that derive robust fdr thresholds such that
a significant change in the fdr threshold yields a small change in the number
of rejected hypotheses. We compare these thresholds with thresholds from other
approaches using two simulated datasets and one cancer microarray dataset. In
the latter, our estimator finds two robust thresholds. Since our fdr estimator
is an extension of the FDR metric, it can be used with many FDR estimation
methods.
  Availability: An R function implementing the proposed estimator is available
at http://www.dbi.tju.edu/dbi/tools/fdr
  Contact: james.schwaber@jefferson.edu
  Supplementary Information: Supplementary figures and code are available at
http://www.dbi.tju.edu/dbi/tools/fdr


  The question of how to validate and interpret correctly the waiting time
probability density functions (WT-PDFs) from single molecule data is addressed.
It is shown by simulation that when a stretched exponential WT-PDF, with a
stretched exponent alfa and a time scale parameter tau, generates the off
periods of a two-state trajectory, a reliable recovery of the input WT-PDF from
the trajectory is obtained even when the bin size used to define the
trajectory, dt, is much larger than the parameter tau. This holds true as long
as the first moment of the WT-PDF is much larger than dt. Our results validate
the results in an earlier study of the activity of single Lipase B molecules
and disprove recent related critique.


  The precise layout of gene expression patterns is a crucial step in
development. Formation of a sharp boundary between high and low expression
domains requires a genetic mechanism which is both sensitive and robust to
fluctuations, a demand that may not be easily achieved by morphogens alone.
Recently it has been demonstrated that small RNAs (and, in particular,
microRNAs) play many roles in embryonic development. While some RNAs are
essential for embryogenesis, others are limited to fine-tuning a predetermined
gene expression pattern. Here we explore the possibility that small RNAs
participate in sharpening a gene expression profile that was crudely
established by a morphogen. To this end we study a model where small RNAs
interact with a target gene and diffusively move from cell to cell. Though
diffusion generally smears spatial expression patterns, we find that
intercellular mobility of small RNAs is actually critical in sharpening the
interface between target expression domains in a robust manner. We discuss the
applicability of our results, as examples, to the case of leaf polarity
establishment in maize and Hox patterning in the early {\it Drosophila} embryo.
Our findings point out the functional significance of some mechanistic
properties, such as mobility of small RNAs and the irreversibility of their
interactions. These properties are yet to be established directly for most
classes of small RNAs. An indirect yet simple experimental test of the proposed
mechanism is suggested in some detail.