  Complex dynamical systems are often modeled as networks, with nodes
representing dynamical units which interact through the network's links. Gene
regulatory networks, responsible for the production of proteins inside a cell,
are an example of system that can be described as a network of interacting
genes. The behavior of a complex dynamical system is characterized by
cooperativity of its units at various scales, leading to emergent dynamics
which is related to the system's function. Among the key problems concerning
complex systems is the issue of stability of their functioning, in relation to
different internal and external interaction parameters.
  In this Thesis we study two-dimensional chaotic maps coupled through
non-directed networks with different topologies. We use a non-symplectic
coupling which involves a time delay in the interaction among the maps. We test
the stability of network topologies through investigation of their collective
motion, done by analyzing the departures from chaotic nature of the isolated
units. The study is done on two network scales: (a) full-size networks (a
computer generated scalefree tree and a tree with addition of cliques); (b)
tree's characteristic sub-graph 4-star, as a tree's typical dynamical motif
which captures its topology in smallest possible number of nodes and is
suitable for time-delayed interaction. We study the dynamical relationship
between these two network structures, examining the emergence of cooperativity
on a large scale (trees) as a consequence of mesoscale dynamical patterns
exhibited by the 4-star. (FULL ABSTRACT INSIDE THE TEXT)


  Cooperative binding of transcription factors (TFs) to cis-regulatory regions
(CRRs) is essential for precision of gene expression in development and other
processes. The classical model of cooperativity requires direct interactions
between TFs, thus constraining the arrangement of TFs sites in a CRR. On the
contrary, genomic and functional studies demonstrate a great deal of
flexibility in such arrangements with variable distances, numbers of sites, and
identities of the involved TFs. Such flexibility is inconsistent with the
cooperativity by direct interactions between TFs. Here we demonstrate that
strong cooperativity among non-interacting TFs can be achieved by their
competition with nucleosomes. We find that the mechanism of nucleosome-mediated
cooperativity is mathematically identical to the Monod-Wyman-Changeux (MWC)
model of cooperativity in hemoglobin. This surprising parallel provides deep
insights, with parallels between heterotropic regulation of hemoglobin (e.g.
Bohr effect) and roles of nucleosome-positioning sequences and chromatin
modifications in gene regulation. Characterized mechanism is consistent with
numerous experimental results, allows substantial flexibility in and modularity
of CRRs, and provides a rationale for a broad range of genomic and evolutionary
observations. Striking parallels between cooperativity in hemoglobin and in
transcription regulation point at a new design principle that may be used in
range of biological systems.


  The aim of this work is to make available to the community a large collection
