  Blood leukocytes can exhibit extensive morphological changes during their
passage through small capillary vessels. The human monocytic THP-1 cell line
was used to explore the metabolic dependence of these shape changes. Cells were
aspirated into micropipettes for determination of the rate of protrusion
formation. They were then released and the kinetics of morphological recovery
was studied. Results were consistent with Evans' model (Blood, 64 : 1028, 1984)
of a viscous liquid droplet surrounded by a tensile membrane. The estimated
values of cytoplasmic viscosity and membrane tension were 162 Pa.s and 0.0142
millinewton/m respectively. The influence of metabolic inhibitors on cell
mechanical behaviour was then studied : results strongly suggested that
deformation involved two sequential phases. The cell elongation rate measured
during the first 30 seconds following the onset of aspiration was unaffected by
azide, an inhibitor of energy production, and it was about doubled by
cytochalasin D, a microfilament inhibitor, and colchicine, a microtubule
inhibitor. However, during the following two minutes, deformation was almost
abolished in cells treated with azide and cytochalasin D, whereas the
protrusion of control cells exhibited about threefold length increase. It is
concluded that, although cells seemed to deform as passive objects, active
metabolic processes were required to allow extensive morphological changes
triggered by external forces.


  Cells are known to utilize biochemical noise to probabilistically switch
between distinct gene expression states. We demonstrate that such noise-driven
switching is dominated by tails of probability distributions and is therefore
exponentially sensitive to changes in physiological parameters such as
transcription and translation rates. However, provided mRNA lifetimes are
short, switching can still be accurately simulated using protein-only models of
gene expression. Exponential sensitivity limits the robustness of noise-driven
switching, suggesting cells may use other mechanisms in order to switch
reliably.


  Learning latent expression themes that best express complex patterns in a
sample is a central problem in data mining and scientific research. For
example, in computational biology we seek a set of salient gene expression
themes that explain a biological process, extracting them from a large pool of
gene expression profiles. In this paper, we introduce probabilistic models to
learn such latent themes in an unsupervised fashion. Our models capture
contagion, i.e., dependence among multiple occurrences of the same feature,
using a hierarchical Bayesian scheme. Contagion is a convenient analytical
formalism to characterize semantic themes underlying observed feature patterns,
such as biological context. We present model variants tailored to different
properties of biological data, and we outline a general variational inference
scheme for approximate posterior inference. We validate our methods on both
simulated data and realistic high-throughput gene expression profiles via SAGE.
Our results show improved predictions of gene functions over existing methods
based on stronger independence assumptions, and demonstrate feasibility of a
promising hierarchical Bayesian formalism for soft clustering and latent
aspects analysis.


  While designing oligonucleotide-based microarrays, cross-hybridization
between surface-bound oligos and non-intended labeled targets is probably the
most difficult parameter to predict. Although literature describes
rules-of-thumb concerning oligo length, overall similarity, and continuous
stretches, the final behavior is difficult to predict. The aim of this study
was to investigate the effect of well-defined mismatches on hybridization
specificity using CodeLink Activated Slides, and to study quantitatively the
relation between hybridization intensity and Gibbs free energy (Delta G),
taking the mismatches into account. Our data clearly showed a correlation
between the hybridization intensity and Delta G of the oligos over three orders
of magnitude for the hybridization intensity, which could be described by the
Langmuir model. As Delta G was calculated according to the nearest-neighbor
model, using values related to DNA hybridizations in solution, this study
clearly shows that target-probe hybridizations on microarrays with a
three-dimensional coating are in quantitative agreement with the corresponding
reaction in solution. These results can be interesting for some practical
applications. The correlation between intensity and Delta G can be used in
quality control of microarray hybridizations by designing probes and
corresponding RNA spikes with a range of Delta G values. Furthermore, this
correlation might be of use to fine-tune oligonucleotide design algorithms in a
way to improve the prediction of the influence of mismatching targets on
microarray hybridizations.


  During bouts of evolutionary diversification, such as adaptive radiations,
the emerging species cluster around different locations in phenotype space, How
such multimodal patterns in phenotype space can emerge from a single ancestral
species is a fundamental question in biology. Frequency-dependent competition
is one potential mechanism for such pattern formation, as has previously been
shown in models based on the theory of adaptive dynamics. Here we demonstrate
that also in models similar to those used in quantitative genetics, phenotype
distributions can split into multiple modes under the force of
frequency-dependent competition. In sexual populations, this requires
assortative mating, and we show that the multimodal splitting of initially
unimodal distributions occurs over a range of assortment parameters. In
addition, assortative mating can be favoured evolutionarily even if it incurs
costs, because it provides a means of alleviating the effects of frequency
dependence. Our results reveal that models at both ends of the spectrum between
essentially monomorphic (adaptive dynamics) and fully polymorphic (quantitative
genetics) yield similar results. This underscores that frequency-dependent
selection is a strong agent of pattern formation in phenotype distributions,
potentially resulting in adaptive speciation.


  Proteins are an important class of biomolecules that serve as essential
building blocks of the cells. Their three-dimensional structures are
responsible for their functions. In this thesis we have investigated the
protein structures using a network theoretical approach. While doing so we used
a coarse-grained method, viz., complex network analysis. We model protein
structures at two length scales as Protein Contact Networks (PCN) and as
Long-range Interaction Networks (LINs). We found that proteins by virtue of
being characterised by high amount of clustering, are small-world networks.
Apart from the small-world nature, we found that proteins have another general
property, viz., assortativity. This is an interesting and exceptional finding
as all other complex networks (except for social networks) are known to be
disassortative. Importantly, we could identify one of the major topological
determinant of assortativity by building appropriate controls.


  We have developed a new extended replica exchange method to study
thermodynamics of a system in the presence of external force. Our idea is based
on the exchange between different force replicas to accelerate the equilibrium
process. We have shown that the refolding pathways of single ubiquitin depend
on which terminus is fixed. If the N-end is fixed then the folding pathways are
different compared to the case when both termini are free, but fixing the
C-terminal does not change them. Surprisingly, we have found that the anchoring
terminal does not affect the pathways of individual secondary structures of
three-domain ubiquitin, indicating the important role of the multi-domain
construction. Therefore, force-clamp experiments, in which one end of a protein
is kept fixed, can probe the refolding pathways of a single free-end ubiquitin
if one uses either the poly-ubiquitin or a single domain with the C-terminus
anchored. However, it is shown that anchoring one end does not affect refolding
pathways of the titin domain I27, and the force-clamp spectroscopy is always
capable to predict folding sequencing of this protein. We have obtained the
reasonable estimate for unfolding barrier of ubiqutin. The linkage between
residue Lys48 and the C-terminal of ubiquitin is found to have the dramatic
effect on the location of the transition state along the end-to-end distance
reaction coordinate, but the multi-domain construction leaves the transition
state almost unchanged. We have found that the maximum force in the
force-extension profile from constant velocity force pulling simulations
depends on temperature nonlinearly. However, for some narrow temperature
interval this dependence becomes linear, as have been observed in recent
experiments.


  Feedback circuits in biochemical networks which underly cellular signaling
pathways are important elements in creating complex behavior. A specific aspect
thereof is how stability of equilibrium points depends on model parameters. For
biochemical networks, which are modelled using many parameters, it is typically
very difficult to estimate the influence of parameters on stability. Finding
parameters which result in a change in stability is a key step for a meaningful
bifurcation analysis. We describe a method based on well known approaches from
control theory, which can locate parameters leading to a change in stability.
The method considers a feedback circuit in the biochemical network and relates
stability properties to the control system obtained by loop--breaking. The
method is applied to a model of a MAPK cascade as an illustrative example.


  Starting from linear chains of amino acids, the spontaneous folding of
proteins into their elaborate three-dimensional structures is one of the
remarkable examples of biological self-organization. We investigated native
state structures of 30 single-domain, two-state proteins, from complex networks
perspective, to understand the role of topological parameters in proteins'
folding kinetics, at two length scales-- as ``Protein Contact Networks (PCNs)''
and their corresponding ``Long-range Interaction Networks (LINs)'' constructed
by ignoring the short-range interactions. Our results show that, both PCNs and
LINs exhibit the exceptional topological property of ``assortative mixing''
that is absent in all other biological and technological networks studied so
far. We show that the degree distribution of these contact networks is partly
responsible for the observed assortativity. The coefficient of assortativity
also shows a positive correlation with the rate of protein folding at both
short and long contact scale, whereas, the clustering coefficients of only the
LINs exhibit a negative correlation. The results indicate that the general
topological parameters of these naturally-evolved protein networks can
effectively represent the structural and functional properties required for
fast information transfer among the residues facilitating biochemical/kinetic
functions, such as, allostery, stability, and the rate of folding.


  Boolean networks have been the object of much attention, especially since S.
Kauffman proposed them in the 1960's as models for gene regulatory networks.
These systems are characterized by being defined on a Boolean state space and
by simultaneous updating at discrete time steps. Of particular importance for
biological applications are networks in which the indegree for each variable is
bounded by a fixed constant, as was stressed by Kauffman in his original
papers.
  An important question is which conditions on the network topology can rule
out exponentially long periodic orbits in the system. In this paper, we
consider systems with positive feedback interconnections among all variables
(known as cooperative systems), which in a continuous setting guarantees a very
stable dynamics. We show that for an arbitrary constant 0<c<2 and sufficiently
large n there exist n-dimensional cooperative Boolean networks in which both
the indegree and outdegree of each variable is bounded by two, and which
nevertheless contain periodic orbits of length at least c^n. In Part II of this
paper we will prove an inverse result showing that any system with such a
dynamic behavior must in a sense be similar to the example described.


  The birth/death process with mutation describes the evolution of a
