  Small, non-coding RNAs (sRNAs) play important roles as genetic regulators in
prokaryotes. sRNAs act post-transcriptionally via complementary pairing with
target mRNAs to regulate protein expression. We use a quantitative approach to
compare and contrast sRNAs with conventional transcription factors (TFs) to
better understand the advantages of each form of regulation. In particular, we
calculate the steady-state behavior, noise properties, frequency-dependent gain
(amplification), and dynamical response to large input signals of both forms of
regulation. While the mean steady-state behavior of sRNA-regulated proteins
exhibits a distinctive tunable threshold-linear behavior, our analysis shows
that transcriptional bursting leads to significantly higher intrinsic noise in
sRNA-based regulation than in TF-based regulation in a large range of
expression levels and limits the ability of sRNAs to perform quantitative
signaling. Nonetheless, we find that sRNAs are better than TFs at filtering
noise in input signals. Additionally, we find that sRNAs allow cells to respond
rapidly to large changes in input signals. These features suggest a niche for
sRNAs in allowing cells to transition quickly yet reliably between distinct
states. This functional niche is consistent with the widespread appearance of
sRNAs in stress-response and quasi-developmental networks in prokaryotes.


  Isolated human buccal epithelium cell were irradiated by microwaves at
frequency f=35 GHz and surface power density E=30 mcW/cm2. The state of
chromatin in human cells was determined by methodsof light and electron
microscopy. The state of cell membranes was evaluated by the method of vital
indigo carmine staining. The microwave-induced condensation of chromatin in
human cells was revealed. Left side circulary polarized waves induced less
effect than linearly polarized radiation. The linearly polarized
electromagnetic waves induced cell membrane damage revealed by the increase of
cell stainability. The data obtained are discussed in connection with the
mechanisms of biologica effect of electromagnetic waves.


  The binding of a ligand molecule to a protein is often accompanied by
conformational changes of the protein. A central question is whether the ligand
induces the conformational change (induced-fit), or rather selects and
stabilizes a complementary conformation from a pre-existing equilibrium of
ground and excited states of the protein (selected-fit). We consider here the
binding kinetics in a simple four-state model of ligand-protein binding. In
this model, the protein has two conformations, which can both bind the ligand.
The first conformation is the ground state of the protein when the ligand is
off, and the second conformation is the ground state when the ligand is bound.
The induced-fit mechanism corresponds to ligand binding in the unbound ground
state, and the selected-fit mechanism to ligand binding in the excited state.
We find a simple, characteristic difference between the on- and off-rates in
the two mechanisms if the conformational relaxation into the ground states is
fast. In the case of selected-fit binding, the on-rate depends on the
conformational equilibrium constant, while the off-rate is independent. In the
case of induced-fit binding, in contrast, the off-rate depends on the
conformational equilibrium, while the on-rate is independent. Whether a protein
binds a ligand via selected-fit or induced-fit thus may be revealed by
mutations far from the protein's binding pocket, or other "perturbations" that
only affect the conformational equilibrium. In the case of selected-fit, such
mutations will only change the on-rate, and in the case of induced-fit, only
the off-rate.


  While ordinary differential equations (ODEs) form the conceptual framework
for modelling many cellular processes, specific situations demand stochastic
models to capture the influence of noise. The most common formulation of
stochastic models for biochemical networks is the chemical master equation
(CME). While stochastic simulations are a practical way to realise the CME,
analytical approximations offer more insight into the influence of noise.
Towards that end, the two-moment approximation (2MA) is a promising addition to
the established analytical approaches including the chemical Langevin equation
(CLE) and the related linear noise approximation (LNA). The 2MA approach
directly tracks the mean and (co)variance which are coupled in general. This
coupling is not obvious in CME and CLE and ignored by LNA and conventional ODE
models. We extend previous derivations of 2MA by allowing a) non-elementary
reactions and b) relative concentrations. Often, several elementary reactions
are approximated by a single step. Furthermore, practical situations often
require the use relative concentrations. We investigate the applicability of
the 2MA approach to the well established fission yeast cell cycle model. Our
analytical model reproduces the clustering of cycle times observed in
experiments. This is explained through multiple resettings of MPF, caused by
the coupling between mean and (co)variance, near the G2/M transition.


  Based on the intensity and rates of various kinds of intense ionizing
radiation events such as supernovae and gamma-ray bursts, it is likely that the
Earth has been subjected to one or extinction level events during the
Phanerozoic. These induce changes in atmospheric chemistry so that the level of
Solar ultraviolet-B radiation reaching the surface and near-surface waters may
be doubled for up to a decade. This UVB level is known from experiment to be
more than enough to kill off many kinds of organisms, particularly
phytoplankton. It could easily induce a crash of the photosynthetic-based food
chain in the oceans. Regularities in the latitudinal distribution of damage are
apparent in simulations of the atmospheric changes. We previously proposed that
the late Ordovician extinction is a plausible candidate for a contribution from
an ionizing radiation event, based on environmental selectivity in trilobites.
To test a null hypothesis based on this proposal, we confront latitudinal
differential extinction rates predicted from the simulations with data from a
published analysis of latitudinal gradients in the Ordovician extinction. The
pattern of UVB damage always shows a strong maximum at some latitude, with
substantially lower intensity to the north and south of this maximum. We find
that the pattern of damage predicted from our simulations is consistent with
the data assuming a burst approximately over the South Pole, and no further
north than -75 degrees. We predict that any land mass (such as parts of north
China, Laurentia, and New Guinea) which then lay north of the equator should be
a refuge from UVB effects, and show a different pattern of extinction in the
first strike of the end-Ordovician extinction, if induced by such a radiation
event.


  We study a matrix model of RNA in which an external perturbation acts on n
nucleotides of the polymer chain. The effect of the perturbation appears in the
exponential generating function of the partition function as a factor
$(1-\frac{n\alpha}{L})$ [where $\alpha$ is the ratio of strengths of the
original to the perturbed term and L is length of the chain]. The asymptotic
behaviour of the genus distribution functions for the extended matrix model are
analyzed numerically when (i) $n=L$ and (ii) $n=1$. In these matrix models of
RNA, as $n\alpha/L$ is increased from 0 to 1, it is found that the universality
of the number of diagrams $a_{L, g}$ at a fixed length L and genus g changes
from $3^{L}$ to $(3-\frac{n\alpha}{L})^{L}$ ($2^{L}$ when $n\alpha/L=1$) and
the asymptotic expression of the total number of diagrams $\cal N$ at a fixed
length L but independent of genus g, changes in the factor $\exp^{\sqrt{L}}$ to
$\exp^{(1-\frac{n\alpha}{L})\sqrt{L}}$ ($exp^{0}=1$ when $n\alpha/L=1$)


  The model of facilitated diffusion describes how DNA-binding proteins, such
as transcription factors (TFs), find their chromosomal targets by combining 3D
diffusion through the cytoplasm and 1D sliding along nonspecific DNA sequences.
The redundant 1D diffusion near the specific binding site extends the target
size and facilitates target location. While this model successfully predicts
the kinetics measured in test tubes, it has not been extended to account for
the highly crowded environment in living cells. Here, we investigate the effect
of other DNA-binding proteins that partially occupy the bacterial chromosome.
We show how they would slow down the search process, mainly through restricted
sliding near the target. This implies that increasing the overall DNA-binding
protein concentration would have a marginal effect in reducing the search time,
because any additional proteins would restrict the search process for each
other. While the presence of other proteins prevents sliding from the flanking
DNA, DNA looping provides an alternative path to transfer from neighboring
sites to the target efficiently. We propose that, when looping is faster than
the initial search process, the auxiliary binding sites further extend the
effective target region and therefore facilitate the target location.


  Interpretation of animal behavior, especially as cooperative or selfish, is a
challenge for evolutionary theory. Strategy of a competition should follow from
corresponding Darwinian payoffs for the available behavioral options. The
payoffs and decision making processes, however, are difficult to observe and
quantify. Here we present a general method for the derivation of evolutionary
payoffs from observable statistics of interactions. The method is applied to
combat of male bowl and doily spiders, to predator inspection by sticklebacks
and to territorial defense by lions, demonstrating animal behavior as a new
type of game theoretical equilibrium. Games animals play may be derived
unequivocally from their observable behavior, the reconstruction, however, can
be subjected to fundamental limitations due to our inability to observe all
information exchange mechanisms (communication).


  In protein folding the term plasticity refers to the number of alternative