of mass-action reaction networks of a given size for further research. The set
is limited to what can be computed on a modern multi-core desktop in reasonable
time (< 20 days). We have currently generated over 47 million unique reaction
networks. All currently generated sets of networks are available and as new
sets are completed they will also be made available. Also provided are programs
for translating them into different formats, along with documentation and
examples. Source code and binaries for all the programs are included. These can
be downloaded from (http://www.sys-bio.org/networkenumeration). This library of
networks will allow for thorough studies of the reaction network space.
Additionally, these methods serve as an example for future work on enumerating
other types of biological networks, such as genetic regulatory networks and
mass-action networks that include regulation.


  Living cells exhibit an important out-of-equilibrium mechanical activity,
mainly due to the forces generated by molecular motors. These motor proteins,
acting individually or collectively on the cytoskeleton, contribute to the
violation of the fluctuation-dissipation theorem in living systems. In this
work we probe the cytoskeletal out-of-equilibrium dynamics by performing
simultaneous active and passive microrheology experiments, using the same
micron-sized probe specifically bound to the actin cortex. The free motion of
the probe exhibits a constrained, subdiffusive behavior at short time scales (t
< 2s), and a directed, superdiffusive behavior at larger time scales, while, in
response to a step force, its creep function presents the usual weak power law
dependence with time. Combining the results of both experiments, we precisely
measure for the first time the power spectrum of the force fluctuations exerted
on this probe, which lies more than one order of magnitude above the spectrum
expected at equilibrium, and greatly depends on frequency. We retrieve an
effective temperature Teff of the system, as an estimate of the departure from
thermal equilibrium. This departure is especially pronounced on long time
scales, where Teff bears the footprint of the cooperative activity of motors
pulling on the actin network. ATP depletion reduces the fluctuating force
amplitude and results in a sharp decrease of Teff towards equilibrium.


  One of the longstanding debates in the history of paleontology focuses on the
issue of whether or not there have been long term cycles (operating over tens
of millions of years) in biodiversity and extinction. Here we consider the
history of this debate by connecting the skein from Grabau up to 2008. We focus
on the evidence for periodicity that has emerged thus far, and conclude that
there is indeed some evidence that periodicity may be real, though of course
more work is needed. We also comment on possible causal mechanisms, focusing
especially on the motion of our solar system in the Galaxy. Moreover, we
consider the reasons why some scientists have opposed periodicity over the
years. Finally, we consider the significance of this for our understanding of
evolution and the history of life.


  We describe an efficient algorithm for determining exactly the minimum number
of sires consistent with the multi-locus genotypes of a mother and her progeny.
We consider cases where a simple exhaustive search through all possible sets of
sires is impossible in practice (because it would take too long to complete).
Our algorithm for solving this combinatorial optimisation problem avoids
visiting large parts of search space which would not improve the solution found
so far (i.e., result in a solution with fewer number of sires). This is of
particular importance when the number of allelic types in the progeny array is
large and when the minimum number of sires is expected to be large. Precisely
in such cases it is important to know the minimum number of sires: this number
gives an exact bound on the most likely number of sires estimated by a random
search algorithm in a parameter region where it may be difficult to determine
whether it has converged. We apply our algorithm to data from the marine snail,
Littorina saxatilis.


  In a population with haploid reproduction any individual has a single parent
in the previous generation. If all genealogical distances among pairs of
individuals (generations from the closest common ancestor) are known it is
possible to exactly reconstruct their genealogical tree. Unfortunately, in most
cases, genealogical distances are unknown and only genetic distances are
available. The genetic distance between two individuals is measurable from
differences in mtDNA (mitochondrial DNA) in the case of humans or other complex
organisms while an analogous distance can be also given for languages where it
is measured from lexical differences. Assuming a constant rate of mutation,
these genetic distances are random and proportional only on average to
genealogical ones. The reconstruction of the genealogical tree from the
available genetic distances is forceful imprecise. In this paper we try to
quantify the error one may commit in the reconstruction of the tree for
different degrees of randomness. The errors may concern both topology of the
tree (the branching hierarchy) and, in case of correct topology, the
proportions of the tree (length of various branches).


  The effects of intrinsic noise on stochastic delay systems is studied within
an expansion in the inverse system size. We show that the stochastic nature of
the underlying dynamics may induce oscillatory behaviour in parameter ranges
where the deterministic system does not sustain cycles, and compute the power
spectra of these stochastic oscillations analytically, in good agreement with
simulations. The theory is developed in the context of a simple one-dimensional
toy model, but is applicable more generally. Gene regulatory systems in
particular often contain only a small number of molecules, leading to
significant fluctuations in mRNA and protein concentrations. As an application
we therefore study a minimalistic model of the expression levels of hes1 mRNA
and Hes1 protein, representing the simple motif of an auto-inhibitory feedback
loop and motivated by its relevance to somite segmentation.


  Analytical analysis of spatially extended autocatalytic and hypercyclic
systems is presented. It is shown that spatially explicit systems in the form
of reaction-diffusion equations with global regulation possess the same major
qualitative features as the corresponding local models. In particular, using
the introduced notion of the stability in the mean integral sense we prove the
competitive exclusion principle for the autocatalytic system and the permanence
for the hypercycle system. Existence and stability of stationary solutions are
studied. For some parameter values it is proved that stable spatially
non-uniform solutions appear.


  We exhibit a stochastic discrete time model that has exactly the Eigen model
as its deterministic continuous limit. Such model can be divided into two
phases: reproduction followed by neutral selection. This result suggests that
Eigen model describes the competition among individuals differing for
reproductive capability but equivalent as survivors. We explicitly write down
the Markov matrix of the stochastic model in the two species case and compute
numerically the master sequence concentration for various values of the total
population. Finally we compare our results with those obtained with Eigen model
and with Nowak and Schuster birth and death model.


  Oxidative DNA single strand breaks arise continuously in cells and defects in
their repair have been implicated in neurological disease. While much progress
has been made in understanding how chromosomal single strand breaks are
repaired little is known about the changes chromatin structure that accompany
this process. Here, we show that nascent recombinant histone H3.1 protein
accumulates and is deposited into chromatin at sites of DNA strand breakage in
quiescent human cells following oxidative stress, and that core components of
the single-strand break repair machinery are required for this process. We show
that the SSBR sensor and scaffold proteins poly (ADP-ribose) polymerase and
XRCC1 facilitate accumulation of chromatin assembly factor-1 (CAF-1) at sites
of oxidative DNA strand breakage, which in turn mediates the deposition of
Histone H3.1. We also demonstrate that depletion of CAF-1 slows global rates of
DNA strand break repair in quiescent cells following oxidative stress,
demonstrating that single-strand break repair and histone deposition are
tightly coordinated processes. These data describe a novel role for the DNA
singlestrand break repair machinery and implicate histone turnover as a core
component of the cellular response of quiescent cells to oxidative damage.


  In biology, information flows from the environment to the genome by the
process of natural selection. But it has not been clear precisely what sort of
information metric properly describes natural selection. Here, I show that
Fisher information arises as the intrinsic metric of natural selection and
evolutionary dynamics. Maximizing the amount of Fisher information about the
environment captured by the population leads to Fisher's fundamental theorem of
natural selection, the most profound statement about how natural selection
influences evolutionary dynamics. I also show a relation between Fisher
information and Shannon information (entropy) that may help to unify the
correspondence between information and dynamics. Finally, I discuss possible
connections between the fundamental role of Fisher information in statistics,
biology, and other fields of science.


  Metabolic pathways describe chains of enzymatic reactions. Their modelling is
a key point to understand living systems. An enzymatic reaction is an
interaction between one or several metabolites (substrates) and an enzyme
(simple protein or enzymatic complex build of several subunits). In our
Mitochondria in Silico Project, MitoScop, we study the metabolism of the
mitochondria, an intra-cellular organelle. Many ordinary differential equation
models are available in the literature. They well fit experimental results on
flux values inside the metabolic pathways, but many parameters are di$\pm$cult
to transcribe with such models: localization of enzymes, rules about the
reactions scheduler, etc Moreover, a model of a significant part of
mitochondrial metabolism could become very complex and contain more than 50
equations. In this context, the multi-agents systems appear as an alternative
to model the metabolic pathways. Firstly, we have looked after membrane design.
The mitochondria is a particular case because the inner mitochondrial space, ie
matricial space, is delimited by two membranes: the inner and the outer one. In
addition to matricial enzymes, other enzymes are located inside the membranes
or in the inter-membrane space. Analysis of mitochondrial metabolism must take
into account this kind of architecture.


  Many studies on animal and human movement patterns report the existence of
scaling laws and power-law distributions. Whereas a number of random walk
models have been proposed to explain observations, in many situations
individuals actually rely on mental maps to explore strongly heterogeneous
environments. In this work we study a model of a deterministic walker, visiting
sites randomly distributed on the plane and with varying weight or
attractiveness. At each step, the walker minimizes a function that depends on
the distance to the next unvisited target (cost) and on the weight of that
target (gain). If the target weight distribution is a power-law, $p(k)\sim
k^{-\beta}$, in some range of the exponent $\beta$, the foraging medium induces
movements that are similar to L\'evy flights and are characterized by
non-trivial exponents. We explore variations of the choice rule in order to
test the robustness of the model and argue that the addition of noise has a
limited impact on the dynamics in strongly disordered media.


  Genomic expression depends critically both on the ability of regulatory
proteins to locate specific target sites on a DNA within seconds and on the
formation of long lived (many minutes) complexes between these proteins and the
DNA. Equilibrium experiments show that indeed regulatory proteins bind tightly
to their target site. However, they also find strong binding to other
non-specific sites which act as traps that can dramatically increase the time
needed to locate the target. This gives rise to a conflict between the speed
and stability requirements. Here we suggest a simple mechanism which can
resolve this long-standing paradox by allowing the target sites to be located
by proteins within short time scales even in the presence of traps. Our
theoretical analysis shows that the mechanism is robust in the presence of
generic disorder in the DNA sequence and does not require a specially designed
target site.


  A novel group of small non-coding RNA, known as microRNA (miRNA) is predicted
to regulate as high as 90% of the coding genes in human. The diversity and
abundance of miRNA targets offer an enormous level of combinatorial
possibilities and suggest that miRNAs and their targets form a complex
regulatory network. In the present study, we analyzed 711 miRNAs and their 34,
525 predicted targets in the miRBase database which generate a complex
bipartite network having numerous numbers of genes forming the hub. Genes at
the hub (total 9877) are significantly over represented in genes with specific
molecular functions, biological processes and biological pathways as revealed
from the analysis using PANTHER. We further construct a miRNA co-target network
by linking every pair of miRNAs which co-target at least one gene. The weight
of the link, which is taken to be the number of co-targets of the pair of
miRNAs vary widely, and we could erase several links while keeping the relevant
features of the network intact. The largest connected sub-graph, thus obtained,
contains 479 miRNAs. More than 75% of the miRNAs deregulated in 15 different
diseases collected from published data are found to be in this largest sub
graph. We further analyze this sub-graph to obtain 70 small clusters containing
total 330 miRNAs of 479. We identified the biological pathways where the
co-targeted genes in the clusters are significantly over- represented in
comparison to that obtained with that are not co-targeted by the miRNAs in the
cluster. Using published data, we identified that specific clusters of miRNAs
are associated with specific diseases by altering particular pathways. We
propose that instead of single miRNA, clusters of miRNA that co-targets the
genes are important for the regulation of miRNA in diseases.


  Boolean networks have been proposed as potentially useful models for genetic
control. An important aspect of these networks is the stability of their
dynamics in response to small perturbations. Previous approaches to stability
have assumed uncorrelated random network structure. Real gene networks
typically have nontrivial topology significantly different from the random
network paradigm. In order to address such situations, we present a general
method for determining the stability of large Boolean networks of any specified
network topology and predicting their steady-state behavior in response to
small perturbations. Additionally, we generalize to the case where individual
genes have a distribution of `expression biases,' and we consider
non-synchronous update, as well as extension of our method to non-Boolean
models in which there are more than two possible gene states. We find that
stability is governed by the maximum eigenvalue of a modified adjacency matrix,
and we test this result by comparison with numerical simulations. We also
discuss the possible application of our work to experimentally inferred gene
networks.


  This paper addresses the problem of finding cycles in the state transition
graphs of synchronous Boolean networks. Synchronous Boolean networks are a
class of deterministic finite state machines which are used for the modeling of
gene regulatory networks. Their state transition graph cycles, called
attractors, represent cell types of the organism being modeled. When the effect
of a disease or a mutation on an organism is studied, attractors have to be
re-computed every time a fault is injected in the model. We present an
algorithm for finding attractors which uses a SAT-based bounded model checking.
Novel features of the algorithm compared to the traditional SAT-based bounded
model checking approaches are: (1) a termination condition which does not
require an explicit computation of the diameter and (2) a technique to reduce
the number of additional clauses which are needed to make paths loop-free. The
presented algorithm uses much less space than existing BDD-based approaches and
has a potential to handle several orders of magnitude larger networks.


  Excitation transfer through interacting systems plays an important role in
many areas of physics, chemistry, and biology. The uncontrollable interaction
of the transmission network with a noisy environment is usually assumed to
deteriorate its transport capacity, especially so when the system is
fundamentally quantum mechanical. Here we identify key mechanisms through which
noise such as dephasing, perhaps counter intuitively, may actually aid
transport through a dissipative network by opening up additional pathways for
excitation transfer. We show that these are processes that lead to the
inhibition of destructive interference and exploitation of line broadening
effects. We illustrate how these mechanisms operate on a fully connected
network by developing a powerful analytical technique that identifies the
invariant (excitation trapping) subspaces of a given Hamiltonian. Finally, we
show how these principles can explain the remarkable efficiency and robustness
of excitation energy transfer from the light-harvesting chlorosomes to the
bacterial reaction center in photosynthetic complexes and present a numerical
analysis of excitation transport across the Fenna-Matthew-Olson (FMO) complex
together with a brief analysis of its entanglement properties. Our results show
that, in general, it is the careful interplay of quantum mechanical features
and the unavoidable environmental noise that will lead to an optimal system
performance.


  Analogous to biological sequence comparison, comparing cellular networks is
an important problem that could provide insight into biological understanding
and therapeutics. For technical reasons, comparing large networks is
computationally infeasible, and thus heuristics such as the degree distribution
have been sought. It is easy to demonstrate that two networks are different by
simply showing a short list of properties in which they differ. It is much
harder to show that two networks are similar, as it requires demonstrating
their similarity in all of their exponentially many properties. Clearly, it is
computationally prohibitive to analyze all network properties, but the larger
the number of constraints we impose in determining network similarity, the more
likely it is that the networks will truly be similar.
  We introduce a new systematic measure of a network's local structure that
imposes a large number of similarity constraints on networks being compared. In
particular, we generalize the degree distribution, which measures the number of
nodes 'touching' k edges, into distributions measuring the number of nodes
'touching' k graphlets, where graphlets are small connected non-isomorphic
subgraphs of a large network. Our new measure of network local structure
consists of 73 graphlet degree distributions (GDDs) of graphlets with 2-5
nodes, but it is easily extendible to a greater number of constraints (i.e.
graphlets). Furthermore, we show a way to combine the 73 GDDs into a network
'agreement' measure. Based on this new network agreement measure, we show that
almost all of the 14 eukaryotic PPI networks, including human, are better
modeled by geometric random graphs than by Erdos-Reny, random scale-free, or
Barabasi-Albert scale-free networks.


  It has been recently argued that the depletion attraction may play an
important role in different aspects of the cellular organization, ranging from
the organization of transcriptional activity in transcription factories to the
formation of the nuclear bodies. In this paper we suggest a new application of
these ideas in the context of the splicing process, a crucial step of messanger
RNA maturation in Eukaryotes. We shall show that entropy effects and the
resulting depletion attraction may explain the relevance of the aspecific
intron length variable in the choice of the splice-site recognition modality.
On top of that, some qualitative features of the genome architecture of higher
Eukaryotes can find an evolutionary realistic motivation in the light of our
model.


  Inspired by problems in biochemical kinetics, we study statistical properties
of an overdamped Langevin process whose friction coefficient depends on the
state of a similar, unobserved process. Integrating out the latter, we derive
the long time behaviour of the mean square displacement. Anomalous diffusion is
found. Since the diffusion exponent can not be predicted using a simple scaling
argument, anomalous scaling appears as well. We also find that the coupling can
lead to ergodic or non-ergodic behaviour of the studied process. We compare our
theoretical predictions with numerical simulations and find an excellent
agreement. The findings caution against treating biochemical systems coupled
with unobserved dynamical degrees of freedom by means of standard, diffusive
Langevin descriptions.


  Translocation through a nanopore is a new experimental technique to probe
physical properties of biomolecules. A bulk of theoretical and computational
work exists on the dependence of the time to translocate a single unstructured
molecule on the length of the molecule. Here, we study the same problem but for
RNA molecules for which the breaking of the secondary structure is the main
barrier for translocation. To this end, we calculate the mean translocation
time of single-stranded RNA through a nanopore of zero thickness and at zero
voltage for many randomly chosen RNA sequences. We find the translocation time
to depend on the length of the RNA molecule with a power law. The exponent
changes as a function of temperature and exceeds the naively expected exponent
of two for purely diffusive transport at all temperatures. We interpret the
power law scaling in terms of diffusion in a one-dimensional energy landscape
with a logarithmic barrier.


  The purpose of the present study is to evaluate changes in chromatin of human
buccal epithelium under the influence of stressing factor - dosed physical
activity. Investigations were performed in a group of students (13 men) of age
19-23. Cells were stained on a slide by a 2% orcein solution in 45% acetic acid
during 1 h. The following physiological indexes were determined: arterial blood
pressure, pulse frequency, and frequency of breathing. The physical stress
produced by the dosed physical activity causes the considerable increase of
degree of heterochromatinization in the cell nuclei of human buccal epithelium.
As a rule, the level of heterochromatinization increases after first stage of
training, but in some donors it increases significantly only after the second
stage of training.


  We investigate the dynamics of the Fisher equation for the spreading of
micro-organisms in one dimension subject to both turbulent convection and
diffusion. We show that for strong enough turbulence, bacteria, for example,
track in a quasilocalized fashion (with remakably long persistance times) sinks
in the turbulent field. An important consequence is a large reduction in the
carrying capacity of the fluid medium. We determine analytically the regimes
where this quasi-localized behavior occurs and test our predictions by
numerical simulations.


  Color centers in diamond nanocrystals are a new class of fluorescence markers
that attract significant interest due to matchless brightness, photostability
and biochemical inertness. Fluorescing diamond nanocrystals containing defects
can be used as markers replacing conventional organic dye molecules, quantum
dots or autofluorescent proteins. They can be applied for tracking and
ultrahigh-resolution localization of the single markers. In addition the spin
properties of diamond defects can be utilized for novel magneto-optical imaging
(MOI) with nanometer resolution. We develop this technique to unravel the
details of the rotary motions and the elastic energy storage mechanism of a
single biological nanomotor FoF1-ATP synthase. FoF1-ATP synthase is the enzyme
that provides the 'chemical energy currency' adenosine triphosphate, ATP, for
living cells. The formation of ATP is accomplished by a stepwise internal
rotation of subunits within the enzyme. Previously subunit rotation has been
monitored by single-molecule fluorescence resonance energy transfer (FRET) and
was limited by the photostability of the fluorophores. Fluorescent nanodiamonds
advance these FRET measurements to long time scales.


  Cellular cargo can be bound to cytoskeletal filaments by one or multiple
active or passive molecular motors. Recent experiments have shown that the
presence of auxiliary, nondriving motors, results in an enhanced processivity
of the cargo, compared to the case of a single active motor alone. We model the
observed cooperative transport process using a stochastic model that describes
the dynamics of two molecular motors, an active one that moves cargo
unidirectionally along a filament track and a passive one that acts as a
tether. Analytical expressions obtained from our analysis are fit to
experimental data to estimate the microscopic kinetic parameters of our model.
Our analysis reveals two qualitatively distinct processivity-enhancing
mechanisms: the passive tether can decrease the typical detachment rate of the
active motor from the filament track or it can increase the corresponding
reattachment rate. Our estimates unambiguously show that in the case of
microtubular transport, a higher average run length arises mainly from the
ability of the passive motor to keep the cargo close to the filament, enhancing
the reattachment rate of an active kinesin motor that has recently detached.
Instead, for myosin-driven transport along actin, the passive motor tightly
tethers the cargo to the filament, suppressing the detachment rate of the
active myosin.


  The conformational kinetics of enzymes can be reliably revealed when they are
governed by Markovian dynamics. Hidden Markov Models (HMMs) are appropriate
especially in the case of conformational states that are hardly
distinguishable. However, the evolution of the conformational states of
proteins mostly shows non-Markovian behavior, recognizable by
non-monoexponential state dwell time histograms. The application of a Hidden
Markov Model technique to a cyclic system demonstrating semi-Markovian dynamics
is presented in this paper and the required extension of the model design is
discussed. As standard ranking criteria of models cannot deal with these
systems properly, a new approach is proposed considering the shape of the dwell
time histograms. We observed the rotational kinetics of a single F1-ATPase
alpha3beta3gamma sub-complex over six orders of magnitude of different ATP to
ADP and Pi concentration ratios, and established a general model describing the
kinetics for the entire range of concentrations. The HMM extension described
here is applicable in general to the accurate analysis of protein dynamics.


  We derive exact solutions of simplified models for the temporal evolution of
the protein concentration within a cell population arbitrarily far from the
stationary state. We show that monitoring the dynamics can assist in modeling
and understanding the nature of the noise and its role in gene expression and
protein production. We introduce a new measure, the cell turnover distribution,
which can be used to probe the phase of transcription of DNA into messenger
RNA.


  We introduce a model of traveling agents ({\it e.g.} frugivorous animals) who
feed on randomly located vegetation patches and disperse their seeds, thus
modifying the spatial distribution of resources in the long term. It is assumed
that the survival probability of a seed increases with the distance to the
parent patch and decreases with the size of the colonized patch. In turn, the
foraging agents use a deterministic strategy with memory, that makes them visit
the largest possible patches accessible within minimal travelling distances.
The combination of these interactions produce complex spatio-temporal patterns.
If the patches have a small initial size, the vegetation total mass (biomass)
increases with time and reaches a maximum corresponding to a self-organized
critical state with power-law distributed patch sizes and L\'evy-like movement
patterns for the foragers. However, this state collapses as the biomass sharply
decreases to reach a noisy stationary regime characterized by corrections to
scaling. In systems with low plant competition, the efficiency of the foraging
rules leads to the formation of heterogeneous vegetation patterns with
$1/f^{\alpha}$ frequency spectra, and contributes, rather counter-intuitively,
to lower the biomass levels.


  We have developed a mathematical model of transcriptional activation by MarA
in Escherichia coli, and used the model to analyze measurements of
MarA-dependent activity of the marRAB, sodA, and micF promoters in mar-rob-
cells. The model rationalizes an unexpected poor correlation between the
mid-point of in vivo promoter activity profiles and in vitro equilibrium
constants for MarA binding to promoter sequences. Analysis of the promoter
activity data using the model yielded the following predictions regarding
activation mechanisms: (1) MarA activation of the marRAB, sodA, and micF
promoters involves a net acceleration of the kinetics of transitions after RNA
polymerase binding, up to and including promoter escape and message elongation;
(2) RNA polymerase binds to these promoters with nearly unit occupancy in the
absence of MarA, making recruitment of polymerase an insignificant factor in
activation of these promoters; and (3) instead of recruitment, activation of
the micF promoter might involve a repulsion of polymerase combined with a large
acceleration of the kinetics of polymerase activity. These predictions are
consistent with published chromatin immunoprecipitation assays of interactions
between polymerase and the E. coli chromosome. A lack of recruitment in
transcriptional activation represents an exception to the textbook description
of activation of bacterial sigma-70 promoters. However, use of accelerated
polymerase kinetics instead of recruitment might confer a competitive advantage
to E. coli by decreasing latency in gene regulation.


  A phylogenetic birth-and-death model is a probabilistic graphical model for a
so-called phylogenetic profile, i.e., the size distribution for a homolog gene
family at the terminal nodes of a phylogeny. Profile datasets are used in
bioinformatics analyses for the inference of evolutionary trees, and of
functional associations between gene families, as well as for the
quantification of various processes guiding genome evolution. Here we describe
the mathematical formalism for phylogenetic birth-and-death models. We also
present an algorithm for computing the likelihood in a gain-loss-duplication
model.


  Quantifying the adaptive landscape in a given dynamical processes has been
one of most important goals in theoretical biology. It can have immediate
implications for many dynamical properties, such as robustness and plasticity.
Based on recent work here I give a nontechnical brief review of this powerful
quantitative concept in biology. This concept was initially proposed by S
Wright 70 years ago, re-introduced by one of the founders of molecular biology
and by others in different biological contexts. It was apparently forgotten by
mainstream modern biologists for many years. Currently, this concept has found
its increasingly important role in the development of systems biology and the
modeling of bionetwork dynamics, from phage lambda genetic switch to endogenous
network of cancer genesis and progression. It is an ideal quantify to describe
the robustness and stability of bionetworks. I will first introduce five
landmark proposals in biology on this concept, to demonstrate the important
common thread in its theoretical biology development. Then I will discuss a few
recent results, focusing on the work showing the logical consistency of
adaptive landscape. From the perspective of a working scientist and of what
needed for a dynamical theory when confronting empirical data, the adaptive
landscape is useful both metaphorically and quantitatively and has captured an
essential aspect of biological dynamical processes. Still, many important open
problems remain to be solved. Having this important problem under control, we
may expect that we are on the right road to quantitatively formulate the
evolutionary dynamics discovered by Darwin and Wallace.


  We introduce a new model of proteins, which extends and enhances the
traditional graphical representation by associating a combinatorial object
called a fatgraph to any protein based upon its intrinsic geometry. Fatgraphs
can easily be stored and manipulated as triples of permutations, and these
methods are therefore amenable to fast computer implementation. Applications
include the refinement of structural protein classifications and the prediction
of geometric and other properties of proteins from their chemical structures.


  Understanding the organization of reaction fluxes in cellular metabolism from
the stoichiometry and the topology of the underlying biochemical network is a
central issue in systems biology. In this task, it is important to devise
reasonable approximation schemes that rely on the stoichiometric data only,
because full-scale kinetic approaches are computationally affordable only for
small networks (e.g. red blood cells, about 50 reactions). Methods commonly
employed are based on finding the stationary flux configurations that satisfy
mass-balance conditions for metabolites, often coupling them to local
optimization rules (e.g. maximization of biomass production) to reduce the size
of the solution space to a single point. Such methods have been widely applied
and have proven able to reproduce experimental findings for relatively simple
organisms in specific conditions. Here we define and study a constraint-based
model of cellular metabolism where neither mass balance nor flux stationarity
are postulated, and where the relevant flux configurations optimize the global
growth of the system. In the case of E. coli, steady flux states are recovered
as solutions, though mass-balance conditions are violated for some metabolites,
implying a non-zero net production of the latter. Such solutions furthermore
turn out to provide the correct statistics of fluxes for the bacterium E. coli
in different environments and compare well with the available experimental
evidence on individual fluxes. Conserved metabolic pools play a key role in
determining growth rate and flux variability. Finally, we are able to connect
phenomenological gene essentiality with `frozen' fluxes (i.e. fluxes with
smaller allowed variability) in E. coli metabolism.


  FoF1-ATP synthase is the enzyme that provides the 'chemical energy currency'