population, and displays rich dynamics including clustering and fluctuations.
We discuss an analytical `field-theoretical' approach to the birth/death
process, using a simple dimensional analysis argument to describe evolution as
a `Super-Brownian Motion' in the infinite population limit. The field theory
technique provides corrections to this for large but finite population, and an
exact description at arbitrary population size. This allows a characterisation
of the difference between the evolution of a phenotype, for which strong local
clustering is observed, and a genotype for which distributions are more
dispersed. We describe the approach with sufficient detail for non-specialists.


  This article reviews quantitative methods to estimate the basic reproduction
number of pandemic influenza, a key threshold quantity to help determine the
intensity of interventions required to control the disease. Although it is
difficult to assess the transmission potential of a probable future pandemic,
historical epidemiologic data is readily available from previous pandemics, and
as a reference quantity for future pandemic planning, mathematical and
statistical analyses of historical data are crucial. In particular, because
many historical records tend to document only the temporal distribution of
cases or deaths (i.e. epidemic curve), our review focuses on methods to
maximize the utility of time-evolution data and to clarify the detailed
mechanisms of the spread of influenza. First, we highlight structured epidemic
models and their parameter estimation method which can quantify the detailed
disease dynamics including those we cannot observe directly.
Duration-structured epidemic systems are subsequently presented, offering firm
understanding of the definition of the basic and effective reproduction
numbers. When the initial growth phase of an epidemic is investigated, the
distribution of the generation time is key statistical information to
appropriately estimate the transmission potential using the intrinsic growth
rate. Applications of stochastic processes are also highlighted to estimate the
transmission potential using the similar data. Critically important
characteristics of influenza data are subsequently summarized, followed by our
conclusions to suggest potential future methodological improvements.


  The replicator dynamics of players choosing either mixed or pure strategies
are usually regarded as equivalent, as long as strategies are played with
identical frequencies. In this paper we show that a population of pure
strategists can be invaded by mixed strategists in any two-player game with
equilibrium coexistence upon the addition of an arbitrarily small amount of
noise in the replication process.


  Emergent pattern formation in self-propelled particle (SPP) systems is
extensively studied because it addresses a range of swarming phenomena which
occur without leadership. Here we present a dynamic SPP model in which a
sensory blind zone is introduced into each particle's zone of interaction.
Using numerical simulations we discovered that the degradation of milling
patterns with increasing blind zone ranges undergoes two distinct transitions,
including a new, spatially nonhomogeneous transition that involves cessation of
particles' motion caused by broken symmetries in their interaction fields. Our
results also show the necessity of nearly complete panoramic sensory ability
for milling behavior to emerge in dynamic SPP models, suggesting a possible
relationship between collective behavior and sensory systems of biological
organisms.


  A hierarchical model for the growth of planar arch structures for RNA
secondary structures is presented, and shown to be equivalent to a tree-growth
model. Both models can be solved analytically, giving access to scaling
functions for large molecules, and corrections to scaling, checked by numerical
simulations of up to 6500 bases. The equivalence of both models should be
helpful in understanding more general tree-growth processes.


  Mechanisms of immunity, and of the host-pathogen interactions in general are
among the most fundamental problems of medicine, ecology, and evolution
studies. Here, we present a microscopic, protein-level, sequence-based model of
immune system, with explicitly defined interactions between host and pathogen
proteins.. Simulations of this model show that possible outcomes of the
infection (extinction of cells, survival with complete elimination of viruses,
or chronic infection with continuous coexistence of cells and viruses)
crucially depend on mutation rates of the viral and immunoglobulin proteins.
Infection is always lethal if the virus mutation rate exceeds a certain
threshold. Potent immunoglobulins are discovered in this model via clonal
selection and affinity maturation. Surviving cells acquire lasting immunity
against subsequent infection by the same virus strain. As a second line of
defense cells develop apoptosis-like behavior by reducing their lifetimes to
eliminate viruses. These results demonstrate the feasibility of microscopic
sequence-based models of immune system, where population dynamics of the
evolving B-cells is explicitly tied to the molecular properties of their
proteins.


  Phylogenetic invariants are not the only constraints on site-pattern
frequency vectors for phylogenetic trees. A mutation matrix, by its definition,
is the exponential of a matrix with non-negative off-diagonal entries; this
positivity requirement implies non-trivial constraints on the site-pattern
frequency vectors. We call these additional constraints ``edge-parameter
inequalities.'' In this paper, we first motivate the edge-parameter
inequalities by considering a pathological site-pattern frequency vector
corresponding to a quartet tree with a negative internal edge. This
site-pattern frequency vector nevertheless satisfies all of the constraints
described up to now in the literature. We next describe two complete sets of
edge-parameter inequalities for the group-based models; these constraints are
square-free monomial inequalities in the Fourier transformed coordinates. These
inequalities, along with the phylogenetic invariants, form a complete
description of the set of site-pattern frequency vectors corresponding to
\emph{bona fide} trees. Said in mathematical language, this paper explicitly
presents two finite lists of inequalities in Fourier coordinates of the form
``monomial $\leq 1$,'' each list characterizing the phylogenetically relevant
semialgebraic subsets of the phylogenetic varieties.


  We explore model based techniques of phylogenetic tree inference exercising
Markov invariants. Markov invariants are group invariant polynomials and are
distinct from what is known in the literature as phylogenetic invariants,
although we establish a commonality in some special cases. We show that the
simplest Markov invariant forms the foundation of the Log-Det distance measure.
We take as our primary tool group representation theory, and show that it
provides a general framework for analysing Markov processes on trees. From this
algebraic perspective, the inherent symmetries of these processes become
apparent, and focusing on plethysms, we are able to define Markov invariants
and give existence proofs. We give an explicit technique for constructing the
invariants, valid for any number of character states and taxa. For phylogenetic
trees with three and four leaves, we demonstrate that the corresponding Markov
invariants can be fruitfully exploited in applied phylogenetic studies.


  I discuss the symmetry of fullerenes, viruses and geodesic domes within a
unified framework of icosadeltahedral representation of these objects. The
icosadeltahedral symmetry is explained in details by examination of all of
these structures. Using Euler's theorem on polyhedra, it is shown how to
calculate the number of vertices, edges, and faces in domes, and number of
atoms, bonds and pentagonal and hexagonal rings in fullerenes. Caspar-Klug
classification of viruses is elaborated as a specific case of icosadeltahedral
geometry.


  The development of a large non-coding fraction in eukaryotic DNA and the
phenomenon of the code-bloat in the field of evolutionary computations show a
striking similarity. This seems to suggest that (in the presence of mechanisms
of code growth) the evolution of a complex code can't be attained without
maintaining a large inactive fraction. To test this hypothesis we performed
computer simulations of an evolutionary toy model for Turing machines, studying
the relations among fitness and coding/non-coding ratio while varying mutation
and code growth rates. The results suggest that, in our model, having a large
reservoir of non-coding states constitutes a great (long term) evolutionary
advantage.


  Phylogenetic networks are a generalization of phylogenetic trees that allow
for the representation of evolutionary events acting at the population level,
like recombination between genes, hybridization between lineages, and lateral
gene transfer. While most phylogenetics tools implement a wide range of
algorithms on phylogenetic trees, there exist only a few applications to work
with phylogenetic networks, and there are no open-source libraries either.
  In order to improve this situation, we have developed a Perl package that
relies on the BioPerl bundle and implements many algorithms on phylogenetic
networks. We have also developed a Java applet that makes use of the
aforementioned Perl package and allows the user to make simple experiments with
phylogenetic networks without having to develop a program or Perl script by
herself.
  The Perl package has been accepted as part of the BioPerl bundle. It can be
downloaded from http://dmi.uib.es/~gcardona/BioInfo/Bio-PhyloNetwork.tgz. The
web-based application is available at http://dmi.uib.es/~gcardona/BioInfo/. The
Perl package includes full documentation of all its features.


  Experiments in recent years have vividly demonstrated that gene expression
can be highly stochastic. How protein concentration fluctuations affect the
growth rate of a population of cells, is, however, a wide open question. We
present a mathematical model that makes it possible to quantify the effect of
protein concentration fluctuations on the growth rate of a population of
genetically identical cells. The model predicts that the population's growth
rate depends on how the growth rate of a single cell varies with protein
concentration, the variance of the protein concentration fluctuations, and the
correlation time of these fluctuations. The model also predicts that when the
average concentration of a protein is close to the value that maximizes the
growth rate, fluctuations in its concentration always reduce the growth rate.
However, when the average protein concentration deviates sufficiently from the
optimal level, fluctuations can enhance the growth rate of the population, even
when the growth rate of a cell depends linearly on the protein concentration.
The model also shows that the ensemble or population average of a quantity,
such as the average protein expression level or its variance, is in general not
equal to its time average as obtained from tracing a single cell and its
descendants. We apply our model to perform a cost-benefit analysis of gene
regulatory control. Our analysis predicts that the optimal expression level of
a gene regulatory protein is determined by the trade-off between the cost of
synthesizing the regulatory protein and the benefit of minimizing the
fluctuations in the expression of its target gene. We discuss possible
experiments that could test our predictions.


  We consider a (sub) critical Galton-Watson process with neutral mutations
(infinite alleles model), and decompose the entire population into clusters of
individuals carrying the same allele. We specify the law of this allelic
partition in terms of the distribution of the number of clone-children and the
number of mutant-children of a typical individual. The approach combines an
extension of Harris representation of Galton-Watson processes and a version of
the ballot theorem. Some limit theorems related to the distribution of the
allelic partition are also given.


  Intracellular transport of DNA carriers is a fundamental step of gene