folding pathways encountered in response to free energy perturbations such as
those induced by mutation. Here we explore the relation between folding
plasticity and a gross, generic feature of the native geometry, namely, the
relative number of local and non-local native contacts. The results from our
study, which is based on Monte Carlo simulations of simple lattice proteins,
show that folding to a structure that is rich in local contacts is considerably
more plastic than folding to a native geometry characterized by having a very
large number of long-range contacts (i.e., contacts between amino acids that
are separated by more than 12 units of backbone distance). The smaller folding
plasticity of `non-local' native geometries is probably a direct consequence of
their higher folding cooperativity that renders the folding reaction more
robust against single- and multiple-point mutations.


  The concept of evolutionary development of structures constituted a
\emph{real} revolution in biology: it was possible to understand how the very
complex structures of life can arise in an out-of-equilibrium system. The
investigation of such systems has shown that indeed, systems under a flux of
energy or matter can self-organize into complex patterns, think for instance to
Rayleigh-Bernard convection, Liesegang rings, patterns formed by granular
systems under shear. Following this line, one could characterize life as a
state of matter, characterized by the slow, continuous process that we call
evolution. In this paper we try to identify the organizational level of life,
that spans several orders of magnitude from the elementary constituents to
whole ecosystems.
  Although similar structures can be found in other contexts like ideas (memes)
in neural systems and self-replicating elements (computer viruses, worms, etc.)
in computer systems, we shall concentrate on biological evolutionary structure,
and try to put into evidence the role and the emergence of network structure in
such systems.


  Copy Number Variations (CNVs) of regions of the human genome are important in
disease association studies.The digital array is a nanofluidic biochip which
utilizes integrated channels and valves that partition mixtures of sample and
reagents into 765 nanovolume reaction chambers. It was recently shown how one
can perform statistical analysis of CNV in a DNA sample the digital array. In
particular, it was shown how one can accurately estimate the true concentration
of the molecules in the DNA sample and then determine the ratios of different
sequences along with statistical confidence intervals on these estimations. In
this paper we perform computation of maximum number of copies which can be
distinguished using the digital array which gives its resolution in terms of
its ability to determine CNV. Then, we demonstrate the usefulness of the
mathematical analysis to solve an important real-world problem of determination
of the copy number of X chromosome as our example application.


  Optimal sequence alignments depend heavily on alignment scoring parameters.
Given input sequences, {\em parametric alignment} is the well-studied problem
that asks for all possible optimal alignment summaries as parameters vary, as
well as the {\em optimality region} of alignment scoring parameters which yield
each optimal alignment. But biologically correct alignments might be {\em
suboptimal} for all parameter choices. Thus we extend parametric alignment to
{\em parametric $k$-best alignment}, which asks for all possible $k$-tuples of
$k$-best alignment summaries $(s_1, s_2, ..., s_k)$, as well as the {\em
$k$-best optimality region} of scoring parameters which make $s_1, s_2, ...,
s_k$ the top $k$ summaries. By exploiting the integer-structure of alignment
summaries, we show that, astonishingly, the complexity of parametric $k$-best
alignment is only polynomial in $k$. Thus parametric $k$-best alignment is
tractable, and can be applied at the whole-genome scale like parametric
alignment.


  The effects of an unconditional move rule in the spatial Prisoner's Dilemma,
Snowdrift and Stag Hunt games are studied. Spatial structure by itself is known
to modify the outcome of many games when compared with a randomly mixed
population, sometimes promoting, sometimes inhibiting cooperation. Here we show
that random dilution and mobility may suppress the inhibiting factors of the
spatial structure in the Snowdrift game, while enhancing the already larger
cooperation found in the Prisoner's dilemma and Stag Hunt games.


  In this article, we present an effective encoding of dendrograms by embedding
them into the Bruhat-Tits trees associated to $p$-adic number fields. As an
application, we show how strings over a finite alphabet can be encoded in
cyclotomic extensions of $\mathbb{Q}_p$ and discuss $p$-adic DNA encoding. The
application leads to fast $p$-adic agglomerative hierarchic algorithms similar
to the ones recently used e.g. by A. Khrennikov and others. From the viewpoint
of $p$-adic geometry, to encode a dendrogram $X$ in a $p$-adic field $K$ means
to fix a set $S$ of $K$-rational punctures on the $p$-adic projective line
$\mathbb{P}^1$. To $\mathbb{P}^1\setminus S$ is associated in a natural way a
subtree inside the Bruhat-Tits tree which recovers $X$, a method first used by
F. Kato in 1999 in the classification of discrete subgroups of
$\textrm{PGL}_2(K)$.
  Next, we show how the $p$-adic moduli space $\mathfrak{M}_{0,n}$ of
$\mathbb{P}^1$ with $n$ punctures can be applied to the study of time series of
dendrograms and those symmetries arising from hyperbolic actions on
$\mathbb{P}^1$. In this way, we can associate to certain classes of dynamical
systems a Mumford curve, i.e. a $p$-adic algebraic curve with totally
degenerate reduction modulo $p$.
  Finally, we indicate some of our results in the study of general discrete
actions on $\mathbb{P}^1$, and their relation to $p$-adic Hurwitz spaces.


  In many epidemiological models a nonlinear transmission function is used in
the form of power law relationship. It is constantly argued that such form
reflects population heterogeneities including differences in the mixing
pattern, susceptibility, and spatial patchiness, although the function itself
is considered phenomenological. Comparison with large-scale simulations show
that models with this transmission function accurately approximate data from
highly heterogeneous sources. In this note we provide a mechanistic derivation
of the power law transmission function, starting with a simple heterogeneous
susceptibles--infectives (SI) model, which is based on a standard mass action
assumption. We also consider the simplest SI model with separable mixing and
compare our results with known results from the literature.


  The G-protein-coupled receptor (GPCR) activated by the neurotransmitter GABA
is made up of two subunits, GABA(B1) and GABA(B2). GABA(B1) binds agonists,
whereas GABA(B2) is required for trafficking GABA(B1) to the cell surface,
increasing agonist affinity to GABA(B1), and activating associated G proteins.
These subunits each comprise two domains, a Venus flytrap domain (VFT) and a
heptahelical transmembrane domain (7TM). How agonist binding to the GABA(B1)
VFT leads to GABA(B2) 7TM activation remains unknown. Here, we used a glycan
wedge scanning approach to investigate how the GABA(B) VFT dimer controls
receptor activity. We first identified the dimerization interface using a
bioinformatics approach and then showed that introducing an N-glycan at this
interface prevents the association of the two subunits and abolishes all
activities of GABA(B2), including agonist activation of the G protein. We also
identified a second region in the VFT where insertion of an N-glycan does not
prevent dimerization, but blocks agonist activation of the receptor. These data
provide new insight into the function of this prototypical GPCR and demonstrate
that a change in the dimerization interface is required for receptor
activation.


  The gamma-aminobutyric acid type B (GABAB) receptor is an allosteric complex
made of two subunits, GABAB1 (GB1) and GABAB2 (GB2). Both subunits are composed
of an extracellular Venus flytrap domain (VFT) and a heptahelical domain (HD).
GB1 binds GABA, and GB2 plays a major role in G-protein activation as well as
in the high agonist affinity state of GB1. How agonist affinity in GB1 is
regulated in the receptor remains unknown. Here, we demonstrate that GB2 VFT is
a major molecular determinant involved in this control. We show that isolated
versions of GB1 and GB2 VFTs in the absence of the HD and C-terminal tail can
form hetero-oligomers as shown by time-resolved fluorescence resonance energy
transfer (based on HTRF technology). GB2 VFT and its association with GB1 VFT
controlled agonist affinity in GB1 in two ways. First, GB2 VFT exerted a direct
action on GB1 VFT, as it slightly increased agonist affinity in isolated GB1
VFT. Second and most importantly, GB2 VFT prevented inhibitory interaction
between the two main domains (VFT and HD) of GB1. According to this model, we
propose that GB1 HD prevents the possible natural closure of GB1 VFT. In
contrast, GB2 VFT facilitates this closure. Finally, such inhibitory contacts
between HD and VFT in GB1 could be similar to those important to maintain the
inactive state of the receptor.


  Here we introduce researchers in algebraic biology to the exciting new field
of cophylogenetics. Cophylogenetics is the study of concomitantly evolving
organisms (or genes), such as host and parasite species. Thus the natural
objects of study in cophylogenetics are tuples of related trees, instead of
individual trees. We review various research topics in algebraic statistics for
phylogenetics, and propose analogs for cophylogenetics. In particular we
propose spaces of cophylogenetic trees, cophylogenetic reconstruction, and
cophylogenetic invariants. We conclude with open problems.


  An impulsive model of augmentative biological control consisting of a general
continuous predator-prey model in ordinary differential equations augmented by
a discrete part describing periodic introductions of predators is considered.
It is shown that there exists an invariant periodic solution that corresponds
to prey eradication and a condition ensuring its global asymptotic stability is
given. An optimisation problem related to the preemptive use of augmentative
biological control is then considered. It is assumed that the per time unit
budget of biological control (i.e. the number of predators to be released) is
fixed and the best deployment of this budget is sought after in terms of
release frequency. The cost function to be minimised is the time taken to
reduce an unforeseen prey (pest) invasion under some harmless level. The
analysis shows that the optimisation problem admits a countable infinite number
of solutions. An argumentation considering the required robustness of the
optimisation result is then conducted and it is shown that the best deployment
is to use as frequent (and thus as small) as possible predator introductions.


  Spatially explicit models have been widely used in today's mathematical
ecology and epidemiology to study persistence and extinction of populations as
well as their spatial patterns. Here we extend the earlier work--static
dispersal between neighbouring individuals to mobility of individuals as well
as multi-patches environment. As is commonly found, the basic reproductive
ratio is maximized for the evolutionary stable strategy (ESS) on diseases'
persistence in mean-field theory. This has important implications, as it
implies that for a wide range of parameters that infection rate will tend
maximum. This is opposite with present results obtained in spatial explicit
models that infection rate is limited by upper bound. We observe the emergence
of trade-offs of extinction and persistence on the parameters of the infection
period and infection rate and show the extinction time having a linear
relationship with respect to system size. We further find that the higher
mobility can pronouncedly promote the persistence of spread of epidemics, i.e.,
the phase transition occurs from extinction domain to persistence domain, and
the spirals' wavelength increases as the mobility increasing and ultimately, it
will saturate at a certain value. Furthermore, for multi-patches case, we find
that the lower coupling strength leads to anti-phase oscillation of infected
fraction, while higher coupling strength corresponds to in-phase oscillation.


  A geometric analysis of protein folding, which complements many of the models
in the literature, is presented. We examine the process from unfolded strand to
the point where the strand becomes self-interacting. A central question is how
it is possible that so many initial configurations proceed to fold to a unique
final configuration. We put energy and dynamical considerations temporarily
aside and focus upon the geometry alone. We parameterize the structure of an
idealized protein using the concept of a ribbon from differential geometry. The
deformation of the ribbon is described by introducing a generic twisting
Ansatz. The folding process in this picture entails a change in shape guided by
the local amino acid geometry. The theory is reparamaterization invariant from
the start, so the final shape is independent of folding time. We develop
differential equations for the changing shape. For some parameter ranges, a
sine-Gordon torsion soliton is found. This purely geometric waveform has
properties similar to dynamical solitons. Namely: A threshold distortion of the
molecule is required to initiate the soliton, after which, small additional
distortions do not change the waveform. In this analysis, the soliton twists
the molecule until bonds form. The analysis reveals a quantitative relationship
between the geometry of the amino acids and the folded form.


  We study the stochastic dynamics of growth and shrinkage of single actin
filaments or microtubules taking into account insertion, removal, and ATP/GTP
hydrolysis of subunits. The resulting phase diagram contains three different
phases: a rapidly growing phase, an intermediate phase and a bound phase. We
analyze all these phases, with an emphasis on the bound phase. We also discuss
how hydrolysis affects force-velocity curves. The bound phase shows features of
dynamic instability, which we characterize in terms of the time needed for the
ATP/GTP cap to disappear as well as the time needed for the filament to reach a
length of zero, i.e., (to collapse) for the first time. We obtain exact
expressions for all these quantities, which we test using Monte Carlo
simulations.


  Significant overweight represents a major health problem in industrialized
countries. Besides its known metabolic origins, this condition may also have an
infectious cause, as recently postulated. Here, it is surmised that the
potentially causative adenovirus 36 contributes to such disorder by
inactivating the retinoblastoma tumor suppressor protein (RB) in a manner
reminiscent of a mechanism employed by both another pathogenic adenoviral agent
and insulin. The present insight additionally suggests novel modes of
interfering with obesity-associated pathology.


  The article continues an analysis of the genetic 8-dimensional
Yin-Yang-algebra. This algebra was revealed in a course of matrix researches of
structures of the genetic code and it was described in the author's articles
arXiv:0803.3330 and arXiv:0805.4692. The article presents data about many kinds
of cyclic permutations of elements of the genetic code in the genetic
(8x8)-matrix [C A; U G](3) of 64 triplets, where C, A, U, G are letters of the
genetic alphabet. These cyclic permutations lead to such reorganizations of the
matrix form of presentation of the initial genetic Yin-Yang-algebra that arisen
matrices serve as matrix forms of presentations of new Yin-Yang-algebras as
well. They are connected algorithmically with Hadamard matrices. The discovered
existence of a hierarchy of the cyclic changes of types of genetic
Yin-Yang-algebras allows thinking about new algebraic-genetic models of cyclic
processes in inherited biological systems including models of cyclic
metamorphoses of animals. These cycles of changes of the genetic 8-dimensional
algebras and of their 8-dimensional numeric systems have many analogies with
famous facts and doctrines of modern and ancient physiology, medicine, etc.
This viewpoint proposes that the famous idea by Pythagoras (about organization
of natural systems in accordance with harmony of numerical systems) should be
combined with the idea of cyclic changes of Yin-Yang-numeric systems in
considered cases. This second idea reminds of the ancient idea of cyclic
changes in nature. From such algebraic-genetic viewpoint, the notion of
biological time can be considered as a factor of coordinating these
hierarchical ensembles of cyclic changes of types of the genetic
multi-dimensional algebras.


  Biological forces govern essential cellular and molecular processes in all
living organisms. Many cellular forces, e.g. those generated in cyclic
conformational changes of biological machines, have repetitive components.
However, little is known about how proteins process repetitive mechanical
stresses. To obtain first insights into dynamic protein mechanics, we probed
the mechanical stability of single and multimeric ubiquitins perturbed by
periodic forces. Using coarse-grained molecular dynamics simulations, we were
able to model repetitive forces with periods about two orders of magnitude
longer than the relaxation time of folded ubiquitins. We found that even a
small periodic force weakened the protein and shifted its unfolding pathways in
a frequency- and amplitude-dependent manner. Our results also showed that the
dynamic response of even a small protein can be complex with transient
refolding of secondary structures and an increasing importance of local
interactions in asymmetric protein stability. These observations were
qualitatively and quantitatively explained using an energy landscape model and
discussed in the light of dynamic single-molecule measurements and
physiological forces. We believe that our approach and results provide first
steps towards a framework to better understand dynamic protein biomechanics and
biological force generation.


  We study evolutionary game dynamics in a well-mixed populations of finite
size, N. A well-mixed population means that any two individuals are equally
likely to interact. In particular we consider the average abundances of two
strategies, A and B, under mutation and selection. The game dynamical
interaction between the two strategies is given by the 2x2 payoff matrix
[(a,b), (c,d)]. It has previously been shown that A is more abundant than B, if
(N-2)a+Nb>Nc+(N-2)d. This result has been derived for particular stochastic
processes that operate either in the limit of asymptotically small mutation
rates or in the limit of weak selection. Here we show that this result holds in
fact for a wide class of stochastic birth-death processes for arbitrary
mutation rate and for any intensity of selection.


  We study the formation of protein-protein encounter complexes with a Langevin
equation approach that considers direct, steric and thermal forces. As three
model systems with distinctly different properties we consider the pairs
barnase:barstar, cytochrome c:cytochrome c peroxidase and p53:MDM2. In each
case, proteins are modeled either as spherical particles, as dipolar spheres or
as collection of several small beads with one dipole. Spherical reaction
patches are placed on the model proteins according to the known experimental
structures of the protein complexes. In the computer simulations, concentration
is varied by changing box size. Encounter is defined as overlap of the reaction
patches and the corresponding first passage times are recorded together with
the number of unsuccessful contacts before encounter. We find that encounter
frequency scales linearly with protein concentration, thus proving that our
microscopic model results in a well-defined macroscopic encounter rate. The
number of unsuccessful contacts before encounter decreases with increasing
encounter rate and ranges from 20-9000. For all three models, encounter rates
are obtained within one order of magnitude of the experimentally measured
association rates. Electrostatic steering enhances association up to 50-fold.
If diffusional encounter is dominant (p53:MDM2) or similarly important as
electrostatic steering (barnase:barstar), then encounter rate decreases with
decreasing patch radius. More detailed modeling of protein shapes decreases
encounter rates by 5-95 percent. Our study shows how generic principles of
protein-protein association are modulated by molecular features of the systems
under consideration. Moreover it allows us to assess different coarse-graining
strategies for the future modelling of the dynamics of large protein complexes.


  Environmental noise can cause an exponential reduction in the mean time to
extinction (MTE) of an isolated population. We study this effect on an example
of a stochastic birth-death process with rates modulated by a colored Gaussian
noise. A path integral formulation yields a transparent way of evaluating the
MTE and finding the optimal realization of the environmental noise that
determines the most probable path to extinction. The population-size dependence
of the MTE changes from exponential in the absence of the environmental noise
to a power law for a short-correlated noise and to no dependence for
long-correlated noise. We also establish the validity domains of the limits of
white noise and adiabatic noise.


  Stochasticity is both exploited and controlled by cells. Although the
intrinsic stochasticity inherent in biochemistry is relatively well understood,
cellular variation, or 'noise', is predominantly generated by interactions of
the system of interest with other stochastic systems in the cell or its
environment. Such extrinsic fluctuations are nonspecific, affecting many system
components, and have a substantial lifetime, comparable to the cell cycle (they
are 'colored'). Here, we extend the standard stochastic simulation algorithm to
include extrinsic fluctuations. We show that these fluctuations affect mean
protein numbers and intrinsic noise, can speed up typical network response
times, and can explain trends in high-throughput measurements of variation. If
extrinsic fluctuations in two components of the network are correlated, they
may combine constructively (amplifying each other) or destructively
(attenuating each other). Consequently, we predict that incoherent feedforward
loops attenuate stochasticity, while coherent feedforwards amplify it. Our
results demonstrate that both the timescales of extrinsic fluctuations and
their nonspecificity substantially affect the function and performance of
biochemical networks.


  The purpose of this article is to show how the isotropy subgroup of leaf
permutations on binary trees can be used to systematically identify
tree-informative invariants relevant to models of phylogenetic evolution. In
the quartet case, we give an explicit construction of the full set of
representations and describe their properties. We apply these results directly
to Markov invariants, thereby extending previous theoretical results by
systematically identifying linear combinations that vanish for a given quartet.
We also note that the theory is fully generalizable to arbitrary trees and is
equally applicable to the related case of phylogenetic invariants. All results
follow from elementary consideration of the representation theory of finite
groups.


  We study the secondary structure of RNA determined by Watson-Crick pairing
without pseudo-knots using Milnor invariants of links. We focus on the first
non-trivial invariant, which we call the Heisenberg invariant. The Heisenberg
invariant, which is an integer, can be interpreted in terms of the Heisenberg
group as well as in terms of lattice paths.
  We show that the Heisenberg invariant gives a lower bound on the number of
unpaired bases in an RNA secondary structure. We also show that the Heisenberg
invariant can predict \emph{allosteric structures} for RNA. Namely, if the
Heisenberg invariant is large, then there are widely separated local maxima
(i.e., allosteric structures) for the number of Watson-Crick pairs found.


  In this paper we investigate mathematical questions concerning the
reliability (reconstruction accuracy) of Fitch's maximum parsimony algorithm
for reconstructing the ancestral state given a phylogenetic tree and a
character. In particular, we consider the question whether the maximum
parsimony method applied to a subset of taxa can reconstruct the ancestral
state of the root more accurately than when applied to all taxa, and we give an
example showing that this indeed is possible. A surprising feature of our
example is that ignoring a taxon closer to the root improves the reliability of
the method. On the other hand, in the case of the two-state symmetric
substitution model, we answer affirmatively a conjecture of Li, Steel and Zhang
which states that under a molecular clock the probability that the state at a
single taxon is a correct guess of the ancestral state is a lower bound on the
reconstruction accuracy of Fitch's method applied to all taxa.


  The past decade has witnessed the development and success of coarse-grained
network models of proteins for predicting many equilibrium properties related
to collective modes of motion. Curiously, the results are usually robust
towards the different methodologies used for constructing the residue networks
from knowledge of the experimental coordinates. We present a systematical study
of network construction strategies, and their effect on the predicted
properties. The analysis is based on the radial distribution function and the
spectral dimensions of a large set of proteins as well as a newly defined
quantity, the angular distribution function. By partitioning the interactions
into an essential and a residual set, we show that the robustness originates
from a large number of long-distance interactions belonging to the latter.
These residuals have a vanishingly small effect on the force vectors on each
residue. The overall force balance then translates into the Hessian as small
shifts in the slow modes of motion and an invariance of the corresponding
eigenvectors. Implications for the study of biologically relevant properties of
proteins are discussed.


  Understanding how stochastic and non-linear deterministic processes interact
is a major challenge in population dynamics theory. After a short review, we
introduce a stochastic individual-centered particle model to describe the
evolution in continuous time of a population with (continuous) age and trait
structures. The individuals reproduce asexually, age, interact and die. The
'trait' is an individual heritable property (d-dimensional vector) that may
influence birth and death rates and interactions between individuals, and vary
by mutation. In a large population limit, the random process converges to the
solution of a Gurtin-McCamy type PDE. We show that the random model has a long
time behavior that differs from its deterministic limit. However, the results
on the limiting PDE and large deviation techniques \textit{\`a la}
Freidlin-Wentzell provide estimates of the extinction time and a better
understanding of the long time behavior of the stochastic process. This has
applications to the theory of adaptive dynamics used in evolutionary biology.
We present simulations for two biological problems involving life-history trait
evolution when body size is plastic and individual growth is taken into
account.


  It is a standard exercise in mechanical engineering to infer the external
forces and torques on a body from its static shape and known elastic
properties. Here we apply this kind of analysis to distorted double-helical DNA
in complexes with proteins. We extract the local mean forces and torques acting
on each base-pair of bound DNA from high-resolution complex structures. Our
method relies on known elastic potentials and a careful choice of coordinates
of the well-established rigid base-pair model of DNA. The results are robust
with respect to parameter and conformation uncertainty. They reveal the complex
nano-mechanical patterns of interaction between proteins and DNA. Being
non-trivially and non-locally related to observed DNA conformations, base-pair
forces and torques provide a new view on DNA-protein binding that complements
structural analysis.


  The binding of a transcription factor (TF) to a DNA operator site can
initiate or repress the expression of a gene. Computational prediction of sites
recognized by a TF has traditionally relied upon knowledge of several cognate
sites, rather than an ab initio approach. Here, we examine the possibility of
using structure-based energy calculations that require no knowledge of bound
sites but rather start with the structure of a protein-DNA complex. We study
the PurR E. coli TF, and explore to which extent atomistic models of
protein-DNA complexes can be used to distinguish between cognate and
non-cognate DNA sites. Particular emphasis is placed on systematic evaluation
of this approach by comparing its performance with bioinformatic methods, by
testing it against random decoys and sites of homologous TFs. We also examine a
set of experimental mutations in both DNA and the protein. Using our explicit
estimates of energy, we show that the specificity for PurR is dominated by
direct protein-DNA interactions, and weakly influenced by bending of DNA.


  Traditionally evolution is seen as a process where from a pool of possible
variations of a population (e.g. biological species or industrial goods) a few
variations get selected which survive and proliferate, whereas the others
vanish. Survival probability is typically associated with the 'fitness' of a
particular variation. In this paper we argue that the notion of fitness is an a
posteriori concept, in the sense that one can assign higher fitness to species
that survive but one can generally not derive or even measure fitness - or
fitness landscapes - per se. For this reason we think that in a 'physical'
theory of evolution such notions should be avoided. In this spirit, here we
propose a random matrix model of evolution where selection mechanisms are
encoded in the interaction matrices of species. We are able to recover some key
facts of evolution dynamics, such as punctuated equilibrium, i.e. the existence
of intrinsic large extinctions events, and, at the same time, periods of
dramatic diversification, as known e.g. from fossil records. Further, we
comment on two fundamental technical problems of a 'physics of evolution', the
non-closedness of its phase space and the problem of co-evolving boundary
conditions, apparent in all systems subject to evolution.


  The response of ecosystems to perturbations is considered from a
thermodynamic perspective by acknowledging that, as for all macroscopic systems
and processes, the dynamics and stability of ecosystems is subject to definite
thermodynamic law. For open ecosystems, exchanging energy, work, and mass with
the environment, the thermodynamic criteria come from non-equilibrium or
irreversible thermodynamics. For ecosystems during periods in which the
boundary conditions may be considered as being constant, it is shown that
criteria from irreversible thermodynamic theory are sufficient to permit a
quantitative prediction of ecosystem response to perturbation. This framework
is shown to provide a new perspective on the population dynamics of real
ecosystems. The formalism is applied to the problem of the population
oscillations of the southern pine beetle.


  Individual mRNAs are translated by multiple ribosomes that initiate
translation with a few seconds interval. The ribosome speed is codon dependant,
and ribosome queuing has been suggested to explain specific data for
translation of some mRNAs in vivo. By modelling the stochastic translation
process as a traffic problem, we here analyze conditions and consequences of
collisions and queuing. The model allowed us to determine the on-rate (0.8 to
1.1 initiations per sec) and the time (1 sec) the preceding ribosome occludes
initiation for Escherichia coli lacZ mRNA in vivo. We find that ribosome
collisions and queues are inevitable consequences of a stochastic translation
mechanism that reduce the translation efficiency substantially on natural
mRNAs. The cells minimize collisions by having its mRNAs being unstable and by
a highly selected codon usage in the start of the mRNA. The cost of mRNA
breakdown is offset by the concomitant increase in translational efficiency.


  We have recently introduced a multistep extension of the greedy algorithm for
modularity optimization. The extension is based on the idea that merging l
pairs of communities (l>1) at each iteration prevents premature condensation
into few large communities. Here, an empirical formula is presented for the
choice of the step width l that generates partitions with (close to) optimal
modularity for 17 real-world and 1100 computer-generated networks. Furthermore,
an in-depth analysis of the communities of two real-world networks (the
metabolic network of the bacterium E. coli and the graph of coappearing words
in the titles of papers coauthored by Martin Karplus) provides evidence that
the partition obtained by the multistep greedy algorithm is superior to the one
generated by the original greedy algorithm not only with respect to modularity
but also according to objective criteria. In other words, the multistep
extension of the greedy algorithm reduces the danger of getting trapped in
local optima of modularity and generates more reasonable partitions.


  The spatio-temporal dynamics of a population present one of the most
fascinating aspects and challenges for ecological modelling. In this article we
review some simple mathematical models, based on one dimensional
reaction-diffusion-advection equations, for the growth of a population on a
heterogeneous habitat. Considering a number of models of increasing complexity
we investigate the often contrary roles of advection and diffusion for the
persistence of the population. When it is possible we demonstrate basic
mathematical techniques and give the critical conditions providing the survival
of a population, in simple systems and in more complex resource-consumer models
which describe the dynamics of phytoplankton in a water column.


  Although most networks in nature exhibit complex topology the origins of such
complexity remains unclear. We introduce a model of a growing network of
interacting agents in which each new agent's membership to the network is
determined by the agent's effect on the network's global stability. It is shown
that out of this stability constraint, scale free networks emerges in a self
organized manner, offering an explanation for the ubiquity of complex
topological properties observed in biological networks.


  After Laskar, the Lyapunov time in the solar system is about five millions
years (5.000.000 [years]). On the other hand, after Kimura, the evolutionary
(phenotypic) rate, for hominids, is 1/5.000.000 [1/years]. Why are these two
quantities so closely related? In this work, following a proposition by
Finlayson and Hutchings et al, I found an inequality, which relates Lyapunov
time and evolution rate. This inequality fits well with some known cases in
biological evolution.


  A general model for the early recognition and colocalization of homologous
DNA sequences is proposed. We show, on a thermodynamic ground, how the distance
between two homologous DNA sequences is spontaneously regulated by the
concentration and affinity of diffusible mediators binding them, which act as a
switch between two phases corresponding to independence or colocalization of
pairing regions.


  A single protein molecule is regarded as a contact network of amino-acid
residues. Some studies have indicated that this network is a small world
network (SWN), while other results have implied that this is a fractal network
(FN). However, SWN and FN are essentially different in the dependence of the
shortest path length on the number of nodes. In this paper, we investigate this
dependence in the residue contact networks of proteins in native structures,
and show that the networks are not SWN but FN. FN is generally characterized by
several dimensions. Among them, we focus on three dimensions; the network
topological dimension $D_c$, the fractal dimension $D_f$, and the spectral
dimension $D_s$. We find that proteins universally yield $D_c \approx 1.9$,
$D_f \approx 2.5$ and $Ds \approx 1.3$. These values are in surprisingly good
coincidence with those in three dimensional critical percolation cluster. Hence
the residue contact networks in the protein native structures belong to the
universality class of three dimensional percolation cluster. The criticality is
relevant to the ambivalent nature of the protein native structures, i.e., the
coexistence of stability and instability, both of which are necessary for a
protein to function as a molecular machine or an allosteric enzyme.


  Until recently, the only examples of inverted biomass pyramids have been in
freshwater and marine planktonic communities. In 2002 and 2008 investigators
documented inverted biomass pyramids for nearly pristine coral reef ecosystems
within the NW Hawaiian islands and the Line Islands, where apex predator
abundance comprises up to 85% of the fish biomass. We build a new refuge based
predator-prey model to study the fish biomass structure at coral reefs and
investigate the effect of fishing on biomass pyramids. Utilizing realistic life
history parameters of coral reef fish, our model exhibits a stable inverted
biomass pyramid. Since the predators and prey are not well mixed, our model
does not incorporate homogeneous mixing and the inverted biomass pyramid is a
consequence of the refuge. Understanding predator-prey dynamics in nearly
pristine conditions provides a more realistic historical framework for
comparison with fished reefs. Finally, we show that fishing transforms the
inverted biomass pyramid to be bottom heavy.


  Recent microbial experiments suggest that enhanced genetic drift at the
frontier of a two-dimensional range expansion can cause genetic sectoring
patterns with fractal domain boundaries. Here, we propose and analyze a simple
model of asexual biological evolution at expanding frontiers to explain these
neutral patterns and predict the effect of natural selection. Our model
attributes the observed gradual decrease in the number of sectors at the
leading edge to an unbiased random walk of sector boundaries. Natural selection
introduces a deterministic bias in the wandering of domain boundaries that
renders beneficial mutations more likely to escape genetic drift and become
established in a sector. We find that the opening angle of those sectors and
the rate at which they become established depend sensitively on the selective
advantage of the mutants. Deleterious mutations, on the other hand, are not
able to establish a sector permanently. They can, however, temporarily "surf"
on the population front, and thereby reach unusual high frequencies. As a
consequence, expanding frontiers are susceptible to deleterious mutations as
revealed by the high fraction of mutants at mutation-selection balance.
Numerically, we also determine the condition at which the wild type is lost in
favor of deleterious mutants (genetic meltdown) at a growing front. Our
prediction for this error threshold differs qualitatively from existing
well-mixed theories, and sets tight constraints on sustainable mutation rates
for populations that undergo frequent range expansions.


  We investigate the mechanical unfolding of the tenth type III domain from
fibronectin, FnIII10, both at constant force and at constant pulling velocity,
by all-atom Monte Carlo simulations. We observe both apparent two-state
unfolding and several unfolding pathways involving one of three major, mutually
exclusive intermediate states. All the three major intermediates lack two of
seven native beta-strands, and share a quite similar extension. The unfolding
behavior is found to depend strongly on the pulling conditions. In particular,
we observe large variations in the relative frequencies of occurrence for the
intermediates. At low constant force or low constant velocity, all the three
major intermediates occur with a significant frequency. At high constant force
or high constant velocity, one of them, with the N- and C-terminal beta-strands
detached, dominates over the other two. Using the extended Jarzynski equality,
we also estimate the equilibrium free-energy landscape, calculated as a
function of chain extension. The application of a constant pulling force leads
to a free-energy profile with three major local minima. Two of these correspond
to the native and fully unfolded states, respectively, whereas the third one
can be associated with the major unfolding intermediates.


  Regulatory gene networks contain generic modules like those involving
feedback loops, which are essential for the regulation of many biological
functions. We consider a class of self-regulated genes which are the building
blocks of many regulatory gene networks, and study the steady state
distributions of the associated Gillespie algorithm by providing efficient
numerical algorithms. We also study a regulatory gene network of interest in
synthetic biology and in gene therapy, using mean-field models with time
delays. Convergence of the related time-nonhomogeneous Markov chain is
established for a class of linear catalytic networks with feedback loops


  We address the problem of diversification in religions by studying selection
on cultural memes that colonize humans hosts. In analogy to studying the
evolution of pathogens or symbionts colonizing animal hosts, we use models for
host-pathogen dynamics known from theoretical epidemiology. In these models,
religious memes colonize individual humans. Rates of transmission of memes
between humans, i.e., transmission of cultural content, and rates of loss of
memes (loss of faith) are determined by the phenotype of the cultural memes,
and by interactions between hosts carrying different memes. In particular,
based on the notion that religion can lead to oppression of lower classes once
a religious society has reached a certain size, we assume that the rate of loss
increases as the number of humans colonized by a particular meme phenotype
increases. This generates frequency-dependent selection on cultural memes, and
we use evolutionary theory to show that this frequency dependence can generate
the emergence of coexisting clusters of different meme types. The different
clusters correspond to different religions, and hence our model describes the
emergence of distinct descendent religions from single ancestral religions.


  The transmission of vector infectious diseases, which produces complex
spatiotemporal patterns, is analyzed by a periodically forced two-dimensional
cellular automata model. The system, which comprises three population levels,
is introduced to describe complex features of the dynamics of the vector
transmitted dengue epidemics, known to be very sensitive to seasonal variables.
The three coupled levels represent the human, the adult and immature vector
populations. The dynamics includes external seasonality forcing (rainfall
intensity data), human and mosquito mobility, and vector control effects. The
model parameters, even if bounded to well defined intervals obtained from
reported data, can be selected to reproduce specific epidemic outbursts. In the
current study, explicit results are obtained by comparison with actual data
retrieved from the time-series of dengue epidemics in two cities in Brazil. The
results show fluctuations that are not captured by mean-field models. It also
reveals the qualitative behavior of the spatiotemporal patterns of the
epidemics. In the extreme situation of absence of external periodic drive, the
model predicts completely distinct long time evolution. The model is robust in
the sense that it is able to reproduce the time series of dengue epidemics of
different cities, provided the forcing term takes into account the local
rainfall modulation. Finally, the dependence between epidemics threshold and
vector control undergoes a transition from power law to stretched exponential
behavior due to human mobility effect.


  Molecular dynamics studies within a coarse-grained structure based model were
used on two similar proteins belonging to the transcarbamylase family to probe
the effects in the native structure of a knot. The first protein,
N-acetylornithine transcarbamylase, contains no knot whereas human ormithine
transcarbamylase contains a trefoil knot located deep within the sequence. In
addition, we also analyzed a modified transferase with the knot removed by the
appropriate change of a knot-making crossing of the protein chain. The studies
of thermally- and mechanically-induced unfolding processes suggest a larger
intrinsic stability of the protein with the knot.


  We present a novel method for identifying a biochemical reaction network
based on multiple sets of estimated reaction rates in the corresponding
reaction rate equations arriving from various (possibly different) experiments.
The current method, unlike some of the graphical approaches proposed in the
literature, uses the values of the experimental measurements only relative to
the geometry of the biochemical reactions under the assumption that the
underlying reaction network is the same for all the experiments.
  The proposed approach utilizes algebraic statistical methods in order to
parametrize the set of possible reactions so as to identify the most likely
network structure, and is easily scalable to very complicated biochemical
systems involving a large number of species and reactions. The method is
illustrated with a numerical example of a hypothetical network arising form a
"mass transfer"-type model.


  Recent studies of cluster distribution in various ecosystems revealed Pareto
statistics for the size of spatial colonies. These results were supported by
cellular automata simulations that yield robust criticality for endogenous
pattern formation based on positive feedback. We show that this self-organized
criticality is a manifestation of the law of proportionate effect, first
discovered in the context of business firm size. Mapping the stochastic model
to a Markov birth-death process, the transition rates are shown to scale
linearly with cluster size. This mapping provides a connection between patch
statistics and the dynamics of the ecosystem; the "first passage time" for
different colonies emerges as a powerful tool that discriminates between
endogenous and exogenous clustering mechanisms. Imminent catastrophic shifts
(like desertification) manifest themselves in a drastic change of the stability
properties of spatial colonies, as the chance of a cluster to disappear depends
logarithmically, rather than linearly, on its size.


  The Robinson-Foulds (RF) distance is by far the most widely used measure of
dissimilarity between trees. Although the distribution of these distances has
been investigated for twenty years, an algorithm that is explicitly polynomial
time has yet to be described for computing this distribution (which is also the
distribution of trees around a given tree under the popular Robinson-Foulds
metric). In this paper we derive a polynomial-time algorithm for this
distribution. We show how the distribution can be approximated by a Poisson
distribution determined by the proportion of leaves that lie in `cherries' of
the given tree. We also describe how our results can be used to derive
normalization constants that are required in a recently-proposed maximum
likelihood approach to supertree construction.


  We analyze the spatial distributions of two groups of benthic foraminifera