adenosine triphosphate, ATP, for living cells. The formation of ATP is
accomplished by a stepwise internal rotation of subunits within the enzyme.
Briefly, proton translocation through the membrane-bound Fo part of ATP
synthase drives a 10-step rotary motion of the ring of c subunits with respect
to the non-rotating subunits a and b. This rotation is transmitted to the gamma
and epsilon subunits of the F1 sector resulting in 120 degree steps. In order
to unravel this symmetry mismatch we monitor subunit rotation by a
single-molecule fluorescence resonance energy transfer (FRET) approach using
three fluorophores specifically attached to the enzyme: one attached to the F1
motor, another one to the Fo motor, and the third one to a non-rotating
subunit. To reduce photophysical artifacts due to spectral fluctuations of the
single fluorophores, a duty cycle-optimized alternating three-laser scheme
(DCO-ALEX) has been developed. Simultaneous observation of the stepsizes for
both motors allows the detection of reversible elastic deformations between the
rotor parts of Fo and F1.


  To what extent must a hydrophobic gate expand for the channel to count as
open? We address this question using the nicotinic acetylcholine receptor
(nAChR) as the exemplar. The nAChR is an integral membrane protein which forms
a cation selective channel gated by neurotransmitter binding to its
extracellular domain. A hydrophobic gating model has been proposed for the
nAChR, whereby the pore is incompletely occluded in the closed state channel,
with a narrow hydrophobic central gate region which presents an energetic
barrier to ion permeation. The nAChR pore is lined by a parallel bundle of five
M2 alpha-helices, with the gate formed by three rings of hydrophobic sidechains
(9', 13', and 17' of M2). A number of models have been proposed to describe the
nature of the conformational change underlying the closed to open transition of
the nAChR. These models involve different degrees of M2 helix displacement,
rotation, and/or kinking. In this study, we use a simple pore expansion method
(previously used to model opening of potassium channels) to generate a series
of progressively wider models of the nAChR transmembrane domain. Continuum
electrostatics calculations are used to assess the change in the barrier height
of the hydrophobic gate as a function of pore expansion. The results suggest
that an increase in radius of Delta r ~ 1.5 angstrom is sufficient to
functionally open the pore without, for example, a requirement for rotation of
the M2 helices. This is evaluated in the context of current mutational and
structural data on the nAChR and its homologues.


  Recent publications have described and applied a novel metric that quantifies
the genetic distance of an individual with respect to two population samples,
and have suggested that the metric makes it possible to infer the presence of
an individual of known genotype in a sample for which only the marginal allele
frequencies are known. However, the assumptions, limitations, and utility of
this metric remained incompletely characterized. Here we present an exploration
of the strengths and limitations of that method. In addition to analytical
investigations of the underlying assumptions, we use both real and simulated
genotypes to test empirically the method's accuracy. The results reveal that,
when used as a means by which to identify individuals as members of a
population sample, the specificity is low in several circumstances. We find
that the misclassifications stem from violations of assumptions that are
crucial to the technique yet hard to control in practice, and we explore the
feasibility of several methods to improve the sensitivity. Additionally, we
find that the specificity may still be lower than expected even in ideal
circumstances. However, despite the metric's inadequacies for identifying the
presence of an individual in a sample, our results suggest potential avenues
for future research on tuning this method to problems of ancestry inference or
disease prediction. By revealing both the strengths and limitations of the
proposed method, we hope to elucidate situations in which this distance metric
may be used in an appropriate manner. We also discuss the implications of our
findings in forensics applications and in the protection of GWAS participant
privacy.


  In this paper, we model the mechanics of a collagen pair in the connective
tissue extracellular matrix that exists in abundance throughout animals,
including the human body. This connective tissue comprises repeated units of
two main structures, namely collagens as well as axial, parallel and regular
anionic glycosaminoglycan between collagens. The collagen fibril can be modeled
by Hooke's law whereas anionic glycosaminoglycan behaves more like a
rubber-band rod and as such can be better modeled by the worm-like chain model.
While both computer simulations and continuum mechanics models have been
investigated the behavior of this connective tissue typically, authors either
assume a simple form of the molecular potential energy or entirely ignore the
microscopic structure of the connective tissue. Here, we apply basic physical
methodologies and simple applied mathematical modeling techniques to describe
the collagen pair quantitatively. We find that the growth of fibrils is
intimately related to the maximum length of the anionic glycosaminoglycan and
the relative displacement of two adjacent fibrils, which in return is closely
related to the effectiveness of anionic glycosaminoglycan in transmitting
forces between fibrils. These reveal the importance of the anionic
glycosaminoglycan in maintaining the structural shape of the connective tissue
extracellular matrix and eventually the shape modulus of human tissues. We also
find that some macroscopic properties, like the maximum molecular energy and
the breaking fraction of the collagen, are also related to the microscopic
characteristics of the anionic glycosaminoglycan.


  The hypothetical impact of aerial spraying on tsetse fly populations is
investigated. Spray cycles are scheduled at intervals two days short of the
first interlarval period and halted once the last of the female flies that
originated from pre-spray-deposited pupae have been sprayed twice. The effect
of temperature on the aerial spraying of tsetse, through its reproductive cycle
and general population dynamics, is of particular interest, given that cooler
weather is preferred for the settling of insecticidal droplets. Spray efficacy
is found to come at a price due to the greater number of cycles necessitated by
cooler weather. The extra cost is argued to be worth while. Pupae, still in the
ground at the end of spraying, are identified as the main threat to a
successful operation. They are slightly more vulnerable at the low temperature
extreme of tsetse habitat (16 $^\circ$C), when the cumulative, natural pupal
mortality is high. One can otherwise base one's expectations on the closeness
with which the time to the third last spray approaches one puparial duration. A
disparity of anything close to the length of a spray cycle advocates caution,
whereas one which comes close to vanishing should be interpreted as being
auspicious. Three such key temperatures, just below which one can anticipate an
improved outcome and just above which caution should be exercised, are 17.146
$^\circ$C, 19.278 $^\circ$C and 23.645 $^\circ$C. A refinement of the existing
formulae for the puparial duration and the first interlarval period might be
prudent in the South African context of a sympatric Glossina brevipalpis-G.
austeni, tsetse population. The resulting aerial spraying strategy would then
be formulated using a G. brevipalpis puparial duration and a G. austeni first
interlarval period.


  The wobble hypothesis does not discriminate between uracil and thymine.
Methylation could favor further stabilization of uracil in the keto form.
Thymine is present in keto form only and can pair up but with adenine. Uracil
can easily construct the enol form; that is why it forms the U-G pair.


  We propose a new distance metric for DNA sequences, which can be defined on
any evolutionary Markov model with infinitesimal generator matrix Q. That is
the new metric can be defined under existing models such as Jukes-Cantor model,
Kimura-2-parameter model, F84 model, GTR model etc. Since our metric does not
depend on the form of the generator matrix Q, it can be defined for very
general models including those with varying nucleotide substitution rates among
lineages. This makes our metric widely applicable. The simulation experiments
carried out shows that the new metric, when defined under classical models such
as the JC, F84 and Kimura-2-parameter models, performs better than these
existing metrics in recovering phylogenetic trees from sequence data. Our
simulation experiments also show that the new metric, under a model that allows
varying nucleotide substitution rates among lineages, performs equally well or
better than its other forms studied.


  We study the functional characteristics of a two-gene motif consisting of a
double positive feedback loop and an autoregulatory negative feedback loop. The
motif appears in the gene regulatory network controlling the functional
activity of pancreatic $\beta$-cells. The model exhibits bistability and
hysteresis in appropriate parameter regions. The two stable steady states
correspond to low (OFF state) and high (ON state) protein levels respectively.
Using a deterministic approach, we show that the region of bistability
increases in extent when the copy number of one of the genes is reduced from
two to one. The negative feedback loop has the effect of reducing the size of
the bistable region. Loss of a gene copy, brought about by mutations, hampers
the normal functioning of the $\beta$-cells giving rise to the genetic
disorder, maturity-onset diabetes of the young (MODY). The diabetic phenotype
makes its appearance when a sizable fraction of the $\beta$-cells is in the OFF
state. Using stochastic simulation techniques, we show that, on reduction of
the gene copy number, there is a transition from the monostable ON to the ON
state in the bistable region of the parameter space. Fluctuations in the
protein levels, arising due to the stochastic nature of gene expression, can
give rise to transitions between the ON and OFF states. We show that as the
strength of autorepression increases, the ON$\to$OFF state transitions become
less probable whereas the reverse transitions are more probable. The
implications of the results in the context of the occurrence of MODY are
pointed out..


  Models of ordering dynamics allow to understand natural systems in which an
initially disordered population homogenizes some traits via local interactions.
The simplest of these models, with wide applications ranging from evolutionary
to social dynamics, are the Voter and Moran processes, usually defined in terms
of static or randomly mixed individuals that interact with a neighbor to copy
or modify a discrete trait. Here we study the effects of diffusion in
Voter/Moran processes by proposing a generalization of ordering dynamics in a
metapopulation framework, in which individuals are endowed with mobility and
diffuse through a spatial structure represented as a graph of patches upon
which interactions take place. We show that diffusion dramatically affects the
time to reach the homogeneous state, independently of the underlying network's
topology, while the final consensus emerges through different local/global
mechanisms, depending on the mobility strength. Our results highlight the
crucial role played by mobility in ordering processes and set up a general
framework that allows to study its effect on a large class of models, with
implications in the understanding of evolutionary and social phenomena.


  In this paper we analyze the vibrational spectra of a large ensemble of
non-homologous protein structures by means of a novel tool, that we coin the
Hierarchical Network Model (HNM). Our coarse-grained scheme accounts for the
intrinsic heterogeneity of force constants displayed by protein arrangements
and also incorporates side-chain degrees of freedom.
  Our analysis shows that vibrational entropy per unit residue correlates with
the content of secondary structure. Furthermore, we assess the individual
contribution to vibrational entropy of the novel features of our scheme as
compared with the predictions of state-of-the-art network models. This analysis
highlights the importance of properly accounting for the intrinsic hierarchy in
force strengths typical of the different atomic bonds that build up and
stabilize protein scaffolds.
  Finally, we discuss possible implications of our findings in the context of
protein aggregation phenomena.


  High-throughput experiments are shedding light on the topology of large
regulatory networks and at the same time their functional states, namely the
states of activation of the nodes (for example transcript or protein levels) in
different conditions, times, environments. We now possess a certain amount of
information about these two levels of description, stored in libraries,
databases and ontologies. A current challenge is to bridge the gap between
topology and function, i.e. developing quantitative models aimed at
characterizing the expression patterns of large sets of genes. However,
approaches that work well for small networks become impossible to master at
large scales, mainly because parameters proliferate. In this review we discuss
the state of the art of large-scale functional network models, addressing the
issue of what can be considered as realistic and what the main limitations may
be. We also show some directions for future work, trying to set the goals that
future models should try to achieve. Finally, we will emphasize the possible
benefits in the understanding of biological mechanisms underlying complex
multifactorial diseases, and in the development of novel strategies for the
description and the treatment of such pathologies.


  Changes in the extent of local concavity along with changes in surface
roughness of binding sites of proteins have long been considered as useful
markers to identify functional sites of proteins. However, an algorithm that
describes the connection between the simultaneous changes of these important
parameters - eludes the students of structural biology. Here a simple yet
general mathematical scheme is proposed that attempts to achieve the same.
Instead of n-dimensional random vector description, protein surface roughness
is described here as a system of algebraic equations. Such description resulted
in the construction of a generalized index that not only describes the
shape-change-vs-surface-roughness-change process but also reduces the
estimation error in local shape characterization. Suitable algorithmic
implementation of it in context-specific macromolecular recognition can be
attempted easily. Contemporary drug discovery studies will be enormously
benefited from this work because it is the first algorithm that can estimate
the change in protein surface roughness as the local shape of the protein is
changing (and vice-versa).


  Theory of multi-dimensional representation of free energy surface of protein
folding is developed by adopting structural order parameters of multiple
regions in protein as multiple coordinates. Various scenarios of folding are
classified in terms of cooperativity within individual regions and interactions
among multiple regions and thus obtained classification is used to analyze the
folding process of several example proteins. Ribosomal protein S6, src-SH3
domain, CheY, barnase, and BBL domain are analyzed with the two-dimensional
representation by using a structure-based Hamiltonian model. Extension to the
higher dimensional representation leads to the finer description of the folding
process. Barnase, NtrC, and an ankyrin repeat protein are examined with the
three-dimensional representation. The multi-dimensional representation allows
us to directly address questions on folding pathways, intermediates, and
transition states.


  The importance of mutator clones in the adaptive evolution of asexual
populations is not fully understood. Here we address this problem by using an
ab initio microscopic model of living cells, whose fitness is derived directly
from their genomes using a biophysically realistic model of protein folding and
interactions in the cytoplasm. The model organisms contain replication
controlling genes (DCGs) and genes modeling the mismatch repair (MMR)
complexes. We find that adaptation occurs through the transient fixation of a
mutator phenotype, regardless of particular perturbations in the fitness
landscape. The microscopic pathway of adaptation follows a well-defined set of
events: stochastic switching to the mutator phenotype first, then mutation in
the MMR complex that hitchhikes with a beneficial mutation in the DCGs, and
finally a compensating mutation in the MMR complex returning the population to
a non-mutator phenotype. Similarity of these results to reported adaptation
events points out to robust universal physical principles of evolutionary
adaptation.


  Synthetic biology is the engineering of cellular networks. It combines
principles of engineering and the knowledge of biological networks to program
the behavior of cells. Computational modeling techniques in conjunction with
molecular biology techniques have been successful in constructing biological
devices such as switches, oscillators, and gates. The ambition of synthetic
biology is to construct complex systems from such fundamental devices, much in
the same way electronic circuits are built from basic parts. As this ambition
becomes a reality, engineering concepts such as interchangeable parts and
encapsulation will find their way into biology. We realize that there is a need
for computational tools that would support such engineering concepts in
biology. As a solution, we have developed the software Athena that allows
biological models to be constructed as modules. Modules can be connected to one
another without altering the modules themselves. In addition, Athena houses
various tools useful for designing synthetic networks including tools to
perform simulations, automatically derive transcription rate expressions, and
view and edit synthetic DNA sequences. New tools can be incorporated into
Athena without modifying existing program via a plugin interface, IronPython
scripts, Systems Biology Workbench interfacing and the R statistical language.
The program is currently for Windows operating systems, and the source code for
Athena is made freely available through CodePlex, www.codeplex.com/athena.


  The structure of DNA Binding Proteins enables a strong interaction with their
specific target site on DNA. However, recent single molecule experiment
reported that proteins can diffuse on DNA. This suggests that the interactions
between proteins and DNA play a role during the target search even far from the
specific site. It is unclear how these non-specific interactions optimize the
search process, and how the protein structure comes into play. Each nucleotide
being negatively charged, one may think that the positive surface of DNA-BPs
should electrostatically collapse onto DNA. Here we show by means of Monte
Carlo simulations and analytical calculations that a counter-intuitive
repulsion between the two oppositely charged macromolecules exists at a
nanometer range. We also show that this repulsion is due to a local increase of
the osmotic pressure exerted by the ions which are trapped at the interface.
For the concave shape of DNA-BPs, and for realistic protein charge densities,
we find that the repulsion pushes the protein in a free energy minimum at a
distance from DNA. As a consequence, a favourable path exists along which
proteins can slide without interacting with the DNA bases. When a protein
encounters its target, the osmotic barrier is completely counter-balanced by
the H-bond interaction, thus enabling the sequence recognition.


  The hair cell's mechanoreceptive organelle, the hair bundle, is highly
sensitive because its transduction channels open over a very narrow range of
displacements. The synchronous gating of transduction channels also underlies
the active hair-bundle motility that amplifies and tunes responsiveness. The
extent to which the gating of independent transduction channels is coordinated
depends on how tightly individual stereocilia are constrained to move as a
unit. Using dual-beam interferometry in the bullfrog's sacculus, we found that
thermal movements of stereocilia located as far apart as a bundle's opposite
edges display high coherence and negligible phase lag. Because the mechanical
degrees of freedom of stereocilia are strongly constrained, a force applied
anywhere in the hair bundle deflects the structure as a unit. This feature
assures the concerted gating of transduction channels that maximizes the
sensitivity of mechanoelectrical transduction and enhances the hair bundle's
capacity to amplify its inputs.


  A bitset is a set that encodes for a binary number. Bitsets are at the basis
of a beautiful theory of recombination with n-loci and here we begin from
scratch and advance to include the derivation of the fundamental results about
the evolution of gamete frequencies and of disequilibrium measures with and
without migration. All techniques have been illustrated and we have invested
moreover a great effort to make the mathematics of this work accessible even
for students in their first year at the university.


  We have developed URDME, a general software for simulation of stochastic
reaction-diffusion processes on unstructured meshes. This allows for a more
flexible handling of complicated geometries and curved boundaries compared to
simulations on structured, cartesian meshes. The underlying algorithm is the
next subvolume method, extended to unstructured meshes by obtaining jump
coefficients from a finite element formulation of the corresponding macroscopic
equation.
  This manual describes version 1.4 of the software. Refer to www.urdme.org for
the latest updates.


  Regulatory interactions between genes show a large amount of cross-species
variability, even when the underlying functions are conserved: There are many
ways to achieve the same function. Here we investigate the ability of
regulatory networks to reproduce given expression levels within a simple model
of gene regulation. We find an exponentially large space of regulatory networks
compatible with a given set of expression levels, giving rise to an extensive
entropy of networks. Typical realisations of regulatory networks are found to
share a bias towards symmetric interactions, in line with empirical evidence.


  We combine the principle of gauge invariance with extrinsic string geometry
to develop a lattice model that can be employed to theoretically describe
properties of chiral, unbranched homopolymers. We find that in its low
temperature phase the model is in the same universality class with proteins
that are deposited in the Protein Data Bank, in the sense of the compactness
index. We apply the model to analyze various statistical aspects of folded
proteins. Curiously we find that it can produce results that are a very good
good match to the data in the Protein Data Bank.


  Proteins are the "work horses" in biological systems. In almost all functions
specific proteins are involved. They control molecular transport processes,
stabilize the cell structure, enzymatically catalyze chemical reactions; others
act as molecular motors in the complex machinery of molecular synthetization
processes. Due to their significance, misfolds and malfunctions of proteins
typically entail disastrous diseases, such as Alzheimer's disease and bovine
spongiform encephalopathy (BSE). Therefore, the understanding of the trinity of
amino acid composition, geometric structure, and biological function is one of
the most essential challenges for the natural sciences. Here, we glance at
conformational transitions accompanying the structure formation in protein
folding processes.


  A rooted acyclic digraph N with labelled leaves displays a tree T when there
exists a way to select a unique parent of each hybrid vertex resulting in the
tree T. Let Tr(N) denote the set of all trees displayed by the network N. In
general, there may be many other networks M such that Tr(M) = Tr(N). A network
is regular if it is isomorphic with its cover digraph. This paper shows that if
N is regular, there is a procedure to reconstruct N given Tr(N). Hence if N and
M are regular networks and Tr(N) = Tr(M), it follows that N = M, proving that a
regular network is uniquely determined by its displayed trees.


  Molecular motors are single macromolecules that generate forces at the
piconewton range and nanometer scale. They convert chemical energy into
mechanical work by moving along filamentous structures. In this paper, we study
the velocity of two-head molecular motors in the framework of a mechanochemical
network theory. The network model, a generalization of the recently work of
Liepelt and Lipowsky (PRL 98, 258102 (2007)), is based on the discrete
mechanochemical states of a molecular motor with multiple cycles. By
generalizing the mathematical method developed by Fisher and Kolomeisky for
single cycle motor (PNAS(2001) 98(14) P7748-7753), we are able to obtain an
explicit formula for the velocity of a molecular motor.


  The biogenesis of lipid droplets (LD) in the yeast Saccharomyces cerevisiae
was theoretically investigated on basis of a biophysical model. In accordance
with the prevailing model of LD formation, we assumed that neutral lipids
oil-out between the membrane leaflets of the endoplasmic reticulum (ER),
resulting in LD that bud-off when a critical size is reached.
  Mathematically, LD were modeled as spherical protuberances in an otherwise
planar ER membrane. We estimated the local phospholipid composition, and
calculated the change in elastic free energy of the membrane caused by nascent
LD. Based on this model calculation, we found a gradual demixing of lipids in
the membrane leaflet that goes along with an increase in surface curvature at
the site of LD formation. During demixing, the phospholipid monolayer was able
to gain energy during LD growth, which suggested that the formation of curved
interfaces was supported by or even driven by lipid demixing. In addition, we
show that demixing is thermodynamically necessary as LD cannot bud-off
otherwise.
  In the case of Saccharomyces cerevisiae our model predicts a LD bud-off
diameter of about 13 nm. This diameter is far below the experimentally
determined size of typical yeast LD. Thus, we concluded that if the standard
model of LD formation is valid, LD biogenesis is a two step process. Small LD
are produced from the ER, which subsequently ripe within the cytosol through a
series of fusions.


  We investigate the thermodynamics and kinetics of DNA hairpins that
fold/unfold under the action of applied mechanical force. We introduce the
concept of the molecular free energy landscape and derive simplified
expressions for the force dependent Kramers-Bell rates. To test the theory we
have designed a specific DNA hairpin sequence that shows two-state cooperative
folding under mechanical tension and carried out pulling experiments using
optical tweezers. We show how we can determine the parameters that characterize
the molecular free energy landscape of such sequence from rupture force kinetic
studies. Finally we combine such kinetic studies with experimental
investigations of the Crooks fluctuation relation to derive the free energy of
formation of the hairpin at zero force.


  We investigate irreversibility and dissipation in single molecules that
cooperatively fold/unfold in a two state manner under the action of mechanical
force. We apply path thermodynamics to derive analytical expressions for the
average dissipated work and the average hopping number in two state systems. It
is shown how these quantities only depend on two parameters that characterize
the folding/unfolding kinetics of the molecule: the fragility and the
coexistence hopping rate. The latter has to be rescaled to take into account
the appropriate experimental setup. Finally we carry out pulling experiments
with optical tweezers in a specifically designed DNA hairpin that shows
two-state cooperative folding. We then use these experimental results to
validate our theoretical predictions.


  An $k$-noncrossing RNA structure can be identified with an $k$-noncrossing
diagram over $[n]$, which in turn corresponds to a vacillating tableaux having
at most $(k-1)$ rows. In this paper we derive the limit distribution of
irreducible substructures via studying their corresponding vacillating
tableaux. Our main result proves, that the limit distribution of the numbers of
irreducible substructures in $k$-noncrossing, $\sigma$-canonical RNA structures
is determined by the density function of a $\Gamma(-\ln\tau_k,2)$-distribution
for some $\tau_k<1$.


  In this paper we study irreducibility in RNA structures. By RNA structure we
mean RNA secondary as well as RNA pseudoknot structures. In our analysis we
shall contrast random and minimum free energy (mfe) configurations. We compute
various distributions: of the numbers of irreducible substructures, their
locations and sizes, parameterized in terms of the maximal number of mutually
crossing arcs, $k-1$, and the minimal size of stacks $\sigma$. In particular,
we analyze the size of the largest irreducible substructure for random and mfe
structures, which is the key factor for the folding time of mfe configurations.


  The general tendency for species number (S) to increase with sampled area (A)
constitutes one of the most robust empirical laws of ecology, quantified by
species-area relationships (SAR). In many ecosystems, SAR curves display a
power-law dependence, $S\propto A^z$. The exponent $z$ is always less than one
but shows significant variation in different ecosystems. We study the multitype
voter model as one of the simplest models able to reproduce SAR similar to
those observed in real ecosystems in terms of basic ecological processes such
as birth, dispersal and speciation. Within the model, the species-area exponent
$z$ depends on the dimensionless speciation rate $\nu$, even though the
detailed dependence is still matter of controversy. We present extensive
numerical simulations in a broad range of speciation rates from $\nu =10^{-3}$
down to $\nu = 10^{-11}$, where the model reproduces values of the exponent
observed in nature. In particular, we show that the inverse of the species-area
exponent linearly depends on the logarithm of $\nu$. Further, we compare the
model outcomes with field data collected from previous studies, for which we
separate the effect of the speciation rate from that of the different species
lifespans. We find a good linear relationship between inverse exponents and
logarithm of species lifespans. However, the slope sets bounds on the
speciation rates that can hardly be justified on evolutionary basis, suggesting
that additional effects should be taken into account to consistently interpret
the observed exponents.


  In this paper I describe a cellular automaton model of a multi-species
ecosystem, suitable for the study of emergent properties of macroevolution.
Unlike majority of ecological models, the number of coexisting species is not
fixed. Starting from one common ancestor they appear by "mutations" of existent
species, and then survive or extinct depending on the balance of local
ecological interactions. Monte-Carlo numerical simulations show that this model
is able to qualitatively reproduce phenomena that have been observed in other
models and in nature.


  We investigate the selective forces that promote the emergence of modularity
in nature. We demonstrate the spontaneous emergence of modularity in a
population of individuals that evolve in a changing environment. We show that
the level of modularity correlates with the rapidity and severity of
environmental change. The modularity arises as a synergistic response to the
noise in the environment in the presence of horizontal gene transfer. We
suggest that the hierarchical structure observed in the natural world may be a
broken symmetry state, which generically results from evolution in a changing
environment. To support our results, we analyze experimental protein
interaction data and show that protein interaction networks became increasingly
modular as evolution proceeded over the last four billion years. We also
discuss a method to determine the divergence time of a protein.


  We investigate the dynamical properties of the transcriptional regulation of
gene expression in the yeast Saccharomyces Cerevisiae within the framework of a
synchronously and deterministically updated Boolean network model. By means of
a dynamically determinant subnetwork, we explore the robustness of
transcriptional regulation as a function of the type of Boolean functions used
in the model that mimic the influence of regulating agents on the transcription
level of a gene. We compare the results obtained for the actual yeast network
with those from two different model networks, one with similar in-degree
distribution as the yeast and random otherwise, and another due to Balcan et
al., where the global topology of the yeast network is reproduced faithfully.
We, surprisingly, find that the first set of model networks better reproduce
the results found with the actual yeast network, even though the Balcan et al.
model networks are structurally more similar to that of yeast.


  Fontanari et al introduced [Phys. Rev. Lett. 91, 218101 (2003)] a model for
studying the Muller's ratchet phenomenon in growing asexual populations. They
studied two situations, either including or not a death probability for each
newborn, but were able to find analytical (recursive) expressions only in the
no-decay case. In this paper a branching process formalism is used to find
recorrence equations that generalize the analytical results of the original
paper besides confirming the interesting effects their simulations revealed.


  We present herein an extension of an algebraic statistical method for
inferring biochemical reaction networks from experimental data, proposed
recently in [3]. This extension allows us to analyze reaction networks that are
not necessarily full-dimensional, i.e., the dimension of their stoichiometric
space is smaller than the number of species. Specifically, we propose to
augment the original algebraic-statistical algorithm for network inference with
a preprocessing step that identifies the subspace spanned by the correct
reaction vectors, within the space spanned by the species. This dimension
reduction step is based on principal component analysis of the input data and
its relationship with various subspaces generated by sets of candidate reaction
vectors. Simulated examples are provided to illustrate the main ideas involved
in implementing this method, and to asses its performance.


  Since the first experimental evidences of active conductances in dendrites,
most neurons have been shown to exhibit dendritic excitability through the
expression of a variety of voltage-gated ion channels. However, despite
experimental and theoretical efforts undertaken in the last decades, the role
of this excitability for some kind of dendritic computation has remained
elusive. Here we show that, owing to very general properties of excitable
media, the average output of a model of active dendritic trees is a highly
non-linear function of their afferent rate, attaining extremely large dynamic
ranges (above 50 dB). Moreover, the model yields double-sigmoid response
functions as experimentally observed in retinal ganglion cells. We claim that
enhancement of dynamic range is the primary functional role of active dendritic
conductances. We predict that neurons with larger dendritic trees should have
larger dynamic range and that blocking of active conductances should lead to a
decrease of dynamic range.


  In a previous work, we gave a metric on the class of semibinary tree-sibling
time consistent phylogenetic networks that is computable in polynomial time; in
particular, the problem of deciding if two networks of this kind are isomorphic
is in P. In this paper, we show that if we remove the semibinarity condition
above, then the problem becomes much harder. More precisely, we proof that the
isomorphism problem for generic tree-sibling time consistent phylogenetic
networks is polynomially equivalent to the graph isomorphism problem. Since the
latter is believed to be neither in P nor NP-complete, the chances are that it
is impossible to define a metric on the class of all tree-sibling time
consistent phylogenetic networks that can be computed in polynomial time.


  We study the evolution of cooperation in the prisoner's dilemma and the
snowdrift game on scale-free networks that are subjected to intentional and
random removal of vertices. We show that, irrespective of the game type,
cooperation on scale-free networks is extremely robust against random deletion
of vertices, but declines fast if vertices with the maximal degree are
targeted. In particular, attack tolerance is lowest if the temptation to defect
is largest, whereby a small fraction of removed vertices suffices to decimate
cooperators. The decline of cooperation can be directly linked to the decrease
of heterogeneity of scale-free networks that sets in due to the removal of high
degree vertices. We conclude that the evolution of cooperation is characterized
by similar attack and error tolerance as was previously reported for
information readiness and spread of viruses on scale-free networks.


  We discuss separate and integrated approaches to building scalable solar
power plants and wildlife sanctuaries. Both solar power plants and wildlife
sanctuaries need a lot of land. We quantify some of the requirements using
various estimates of the rate of solar power production as well as the rate of
adding wildlife to a sanctuary over the time range 2010-2050. We use population
dynamics equations to study the evolution of solar energy and tiger populations
up to and beyond 2050.


  This paper considers the relative frequencies of distinct types of
individuals in multitype branching processes. We prove that the frequencies are
asymptotically multivariate normal when the initial number of ancestors is
large and the time of observation is fixed. The result is valid for any
branching process with a finite number of types; the only assumption required
is that of independent individual evolutions. The problem under consideration
is motivated by applications in the area of cell biology. Specifically, the
reported limiting results are of advantage in cell kinetics studies where the
relative frequencies but not the absolute cell counts are accessible to
measurement. Relevant statistical applications are discussed in the context of
asymptotic maximum likelihood inference for multitype branching processes.


  New approach to design a dynamic model of genes with multiple autonomous
regulatory modules by evolution in silico is proposed. The approach is based on
Genetic Algorithms, enforced by new crossover operators, especially worked out
for these purposes. The approach exploits the subbasin-portal architecture of
the fitness functions suitable for this kind of evolutionary modeling. The
effectiveness of the approach is demonstrated on a series of benchmark tests.


  In this paper we denonstrate that the use of the system of 2-adic numbers
provides a new insight to some problems of genetics, in particular, generacy of
the genetic code and the structure of the PAM matrix in bioinformatics. The
2-adic distance is an ultrametric and applications of ultrametrics in
bioinformatics are not surprising. However, by using the 2-adic numbers we
match ultrametric with a number theoretic structure. In this way we find new
applications of an ultrametric which differ from known up to now in
bioinformatics.
  We obtain the following results. We show that the PAM matrix A allows the
expansion into the sum of the two matrices A=A^{(2)}+A^{(\infty)}, where the
matrix A^{(2)} is 2-adically regular (i.e. matrix elements of this matrix are
close to locally constant with respect to the discussed earlier by the authors
2-adic parametrization of the genetic code), and the matrix A^{(\infty)} is
sparse. We discuss the structure of the matrix A^{(\infty)} in relation to the
side chain properties of the corresponding amino acids.


  Synthesis of the biological "energy currency molecule" adenosine triphosphate
ATP is accomplished by FoF1-ATP synthase. In the plasma membrane of Escherichia
coli, proton-driven rotation of a ring of 10 c subunits in the Fo motor powers
catalysis in the F1 motor. While F1 uses 120 degree stepping, Fo models predict
a step-by-step rotation of c subunits 36 degree at a time, which is here
demonstrated by single-molecule fluorescence resonance energy transfer.


  The dynamics of two groups of molecular motors pulling in opposite directions
on a rigid filament is studied theoretically. To this end we first consider the
behavior of one set of motors pulling in a single direction against an external
force using a new mean-field approach. Based on these results we analyze a
similar setup with two sets of motors pulling in opposite directions in a
tug-of-war in the presence of an external force. In both cases we find that the
interplay of fluid friction and protein friction leads to a complex phase
diagram where the force-velocity relations can exhibit regions of bistability
and spontaneous symmetry breaking. Finally, motivated by recent work, we turn
to the case of motility assay experiments where motors bound to a surface push
on a bundle of filaments. We find that, depending on the absence or the
presence of a bistability in the force-velocity curve at zero force, the bundle
exhibits anomalous or biased diffusion on long-time and large-length scales.


  We consider a charged Brownian particle in an asymmetric bistable
electrostatic potential biased by an externally applied or induced time
periodic electric field. While the amplitude of the applied field is
independent of frequency, that of the one induced by a magnetic field is.
Borrowing from protein channel terminology, we define the open probability as
the relative time the Brownian particle spends on a prescribed side of the
potential barrier. We show that while there is no peak in the open probability
as the frequency of the applied field and the bias (depolarization) of the
potential are varied, there is a narrow range of low frequencies of the induced
field and a narrow range of the low bias of the potential where the open
probability peaks. This manifestation of stochastic resonance is consistent
with experimental results on the voltage gated Iks and KCNQ1 potassium channels
of biological membranes and on cardiac myocytes.


  CLeFAPS, a fast and flexible pairwise structural alignment algorithm based on
a rigid-body framework, namely CLePAPS, is proposed. Instead of allowing twists
(or bends), the flexible in CLeFAPS means: (a) flexibilization of the
algorithm's parameters through self-adapting with the input structures' size,
(b) flexibilization of adding the aligned fragment pairs (AFPs) into an
one-to-multi correspondence set instead of checking their position conflict,
(c) flexible fragment may be found through an elongation procedure rooted in a
vector-based score instead of a distance-based score. We perform a comparison
between CLeFAPS and other popular algorithms including rigid-body and flexible
on a closely-related protein benchmark (HOMSTRAD) and a distantly-related
protein benchmark (SABmark) while the latter is also for the discrimination
test, the result shows that CLeFAPS is competitive with or even outperforms
other algorithms while the running time is only 1/150 to 1/50 of them.


  Effective therapy of complex diseases requires control of highly non-linear
complex networks that remain incompletely characterized. In particular, drug
intervention can be seen as control of signaling in cellular networks.
Identification of control parameters presents an extreme challenge due to the
combinatorial explosion of control possibilities in combination therapy and to
the incomplete knowledge of the systems biology of cells. In this review paper
we describe the main current and proposed approaches to the design of
combinatorial therapies, including the empirical methods used now by clinicians
and alternative approaches suggested recently by several authors. New
approaches for designing combinations arising from systems biology are
described. We discuss in special detail the design of algorithms that identify
optimal control parameters in cellular networks based on a quantitative
characterization of control landscapes, maximizing utilization of incomplete
knowledge of the state and structure of intracellular networks. The use of new
technology for high-throughput measurements is key to these new approaches to
combination therapy and essential for the characterization of control
landscapes and implementation of the algorithms. Combinatorial optimization in
medical therapy is also compared with the combinatorial optimization of
engineering and materials science and similarities and differences are
delineated.


  We present two algorithms for computing the geodesic distance between
phylogenetic trees in tree space, as introduced by Billera, Holmes, and
Vogtmann (2001). We show that the possible combinatorial types of shortest
paths between two trees can be compactly represented by a partially ordered
set. We calculate the shortest distance along each candidate path by converting
the problem into one of finding the shortest path through a certain region of
Euclidean space. In particular, we show there is a linear time algorithm for
finding the shortest path between a point in the all positive orthant and a
point in the all negative orthant of R^k contained in the subspace of R^k
consisting of all orthants with the first i coordinates non-positive and the
remaining coordinates non-negative for 0 <= i <= k.


  Unshared consensus decision-making processes, in which one or a small number
of individuals make the decision for the rest of a group, are rarely
documented. However, this mechanism can be beneficial for all group members
when one individual has greater knowledge about the benefits of the decision
than other group members. Such decisions are reached during certain activity
shifts within the population of bottlenose dolphins residing in Doubtful Sound,
New Zealand. Behavioral signals are performed by one individual and seem to
precipitate shifts in the behavior of the entire group: side flops are
performed by males and initiate traveling bouts while upside-down lobtails are
performed by females and terminate traveling bouts. However, these signals are
not observed at all activity shifts. We find that while side flops were
performed by males that have greater knowledge than other male group members,
this was not the case for females performing upside-down lobtails. The reason
for this could have been that a generally high knowledge about the optimal
timing of travel terminations rendered it less important which individual
female made the decision.


  In this paper, we investigate the epidemic spreading for SIR model in
weighted scale-free networks with nonlinear infectivity, where the transmission
rate in our analytical model is weighted. Concretely, we introduce the
infectivity exponent $\alpha$ and the weight exponent $\beta$ into the
analytical SIR model, then examine the combination effects of $\alpha$ and
$\beta$ on the epidemic threshold and phase transition. We show that one can
adjust the values of $\alpha$ and $\beta$ to rebuild the epidemic threshold to
a finite value, and it is observed that the steady epidemic prevalence $R$
grows in an exponential form in the early stage, then follows hierarchical
dynamics. Furthermore, we find $\alpha$ is more sensitive than $\beta$ in the
transformation of the epidemic threshold and epidemic prevalence, which might
deliver some useful information or new insights in the epidemic spreading and
the correlative immunization schemes.


  The paper aims to discuss statistical properties of the multi-agent based
model of competitive growth. Each of the agents is described by growth (or
decay) rule of its virtual "mass" with the rate affected by the interaction
with other agents. The interaction depends on the strategy vector and mutual
distance between agents and both are subjected to the agent's individual
optimization process. Steady-state simulations yield phase diagrams with the
high and low competition phases (HCP and LCP, respectively) separated by
critical point. Particular focus has been made on the indicators of the
power-law behavior of the mass distributions with respect to the critical
regime. In this regime the study has revealed remarkable anomaly in the
optimization efficiency.


  Fluorescence recovery after photobleaching (FRAP) measurements offer an
important tool for analyzing diffusion and binding processes. Confocal scanning
laser microscopes that are used in FRAP experiments bleach regions with a
radially Gaussian distributed profile. Previous attempts to derive analytical
expressions in the case of processes governed by fast diffusion have overlooked
the characteristics of the instruments used to perform FRAP measurements and
therefore led to approximating solutions. In the present paper, bleaching laser
beam characteristics are incorporated into an improved model to provide a more
rigorous and accurate method. The proposed model simulates binding inside
bounded regions, and it leads to FRAP curves that depend on the on and off
rates that can be employed to determine the rate constants. It can be used in
conjunction with experimental data acquired with confocal scanning laser
microscopes to investigate the biophysical properties of proteins in living
cells. The model aims to improve the accuracy when determining rate constants
by taking into account amore realistic scenario of the light-matter
interaction.


  Here we study how mutations which change physical properties of cell proteins
(stability) impact population survival and growth. In our model the genotype is
presented as a set of N numbers, folding free energies of cells N proteins.
Mutations occur upon replications so that stabilities of some proteins in
daughter cells differ from those in parent cell by random amounts drawn from
experimental distribution of mutational effects on protein stability. The
genotype-phenotype relationship posits that unstable proteins confer lethal
phenotype to a cell and in addition the cells fitness (duplication rate) is
proportional to the concentration of its folded proteins. Simulations reveal
that lethal mutagenesis occurs at mutation rates close to 7 mutations per
genome per replications for RNA viruses and about half of that for DNA based
organisms, in accord with earlier predictions from analytical theory and
experiment. This number appears somewhat dependent on the number of genes in
the organisms and natural death rate. Further, our model reproduces the
distribution of stabilities of natural proteins in excellent agreement with
experiment. Our model predicts that species with high mutation rates, tend to
have less stable proteins compared to species with low mutation rate.


  Over the past decade network theory has been applied successfully to the
study of a variety of complex adaptive systems. However, the application of
these techniques to non-human social networks has several shortfalls. Firstly,
in most cases the strength of associations between individuals is disregarded.
Secondly, present techniques assume that observed interactions are invariant
values and not statistical samples taken from a population. These two
simplifications have weakened the value of these techniques when applied to the
study of animal social systems. Here we introduce a set of behaviorally
meaningful weighted network statistics that can be readily applied to matrices
of association indices between pairs of individual animals. We also introduce
bootstrapping techniques that estimate the effects of sampling uncertainty on
the network statistics and structure. Finally, we discuss the use of
randomisation tests to detect the departure of observed network statistics from
expected values under null hypotheses of random association given the sampling
structure of the data. We use two case studies to show that these techniques
provide invaluable insight in the dynamics of interactions within social units
and in the community structure of societies.


  Because of the double-helical structure of DNA, in which two strands of
complementary nucleotides intertwine around each other, a covalently closed DNA
molecule with no interruptions in either strand can be viewed as two
interlocked single-stranded rings. Two closed space curves have long been known
by mathematicians to exhibit a property called the linking number, a
topologically invariant integer, expressible as the sum of two other
quantities, the twist of one of the curves about the other, and the writhing
number, or writhe, a measure of the chiral distortion from planarity of one of
the two closed curves. We here derive expressions for the twist of supercoiled
DNA and the writhe of a closed molecule consistent with the modern view of DNA
as a sequence of base-pair steps. Structural biologists commonly characterize
the spatial disposition of each step in terms of six rigid-body parameters, one
of which, coincidentally, is also called the twist. Of interest is the
difference in the mathematical properties between this step-parameter twist and
the twist of supercoiling associated with a given base-pair step. For example,
it turns out that the latter twist, unlike the former, is sensitive to certain
translational shearing distortions of the molecule that are chiral in nature.
Thus, by comparing the values for the two twists for each step of a
high-resolution structure of a protein-DNA complex, the nucleosome considered
here, for example, we may be able to determine how the binding of various
proteins contributes to chiral structural changes of the DNA.


  We present a non-neutral stochastic model for the dynamics taking place in a
meta-community ecosystems in presence of migration. The model provides a
framework for describing the emergence of multiple ecological scenarios and
behaves in two extreme limits either as the unified neutral theory of
biodiversity or as the Bak-Sneppen model. Interestingly, the model shows a
condensation phase transition where one species becomes the dominant one, the
diversity in the ecosystems is strongly reduced and the ecosystem is
non-stationary. This phase transition extend the principle of competitive
exclusion to open ecosystems and might be relevant for the study of the impact
of invasive species in native ecologies.


  In this paper we consider piecewise affine differential equations modeling
gene networks. We work with arbitrary decay rates, and under a local hypothesis
expressed as an alignment condition of successive focal points. The interaction
graph of the system may be rather complex (multiple intricate loops of any
sign, multiple thresholds...). Our main result is an alternative theorem
showing that, if a sequence of region is periodically visited by trajectories,
then under our hypotheses, there exists either a unique stable periodic
solution, or the origin attracts all trajectories in this sequence of regions.
This result extends greatly our previous work on a single negative feedback
loop. We give several examples and simulations illustrating different cases.


  To represent the evolution of nucleic acid and protein sequence, we express
the parallel and Eigen models for molecular evolution in terms of a functional
integral representation with an $h$-letter alphabet, lifting the two-state,
purine/pyrimidine assumption often made in quasi-species theory. For arbitrary
$h$ and a general mutation scheme, we obtain the solution of this model in
terms of a maximum principle. Euler's theorem for homogeneous functions is used
to derive this `thermodynamic' formulation of evolution. The general result for
the parallel model reduces to known results for the purine/pyrimidine $h=2$
alphabet and the nucleic acid $h=4$ alphabet for the Kimura 3 ST mutation
scheme. Examples are presented for the $h=4$ and $h=20$ cases. We derive the
maximum principle for the Eigen model for general $h$. The general result for
the Eigen model reduces to a known result for $h=2$. Examples are presented for
the nucleic acid $h=4$ and the amino acid $h=20$ alphabet. An error catastrophe
phase transition occurs in these models, and the order of the phase transition
changes from second to first order for smooth fitness functions when the
alphabet size is increased beyond two letters to the generic case. As examples,
we analyze the general analytic solution for sharp peak, linear, quadratic, and
quartic fitness functions.


  Large protein assemblies, such as virus capsids, may be coarse-grained as a
set of rigid domains linked by generalized (rotational and stretching) harmonic
springs. We present a method to obtain the elastic parameters and overdamped
dynamics for these springs from all-atom molecular dynamics simulations of one
pair of domains at a time. The computed relaxation times of this pair give a
consistency check for the simulation, and (using a fluctuation-dissipation
relationship) we find the corrective force needed to null systematic drifts. As
a first application we predict the stiffness of an HIV capsid layer and the
relaxation time for its breathing mode.


  In the course of Darwinian evolution of a population, punctualism is an
important phenomenon whereby long periods of genetic stasis alternate with
short periods of rapid evolutionary change. This paper provides a mathematical
interpretation of punctualism as a sequence of change of basin of attraction
for a diffusion model of the theory of adaptive dynamics. Such results rely on
large deviation estimates for the diffusion process. The main difficulty lies
in the fact that this diffusion process has degenerate and non-Lipschitz
diffusion part at isolated points of the space and non-continuous drift part at
the same points. Nevertheless, we are able to prove strong existence and the
strong Markov property for these diffusions, and to give conditions under which
pathwise uniqueness holds. Next, we prove a large deviation principle involving
a rate function which has not the standard form of diffusions with small noise,
due to the specific singularities of the model. Finally, this result is used to
obtain asymptotic estimates for the time needed to exit an attracting domain,
and to identify the points where this exit is more likely to occur.


  Synthesis of protein molecules in a cell are carried out by ribosomes. A
ribosome can be regarded as a molecular motor which utilizes the input chemical
energy to move on a messenger RNA (mRNA) track that also serves as a template
for the polymerization of the corresponding protein. The forward movement,
however, is characterized by an alternating sequence of translocation and
pause. Using a quantitative model, which captures the mechanochemical cycle of
an individual ribosome, we derive an {\it exact} analytical expression for the
distribution of its dwell times at the successive positions on the mRNA track.
Inverse of the average dwell time satisfies a ``Michaelis-Menten-like''
equation and is consistent with the general formula for the average velocity of
a molecular motor with an unbranched mechano-chemical cycle. Extending this
formula appropriately, we also derive the exact force-velocity relation for a
ribosome. Often many ribosomes simultaneously move on the same mRNA track,
while each synthesizes a copy of the same protein. We extend the model of a
single ribosome by incorporating steric exclusion of different individuals on
the same track. We draw the phase diagram of this model of ribosome traffic in
3-dimensional spaces spanned by experimentally controllable parameters. We
suggest new experimental tests of our theoretical predictions.