delivery. We present here a theoretical approach to study generically a single
virus or DNA particle trafficking in a cell cytoplasm. Cellular trafficking has
been studied experimentally mostly at the macroscopic level, but very little
has been done so far at the microscopic level. We present here a physical model
to account for certain aspects of cellular organization, starting with the
observation that a viral particle trajectory consists of epochs of pure
diffusion and epochs of active transport along microtubules. We define a
general degradation rate to describe the limitations of the delivery of plasmid
or viral particles to the nucleus imposed by various types of direct and
indirect hydrolysis activity inside the cytoplasm. Following a homogenization
procedure, which consists of replacing the switching dynamics by a single
steady state stochastic description, not only can we study the spatio-temporal
dynamics of moving objects in the cytosol, but also estimate the probability
and the mean time to go from the cell membrane to a nuclear pore. Computational
simulations confirm that our model can be used to analyze and interpret viral
trajectories and estimate quantitatively the success of nuclear delivery.


  A construction method for duplex cage structures with icosahedral sym- metry
made out of single-stranded DNA molecules is presented and applied to an
icosidodecahedral cage. It is shown via a mixture of analytic and computer
techniques that there exist realisations of this graph in terms of two circular
DNA molecules. These blueprints for the organisation of a cage structure with a
noncrystallographic symmetry may assist in the design of containers made from
DNA for applications in nanotechnology.


  We present a detailed analysis, based on the Forward Flux Sampling (FFS)
simulation method, of the switching dynamics and stability of two models of
genetic toggle switches, consisting of two mutually-repressing genes encoding
transcription factors (TFs); in one model (the exclusive switch), they mutually
exclude each other's binding, while in the other model (general switch) the two
transcription factors can bind simultaneously to the shared operator region. We
assess the role of two pairs of reactions that influence the stability of these
switches: TF-TF homodimerisation and TF-DNA association/dissociation. We
factorise the flipping rate k into the product of the probability rho(q*) of
finding the system at the dividing surface (separatrix) between the two stable
states, and a kinetic prefactor R. In the case of the exclusive switch, the
rate of TF-operator binding affects both rho(q*) and R, while the rate of TF
dimerisation affects only R. In the case of the general switch both TF-operator
binding and TF dimerisation affect k, R and rho(q*). To elucidate this, we
analyse the transition state ensemble (TSE). For the exclusive switch, varying
the rate of TF-operator binding can drastically change the pathway of
switching, while changing the rate of dimerisation changes the switching rate
without altering the mechanism. The switching pathways of the general switch
are highly robust to changes in the rate constants of both TF-operator and
TF-TF binding, even though these rate constants do affect the flipping rate;
this feature is unique for non-equilibrium systems.


  It has been found in in vitro experiments that cytoskeletal filaments driven
by molecular motors show finite diffusion in sliding motion even in the long
filament limit [Y. Imafuku et al., Biophys. J. 70 (1996) 878-886; N. Noda et
al., Biophys. 1 (2005) 45-53]. This anomalous fluctuation can be an evidence
for cooperativity among the motors in action because fluctuation should be
averaged out for a long filament if the action of each motor is independent. In
order to understand the nature of the fluctuation in molecular motors, we
perform numerical simulations and analyse velocity correlation in three
existing models that are known to show some kind of cooperativity and/or large
diffusion coefficient, i.e. Sekimoto-Tawada model [K. Sekimoto and K. Tawada,
Phys. Rev. Lett. 75 (1995) 180], Prost model [J. Prost et al., Phys. Rev. Lett.
72 (1994) 2652], and Duke model [T. Duke, Proc. Natl. Acad. Sci. USA, 96 (1999)
2770]. It is shown that Prost model and Duke model do not give a finite
diffusion in the long filament limit in spite of collective action of motors.
On the other hand, Sekimoto-Tawada model has been shown to give the diffusion
coefficient that is independent of filament length, but it comes from the long
time correlation whose time scale is proportional to filament length, and our
simulations show that such a long correlation time conflicts with the
experimental time scales. We conclude that none of the three models do not
represent experimental findings. In order to explain the observed anomalous
diffusion, we have to seek for the mechanism that should allow both the
amplitude and the time scale of the velocity correlation to be independent of
the filament length.


  In many stochastic simulations of biochemical reaction networks, it is
desirable to ``coarse-grain'' the reaction set, removing fast reactions while
retaining the correct system dynamics. Various coarse-graining methods have
been proposed, but it remains unclear which methods are reliable and which
reactions can safely be eliminated. We address these issues for a model gene
regulatory network that is particularly sensitive to dynamical fluctuations: a
bistable genetic switch. We remove protein-DNA and/or protein-protein
association-dissociation reactions from the reaction set, using various
coarse-graining strategies. We determine the effects on the steady-state
probability distribution function and on the rate of fluctuation-driven switch
flipping transitions. We find that protein-protein interactions may be safely
eliminated from the reaction set, but protein-DNA interactions may not. We also
find that it is important to use the chemical master equation rather than
macroscopic rate equations to compute effective propensity functions for the
coarse-grained reactions.


  It has recently been proposed, with the help of numerical investigations,
that fluorescence correlation spectroscopy at variable observation area can
reveal the existence of a meshgrid of semi-permeable barriers hindering the
two-dimensional diffusion of tagged particles, such as plasmic membrane
constituents. We present a complete theory confirming and accounting for these
findings. It enables a reliable, quantitative exploitation of experimental data
from which the sub-wavelength mesh size can be extracted. Time scales at which
fluorescence correlation spectroscopy must be performed experimentally are
discussed in detail.


  Background : The carnivorous plants of the genus Nepenthes, widely
distributed in the Asian tropics, rely mostly on nutrients derived from
arthropods trapped in their pitcher-shaped leaves and digested by their
enzymatic fluid. The genus exhibits a great diversity of prey and pitcher forms
and its mechanism of trapping has long intrigued scientists. The slippery inner
surfaces of the pitchers, which can be waxy or highly wettable, have so far
been considered as the key trapping devices. However, the occurrence of species
lacking such epidermal specializations but still effective at trapping insects
suggests the possible implication of other mechanisms. Methodology/Principal
Findings : Using a combination of insect bioassays, high-speed video and
rheological measurements, we show that the digestive fluid of Nepenthes
rafflesiana is highly viscoelastic and that this physical property is crucial
for the retention of insects in its traps. Trapping efficiency is shown to
remain strong even when the fluid is highly diluted by water, as long as the
elastic relaxation time of the fluid is higher than the typical time scale of
insect movements. Conclusions/Significance : This finding challenges the common
classification of Nepenthes pitchers as simple passive traps and is of great
adaptive significance for these tropical plants, which are often submitted to
high rainfalls and variations in fluid concentration. The viscoelastic trap
constitutes a cryptic but potentially widespread adaptation of Nepenthes
species and could be a homologous trait shared through common ancestry with the
sundew (Drosera) flypaper plants. Such large production of a highly
viscoelastic biopolymer fluid in permanent pools is nevertheless unique in the
plant kingdom and suggests novel applications for pest control.


  Miniaturization is a hallmark of modern technologies. Notably, this feature
has not spared molecular biology and its potential applications. Towards
developing more effective therapeutics against cancer, studies began to explore
more than a decade ago how natural tumor suppression could be translated into
antineoplastic drugs. To this end, investigators focused on major constituents
of a central pathway that protects cells against neoplastic transformation: the
nuclear retinoblastoma protein (RB) pathway. As such, peptide mimetics of RB,
p16 and p21 were developed. Likewise, the p53 and von Hippel-Lindau gene
products which affect indirectly the RB pathway provided additional templates
for the development of anti-proliferative peptides. Each of the peptides
derived from these distinct tumor suppressors was made cell-permeable by its
ligation to an amino acid sequence conferring cellular internalization. Details
reviewed here reveal that through the application of such anti-cancer peptide
therapeutics alone or in conjunction whenever synergy is to expect, the dark
era of chemotherapy will likely be overcome, at last.


  As a model for gene and protein interactions we study a set for molecular
catalytic reactions. The model is based on experimentally motivated interaction
network topologies, and is designed to capture some key statistics of gene
expression statistics. We impose a non-linearity to the system by a boundary
condition which guarantees non-negative concentrations of chemical
concentrations and study the system stability quantified by maximum Lyapunov
exponents. We find that the non-negativity constraint leads to a drastic
inflation of those regions in parameter space where the Lyapunov exponent
exactly vanishes. We explain the finding as a self-organized critical
phenomenon. The robustness of this finding with respect to different network
topologies and the role of intrinsic molecular- and external noise is
discussed. We argue that systems with inflated 'edges of chaos' could be much
more easily favored by natural selection than systems where the Lyapunov
exponent vanishes only on a parameter set of measure zero.


  We investigate the evolution of cooperation on a non - growth network model
with death/birth dynamics. Nodes reproduce under selection for higher payoffs
in a prisoners dilemma game played between network neighbours. The mean field
characteristics of the model are explored and an attempt is made to understand
the size dependent behaviour of the model in terms of fluctuations in the
strategy densities. We also briefly comment on the role of strategy mutation in
regulating the strategy densties.


  Representation of intracellular signaling networks as directed graphs allows