(Adelosina spp. + Quinqueloculina spp. and Elphidium spp.), along Sicilian
coast, and their correlation with six different heavy metals, responsible for
the pollution. Samples were collected inside the Gulf of Palermo, which has a
high level of pollution due to heavy metals, and along the coast of Lampedusa
island (Sicily Channel, Southern Mediterranean), which is characterized by
unpolluted sea waters. Because of the environmental pollution we find: (i) an
anticorrelated spatial behaviour between the two groups of benthic foraminifera
analyzed; (ii) an anticorrelated (correlated) spatial behaviour between the
first (second) group of benthic foraminifera with metal concentrations; (iii)
an almost uncorrelated spatial behaviour between low concentrations of metals
and the first group of foraminifera in clean sea water sites. We introduce a
two-species model based on the generalized Lotka-Volterra equations in the
presence of a multiplicative noise, which models the interaction between
species and environmental pollution due to the presence in top-soft sediments
of heavy metals. The interaction coefficients between the two species are kept
constant with values in the coexistence regime. Using proper values for the
initial conditions and the model parameters, we find for the two species a
theoretical spatial distribution behaviour in a good agreement with the data
obtained from the 63 sites analyzed in our study.


  Molecular motors convert chemical energy into mechanical work while operating
in an environment dominated by Brownian motion. The aim of this paper is to
explore the flow of energy between the molecular motors and its surroundings,
in particular, its efficiency. Based on the Fokker-Planck equation with either
$N$ or infinite chemical states, we find that the energy efficiency of the
molecular motors, whether the Stokes efficiency or the usual thermodynamic
efficiency, is strictly bounded by 1, because of the dissipation of the energy
in both the overdamped surroundings and in the process of the chemical
reaction.


  Main aim of this topical issue is to report recent advances in noisy