for the identification of regulatory motifs. Regulatory motifs are groups of
nodes with the same connectivity structure, capable of processing information.
The bifan motif, made of two source nodes directly cross-regulating two target
nodes, is an over-represented motif in a mammalian cell signaling network and
in transcriptional networks. One example of a bifan is the two MAP-kinases, p38
and JNK that phosphorylate and activate the two transcription factors ATF2 and
Elk-1. We have used a system of coupled ordinary differential equations to
analyze the regulatory capability of this bifan motif by itself, and when it
interacts with other motifs such as positive and negative feedback loops. Our
results indicate that bifans provide temporal regulation of signal propagation
and act as signal sorters, filters, and synchronizers. Bifans that have OR gate
configurations show rapid responses while AND gate bifans can introduce delays
and allow prolongation of signal outputs. Bifans that are AND gates can filter
noisy signal inputs. The p38/JNK-ATF2/Elk-1bifan synchronizes the output of
activated transcription factors. Synchronization is a robust property of bifans
and is exhibited even when the bifan is adjacent to a positive feedback loop.
The presence of the bifan promotes the transcription and translation of the
dual specificity protein phosphatase MKP-1 that inhibits p38 and JNK thus
enabling a negative feedback loop. These results indicate that bifan motifs in
cell signaling networks can contribute to signal processing capability both
intrinsically and by enabling the functions of other regulatory motifs.


  The natural evolution of life seems to proceed through steps characterized by
phases of relatively rapid changes, followed by longer, more stable periods. In
the light of the string-theory derived physical scenario proposed in [1], we
discuss how this behaviour can be related to a sequence of resonances of the
energy of natural sources of radiation and absorption energies of the DNA,
responsible for mutagenesis. In a scenario in which these energy scales run
independently as functions of the age of the Universe, the conditions for
evolutionary mutagenesis are satisfied only at discrete points of the time
axis, and for a short period, corresponding to the width of the resonance. We
consider in particular the evolution of the primates through subsequent steps
of increasing cranio-facial contraction, and the great Eras of life (Paleozoic,
Mesozoic, Cenozoic), showing that the transitions occur at the predicted times
of resonance.


  In order to express specific genes at the right time, the transcription of
genes is regulated by the presence and absence of transcription factor
molecules. With transcription factor concentrations undergoing constant
changes, gene transcription takes place out of equilibrium. In this paper we
discuss a simple mapping between dynamic models of gene expression and
stochastic systems driven out of equilibrium. Using this mapping, results of
nonequilibrium statistical mechanics such as the Jarzynski equality and the
fluctuation theorem are demonstrated for gene expression dynamics. Applications
of this approach include the determination of regulatory interactions between
genes from experimental gene expression data.


  A comparison is made between conventional Michaelis-Menten kinetics and two
models that take into account the duration of the conformational changes that
take place at the molecular level during the catalytic cycle of a monomer. The
models consider the time that elapses from the moment an enzyme-substrate
complex forms until the moment a product molecule is released, as well as the
recovery time needed to reset the conformational change that took place. In the
first model the dynamics is described by a set of delayed differential
equations, instead of the ordinary differential equations associated to
Michaelis-Menten kinetics. In the second model the delay, the discretization
inherent to enzyme reactions and the stochastic binding of substrates to
enzimes at the molecular level is considered. All three models agree at
equilibrium, as expected; however, out-of-equilibrium dynamics can differ
substantially. In particular, both delayed models show oscillations at low
values of the Michaelis constant which are not reproduced by the
Michaelis-Menten model. Additionally, in certain cases, the dynamics shown by
the continuous delayed model differs from the dynamics of the discrete delayed
model when some reactant become scarce.


  Both external environmental selection and internal lower-level evolution are
essential for an integral picture of evolution. This paper proposes that the
division of internal evolution into DNA/RNA pattern formation (genotype) and
protein functional action (phenotype) resolves a universal conflict between
fitness and evolvability. Specifically, this paper explains how this universal
conflict drove the emergence of genotype-phenotype division, why this labor
division is responsible for the extraordinary complexity of life, and how the
specific ways of genotype-phenotype mapping in the labor division determine the
paths and forms of evolution and development.


  Complex networks provide us a new view for investigation of immune systems.
In this paper we collect data through STRING database and present a model with
cooperation network theory. The cytokine-protein network model we consider is
constituted by two kinds of nodes, one is immune cytokine types which can act
as acts, other one is protein type which can act as actors. From act degree
distribution that can be well described by typical SPL -shifted power law
functions, we find that
HRAS.TNFRSF13C.S100A8.S100A1.MAPK8.S100A7.LIF.CCL4.CXCL13 are highly
collaborated with other proteins. It reveals that these mediators are important
in cytokine-protein network to regulate immune activity. Dyad act degree
distribution is another important property to generalized collaboration
network. Dyad is two proteins and they appear in one cytokine collaboration
relationship. The dyad act degree distribution can be well described by typical
SPL functions. The length of the average shortest path is 1.29. These results
show that this model could describe the cytokine-protein collaboration
preferably


  Chemotaxis can be defined as an innate behavioural response by an organism to
a directional stimulus, in which bacteria, and other single-cell or
multicellular organisms direct their movements according to certain chemicals
in their environment. This is important for bacteria to find food (e.g.,
glucose) by swimming towards the highest concentration of food molecules, or to
flee from poisons. Based on self-organized computational approaches and similar
stigmergic concepts we derive a novel swarm intelligent algorithm. What strikes
from these observations is that both eusocial insects as ant colonies and
bacteria have similar natural mechanisms based on stigmergy in order to emerge
coherent and sophisticated patterns of global collective behaviour. Keeping in
mind the above characteristics we will present a simple model to tackle the
collective adaptation of a social swarm based on real ant colony behaviors (SSA
algorithm) for tracking extrema in dynamic environments and highly multimodal
complex functions described in the well-know De Jong test suite. Later, for the
purpose of comparison, a recent model of artificial bacterial foraging (BFOA
algorithm) based on similar stigmergic features is described and analyzed.
Final results indicate that the SSA collective intelligence is able to cope and
quickly adapt to unforeseen situations even when over the same cooperative
foraging period, the community is requested to deal with two different and
contradictory purposes, while outperforming BFOA in adaptive speed. Results
indicate that the present approach deals well in severe Dynamic Optimization
problems.


  The Affymetrix U95 and U133 Latin Square spike-in datasets are reanalysed,
together with a dataset from a version of the U95 spike-in experiment without a
complex non-specific background. The approach uses a physico-chemical model
which includes the effects the specific and non-specific hybridisation and
probe folding at the microarray surface, target folding and hybridisation in
the bulk RNA target solution, and duplex dissociation during the
post-hybridisatoin washing phase. The model predicts a three parameter
hyperbolic response function that fits well with fluorescence intensity data
from all three datasets. The importance of the various hybridisation and
washing effects in determining each of the three parameters is examined, and
some guidance is given as to how a practical algorithm for determining specific
target concentrations might be developed.


  We present a statistical-mechanical analysis of the positioning of
nucleosomes along one of the chromosomes of yeast DNA as a function of the
strength of the binding potential and of the chemical potential of the
nucleosomes. We find a significant density of two-level nucleosome switching
regions where, as a function of the chemical potential, the nucleosome
distribution undergoes a "micro" first-order transition. The location of these
nucleosome switches shows a strong correlation with the location of
transcription-factor binding sites.


  Population stratification is a problem encountered in several areas of
biology and public health. We tackle this problem by mapping a population and
its elements attributes into a hypergraph, a natural extension of the concept
of graph or network to encode associations among any number of elements. On
this hypergraph, we construct a statistical model reflecting our intuition
about how the elements attributes can emerge from a postulated population
structure. Finally, we introduce the concept of stratification
representativeness as a mean to identify the simplest stratification already
containing most of the information about the population structure. We
demonstrate the power of this framework stratifying an animal and a human
population based on phenotypic and genotypic properties, respectively.


  Thermal folding Molecular Dynamics simulations of the domain C5 from Myosin
Binding Protein C were performed using a native-centric model to study the role
of three mutations related to Familial Hypertrophic Cardiomyopathy. Mutation of
Asn755 causes the largest shift of the folding temperature, and the residue is
located in the CFGA' beta-sheet featuring the highest Phi-values. The mutation
thus appears to reduce the thermodynamic stability in agreement with
experimental data. The mutations on Arg654 and Arg668, conversely, cause a
little change in the folding temperature and they reside in the low Phi-value
BDE beta-sheet, so that their pathologic role cannot be related to impairment
of the folding process but possibly to the binding with target molecules. As
the typical signature of Domain C5 is the presence of a longer and
destabilizing CD-loop with respect to the other Ig-like domains we completed
the work with a bioinformatic analysis of this loop showing a high density of
negative charge and low hydrophobicity. This indicates the CD-loop as a
natively unfolded sequence with a likely coupling between folding and ligand
binding.


  The precision of biochemical signaling is limited by randomness in the
diffusive arrival of molecules at their targets. For proteins binding to the
specific sites on the DNA and regulating transcription, the ability of the
proteins to diffuse in one dimension by sliding along the length of the DNA, in
addition to their diffusion in bulk solution, would seem to generate a larger
target for DNA binding, consequently reducing the noise in the occupancy of the
regulatory site. Here we show that this effect is largely cancelled by the
enhanced temporal correlations in one dimensional diffusion. With realistic
parameters, sliding along DNA has surprisingly little effect on the physical
limits to the precision of transcriptional regulation.


  Natively unstructured proteins defy the classical "one sequence-one