nonequilibrium processes useful to describe the dynamics of ecological systems
and to address the mechanisms of spatio-temporal pattern formation in ecology
both from the experimental and theoretical points of view. This is in order to
understand the dynamical behaviour of ecological complex systems through the
interplay between nonlinearity, noise, random and periodic environmental
interactions. Discovering the microscopic rules and the local interactions
which lead to the emergence of specific global patterns or global dynamical
behaviour and the noises role in the nonlinear dynamics is an important, key
aspect to understand and then to model ecological complex systems.


  We consider the dynamics of cargo driven by a collection of interacting
molecular motors in the context of an asymmetric simple exclusion processes
(ASEP). The model is formulated to account for i) excluded volume interactions,
ii) the observed asymmetry of the stochastic movement of individual motors and
iii) interactions between motors and cargo. Items (i) and (ii) form the basis
of ASEP models and have already been considered in the literature to study the
behavior of motor density profile [Parmeggiani 03]. Item (iii) is new. It is
introduced here as an attempt to describe explicitly the dependence of cargo
movement on the dynamics of motors. The steady-state solutions of the model
indicate that the system undergoes a phase transition of condensation type as
the motor density varies. We study the consequences of this transition to the
properties of cargo velocity.


  Comprehensive knowledge of protein-ligand interactions should provide a
useful basis for annotating protein functions, studying protein evolution,
engineering enzymatic activity, and designing drugs. To investigate the
diversity and universality of ligand binding sites in protein structures, we
conducted the all-against-all atomic-level structural comparison of over
180,000 ligand binding sites found in all the known structures in the Protein
Data Bank by using a recently developed database search and alignment
algorithm. By applying a hybrid top-down-bottom-up clustering analysis to the
comparison results, we determined approximately 3000 well-defined structural
motifs of ligand binding sites. Apart from a handful of exceptions, most
structural motifs were found to be confined within single families or
superfamilies, and to be associated with particular ligands. Furthermore, we
analyzed the components of the similarity network and enumerated more than 4000
pairs of ligand binding sites that were shared across different protein folds.


  We discuss effects of the stepping kinetics of molecular motors on their
traffic behavior on crowded filaments using a simple two-state chemomechanical
cycle. While the general traffic behavior is quite robust with respect to the
detailed kinetics of the step, a few observable parameters exhibit a strong
dependence on these parameters. Most strikingly, the effective unbinding rate
of the motors may both increase and decrease with increasing traffic density,
depending on the details of the motor step. Likewise the run length either
exhibits a strong decrease or almost no dependence on the traffic density. We
compare our theoretical results with recent experimental observations on motor
traffic.


  We study a simple reaction-diffusion population model [proposed by A. Windus
and H. J. Jensen, J. Phys. A: Math. Theor. 40, 2287 (2007)] on scale-free
networks. In the case of fully random diffusion, the network topology cannot
affect the critical death rate, whereas the heterogeneous connectivity can
cause smaller steady population density and critical population density. In the
case of modified diffusion, we obtain a larger critical death rate and steady
population density, at the meanwhile, lower critical population density, which
is good for the survival of species. The results were obtained using a
mean-field-like framework and were confirmed by computer simulations.


  Aim: To assess the frequency of secondary transmissions of primary pneumonic