structure" paradigm of protein science. Monomers of these proteins in
pathological conditions can aggregate in the cell, a process that underlies
socially relevant neurodegenerative diseases such as Alzheimer and Parkinson. A
full comprehension of the formation and structure of the so-called misfolded
intermediates from which the aggregated states ensue is still lacking. We
characterized the folding and the conformational diversity of alpha-synuclein
(aSyn), a natively unstructured protein involved in Parkinson disease, by
mechanically stretching single molecules of this protein and recording their
mechanical properties. These experiments permitted us to directly observe
directly and quantify three main classes of conformations that, under in vitro
physiological conditions, exist simultaneously in the aSyn sample, including
disordered and "beta-like" structures. We found that this class of "beta-like"
structures is directly related to aSyn aggregation. In fact, their relative
abundance increases drastically in three different conditions known to promote
the formation of aSyn fibrils: the presence of Cu2+, the occurrence of the
pathogenic A30P mutation, and high ionic strength. We expect that a critical
concentration of aSyn with a "beta-like" structure must be reached to trigger
fibril formation. This critical concentration is therefore controlled by a
chemical equilibrium. Novel pharmacological strategies can now be tailored to
act upstream, before the aggregation process ensues, by targeting this
equilibrium. To this end, Single Molecule Force Spectroscopy can be an
effective tool to tailor and test new pharmacological agents.


  This work applies the methods of signal processing and the concepts of
control system design to model the maintenance and modulation of reading frame
in the process of protein synthesis. The model shows how translational speed
can modulate translational accuracy to accomplish programmed +1 frameshifts and
could have implications for the regulation of translational efficiency. A
series of free energy estimates were calculated from the ribosome's interaction
with mRNA sequences during the process of translation elongation in eubacteria.
A sinusoidal pattern of roughly constant phase was detected in these free
energy signals. Signal phase was identified as a useful parameter for locating
programmed +1 frameshifts encoded in bacterial genes for release factor 2. A
displacement model was developed that captures the mechanism of frameshift
based on the information content of the signal parameters and the relative
abundance of tRNA in the bacterial cell. Results are presented using
experimentally verified frameshift genes across eubacteria.


  We analyzed folding routes predicted by a variational model in terms of a
generalized formalism of the capillarity scaling theory for 28 two-state
proteins. The scaling exponent ranged from 0.2 to 0.45 with an average of 0.33.
This average value corresponds to packing of rigid objects.That is, on average
the folded core of the nucleus is found to be relatively diffuse. We also
studied the growth of the folding nucleus and interface along the folding route
in terms of the density or packing fraction. The evolution of the folded core
and interface regions can be classified into three patterns of growth depending
on how the growth of the folded core is balanced by changes in density of the
interface. Finally, we quantified the diffuse versus polarized structure of the
critical nucleus through direct calculation of the packing fraction of the
folded core and interface regions. Our results support the general picture of
describing protein folding as the capillarity-like growth of folding nuclei.


  Myosin V is a two-headed processive motor protein that walks in a
hand-over-hand fashion along actin filaments. When it encounters a filament
branch, formed by the Arp2/3 complex, it can either stay on the straight mother
filament, or switch to the daughter filament. We study both probabilities using
the elastic lever arm model for myosin V. We calculate the shapes and bending
energies of all relevant configurations in which the trail head is bound to the
actin filament before Arp2/3 and the lead head is bound either to the mother or
to the daughter filament. Based on the assumption that the probability for a
head to bind to a certain actin subunit is proportional to the Boltzmann factor
obtained from the elastic energy, we calculate the mother/daughter filament
branching ratio. Our model predicts a value of 27% for the daughter and 73% for
the mother filament. This result is in good agreement with recent experimental
data.


  Biomolecular structures are assemblies of emergent anisotropic building
modules such as uniaxial helices or biaxial strands. We provide an approach to
understanding a marginally compact phase of matter that is occupied by proteins
and DNA. This phase, which is in some respects analogous to the liquid crystal
phase for chain molecules, stabilizes a range of shapes that can be obtained by
sequence-independent interactions occurring intra- and intermolecularly between
polymeric molecules. We present a singularityfree self-interaction for a tube
in the continuum limit and show that this results in the tube being positioned
in the marginally compact phase. Our work provides a unified framework for
understanding the building blocks of biomolecules.


  Specific binding of proteins to DNA is one of the most common ways in which
gene expression is controlled. Although general rules for the DNA-protein
recognition can be derived, the ambiguous and complex nature of this mechanism
precludes a simple recognition code, therefore the prediction of DNA target
sequences is not straightforward. DNA-protein interactions can be studied using
computational methods which can complement the current experimental methods and
offer some advantages. In the present work we use physical effective potentials
to evaluate the DNA-protein binding affinities for the lambda repressor-DNA
complex for which structural and thermodynamic experimental data are available.
The effect of conformational sampling by Molecular Dynamics simulations on the
computed binding energy is assessed; results show that this effect is in
general negative and the reproducibility of the experimental values decreases
with the increase of simulation time considered. The free energy of binding for
non-specific complexes agrees with earlier theoretical suggestions. Moreover,
as a results of these analyses, we propose a protocol for the prediction of
DNA-binding target sequences. The possibility of searching regulatory elements
within the bacteriophage-lambda genome using this protocol is explored. Our
analysis shows good prediction capabilities, even in the absence of any
thermodynamic data and information on the naturally recognized sequence. This
study supports the conclusion that physics-based methods can offer a completely
complementary methodology to sequence-based methods for the identification of
DNA-binding protein target sequences.


  We investigate in detail the model of a trophic web proposed by Amaral and
Meyer [Phys. Rev. Lett. 82, 652 (1999)]. We focused on small-size systems that
are relevant for real biological food webs and for which the fluctuations are
playing an important role. We show, using Monte Carlo simulations, that such
webs can be non-viable, leading to extinction of all species in small and/or
weakly coupled systems. Estimations of the extinction times and survival
chances are also given. We show that before the extinction the fraction of
highly-connected species ("omnivores") is increasing. Viable food webs exhibit
a pyramidal structure, where the density of occupied niches is higher at lower
trophic levels, and moreover the occupations of adjacent levels are closely
correlated. We also demonstrate that the distribution of the lengths of food
chains has an exponential character and changes weakly with the parameters of
the model. On the contrary, the distribution of avalanche sizes of the extinct
species depends strongly on the connectedness of the web. For rather loosely
connected systems we recover the power-law type of behavior with the same
exponent as found in earlier studies, while for densely-connected webs the
distribution is not of a power-law type.


  Phylogenetic networks provide a way to describe and visualize evolutionary
histories that have undergone so-called reticulate evolutionary events such as
recombination, hybridization or horizontal gene transfer. The level k of a
network determines how non-treelike the evolution can be, with level-0 networks
being trees. We study the problem of constructing level-k phylogenetic networks
from triplets, i.e. phylogenetic trees for three leaves (taxa). We give, for
each k, a level-k network that is uniquely defined by its triplets. We
demonstrate the applicability of this result by using it to prove that (1) for
all k of at least one it is NP-hard to construct a level-k network consistent
with all input triplets, and (2) for all k it is NP-hard to construct a level-k
network consistent with a maximum number of input triplets, even when the input
is dense. As a response to this intractability we give an exact algorithm for
constructing level-1 networks consistent with a maximum number of input
triplets.


  As the past decade barely dawned, a fundamentally novel view of cancer
relating to signal transduction through intracellular hormones/growth factors
and their subunits began to unfold. Further along, it gained additional
substance with the advent of the interdisciplinary fields of particle biology
and peptide strings which explain (onco)protein dynamics in spacetime, for
instance insulin-driven sub- and trans-cellular carcinogenesis, by physical
principles. Here, this new understanding is expanded to introduce the concept
of "oncoprotein metastasis" preceding cancer cell spread and, thereby, a
particular emphasis is placed on its potential role in the emergence of the
pre-metastatic niche. Consistent with this perception, yet unlike currently
advocated treatments that target cancer cells only, future antineoplastic
strategies should aim to mimic natural tumor suppressors as well as involve
both (morphologically) normal and malignant cells. If validated in human
patients with advanced cancer disease, its otherwise frequently lethal course
may be halted and reversed just in time.


  A ubiquitous building block of signaling pathways is a cycle of covalent
modification (e.g., phosphorylation and dephosphorylation in MAPK cascades).
Our paper explores the kind of information processing and filtering that can be
accomplished by this simple biochemical circuit.
  Signaling cycles are particularly known for exhibiting a highly sigmoidal
(ultrasensitive) input-output characteristic in a certain steady-state regime.
Here we systematically study the cycle's steady-state behavior and its response
to time-varying stimuli. We demonstrate that the cycle can actually operate in
four different regimes, each with its specific input-output characteristics.
These results are obtained using the total quasi-steady-state approximation,
which is more generally valid than the typically used Michaelis-Menten
approximation for enzymatic reactions. We invoke experimental data that
suggests the possibility of signaling cycles operating in one of the new
regimes.
  We then consider the cycle's dynamic behavior, which has so far been
relatively neglected. We demonstrate that the intrinsic architecture of the
cycles makes them act - in all four regimes - as tunable low-pass filters,
filtering out high-frequency fluctuations or noise in signals and environmental
cues. Moreover, the cutoff frequency can be adjusted by the cell. Numerical
simulations show that our analytical results hold well even for noise of large
amplitude. We suggest that noise filtering and tunability make signaling cycles
versatile components of more elaborate cell signaling pathways.


  We probe the dynamic strength of multiple biotin-streptavidin adhesion bonds
under linear loading using the biomembrane force probe setup for dynamic force
spectroscopy. Measured rupture force histograms are compared to results from a
master equation model for the stochastic dynamics of bond rupture under load.
This allows us to extract the distribution of the number of initially closed
bonds. We also extract the molecular parameters of the adhesion bonds, in good
agreement with earlier results from single bond experiments. Our analysis shows
that the peaks in the measured histograms are not simple multiples of the
single bond values, but follow from a superposition procedure which generates
different peak positions.


  The distribution of fitness effects of adaptive mutations remains poorly