plague relative to the onset of fever. Methods: A simple backcalculation method
was employed to estimate the frequency of secondary transmissions relative to
disease-age. A likelihood-based procedure was taken using observed
distributions of the serial interval (n = 177) and incubation period (n = 126).
Furthermore, an extended model was developed to account for the survival
probability of cases. Results: The simple backcalculation suggested that 31.0%
(95% confidence intervals (CI): 11.6, 50.4) and 28.0 % (95% CI: 10.2, 45.8) of
the total number of secondary transmissions had occurred at second and third
days of the disease, respectively, and more than four-fifths of the secondary
transmission occurred before the end of third day of disease. The
survivorship-adjusted frequency of secondary transmissions was obtained,
demonstrating that the infectiousness in later stages of illness was not
insignificant and indicates that the obtained frequencies were likely biased on
underlying factors including isolation measures. Conclusion: The simple
exercise suggests a need to implement countermeasures during pre-clinical stage
or immediately after onset. Further information is needed to elucidate the
finer details of the disease-age specific infectiousness.


  Protein electrostatic states have been demonstrated to play crucial roles in
catalysis, ligand binding, protein stability, and in the modulation of
allosteric effects. Electrostatic states are demonstrated to appear conserved
among DEAD-box motifs and evidence is presented that the structural changes
that occur to DEAD box proteins upon ligand binding alter the DEAD-box motif
electrostatics in a way the facilitates the catalytic role of the DEAD-box
glutatmate.


  ESPSim is an open source JAVA program that enables the comparisons of protein
electrostatic potential maps via the computation of an electrostatic similarity
measure. This program has been utilized to demonstrate a high degree of
electrostatic similarity among the potential maps of lysozyme proteins,
suggesting that protein electrostatic states are conserved within lysozyme
proteins. ESPSim is freely available under the AGPL License from
http://www.bioinformatics.org/project/?group_id=830


  Which factors govern the evolution of mutation rates and emergence of
species? Here, we address this question using a first principles model of life
where population dynamics of asexual organisms is coupled to molecular
properties and interactions of proteins encoded in their genomes. Simulating
evolution of populations, we found that fitness increases in punctuated steps
via epistatic events, leading to formation of stable and functionally
interacting proteins. At low mutation rates, species - populations of organisms
with identical genotypes - form, while at higher mutation rates, species are
lost through delocalization in sequence space without an apparent loss of
fitness. However, when mutation rate was a selectable trait, the population
initially maintained high mutation rate until a high fitness level is reached,
after which organisms with low mutation rates are gradually selected, with the
population eventually reaching mutation rates comparable to those of modern
DNA-based organisms. These results provide microscopic insights into the
dynamic fitness landscape of asexual populations of unicellular organisms.


  Rigidity analysis using the "pebble game" can usefully be applied to protein
crystal structures to obtain information on protein folding, assembly and the
structure-function relationship. However, previous work using this technique
has not made clear how sensitive rigidity analysis is to small structural
variations. We present a comparative study in which rigidity analysis is
applied to multiple structures, derived from different organisms and different
conditions of crystallisation, for each of several different proteins. We find
that rigidity analysis is best used as a comparative tool to highlight the
effects of structural variation. Our use of multiple protein structures brings
out a previously unnoticed peculiarity in the rigidity of trypsin.


  Protein electrostatics have been demonstrated to play a vital role in protein
functionality, with many functionally important amino acid residues exhibiting
an electrostatic state that is altered from that of a normal amino acid
residue. Residues with altered electrostatic states can be identified by the
presence of a pKa value that is perturbed by 2 or more pK units, and such
residues have been demonstrated to play critical roles in catalysis, ligand
binding, and protein stability. Within the HCV helicase and polymerase, as well
as the HIV reverse transcriptase, highly conserved regions were demonstrated to
possess a greater number and magnitude of perturbations than lesser conserved
regions, suggesting that there is an interrelationship present between protein
electrostatics and evolution.


  Protein dynamics in cells may be different from that in dilute solutions in
vitro since the environment in cells is highly concentrated with other
macromolecules. This volume exclusion due to macromolecular crowding is
predicted to affect both equilibrium and kinetic processes involving protein
conformational changes. To quantify macromolecular crowding effects on protein
folding mechanisms, here we have investigated the folding energy landscape of
an alpha/beta protein, apoflavodoxin, in the presence of inert macromolecular
crowding agents using in silico and in vitro approaches. By coarse-grained
molecular simulations and topology-based potential interactions, we probed the
effects of increased volume fraction of crowding agents (phi_c) as well as of
crowding agent geometry (sphere or spherocylinder) at high phi_c. Parallel
kinetic folding experiments with purified Desulfovibro desulfuricans
apoflavodoxin in vitro were performed in the presence of Ficoll (sphere) and
Dextran (spherocylinder) synthetic crowding agents. In conclusion, we have
identified in silico crowding conditions that best enhance protein stability
and discovered that upon manipulation of the crowding conditions, folding
routes experiencing topological frustrations can be either enhanced or
relieved. The test-tube experiments confirmed that apoflavodoxin's
time-resolved folding path is modulated by crowding agent geometry. We propose
that macromolecular crowding effects may be a tool for manipulation of protein
folding and function in living cells.


  Noise in the expression of a gene produces fluctuations in the concentration
of the gene product. These fluctuations can interfere with optimal function or
can be exploited to generate beneficial diversity between cells; gene
expression noise is therefore expected to be subject to evolutionary pressure.
Shifts between modes of high and low rates of transcription initiation at a
promoter appear to contribute to this noise both in eukaryotes and prokaryotes.
However, models invoked for eukaryotic promoter noise such as stable activation
scaffolds or persistent nucleosome alterations seem unlikely to apply to
prokaryotic promoters. We consider the relative importance of the steps
required for transcription initiation. The 3-step transcription initiation
model of McClure is extended into a mathematical model that can be used to
predict consequences of additional promoter properties. We show in principle
that the transcriptional bursting observed at an E. coli promoter by Golding et
al. (2005) can be explained by stimulation of initiation by the negative
supercoiling behind a transcribing RNA polymerase (RNAP) or by the formation of
moribund or dead-end RNAP-promoter complexes. Both mechanisms are tunable by
the alteration of promoter kinetics and therefore allow the optimization of
promoter mediated noise.


  The Tethered Particle Motion (TPM) method has been used to observe and
characterize a variety of protein-DNA interactions including DNA looping and
transcription. TPM experiments exploit the Brownian motion of a DNA-tethered
bead to probe biologically relevant conformational changes of the tether. In
these experiments, a change in the extent of the bead's random motion is used
as a reporter of the underlying macromolecular dynamics and is often deemed
sufficient for TPM analysis. However, a complete understanding of how the
motion depends on the physical properties of the tethered particle complex
would permit more quantitative and accurate evaluation of TPM data. For
instance, such understanding can help extract details about a looped complex
geometry (or multiple coexisting geometries) from TPM data. To better
characterize the measurement capabilities of TPM experiments involving DNA
tethers, we have carried out a detailed calibration of TPM magnitude as a
function of DNA length and particle size. We also explore how experimental
parameters such as acquisition time and exposure time affect the apparent
motion of the tethered particle. We vary the DNA length from 200bp to 2.6kbp
and consider particle diameters of 200, 490 and 970nm. We also present a
systematic comparison between measured particle excursions and theoretical
expectations, which helps clarify both the experiments and models of DNA
conformation.


  Certain short polycations, such as TAT and polyarginine, rapidly pass through