understood, both empirically and theoretically. We study this distribution
using a version of Fisher's geometrical model without pleiotropy, such that
each mutation affects only a single trait. We are motivated by the notion of an
organism's chemotype, the set of biochemical reaction constants that govern its
molecular constituents. From physical considerations, we expect the chemotype
to be of high dimension and to exhibit very little pleiotropy. Our model
generically predicts striking cusps in the distribution of the fitness effects
of arising and fixed mutations. It further predicts that a single element of
the chemotype should comprise all mutations at the high-fitness ends of these
distributions. Using extreme value theory, we show that the two cusps with the
highest fitnesses are typically well-separated, even when the chemotype
possesses thousands of elements; this suggests a means to observe these cusps
experimentally. More broadly, our work demonstrates that new insights into
evolution can arise from the chemotype perspective, a perspective between the
genotype and the phenotype.


  This paper develops simplified mathematical models describing the
mutation-selection balance for the asexual and sexual replication pathways in
{\it Saccharomyces cerevisiae}. We assume diploid genomes consisting of two
chromosomes, and we assume that each chromosome is functional if and only if
its base sequence is identical to some master sequence. The growth and
replication of the yeast cells is modeled as a first-order process, with
first-order growth rate constants that are determined by whether a given genome
consists of zero, one, or two functional chromosomes. In the asexual pathway,
we assume that a given diploid cell divides into two diploids. In the sexual
pathway, we assume that a given diploid cell divides into two diploids, each of
which then divide into two haploids. The resulting four haploids enter a
haploid pool, where they grow and replicate until they meet another haploid
with which to fuse. When the cost for sex is low, we find that the selective
mating strategy leads to the highest mean fitness of the population, when
compared to all of the other strategies. We also show that, at low to
intermediate replication fidelities, sexual replication with random mating has
a higher mean fitness than asexual replication, as long as the cost for sex is
low. This is consistent with previous work suggesting that sexual replication
is advantageous at high population densities, low replication rates, and
intermediate replication fidelities. The results of this paper also suggest
that {\it S. cerevisiae} switches from asexual to sexual replication when
stressed, because stressful growth conditions provide an opportunity for the
yeast to clear out deleterious mutations from their genomes.


  We examine the dynamics of biological nanomotors within a simple model of a
rotor having three ion-binding sites. It is shown that in the presence of an
external dc electric field in the plane of the rotor, the loading of the ion
from the positive side of a membrane (rotor charging) provides a torque leading
to the motor rotation. We derive equations for the proton populations of the
sites and solve these equations numerically jointly with the Langevin-type
equation for the rotor angle. Using parameters for biological systems, we
demonstrate that the sequential loading and unloading of the sites lead to the
unidirectional rotation of the motor. The previously unexplained phenomenon of
fast direction-switching in the rotation of a bacterial flagellar motor can
also be understood within our model.


  A molecular dynamics calculation of the amino acid polar requirement is
presented and used to score the canonical genetic code. Monte Carlo simulation
shows that this computational polar requirement has been optimized by the
canonical genetic code more than any previously-known measure. These results
strongly support the idea that the genetic code evolved from a communal state
of life prior to the root of the modern ribosomal tree of life.


  Virus trafficking is fundamental for infection success and plasmid cytosolic
trafficking is a key step of gene delivery. Based on the main physical
properties of the cellular transport machinery such as microtubules, motor
proteins, our goal here is to derive a mathematical model to study cytoplasmic
trafficking. Because experimental results reveal that both active and passive
movement are necessary for a virus to reach the cell nucleus, by taking into
account the complex interactions of the virus with the microtubules, we derive
here an estimate of the mean time a virus reaches the nucleus. In particular,
we present a mathematical procedure in which the complex viral movement,
oscillating between pure diffusion and a deterministic movement along
microtubules, can be approximated by a steady state stochastic equation with a
constant effective drift. An explicit expression for the drift amplitude is
given as a function of the real drift, the density of microtubules and other
physical parameters. The present approach can be used to model viral
trafficking inside the cytoplasm, which is a fundamental step of viral
infection, leading to viral replication and in some cases to cell damage.


  Proteus mirabilis are bacteria that make strikingly regular spatial-temporal
patterns on agar surfaces. In this paper we investigate a mathematical model
that has been shown to display these structures when solved numerically. The
model consists of an ordinary differential equation coupled with a partial
differential equation involving a first-order hyperbolic aging term together
with nonlinear degenerate diffusion. The system is shown to admit global weak
solutions.


  In this paper we explain the biophysical principles which, according to us,
govern the Proteus mirabilis swarm phenomenon. Then, we translate this
explanation into a mathematical model, essentially based on partial
differential equations. This model is then implemented using numerical methods
of the finite volume type in order to make simulations. The simulations show
most of the characteristics which are observed in situ and in particular the
terrace generation.


  The mean time required by a transcription factor (TF) or an enzyme to find a
target in the nucleus is of prime importance for the initialization of
transcription, gene activation or the start of DNA repair. We obtain new
estimates for the mean search time when the TF or enzyme, confined to the cell
nucleus, can switch from a one dimensional motion along the DNA and a free
Brownian regime inside the crowded nucleus. We give analytical expressions for
the mean time the particle stays bound to the DNA, $\tau_{DNA}$, and the mean
time it diffuses freely, $\tau_{free}$. Contrary to previous results but in
agreement with experimental data, we find a factor $\tau_{DNA} \approx 3.7
\tau_{free}$ for the Lac-I TF. The formula obtained for the time required to
bind to a target site is found to be coherent with observed data. We also
conclude that a higher DNA density leads to a more efficient search process.


  DNA microarrays are devices that are able, in principle, to detect and
quantify the presence of specific nucleic acid sequences in complex biological
mixtures. The measurement consists in detecting fluorescence signals from
several spots on the microarray surface onto which different probe sequences
are grafted. One of the problems of the data analysis is that the signal
contains a noisy background component due to non-specific binding. This paper
presents a physical model for background estimation in Affymetrix Genechips. It
combines two different approaches. The first is based on the sequence
composition, specifically its sequence dependent hybridization affinity. The
second is based on the strong correlation of intensities from locations which
are the physical neighbors of a specific spot on the chip. Both effects are
incorporated in a background functional which contains 24 free parameters,
fixed by minimization on a training data set. In all data analyzed the sequence
specific parameters, obtained by minimization, are found to strongly correlate
with empirically determined stacking free energies for RNA/DNA hybridization in
solution. Moreover, there is an overall agreement with experimental background
data and we show that the physics-based model proposed in this paper performs
on average better than purely statistical approaches for background
calculations. The model thus provides an interesting alternative method for
background subtraction schemes in Affymetrix Genechips.


  Recently, microRNAs (miRNAs) have emerged as central posttranscriptional
regulators of gene expression. miRNAs regulate many key biological processes,
including cell growth, death, development and differentiation. This discovery
is challenging the central dogma of molecular biology. Genes are working
together by forming cellular networks. It has become an emerging concept that
miRNAs could intertwine with cellular networks to exert their function. Thus,
it is essential to understand how miRNAs take part in cellular processes at a
systems-level. In this review, I will first introduce basic knowledge of miRNAs
and their relations to heart disaeses and cancer, highlight recently dicovered
functions such as filtering out gene expression noise by miRNAs. I will aslo
introduce basic concepts of cellular networks and interpret their biological
meaning in such a way that the network concepts are digested in a biological
context and are understandable for biologists. Finally, I will summarize the
most recent progress in understanding of miRNA biology at a systems-level: the
principles of miRNA regulation of the major cellular networks including
signaling, metabolic, protein interaction and gene regulatory networks. A
common miRNA regulatory principle is emerging: miRNAs preferentially regulated
the genes that have high regulation complexity. In addition, miRNAs
preferentially regulate positive regulatory motifs, highly connected scaffolds
and the most network downstream components of cellular signaling networks,
while miRNAs selectively regulate the genes which have specific network
structural features on metabolic networks.


  We present a comprehensive computational study of some 900 possible
"lambda-lac" mutants of the lysogeny maintenance switch in phage lambda, of
which up to date 19 have been studied experimentally (Atsumi & Little, PNAS
103: 4558-4563, (2006)). We clarify that these mutants realise regulatory
schemes quite different from wild-type lambda, and can therefore be expected to
behave differently, within the conventional mechanistic setting in which this
problem has often been framed. We verify that indeed, within this framework,
across this wide selection of mutants the lambda-lac mutants for the most part
either have no stable lytic states, or should only be inducible with
difficulty. In particular, the computational results contradicts the
experimental finding that four lambda-lac mutants both show stable lysogeny and
are inducible. This work hence suggests either that the four out of 900 mutants
are special, or that lambda lysogeny and inducibility are holistic effects
involving other molecular players or other mechanisms, or both. The approach
illustrates the power and versatility of computational systems biology to
systematically and quickly test a wide variety of examples and alternative
hypotheses for future closer experimental studies.


  Genomic alterations lead to cancer complexity and form a major hurdle for a
comprehensive understanding of the molecular mechanisms underlying oncogenesis.
In this review, we describe the recent advances in studying cancer-associated
genes from a systems biological point of view. The integration of known cancer
genes onto protein and signaling networks reveals the characteristics of cancer
genes within networks. This approach shows that cancer genes often function as
network hub proteins which are involved in many cellular processes and form
focal nodes in the information exchange between many signaling pathways.
Literature mining allows constructing gene-gene networks, in which new cancer
genes can be identified. The gene expression profiles of cancer cells are used
for reconstructing gene regulatory networks. By doing so, the genes, which are
involved in the regulation of cancer progression, can be picked up from these
networks after which their functions can be further confirmed in the
laboratory.


  We present a kinetic Monte Carlo method for simulating chemical