the plasma membranes of mammalian cells by an unknown mechanism called
transduction as well as by endocytosis and macropinocytosis. These
cell-penetrating peptides (CPPs) promise to be medically useful when fused to
biologically active peptides. I offer a simple model in which one or more CPPs
and the phosphatidylserines of the inner leaflet form a kind of capacitor with
a voltage in excess of 180 mV, high enough to create a molecular electropore.
The model is consistent with an empirical upper limit on the cargo peptide of
40--60 amino acids and with experimental data on how the transduction of a
polyarginine-fluorophore into mouse C2C12 myoblasts depends on the number of
arginines in the CPP and on the CPP concentration. The model makes three
testable predictions.


  We investigate the dynamics of two interacting diffusing particles in an
infinite effectively one dimensional system; the particles interact through a
step-like potential of width b and height phi_0 and are allowed to pass one
another. By solving the corresponding 2+1-variate Fokker-Planck equation an
exact result for the two particle conditional probability density function
(PDF) is obtained for arbitrary initial particle positions. From the
two-particle PDF we obtain the overtake probability, i.e. the probability that
the two particles has exchanged positions at time t compared to the initial
configuration. In addition, we calculate the trapping probability, i.e. the
probability that the two particles are trapped close to each other (within the
barrier width b) at time t, which is mainly of interest for an attractive
potential, phi_0<0. We also investigate the tagged particle PDF, relevant for
describing the dynamics of one particle which is fluorescently labeled. Our
analytic results are in excellent agreement with the results of stochastic
simulations, which are performed using the Gillespie algorithm.


  Due to inherent complexity active transport presents a landmark hurdle for
oral absorption properties prediction. We present a novel approach
carrier-mediated drug absorption parameters calculation based on entirely
different paradigm than QSPR. We capitalize on recently emerged ideas that
molecule activities against a large protein set can be used for prediction of
biological effects and performed a large scale numerical docking of drug-like
compounds to a large diversified set of proteins. As a result we identified for
the first time a protein, binding to which correlates well with the intestinal
permeability of many actively absorbed compounds. Although the protein is not a
transporter, we speculate that it has the binding site force field similar to
that of an important intestinal transporter. The observation helped us to
improve the passive absorption model by adding non-liner flux associated with
the transporting protein to obtain a quantitative model of active transport.
This study demonstrates that binding data to a sufficiently representative set
of proteins can serve as a basis for active absorption prediction for a given
compound.


  Genetic robustness, the preservation of an optimal phenotype in the face of
mutations, is critical to the understanding of evolution as phenotypically
expressed genetic variation is the fuel of natural selection. The origin of
genetic robustness, whether it evolves directly by natural selection or it is a
correlated byproduct of other phenotypic traits, is, however, unresolved.
Examining microRNA (miRNA) genes of several eukaryotic species, Borenstein and
Ruppin (Borenstein et al. 2006, PNAS 103: 6593), showed that the structure of
miRNA precursor stem-loops exhibits significantly increased mutational
robustness in comparison with a sample of random RNA sequences with the same
stem-loop structure. The observed robustness was found to be uncorrelated with
traditional measures of environmental robustness -- implying that miRNA
sequences show evidence of the direct evolution of genetic robustness. These
findings are surprising as theoretical results indicate that the direct
evolution of robustness requires high mutation rates and/or large effective
population sizes only found among RNA viruses, not multicellular eukaryotes.
Introducing a novel measure of environmental robustness based on the
equilibrium thermodynamic ensemble of secondary structures of the miRNA
precursor sequences we demonstrate that the biophysics of RNA folding, induces
a high level of correlation between genetic (mutational) and environmental
(thermodynamic) robustness, as expected from the theory of plastogenetic
congruence introduced by Ancel and Fontana (Ancel et al. 2000, J. Exp. Zool.
288: 242). In light of theoretical considerations we believe that this
correlation strongly suggests that genetic robustness observed in miRNA
sequences is the byproduct of selection for environmental robustness.


  We study the population dynamics of lytic viruses which replicate slowly in
dividing host cells within an organism or cell culture, and find a range of
viral replication rates that allows viruses to persist, avoiding extinction of
host cells or dilution of viruses at too rapid or too slow viral replication.
For the within-host competition between multiple viral strains, a strain with a
"stable" replication rate could outcompete another strain with a higher or
lower replication rate, therefore natural selection of viruses stabilizes the
viral persistence. However, when strains with higher and lower than the
"stable" value replication rates are both present, competition between strains
does not result in dominance of one strain, but in their coexistence.


  We show how to construct a reduced description of interacting genes in noisy,
small regulatory networks using coupled binary "spin" variables. Treating both
the protein number and gene expression state variables stochastically and on
equal footing we propose a mapping which connects the molecular level
description of networks to the binary representation. We construct a phase
diagram indicating when genes can be considered to be independent and when the
coupling between them cannot be neglected leading to synchrony or correlations.
We find that an appropriately mapped boolean description reproduces the
probabilities of gene expression states of the full stochastic system very well
and can be transfered to examples of self-regulatory systems with a larger
number of gene copies.


  We introduce a generalization of the parallel, or Crow-Kimura, and Eigen
models of molecular evolution to represent the exchange of genetic information
between individuals in a population. We study the effect of different schemes
of genetic recombination on the steady-state mean fitness and distribution of
individuals in the population, through an analytic field theoretic mapping. We
investigate both horizontal gene transfer from a population and recombination
between pairs of individuals. Somewhat surprisingly, these nonlinear
generalizations of quasi-species theory to modern biology are analytically
solvable. For two-parent recombination, we find two selected phases, one of
which is spectrally rigid. We present exact analytical formulas for the
equilibrium mean fitness of the population, in terms of a maximum principle,
which are generally applicable to any permutation invariant replication rate
function. For smooth fitness landscapes, we show that when positive epistatic
interactions are present, recombination or horizontal gene transfer introduces
a mild load against selection. Conversely, if the fitness landscape exhibits
negative epistasis, horizontal gene transfer or recombination introduce an
advantage by enhancing selection towards the fittest genotypes. These results
prove that the mutational deterministic hypothesis holds for quasi-species
models. For the discontinuous single sharp peak fitness landscape, we show that
horizontal gene transfer has no effect on the fitness, while recombination
decreases the fitness, for both the parallel and the Eigen models. We present
numerical and analytical results as well as phase diagrams for the different
cases.


  Sequence comparison and alignment has had an enormous impact on our
understanding of evolution, biology, and disease. Comparison and alignment of
biological networks will likely have a similar impact. Existing network
alignments use information external to the networks, such as sequence, because
no good algorithm for purely topological alignment has yet been devised. In
this paper, we present a novel algorithm based solely on network topology, that
can be used to align any two networks. We apply it to biological networks to
produce by far the most complete topological alignments of biological networks
to date. We demonstrate that both species phylogeny and detailed biological
function of individual proteins can be extracted from our alignments.
Topology-based alignments have the potential to provide a completely new,
independent source of phylogenetic information. Our alignment of the
protein-protein interaction networks of two very different species--yeast and
human--indicate that even distant species share a surprising amount of network
topology with each other, suggesting broad similarities in internal cellular
wiring across all life on Earth.


  The identification of key populations shaping the structure and connectivity
of metapopulation systems is a major challenge in population ecology. The use
of molecular markers in the theoretical framework of population genetics has
allowed great advances in this field, but the prime question of quantifying the
role of each population in the system remains unresolved. Furthermore, the use
and interpretation of classical methods are still bounded by the need for a
priori information and underlying assumptions that are seldom respected in
natural systems. Network theory was applied to map the genetic structure in a
metapopulation system using microsatellite data from populations of a
threatened seagrass, Posidonia oceanica, across its whole geographical range.
The network approach, free from a priori assumptions and of usual underlying
hypothesis required for the interpretation of classical analysis, allows both
the straightforward characterization of hierarchical population structure and
the detection of populations acting as hubs critical for relaying gene flow or
sustaining the metapopulation system. This development opens major perspectives
in ecology and evolution in general, particularly in areas such as conservation
biology and epidemiology, where targeting specific populations is crucial.


  Many processes and models --in biological, physical, social, and other
contexts-- produce trees whose depth scales logarithmically with the number of
leaves. Phylogenetic trees, describing the evolutionary relationships between
biological species, are examples of trees for which such scaling is not
observed. With this motivation, we analyze numerically two branching models
leading to non-logarithmic scaling of the depth with the number of leaves. For
Ford's alpha model, although a power-law scaling of the depth with tree size
was established analytically, our numerical results illustrate that the
asymptotic regime is approached only at very large tree sizes. We introduce
here a new model, the activity model, showing analytically and numerically that
it also displays a power-law scaling of the depth with tree size at a critical
parameter value.


  The force generated between actin and myosin acts predominantly along the
direction of the actin filament, resulting in relative sliding of the thick and
thin filaments in muscle or transport of myosin cargos along actin tracks.
Previous studies have also detected lateral forces or torques that are
generated between actin and myosin, but the origin and biological role of these
sideways forces is not known. Here we adapt an actin gliding filament assay in
order to measure the rotation of an actin filament about its axis (twirling) as
it is translocated by myosin. We quantify the rotation by determining the
orientation of sparsely incorporated rhodamine-labeled actin monomers, using
polarized total internal reflection (polTIRF) microscopy. In order to determine
the handedness of the filament rotation, linear incident polarizations in
between the standard s- and p-polarizations were generated, decreasing the
ambiguity of our probe orientation measurement four-fold. We found that whole
myosin II and myosin V both twirl actin with a relatively long (micron),
left-handed pitch that is insensitive to myosin concentration, filament length
and filament velocity.


  In this paper we propose a microscopic model to study the polymerization of