transformations specified by reaction rules, which can be viewed as generators
of chemical reactions, or equivalently, definitions of reaction classes. A rule
identifies the molecular components involved in a transformation, how these
components change, conditions that affect whether a transformation occurs, and
a rate law. The computational cost of the method, unlike conventional
simulation approaches, is independent of the number of possible reactions,
which need not be specified in advance or explicitly generated in a simulation.
To demonstrate the method, we apply it to study the kinetics of multivalent
ligand-receptor interactions. We expect the method will be useful for studying
cellular signaling systems and other physical systems involving aggregation
phenomena.


  We develop coarse-grained models that describe the dynamic encapsidation of
functionalized nanoparticles by viral capsid proteins. We find that some forms
of cooperative interactions between protein subunits and nanoparticles can
dramatically enhance rates and robustness of assembly, as compared to the
spontaneous assembly of subunits into empty capsids. For large core-subunit
interactions, subunits adsorb onto core surfaces en masse in a disordered
manner, and then undergo a cooperative rearrangement into an ordered capsid
structure. These assembly pathways are unlike any identified for empty capsid
formation. Our models can be directly applied to recent experiments in which
viral capsid proteins assemble around the functionalized inorganic
nanoparticles [Sun et al., Proc. Natl. Acad. Sci (2007) 104, 1354]. In
addition, we discuss broader implications for understanding the dynamic
encapsidation of single-stranded genomic molecules during viral replication and
for developing multicomponent nanostructured materials.


  From the response to external stimuli to cell division and death, the
dynamics of living cells is based on the expression of specific genes at
specific times. The decision when to express a gene is implemented by the
binding and unbinding of transcription factor molecules to regulatory DNA.
Here, we construct stochastic models of gene expression dynamics and test them
on experimental time-series data of messenger-RNA concentrations. The models
are used to infer biophysical parameters of gene transcription, including the
statistics of transcription factor-DNA binding and the target genes controlled
by a given transcription factor.


  The genetic code markup is the assignment of stop codons. The standard
genetic code markup ensures the maximum possible stability of genetic
information with respect to two fault classes: frameshift and nonsense
mutations. There are only 528 (about 1,3% of total number) optimal markups in
the set of markups having 3 stop codons. Among the sets of markups with
1,2,...,8 stop codons, the standard case having 3 stop codons has maximum
absolute number of optimal markups.


  The coevolutionary dynamics in finite populations currently is investigated
in a wide range of disciplines, as chemical catalysis, biological evolution,
social and economic systems. The dynamics of those systems can be formulated
within the unifying framework of evolutionary game theory. However it is not a
priori clear which mathematical description is appropriate when populations are
not infinitely large. Whereas the replicator equation approach describes the
infinite population size limit by deterministic differential equations, in
finite populations the dynamics is inherently stochastic which can lead to new
effects. Recently, an explicit mean-field description in the form of a
Fokker-Planck equation was derived for frequency-dependent selection in finite
populations based on microscopic processes. In asymmetric conflicts between two
populations with a cyclic dominance, a finite-size dependent drift reversal was
demonstrated, depending on the underlying microscopic process of the
evolutionary update. Cyclic dynamics appears widely in biological coevolution,
be it within a homogeneous population, or be it between disjunct populations as
female and male. Here explicit analytic address is given and the average drift
is calculated for the frequency-dependent Moran process and for different
pairwise comparison processes. It is explicitely shown that the drift reversal
cannot occur if the process relies on payoff differences between pairs of
individuals. Further, also a linear comparison with the average payoff does not
lead to a drift towards the internal fixed point. Hence the nonlinear
comparison function of the frequency-dependent Moran process, together with its
usage of nonlocal information via the average payoff, is the essential part of
the mechanism.


  Systems biology focuses on the study of entire biological systems rather than
on their individual components. With the emergence of high-throughput data
generation technologies for molecular biology and the development of advanced
mathematical modeling techniques, this field promises to provide important new
insights. At the same time, with the availability of increasingly powerful
computers, computer algebra has developed into a useful tool for many
applications. This article illustrates the use of computer algebra in systems
biology by way of a well-known gene regulatory network, the Lac Operon in the
bacterium E. coli.


  It has been shown that one can accommodate data (Bayes) and constraints
(MaxEnt) in one method, the method of Maximum (relative) Entropy (ME) (Giffin
2007). In this paper we show a complex agent based example of inference with
two different forms of information; moments and data. In this example, several
agents each receive partial information about a system in the form of data. In
addition, each agent agrees or is informed that there are certain global
constraints on the system that are always true. The agents are then asked to
make inferences about the entire system. The system becomes more complex as we
add agents and allow them to share information. This system can have a
geometrical form, such as a crystal structure. The shape may dictate how the
agents are able to share information, such as sharing with nearest neighbors.
This method can be used to model many systems where the agents or cells have
local or partial information but must adhere to some global rules. This could
also illustrate how the agents evolve and could illuminate emergent behavior of
the system.


  KIF1A kinesins are single-headed motor proteins which move on cylindrical
nano-tubes called microtubules (MT). A normal MT consists of 13 protofilaments
on which the equispaced motor binding sites form a periodic array. The
collective movement of the kinesins on a MT is, therefore, analogous to
vehicular traffic on multi-lane highways where each protofilament is the
analogue of a single lane. Does lane-changing increase or decrease the motor
flux per lane? We address this fundamental question here by appropriately
extending a recent model [{\it Phys. Rev. E {\bf 75}, 041905 (2007)}]. By
carrying out analytical calculations and computer simulations of this extended
model, we predict that the flux per lane can increase or decrease with the
increasing rate of lane changing, depending on the concentrations of motors and
the rate of hydrolysis of ATP, the ``fuel'' molecules. Our predictions can be
tested, in principle, by carrying out {\it in-vitro} experiments with
fluorescently labelled KIF1A molecules.


  Estimation of molecular evolutionary divergence times requires models of rate
change. These vary with regard to the assumption of what quantity is penalized.
The possibilities considered are the rate of evolution, the log of the rate of
evolution and the inverse of the rate of evolution. These models also vary with
regard to how time affects the expected variance of rate change. Here the
alternatives are not at all, linearly with time and as the product of rate and
time. This results in a set of nine models, both random walks and Brownian
motion. A priori any of these models could be correct, yet different
researchers may well prefer, or simply use, one rather than the others. Another
variable is whether to use a scaling factor to take account of the variance of
the process of rate change being unknown and therefore avoid minimizing the
penalty function with unrealistically large times. Here the difference these
models and assumptions make on a tree of mammals, with the root fixed and with
a single internal node fixed, is measured. The similarity of models is measured
as the correlation of their time estimates and visualized with a least squares
tree. The fit of model to data is measured and Q-Q plots are shown. Comparing
model estimates with each other, the age of clades within Laurasiatheria are
seen to vary far more across models than those within Supraprimates (informally
called Euarchontoglires). Especially problematic are the often-used fossil
calibrated nodes of horse/rhino and whale/hippo clashing with times within
Supraprimates and in particular no fossil rodent teeth older than ~60 mybp. A
scaling factor in addition to penalizing rate change is seen to yield
consistent relative time estimates irrespective of exactly where the
calibration point is placed.


  The growth rate of organisms depends both on external conditions and on
internal states, such as the expression levels of various genes. We show that
to achieve a criterion mean growth rate over an ensemble of conditions, the
internal variables must carry a minimum number of bits of information about
those conditions. Evolutionary competition thus can select for cellular
mechanisms that are more efficient in an abstract, information theoretic sense.
Estimates based on recent experiments suggest that the minimum information
required for reasonable growth rates is close to the maximum information that
can be conveyed through biologically realistic regulatory mechanisms. These
ideas are applicable most directly to unicellular organisms, but there are
analogies to problems in higher organisms, and we suggest new experiments for
both cases.


  This review was written for the Encyclopedia of Complexity and System Science
(Springer-Verlag, Berlin, 2008), and is intended as a guide to the growing
literature which approaches the phenomena of cell biology from a more
theoretical point of view. We begin with the building blocks of cellular
networks, and proceed toward the different classes of models being explored,
finally discussing the "design principles" which have been suggested for these
systems. Although largely a dispassionate review, we do draw attention to areas
where there seems to be general consensus on ideas that have not been tested
very thoroughly and, more optimistically, to areas where we feel promising
ideas deserve to be more fully explored.


  When simulating biological populations under different evolutionary genetic
models, backward or forward strategies can be followed. Backward simulations,
also called coalescent-based simulations, are computationally very efficient.
However, this framework imposes several limitations that forward simulation
does not. In this work, a new simple and efficient model to perform forward
simulation of populations and/or genomes is proposed. The basic idea considers
an individual as the differences (mutations) between this individual and a
reference or consensus genotype. Thus, this individual is no longer represented
by its complete sequence or genotype. An example of the efficiency of the new
model with respect to a more classical forward one is demonstrated. This
example models the evolution of HIV resistance using the B_FR.HXB2 reference
sequence to study the emergence of known resistance mutants to Zidovudine and
Didanosine drugs


  Design of experiments is a branch of statistics that aims to identify
efficient procedures for planning experiments in order to optimize knowledge
discovery. Network inference is a subfield of systems biology devoted to the
identification of biochemical networks from experimental data. Common to both
areas of research is their focus on the maximization of information gathered
from experimentation. The goal of this paper is to establish a connection
between these two areas coming from the common use of polynomial models and
techniques from computational algebra.


  A spatiotemporal noise is assumed to reflect the environmental fluctuation in
a spatially extended tumor system. We introduce firstly the structure factor to
reveal the invasive tumor growth quantitatively. The homogenous environment can
lead to an expansive growth of the tumor cells, while the inhomogenous
environment underlies an infiltrative growth. The different responses of above
two cases are separated by a characteristic critical intensity of the
spatiotemporal noise. Theoretical and numerical results present a close
annotation to the clinical images.


  The total conformational energy is assumed to consist of pairwise interaction
energies between atoms or residues, each of which is expressed as a product of
a conformation-dependent function (an element of a contact matrix, C-matrix)
and a sequence-dependent energy parameter (an element of a contact energy
matrix, E-matrix). Such pairwise interactions in proteins force native
C-matrices to be in a relationship as if the interactions are a Go-like
potential [N. Go, Annu. Rev. Biophys. Bioeng. 12. 183 (1983)] for the native
C-matrix, because the lowest bound of the total energy function is equal to the
total energy of the native conformation interacting in a Go-like pairwise
potential. This relationship between C- and E-matrices corresponds to (a) a
parallel relationship between the eigenvectors of the C- and E-matrices and a
linear relationship between their eigenvalues, and (b) a parallel relationship
between a contact number vector and the principal eigenvectors of the C- and
E-matrices; the E-matrix is expanded in a series of eigenspaces with an
additional constant term, which corresponds to a threshold of contact energy
that approximately separates native contacts from non-native ones. These
relationships are confirmed in 182 representatives from each family of the SCOP
database by examining inner products between the principal eigenvector of the
C-matrix, that of the E-matrix evaluated with a statistical contact potential,
and a contact number vector. In addition, the spectral representation of C- and
E-matrices reveals that pairwise residue-residue interactions, which depends
only on the types of interacting amino acids but not on other residues in a
protein, are insufficient and other interactions including residue
connectivities and steric hindrance are needed to make native structures the
unique lowest energy conformations.


  Using Monte Carlo model of biological evolution we have discovered that
populations can switch between two different strategies of their genomes'
evolution; Darwinian purifying selection and complementing the haplotypes. The
first one is exploited in the large panmictic populations while the second one
in the small highly inbred populations. The choice depends on the crossover
frequency. There is a power law relation between the critical value of
crossover frequency and the size of panmictic population. Under the constant
inbreeding this critical value of crossover does not depend on the population
size and has a character of phase transition. Close to this value sympatric
speciation is observed.


  A model for epidemics on an adaptive network is considered. Nodes follow an
SIRS (susceptible-infective-recovered-susceptible) pattern. Connections are
rewired to break links from non-infected nodes to infected nodes and are
reformed to connect to other non-infected nodes, as the nodes that are not
infected try to avoid the infection. Monte Carlo simulation and numerical
solution of a mean field model are employed. The introduction of rewiring
affects both the network structure and the epidemic dynamics. Degree
distributions are altered, and the average distance from a node to the nearest
infective increases. The rewiring leads to regions of bistability where either
an endemic or a disease-free steady state can exist. Fluctuations around the
endemic state and the lifetime of the endemic state are considered. The
fluctuations are found to exhibit power law behavior.


  Temporally and spatially resolved measurements of protein transport inside
cells provide important clues to the functional architecture and dynamics of
biological systems. Fluorescence Recovery After Photobleaching (FRAP) technique
has been used over the past three decades to measure the mobility of
macromolecules and protein transport and interaction with immobile structures
inside the cell nucleus. A theoretical model is presented that aims to describe
protein transport inside the nucleus, a process which is influenced by the
presence of a boundary (i.e. membrane). A set of reaction-diffusion equations
is employed to model both the diffusion of proteins and their interaction with
immobile binding sites. The proposed model has been designed to be applied to
biological samples with a Confocal Laser Scanning Microscope (CLSM) equipped
with the feature to bleach regions characterised by a scanning beam that has a
radially Gaussian distributed profile. The proposed model leads to FRAP curves
that depend on the on- and off-rates. Semi-analytical expressions are used to
define the boundaries of on- (off-) rate parameter space in simplified cases
when molecules move within a bounded domain. The theoretical model can be used
in conjunction to experimental data acquired by CLSM to investigate the
biophysical properties of proteins in living cells.


  We have used the Penna ageing model to analyze how the differences in
evolution of sex chromosomes depend on the strategy of reproduction. In
panmictic populations, when females (XX) can freely choose the male partner
(XY) for reproduction from the whole population, the Y chromosome accumulates
defects and eventually the only information it brings is a male sex
determination. As a result of shrinking Y chromosome the males become
hemizygous in respect to the X chromosome content and are characterized by
higher mortality, observed also in the human populations. If it is assumed in
the model that the presence of the male is indispensable at least during the
pregnancy of his female partner and he cannot be seduced by another female at
least during the one reproduction cycle - the Y chromosome preserves its
content, does not shrink and the lifespan of females and males is the same.
Thus, Y chromosome shrinks not because of existing in one copy, without the
possibility of recombination, but because it stays under weaker selection
pressure; in panmictic populations without the necessity of being faithful, a
considerable fraction of males is dispensable and they can be eliminated from
the population without reducing its reproduction potential.


  In this paper, we investigate the emergence of a predator-prey system with
Ivlev-type functional response and reaction-diffusion. We study how diffusion
affects the stability of predator-prey coexistence equilibrium and derive the
conditions for Hopf and Turing bifurcation in the spatial domain. Based on the
bifurcation analysis, we give the spatial pattern formation, the evolution
process of the system near the coexistence equilibrium point, via numerical
simulation. We find that pure Hopf instability leads to the formation of spiral
patterns and pure Turing instability destroys the spiral pattern and leads to
the formation of chaotic spatial pattern. Furthermore, we perform three
categories of initial perturbations which predators are introduced in a small
domain to the coexistence equilibrium point to illustrate the emergence of
spatiotemporal patterns, we also find that in the beginning of evolution of the
spatial pattern, the special initial conditions have an effect on the formation
of spatial patterns, though the effect is less and less with the more and more
iterations. This indicates that for prey-dependent type predator-prey model,
pattern formations do depend on the initial conditions, while for
predator-dependent type they do not. Our results show that modeling by
reaction-diffusion equations is an appropriate tool for investigating
fundamental mechanisms of complex spatiotemporal dynamics.


  In this paper, we investigate the emergence of a predator-prey model with
Beddington-DeAngelis-type functional response and reaction-diffusion. We derive
the conditions for Hopf and Turing bifurcation on the spatial domain. Based on
the stability and bifurcation analysis, we give the spatial pattern formation
via numerical simulation, i.e., the evolution process of the model near the
coexistence equilibrium point. We find that for the model we consider, pure
Turing instability gives birth to the spotted pattern, pure Hopf instability
gives birth to the spiral wave pattern, and both Hopf and Turing instability
give birth to stripe-like pattern. Our results show that reaction-diffusion
model is an appropriate tool for investigating fundamental mechanism of complex
spatiotemporal dynamics. It will be useful for studying the dynamic complexity
of ecosystems.


  In this paper we argue that biological fitness is a multi-objective concept
hence the statement "fittest" is inappropriate. The following statement is
proposed "Survival is mostly for those with non-dominated fitness". Also we use
some TV games to show that under the following conditions: i) There are no
dominant players. ii) At each time step successful players may eliminate some
of their less successful competitors, Then the ultimate winner may not be the
fittest (but close).


  We argue that immune system is an adaptive complex system. It is shown that
it has emergent properties. Its network structure is of the small world network
type. The network is of the threshold type, which helps in avoiding
autoimmunity. It has the property that every antigen (e.g.virus or bacteria) is
typically attacked by more than one effector. This stabilizes the equilibrium
state. Modelling complex systems is discussed. Cellular automata (CA) type
models are successful but there are much less analytic results about CA than
about other less successful models e.g. partial differential equations (PDE). A
compromise is proposed


  It is proposed that using both self-non-self and danger theories give a
better understanding of how the immune system works. It is proposed that
comparing immune system to police force is useful in this case since police
responds both to danger or damage signals and to foreign or suspicious behavior
even if no danger signals existed. We also propose that due to low zone
tolerance immunotherapy needs to be combined with another treatment method for
cancer e.g. chemotherapy or/and radiotherapy to get a sufficient eradication of
tumors. Finally we propose that fractional order differential equations are
more suitable than the familiar integer order differential equations. A
fractional order example of two immune effectors attacking an antigen is given.


  The nematode Caenorhabditis elegans responds to an impressive range of
chemical, mechanical and thermal stimuli and is extensively used to investigate
the molecular mechanisms that mediate chemosensation, mechanotransduction and
thermosensation. The main behavioral output of these responses is manifested as
alterations in animal locomotion. Monitoring and examination of such
alterations requires tools to capture and quantify features of nematode
movement. In this paper, we introduce Nemo (nematode movement), a
computationally efficient and robust two-dimensional object tracking algorithm
for automated detection and analysis of C. elegans locomotion. This algorithm
enables precise measurement and feature extraction of nematode movement
components. In addition, we develop a Graphical User Interface designed to
facilitate processing and interpretation of movement data. While, in this
study, we focus on the simple sinusoidal locomotion of C. elegans, our approach
can be readily adapted to handle complicated locomotory behaviour patterns by
including additional movement characteristics and parameters subject to
quantification. Our software tool offers the capacity to extract, analyze and
measure nematode locomotion features by processing simple video files. By
allowing precise and quantitative assessment of behavioral traits, this tool
will assist the genetic dissection and elucidation of the molecular mechanisms
underlying specific behavioral responses.