microtubules (MTs). Starting from fundamental reactions during MT's assembly
and disassembly processes, we systematically derive a nonlinear system of
equations that determines the dynamics of microtubules in 3D. %coexistence with
tubulin dimers in a solution. We found that the dynamics of a MT is
mathematically expressed via a cubic-quintic nonlinear Schrodinger (NLS)
equation. Interestingly, the generic 3D solution of the NLS equation exhibits
linear growing and shortening in time as well as temporal fluctuations about a
mean value which are qualitatively similar to the dynamic instability of MTs
observed experimentally. By solving equations numerically, we have found
spatio-temporal patterns consistent with experimental observations.


  Biologists rely heavily on the language of information, coding, and
transmission that is commonplace in the field of information theory as
developed by Claude Shannon, but there is open debate about whether such
language is anything more than facile metaphor. Philosophers of biology have
argued that when biologists talk about information in genes and in evolution,
they are not talking about the sort of information that Shannon's theory
addresses. First, philosophers have suggested that Shannon theory is only
useful for developing a shallow notion of correlation, the so-called "causal
sense" of information. Second they typically argue that in genetics and
evolutionary biology, information language is used in a "semantic sense,"
whereas semantics are deliberately omitted from Shannon theory. Neither
critique is well-founded. Here we propose an alternative to the causal and
semantic senses of information: a transmission sense of information, in which
an object X conveys information if the function of X is to reduce, by virtue of
its sequence properties, uncertainty on the part of an agent who observes X.
The transmission sense not only captures much of what biologists intend when
they talk about information in genes, but also brings Shannon's theory back to
the fore. By taking the viewpoint of a communications engineer and focusing on
the decision problem of how information is to be packaged for transport, this
approach resolves several problems that have plagued the information concept in
biology, and highlights a number of important features of the way that
information is encoded, stored, and transmitted as genetic sequence.


  We measured Young's modulus at temperatures ranging from 20 to 100 ^{\circ}$C
for a collagen fibril taken from rat's tendon. The hydration change under
heating and the damping decrement were measured as well. At physiological
temperatures $25-45^{\circ}$C Young's modulus decreases, which can be
interpreted as instability of collagen. For temperatures between
$45-80^{\circ}$C Young's modulus first stabilizes and then increases with
decreasing the temperature. The hydrated water content and the damping
decrement have strong maxima in the interval $70-80^{\circ}$C indicating on
complex inter-molecular structural changes in the fibril. All these effects
disappear after heat-denaturating the sample at $120^\circ$C. Our main result
is a five-stage mechanism by which the instability of a single collagen at
physiological temperatures is compensated by the interaction between collagen
molecules within the fibril.


  Background: Alignment of biological sequences such as DNA, RNA or proteins is
one of the most widely used tools in computational bioscience. All existing
alignment algorithms rely on heuristic scoring schemes based on biological
expertise. Therefore, these algorithms do not provide model independent and
objective measures for how similar two (or more) sequences actually are.
Although information theory provides such a similarity measure -- the mutual
information (MI) -- previous attempts to connect sequence alignment and
information theory have not produced realistic estimates for the MI from a
given alignment.
  Results: Here we describe a simple and flexible approach to get robust
estimates of MI from {\it global} alignments. For mammalian mitochondrial DNA,
our approach gives pairwise MI estimates for commonly used global alignment
algorithms that are strikingly close to estimates obtained by an entirely
unrelated approach -- concatenating and zipping the sequences.
  Conclusions: This remarkable consistency may help establish MI as a reliable
tool for evaluating the quality of global alignments, judging the relative
merits of different alignment algorithms, and estimating the significance of
specific alignments. We expect that our approach can be extended to establish
further connections between information theory and sequence alignment,
including applications to local and multiple alignment procedures.


  We study the global stability of a class of models for in-vivo virus
dynamics, that take into account the CTL immune response and display antigenic
variation. This class includes a number of models that have been extensively
used to model HIV dynamics. We show that models in this class are globally
asymptotically stable, under mild hypothesis, by using appropriate Lyapunov
functions. We also characterise the stable equilibrium points for the entire
biologically relevant parameter range. As a byproduct, we are able to determine
what is the diversity of the persistent strains.


  In this paper, a model for understanding the effects of selection using
systems- level computational approaches is introduced. A number of concepts and
principles essential for understanding the motivation for constructing the
model will be introduced first. This will be followed by a description of
parameters, measurements, and graphical representations used in the model. Four
possible outcomes for this model are then introduced and described. In
addition, the relationship of relative fitness to selection is described.
Finally, the consequences and potential lessons learned from the model are
discussed.


  We show that the simplest stochastic epidemiological models with spatial
correlations exhibit two types of oscillatory behaviour in the endemic phase.
In a large parameter range, the oscillations are due to resonant amplification
of stochastic fluctuations, a general mechanism first reported for
predator-prey dynamics. In a narrow range of parameters that includes many
infectious diseases which confer long lasting immunity the oscillations persist
for infinite populations. This effect is apparent in simulations of the
stochastic process in systems of variable size, and can be understood from the
phase diagram of the deterministic pair approximation equations. The two
mechanisms combined play a central role in explaining the ubiquity of
oscillatory behaviour in real data and in simulation results of epidemic and
other related models.


  Tumour cells have to acquire a number of capabilities if a neoplasm is to
become a cancer. One of these key capabilities is increased motility which is
needed for invasion of other tissues and metastasis. This paper presents a
qualitative mathematical model based on game theory and computer simulations
using cellular automata. With this model we study the circumstances under which
mutations that confer increased motility to cells can spread through a tumour
made of rapidly proliferating cells. The analysis suggests therapies that could
help prevent the progression towards malignancy and invasiveness of benign
tumours.


  Tumour progression has been described as a sequence of traits or phenotypes
that cells have to acquire if the neoplasm is to become an invasive and
malignant cancer. Although the genetic mutations that lead to these phenotypes
are random, the process by which some of these mutations become successful and
spread is influenced by the tumour microenvironment and the presence of other
phenotypes. It is thus likely that some phenotypes that are essential in tumour
progression will emerge in the tumour population only with the prior presence
of other different phenotypes. In this paper we use evolutionary game theory to
analyse the interactions between three different tumour cell phenotypes defined
by autonomous growth, anaerobic glycolysis, and cancer cell invasion. The model
allows to understand certain specific aspects of glioma progression such as the
emergence of diffuse tumour cell invasion in low-grade tumours. We find that
the invasive phenotype is more likely to evolve after the appearance of the
glycolytic phenotype which would explain the ubiquitous presence of invasive
growth in malignant tumours. The result suggests that therapies which increase
the fitness cost of switching to anaerobic glycolysis might decrease the
probability of the emergence of more invasive phenotypes


  Environmental and genetic mutations can transform the cells in a co-operating
healthy tissue into an ecosystem of individualistic tumour cells that compete
for space and resources. Various selection forces are responsible for driving
the evolution of cells in a tumour towards more malignant and aggressive
phenotypes that tend to have a fitness advantage over the older populations.
Although the evolutionary nature of cancer has been recognised for more than
three decades (ever since the seminal work of Nowell) it has been only recently
that tools traditionally used by ecological and evolutionary researchers have
been adopted to study the evolution of cancer phenotypes in populations of
individuals capable of co-operation and competition. In this chapter we will
describe game theory as an important tool to study the emergence of cell
phenotypes in a tumour and will critically review some of its applications in
cancer research. These applications demonstrate that game theory can be used to
understand the dynamics of somatic cancer evolution and suggest new therapies
in which this knowledge could be applied to gain some control over the
evolution of the tumour.


  Normal mode analysis is a widely used technique for reconstructing
conformational changes of proteins from the knowledge of native structures. In
this Letter, we investigate to what extent normal modes capture the salient
features of the dynamics over a range of temperatures from close to T = 0 to
above unfolding. We show that on the one hand, the use of normal modes at
physiological temperatures is justified provided proteins are cooperative. On
the other hand, it is imperative to consider several modes in order to
eliminate the unpredictable temperature dependence of single- mode
contributions to the protein fluctuations.


  Interacting populations often create complicated spatiotemporal behavior, and
understanding it is a basic problem in the dynamics of spatial systems. We
study the two-species case by simulations of a host--parasitoid model. In the
case of co-existence, there are spatial patterns leading to noise-sustained
oscillations. We introduce a new measure for the patterns, and explain the
oscillations as a consequence of a timescale separation and noise. They are
linked together with the patterns by letting the spreading rates depend on
instantaneous population densities. Applications are discussed.


  Recently, computational modelling became a very important research tool that
enables us to study problems that for decades evaded scientific analysis.
Evolutionary systems are certainly examples of such problems: they are composed
of many units that might reproduce, diffuse, mutate, die, or in some cases for
example communicate. These processes might be of some adaptive value, they
influence each other and occur on various time scales. That is why such systems
are so difficult to study. In this paper we briefly review some computational
approaches, as well as our contributions, to the evolution of ecosystems and
language. We start from Lotka-Volterra equations and the modelling of simple
two-species prey-predator systems. Such systems are canonical example for
studying oscillatory behaviour in competitive populations. Then we describe
various approaches to study long-term evolution of multi-species ecosystems. We
emphasize the need to use models that take into account both ecological and
evolutionary processes. Finally, we address the problem of the emergence and
development of language. It is becoming more and more evident that any theory
of language origin and development must be consistent with darwinian principles
of evolution. Consequently, a number of techniques developed for modelling
evolution of complex ecosystems are being applied to the problem of language.
We briefly review some of these approaches.