  Metrics for indirectly predicting the folding rates of RNA sequences are of
interest. In this letter, we introduce a simple metric of RNA structural
complexity, which accounts for differences in the energetic contributions of
RNA base contacts toward RNA structure formation. We apply the metric to RNA
sequences whose folding rates were previously determined experimentally. We
find that the metric has good correlation (correlation coefficient: -0.95, p <<
0.01) with the logarithmically transformed folding rates of those RNA
sequences. This suggests that the metric can be useful for predicting RNA
folding rates. We use the metric to predict the folding rates of bacterial and
eukaryotic group II introns. Future applications of the metric (e.g., to
predict structural RNAs) could prove fruitful.


  Consider a Markov chain $(\xi_v)_{v \in V} \in [k]^V$ on the infinite $b$-ary
tree $T = (V,E)$ with irreducible edge transition matrix $M$, where $b \geq 2$,
$k \geq 2$ and $[k] = \{1,...,k\}$. We denote by $L_n$ the level-$n$ vertices
of $T$. Assume $M$ has a real second-largest (in absolute value) eigenvalue
$\lambda$ with corresponding real eigenvector $\nu \neq 0$. Letting $\sigma_v =
\nu_{\xi_v}$, we consider the following root-state estimator, which was
introduced by Mossel and Peres (2003) in the context of the "recontruction
problem" on trees: \begin{equation*} S_n = (b\lambda)^{-n} \sum_{x\in L_n}
\sigma_x. \end{equation*} As noted by Mossel and Peres, when $b\lambda^2 > 1$
(the so-called Kesten-Stigum reconstruction phase) the quantity $S_n$ has
uniformly bounded variance. Here, we give bounds on the moment-generating
functions of $S_n$ and $S_n^2$ when $b\lambda^2 > 1$. Our results have
implications for the inference of evolutionary trees.


  We introduce a new distance-based phylogeny reconstruction technique which
provably achieves, at sufficiently short branch lengths, a polylogarithmic
sequence-length requirement -- improving significantly over previous polynomial
bounds for distance-based methods. The technique is based on an averaging
procedure that implicitly reconstructs ancestral sequences.
  In the same token, we extend previous results on phase transitions in
phylogeny reconstruction to general time-reversible models. More precisely, we
show that in the so-called Kesten-Stigum zone (roughly, a region of the
parameter space where ancestral sequences are well approximated by ``linear
combinations'' of the observed sequences) sequences of length $\poly(\log n)$
suffice for reconstruction when branch lengths are discretized. Here $n$ is the
number of extant species.
  Our results challenge, to some extent, the conventional wisdom that estimates
of evolutionary distances alone carry significantly less information about
phylogenies than full sequence datasets.


  This study examines a new formulation of non-equilibrium thermodynamics,
which gives a conditional derivation of the ``maximum entropy production''
(MEP) principle for flow and/or chemical reaction systems at steady state. The
analysis uses a dimensionless potential function $\phi_{st}$ for
non-equilibrium systems, analogous to the free energy concept of equilibrium
thermodynamics. Spontaneous reductions in $\phi_{st}$ arise from increases in
the ``flux entropy'' of the system - a measure of the variability of the fluxes
- or in the local entropy production; conditionally, depending on the behaviour
of the flux entropy, the formulation reduces to the MEP principle. The inferred
steady state is also shown to exhibit high variability in its instantaneous
fluxes and rates, consistent with the observed behaviour of turbulent fluid
flow, heat convection and biological systems; one consequence is the
coexistence of energy producers and consumers in ecological systems. The
different paths for attaining steady state are also classified.


  We consider the optimal dynamics in the infinite population evolution models
with general symmetric fitness landscape. The search of optimal evolution
trajectories are complicated due to sharp transitions (like shock waves) in
evolution dynamics with smooth fitness landscapes, which exist even in case of
popular quadratic fitness. We found exact analytical solutions for
discontinuous dynamics at the large genome length limit.
  We found the optimal mutation rates for the fixed fitness landscape. The
single peak fitness landscape gives the fastest dynamics to send the vast
majority of the population from the initial sequence to the neighborhood of the
final sequence.


  Proteins are large and complex molecular machines. In order to perform their
function, most of them need energy, e.g. either in the form of a photon, like
in the case of the visual pigment rhodopsin, or through the breaking of a
chemical bond, as in the presence of adenosine triphosphate (ATP). Such energy,
in turn, has to be transmitted to specific locations, often several tens of
Angstroms away from where it is initially released. Here we show, within the
framework of a coarse-grained nonlinear network model, that energy in a protein
can jump from site to site with high yields, covering in many instances
remarkably large distances. Following single-site excitations, few specific
sites are targeted, systematically within the stiffest regions. Such energy
transfers mark the spontaneous formation of a localized mode of nonlinear
origin at the destination site, which acts as an efficient energy-accumulating
centre. Interestingly, yields are found to be optimum for excitation energies
in the range of biologically relevant ones.


  Aerial displays of starlings (Sturnus vulgaris) at their communal roosts are
complex: thousands of individuals form multiple flocks which are continually
changing shape and density, while splitting and merging. To understand these
complex displays both empirical data and models are needed. Whereas detailed
empirical data were recently collected through video recordings and position
measurements by stereo photography of flocks of thousands of starlings, there
are as yet no models that generate these complex patterns. Numerous computer
models in biology, however, suggest that patterns of single groups of moving
animals may emerge by self-organisation from movement and local coordination
(through attraction, alignment and avoidance of collision). In this paper, we
investigated whether this approach can be extended to generate patterns
resembling these aerial displays of starlings. We show in a model that to
generate many of the patterns measured empirically in real starlings we have to
extend the usual rules of local coordination with specifics of starling
behaviour, mainly 1) their aerial locomotion, 2) a low and constant number of
interaction-partners and 3) preferential movement above a roosting area. Our
model can be used as a tool for the study of these displays, because it
provides new integrative hypotheses about the mechanisms underlying these
displays and of swarming patterns in biological systems in general.


  Linyphiid spiders have evolved the ability to disperse long distances by a
process known as ballooning. It has been hypothesized that ballooning may allow
populations to persist in the highly disturbed agricultural areas that the
spiders prefer. In this study, I develop a stochastic population model to
explore how the propensity for this type of long distance dispersal influences
long term population persistence in a heterogeneous landscape where
catastrophic mortality events are common. Analysis of this model indicates that
although some dispersal does indeed decrease the probability of extinction of
the population, the frequency of dispersal is only important in certain
extremes. Instead, both the mean population birth and death rates, and the
landscape composition, are much more important in determining the probability
of extinction than the dispersal process. Thus, in order to develop effective
conservation strategies for these spiders, better understanding of life history
processes should be prioritized over an understanding of dispersal strategies.


  We present a stochastic computational model to study the mechanism of
signalling between a source and a target ionic transporter, both localized on
the plasma membrane (PM) and in intracellular nanometre-scale subplasmalemmal
signalling compartments comprising the PM, the sarcoplasmic reticulum (SR),
Ca2+ and Na+ transporters, and the intervening cytosol. We refer to these
compartments, sometimes called junctions, as cytoplasmic nanospaces or
nanodomains. In the chain of events leading to Ca2+ influx for SR reloading
during asynchronous Ca2+ waves in vascular smooth muscle (VSM), the physical
and functional link between non-selective cation channels (NSCC) and Na+/Ca2+
exchangers (NCX) needs to be elucidated in view of two recent findings: the
identification of the transient receptor potential canonical channel 6 (TRPC6)
as a crucial NSCC in VSM cells and the observation of localized cytosolic [Na+]
transients following purinergic stimulation of these cells. Having previously
helped clarify the Ca2+ signalling step between NCX and SERCA behind SR Ca2+
refilling, this quantitative approach now allows us to model the upstream
linkage of NSCC and NCX. We have implemented a random walk (RW) Monte Carlo
(MC) model with simulations mimicking a Na+ diffusion process originating at
the NSCC within PM-SR junctions. The model calculates the average [Na+] in the
nanospace and also produces [Na+] profiles as a function of distance from the
Na+ source. Our results highlight the necessity of a strategic juxtaposition of
the relevant signalling channels as well as other physical structures within
the nanospaces to permit adequate [Na+] build-up to provoke NCX reversal and
Ca2+ influx to refill the SR.


  A hierarchical structure describing the inter-relationships of species has
long been a fundamental concept in systematic biology, from Linnean
classification through to the more recent quest for a 'Tree of Life.' In this
paper we use an approach based on discrete mathematics to address a basic
question: Could one delineate this hierarchical structure in nature purely by
reference to the 'genealogy' of present-day individuals, which describes how
they are related with one another by ancestry through a continuous line of
descent? We describe several mathematically precise ways by which one can
naturally define collections of subsets of present day individuals so that
these subsets are nested (and so form a tree) based purely on the directed
graph that describes the ancestry of these individuals. We also explore the
relationship between these and related clustering constructions.


  In spatially distributed cellular systems, it is often convenient to
represent complicated auxiliary pathways and spatial transport by time-delayed
reaction rates. Furthermore, many of the reactants appear in low numbers
necessitating a probabilistic description. The coupling of delayed rates with
stochastic dynamics leads to a probability conservation equation characterizing
a non-Markovian process. A systematic approximation is derived that
incorporates the effect of delayed rates on the characterization of molecular
noise, valid in the limit of long delay time. By way of a simple example, we
show that delayed reaction dynamics can only increase intrinsic fluctuations
about the steady-state. The method is general enough to accommodate nonlinear
transition rates, allowing characterization of fluctuations around a
delay-induced limit cycle.


  Boolean networks have long been used as models of molecular networks and play
an increasingly important role in systems biology. This paper describes a
software package, Polynome, offered as a web service, that helps users
construct Boolean network models based on experimental data and biological
input. The key feature is a discrete analog of parameter estimation for
continuous models. With only experimental data as input, the software can be
used as a tool for reverse-engineering of Boolean network models from
experimental time course data.


  A quantitative theory on the construction and the evolution of the genetic
code is proposed. Through introducing the concept of mutational deterioration
(MD) and developing a theoretical formalism on MD minimization we have proved:
1, the redundancy distribution of codons in the genetic code obeys MD
minimization principle; 2, the hydrophilic-hydrophobic distribution of amino
acids on the code table is global MD (GMD) minimal; 3, the standard genetic
code can be deduced from the adaptive minimization of GMD; 4, the variants of
the standard genetic code can be explained quantitatively by use of GMD
formalism and the general trend of the evolution is GMD non-increasing which
reflects the selection on the code. We have demonstrated that the redundancy
distribution of codons and the hydrophobic-hydrophilic (H-P) distribution of
amino acids are robust in the code relative to the mutational parameter, and
indicated that the GMD can be looked as a non-fitness function on the adaptive
landscape. Finally, an important aspect on the symmetry of the code
construction, the Yin-Yang duality is investigated. The Yin-Yang duality among
codons affords a sound basis for understanding the H-P structure in the genetic
code.


  Recent discoveries show steady improvements in life expectancy during modern
decades. Does this support that humans continue to live longer in future? We
recently put forward the maximum survival tendency, as found in survival curves
of industrialized countries, which is described by extended Weibull model with
age-dependent stretched exponent. The maximum survival tendency suggests that
human survival dynamics may possess its intrinsic limit, beyond which survival
is inevitably forbidden. Based on such tendency, we develop the model and
explore the patterns in the maximum lifespan limits from industrialized
countries during recent three decades. This analysis strategy is simple and
useful to interpret the complicated human survival dynamics.


  We investigated the error-minimization properties of putative primordial
codes that consisted of 16 supercodons, with the third base being completely
redundant, using a previously derived cost function and the error minimization
percentage as the measure of a code's robustness to mistranslation. It is shown
that, when the 16-supercodon table is populated with 10 putative primordial
amino acids, inferred from the results of abiotic synthesis experiments and
other evidence independent of the code evolution, and with minimal assumptions
used to assign the remaining supercodons, the resulting 2-letter codes are
nearly optimal in terms of the error minimization level. The results of the
computational experiments with putative primordial genetic codes that contained
only two meaningful letters in all codons and encoded 10 to 16 amino acids
indicate that such codes are likely to have been nearly optimal with respect to
the minimization of translation errors. This near-optimality could be the
outcome of extensive early selection during the co-evolution of the code with
the primordial, error-prone translation system, or a result of a unique,
accidental event. Under this hypothesis, the subsequent expansion of the code
resulted in a decrease of the error minimization level that became sustainable
owing to the evolution of a high-fidelity translation system.


  Based on a non-equilibrium mechanism for spatial pattern formation we study
how position information can be controlled by locally coupled discrete
dynamical networks, similar to gene regulation networks of cells in a
developing multicellular organism. As an example we study the developmental
problems of domain formation and proportion regulation in the presence of
noise, as well as in the presence of cell flow. We find that networks that
solve this task exhibit a hierarchical structure of information processing and
are of similar complexity as developmental circuits of living cells. Proportion
regulation is scalable with system size and leads to sharp, precisely localized
boundaries of gene expression domains, even for large numbers of cells. A
detailed analysis of noise-induced dynamics, using a mean-field approximation,
shows that noise in gene expression states stabilizes (rather than disrupts)
the spatial pattern in the presence of cell movements, both for stationary as
well as growing systems. Finally, we discuss how this mechanism could be
realized in the highly dynamic environment of growing tissues in multi-cellular
organisms.


  Genetic regulation is a key component in development, but a clear
understanding of the structure and dynamics of genetic networks is not yet at
hand. In this paper we investigate these properties within an artificial genome
model originally introduced by Reil (1999). We analyze statistical properties
of randomly generated genomes both on the sequence- and network level, and show
that this model correctly predicts the frequency of genes in genomes as found
in experimental data. Using an evolutionary algorithm based on stabilizing
selection for a phenotype, we show that dynamical robustness against single
base mutations, as well as against random changes in initial states of
regulatory dynamics that mimic stochastic fluctuations in environmental
conditions, can emerge in parallel. Point mutations at the sequence level have
strongly non-linear effects on network wiring, including as well structurally
neutral mutations and simultaneous rewiring of multiple connections, which
occasionally lead to strong reorganization of the attractor landscape and
metastability of evolutionary dynamics. Evolved genomes exhibit characteristic
patterns on both sequence and network level.


  We study the evolution of cooperation in public goods games on different
regular graphs as a function of the noise level underlying strategy adoptions.
We focus on the effects that are brought about by different group sizes of
public goods games in which individuals participate, revealing that larger
groups of players may induce qualitatively different behavior when approaching
the deterministic limit of strategy adoption. While by pairwise interactions an
intermediate uncertainty by strategy adoptions may ensure optimal conditions
for the survival of cooperators at a specific graph topology, larger groups
warrant this only in the vicinity of the deterministic limit independently from
the underlying graph. These discrepancies are attributed to the indirect
linkage of otherwise not directly connected players, which is brought about by
joint memberships within the larger groups. Thus, we show that increasing the
group size may introduce an effective transition of the interaction topology,
and that the latter shapes the noise dependence of the evolution of cooperation
in case of pairwise interactions only.


  It is a classical result of Stein and Waterman that the asymptotic number of
RNA secondary structures is $1.104366 n^{-3/2} 2.618034^n$. In this paper, we
study combinatorial asymptotics for two special subclasses of RNA secondary
structures - canonical and saturated structures. Canonical secondary structures
were introduced by Bompf\"unewerer et al., who noted that the run time of
Vienna RNA Package is substantially reduced when restricting computations to
canonical structures. Here we provide an explanation for the speed-up.
Saturated secondary structures have the property that no base pairs can be
added without violating the definition of secondary structure (i.e. introducing
a pseudoknot or base triple). Here we compute the asymptotic number of
saturated structures, we show that the asymptotic expected number of base pairs
is $0.337361 n$, and the asymptotic number of saturated stem-loop structures is
$0.323954 1.69562^n$, in contrast to the number $2^{n-2}$ of (arbitrary)
stem-loop structures as classically computed by Stein and Waterman. Finally, we
show that the density of states for [all resp. canonical resp. saturated]
secondary structures is asymptotically Gaussian. We introduce a stochastic
greedy method to sample random saturated structures, called quasi-random
saturated structures, and show that the expected number of base pairs of is
$0.340633 n$.


  The lytic-lysogenic switch of the temperate lactococcal phage TP901-1 is
fundamentally different from that of phage lambda. In phage TP901-1, the lytic
promoter PL is repressed by CI whereas repression of the lysogenic promoter PR
requires the presence of both of the antagonistic regulator proteins, MOR and
CI. We model the central part of the switch and compare the two cases for PR
repression: the one where the two regulators interact only on the DNA, and the
other where the two regulators form a heteromer complex in the cytoplasm prior
to DNA binding. The models are analyzed for bistability, and the predicted
promoter repression folds are compared to experimental data. We conclude that
the experimental data are best reproduced the latter case, where a heteromer
complex forms in solution. We further find that CI sequestration by the
formation of MOR:CI complexes in cytoplasm makes the genetic switch robust.


  Under the conditions prevalent in the late Pleistocene (small hunter-gatherer
groups and frequent inter-group conflicts), coevolution of gene-related
behavior and culturally transmitted group-level institutions provides a
plausible explanation for the parochial altruistic and reciprocator traits of
most modern humans. When, with the agricultural revolution, societies became
larger and more complex, the collective nature of the monitoring and punishment
of norm violators was no longer effective. This led to the emergence of new
institutions of governance and social hierarchies. Likely, the smooth
acceptance of the new institutions was possible only because, in the majority
of the population, the reciprocator trait had become an internalized norm.
However the new ruling class has its own dynamics which in turn may lead to new
social crisis. Using a simple model, inspired on previous work by Bowles and
Gintis, these effects are studied here.


  Recent studies have demonstrated that network approaches are highly
appropriate tools to understand the extreme complexity of the aging process.
The generality of the network concept helps to define and study the aging of
technological, social networks and ecosystems, which may give novel concepts to
cure age-related diseases. The current review focuses on the role of
protein-protein interaction networks (interactomes) in aging. Hubs and
inter-modular elements of both interactomes and signaling networks are key
regulators of the aging process. Aging induces an increase in the permeability
of several cellular compartments, such as the cell nucleus, introducing gross
changes in the representation of network structures. The large overlap between
aging genes and genes of age-related major diseases makes drugs which aid
healthy aging promising candidates for the prevention and treatment of
age-related diseases, such as cancer, atherosclerosis, diabetes and
neurodegenerative disorders. We also discuss a number of possible research
options to further explore the potential of the network concept in this
important field, and show that multi-target drugs (representing
"magic-buckshots" instead of the traditional "magic bullets") may become an
especially useful class of age-related future drugs.


  Complex systems are successfully reduced to interacting elements via the
network concept. Transport plays a key role in the survival of networks. For
example the specialized signaling cascades of cellular networks filter noise
and efficiently adapt the network structure to new stimuli. However, our
general understanding of transport mechanisms and signaling pathways in complex
systems is yet limited. Here we summarize the key network structures involved
in transport, list the solutions available to overloaded systems for relaxing
their load and outline a possible method for the computational determination of
signaling pathways. We highlight that in addition to hubs, bridges and the
network skeleton, the overlapping modular structure is also essential in
network transport. Moreover, by locating network elements in the space of
overlapping network modules and evaluating their distance in this "module
space", it may be possible to approximate signaling pathways computationally,
which, in turn could serve the identification of signaling pathways of complex
systems. Our model may be applicable in a wide range of fields including
traffic control or drug design.


  During an infection, HIV experiences strong selection by immune system T
cells. Recent experimental work has shown that MHC escape mutations form an
important pathway for HIV to avoid such selection. In this paper, we study a
model of MHC escape mutation. The model is a predator-prey model with two prey,
composed of two HIV variants, and one predator, the immune system CD8 cells. We
assume that one HIV variant is visible to CD8 cells and one is not. The model
takes the form of a system of stochastic differential equations. Motivated by
well-known results concerning the short life-cycle of HIV intrahost, we assume
that HIV population dynamics occur on a faster time scale then CD8 population
dynamics. This separation of time scales allows us to analyze our model using
an asymptotic approach.
  Using this model we study the impact of an MHC escape mutation on the
population dynamics and genetic evolution of the intrahost HIV population. From
the perspective of population dynamics, we show that the competition between
the visible and invisible HIV variants can reach steady states in which either
a single variant exists or in which coexistence occurs depending on the
parameter regime. We show that in some parameter regimes the end state of the
system is stochastic. From a genetics perspective, we study the impact of the
population dynamics on the lineages of HIV samples taken after an escape
mutation occurs. We show that the lineages go through severe bottlenecks and
that the lineage distribution can be characterized by a Kingman coalescent.


  Eastern Boundary Upwelling Systems (EBUS) are characterized by a high
productivity of plankton associated with large commercial fisheries, thus
playing key biological and socio-economical roles. The aim of this work is to
make a comparative study of these four upwelling systems focussing on their
surface stirring, using the Finite Size Lyapunov Exponents (FSLEs), and their
biological activity, based on satellite data. First, the spatial distribution
of horizontal mixing is analysed from time averages and from probability
density functions of FSLEs. Then we studied the temporal variability of surface
stirring focussing on the annual and seasonal cycle. There is a global negative
correlation between surface horizontal mixing and chlorophyll standing stocks
over the four areas. To try to better understand this inverse relationship, we
consider the vertical dimension by looking at the Ekman-transport and vertical
velocities. We suggest the possibility of a changing response of the
phytoplankton to sub/mesoscale turbulence, from a negative effect in the very
productive coastal areas to a positive one in the open ocean.


  A protein undergoes conformational dynamics with multiple time scales, which
results in fluctuating enzyme activities. Recent studies in single molecule
enzymology have observe this "age-old" dynamic disorder phenomenon directly.
However, the single molecule technique has its limitation. To be able to
observe this molecular effect with real biochemical functions {\it in situ}, we
propose to couple the fluctuations in enzymatic activity to noise propagations
in small protein interaction networks such as zeroth order ultra-sensitive
phosphorylation-dephosphorylation cycle. We showed that enzyme fluctuations
could indeed be amplified by orders of magnitude into fluctuations in the level
of substrate phosphorylation | a quantity widely interested in cellular
biology. Enzyme conformational fluctuations sufficiently slower than the
catalytic reaction turn over rate result in a bimodal concentration
distribution of the phosphorylated substrate. In return, this network amplified
single enzyme fluctuation can be used as a novel biochemical "reporter" for
measuring single enzyme conformational fluctuation rates.


  The microtubule assembly process has been extensively studied, but the
underlying molecular mechanism remains poorly understood. The structure of an
artificially generated sheet polymer that alternates two types of lateral
contacts and that directly converts into microtubules, has been proposed to
correspond to the intermediate sheet structure observed during microtubule
assembly. We have studied the self-assembly process of GMPCPP tubulins into
sheet and microtubule structures using thermodynamic analysis and stochastic
simulations. With the novel assumptions that tubulins can laterally interact in
two different forms, and allosterically affect neighboring lateral
interactions, we can explain existing experimental observations. At low
temperature, the allosteric effect results in the observed sheet structure with
alternating lateral interactions as the thermodynamically most stable form. At
normal microtubule assembly temperature, our work indicates that a class of
sheet structures resembling those observed at low temperature is transiently
trapped as an intermediate during the assembly process. This work may shed
light on the tubulin molecular interactions, and the role of sheet formation
during microtubule assembly.


  Detection of rare variants by resequencing is important for the
identification of individuals carrying disease variants. Rapid sequencing by
new technologies enables low-cost resequencing of target regions, although it
is still prohibitive to test more than a few individuals. In order to improve
cost trade-offs, it has recently been suggested to apply pooling designs which
enable the detection of carriers of rare alleles in groups of individuals.
However, this was shown to hold only for a relatively low number of individuals
in a pool, and requires the design of pooling schemes for particular cases.
  We propose a novel pooling design, based on a compressed sensing approach,
which is both general, simple and efficient. We model the experimental
procedure and show via computer simulations that it enables the recovery of
rare allele carriers out of larger groups than were possible before, especially
in situations where high coverage is obtained for each individual.
  Our approach can also be combined with barcoding techniques to enhance
performance and provide a feasible solution based on current resequencing
costs. For example, when targeting a small enough genomic region (~100
base-pairs) and using only ~10 sequencing lanes and ~10 distinct barcodes, one
can recover the identity of 4 rare allele carriers out of a population of over
4000 individuals.


  Determining the mechanism by which transfer RNAs (tRNAs) rapidly and
precisely transit through the ribosomal A, P and E sites during translation
remains a major goal in the study of protein synthesis. Here, we report the
real-time dynamics of the L1 stalk, a structural element of the large ribosomal
subunit that is implicated in directing tRNA movements during translation.
Within pre-translocation ribosomal complexes, the L1 stalk exists in a dynamic
equilibrium between open and closed conformations. Binding of elongation factor
G (EF-G) shifts this equilibrium towards the closed conformation through one of
at least two distinct kinetic mechanisms, where the identity of the P-site tRNA
dictates the kinetic route that is taken. Within post-translocation complexes,
L1 stalk dynamics are dependent on the presence and identity of the E-site
tRNA. Collectively, our data demonstrate that EF-G and the L1 stalk
allosterically collaborate to direct tRNA translocation from the P to the E
sites, and suggest a model for the release of E-site tRNA.


  In this work we deal with parameter estimation in a latent variable model,
namely the multiple-hidden i.i.d. model, which is derived from multiple
alignment algorithms. We first provide a rigorous formalism for the homology
structure of k sequences related by a star-shaped phylogenetic tree in the
context of multiple alignment based on indel evolution models. We discuss
possible definitions of likelihoods and compare them to the criterion used in
multiple alignment algorithms. Existence of two different Information
divergence rates is established and a divergence property is shown under
additional assumptions. This would yield consistency for the parameter in
parametrization schemes for which the divergence property holds. We finally
extend the definition of the multiple-hidden i.i.d. model and the results
obtained to the case in which the sequences are related by an arbitrary
phylogenetic tree. Simulations illustrate different cases which are not covered
by our results.


  Background: Hidden Markov models are widely employed by numerous
bioinformatics programs used today. Applications range widely from comparative
gene prediction to time-series analyses of micro-array data. The parameters of
the underlying models need to be adjusted for specific data sets, for example
the genome of a particular species, in order to maximize the prediction
accuracy. Computationally efficient algorithms for parameter training are thus
key to maximizing the usability of a wide range of bioinformatics applications.
  Results: We introduce two computationally efficient training algorithms, one
for Viterbi training and one for stochastic expectation maximization (EM)
training, which render the memory requirements independent of the sequence
length. Unlike the existing algorithms for Viterbi and stochastic EM training
which require a two-step procedure, our two new algorithms require only one
step and scan the input sequence in only one direction. We also implement these
two new algorithms and the already published linear-memory algorithm for EM
training into the hidden Markov model compiler HMM-Converter and examine their
respective practical merits for three small example models.
  Conclusions: Bioinformatics applications employing hidden Markov models can
use the two algorithms in order to make Viterbi training and stochastic EM
training more computationally efficient. Using these algorithms, parameter
training can thus be attempted for more complex models and longer training
sequences. The two new algorithms have the added advantage of being easier to
implement than the corresponding default algorithms for Viterbi training and
stochastic EM training.


  The survival probability of immobile targets, annihilated by a population of
random walkers on inhomogeneous discrete structures, such as disordered solids,
glasses, fractals, polymer networks and gels, is analytically investigated. It
is shown that, while it cannot in general be related to the number of distinct
visited points, as in the case of homogeneous lattices, in the case of bounded
coordination numbers its asymptotic behaviour at large times can still be
expressed in terms of the spectral dimension $\widetilde {d}$, and its exact
analytical expression is given. The results show that the asymptotic survival
probability is site independent on recurrent structures ($\widetilde{d}\leq2$),
while on transient structures ($\widetilde{d}>2$) it can strongly depend on the
target position, and such a dependence is explicitly calculated.


  Demographic models built from genetic data play important roles in
illuminating prehistorical events and serving as null models in genome scans
for selection. We introduce an inference method based on the joint frequency
spectrum of genetic variants within and between populations. For candidate
models we numerically compute the expected spectrum using a diffusion
approximation to the one-locus two-allele Wright-Fisher process, involving up
to three simultaneous populations. Our approach is a composite likelihood
scheme, since linkage between neutral loci alters the variance but not the
expectation of the frequency spectrum. We thus use bootstraps incorporating
linkage to estimate uncertainties for parameters and significance values for
hypothesis tests. Our method can also incorporate selection on single sites,
predicting the joint distribution of selected alleles among populations
experiencing a bevy of evolutionary forces, including expansions, contractions,
migrations, and admixture. As applications, we model human expansion out of
Africa and the settlement of the New World, using 5 Mb of noncoding DNA
resequenced in 68 individuals from 4 populations (YRI, CHB, CEU, and MXL) by
the Environmental Genome Project. We also combine our demographic model with a
previously estimated distribution of selective effects among newly arising
amino acid mutations to accurately predict the frequency spectrum of
nonsynonymous variants across three continental populations (YRI, CHB, CEU).


  Molecular dynamics simulations were used to determine the binding affinities
between between the hormone 17 beta-estradiol (E2) and different estrogen
receptor (ER) isoforms in the rainbow trout, Oncorhynchus mykiss. Previous
phylogenetic analysis indicates that a whole genome duplication prior to the
divergence of ray-finned fish led to two distinct ER beta isoforms, ER beta 1
and ER beta 2, and the recent whole genome duplication in the ancestral
salmonid created two ER alpha isoforms, ER alpha 1 and ER alpha 2. The
objective of our computational studies is to provide insight into the
underlying evolutionary pressures on these isoforms. For the ER alpha subtype
our results show that E2 binds preferentially to ER alpha 1 over ER alpha 2.
Tests of lineage specific dN/dS ratios indicate that the ligand binding domain
of the ER alpha 2 gene is evolving under relaxed selection relative to all
other ER alpha genes. Comparison with the highly conserved DNA binding domain
suggests that ER alpha 2 may be undergoing neofunctionalization possibly by
binding to another ligand. By contrast, both ER beta 1 and ER beta 2 bind
similarly to E2 and the best fitting model of selection indicates that the
ligand binding domain of all ER beta genes are evolving under the same level of
purifying selection, comparable to ER alpha 1.


  We study the kinetics for the search of an immobile target by randomly moving
searchers that detect it only upon encounter. The searchers perform
intermittent random walks on a one-dimensional lattice. Each searcher can step
on a nearest neighbor site with probability "alpha", or go off lattice with
probability "1 - \alpha" to move in a random direction until it lands back on
the lattice at a fixed distance L away from the departure point. Considering
"alpha" and L as optimization parameters, we seek to enhance the chances of
successful detection by minimizing the probability P_N that the target remains
undetected up to the maximal search time N. We show that even in this simple
model a number of very efficient search strategies can lead to a decrease of
P_N by orders of magnitude upon appropriate choices of "alpha" and L. We
demonstrate that, in general, such optimal intermittent strategies are much
more efficient than Brownian searches and are as efficient as search algorithms
based on random walks with heavy-tailed Cauchy jump-length distributions. In
addition, such intermittent strategies appear to be more advantageous than
Levy-based ones in that they lead to more thorough exploration of visited
regions in space and thus lend themselves to parallelization of the search
processes.


  Assume two different communities each of which maintain their respective
opinions mainly because of the weak interaction between the two communities. In
such a case, it is an interesting problem to find the necessary strength of
inter-community interaction in order for the two communities to reach a
consensus. In this paper, the Information Accumulation System (IAS) model is
applied to investigate the problem. With the application of the IAS model, the
opinion dynamics of the two-community problem is found to belong to a wider
class of two-species problems appearing in population dynamics or in
competition of two languages, for all of which the governing equations can be
described in terms of coupled logistic maps. Tipping diffusivity is defined as
the maximal inter-community interaction that the two communities can stay in
different opinions. For a problem with simple community structure and
homogeneous individuals, the tipping diffusivity is calculated theoretically.
As a conclusion of the paper, the convergence of the two communities to the
same value is less possible the more overall interaction, intra-community and
inter-community, takes place. This implies, for example, that the increase in
the interaction between individuals caused by the development of modern
communication tools, such as facebook and twitter, does not necessarily improve
the tendency towards global convergence between different communities. If the
number of internal links increases by a factor, the number of inter-community
links must be increased by an even higher factor, in order for consensus to be
the only stable attractor.


  Word matches are often used in sequence comparison methods, either as a
measure of sequence similarity or in the first search steps of algorithms such
as BLAST or BLAT. The D2 statistic is the number of matches of words of k
letters between two sequences. Recent advances have been made in the
characterisation of this statistic and in the approximation of its
distribution. Here, these results are extended to the case of approximate word
matches.
  We compute the exact value of the variance of the D2 statistic for the case
of a uniform letter distribution, and introduce a method to provide accurate
approximations of the variance in the remaining cases. This enables the
distribution of D2 to be approximated for typical situations arising in
biological research. We apply these results to the identification of
cis-regulatory modules, and show that this method detects such sequences with a
high accuracy.
  The ability to approximate the distribution of D2 for both exact and
approximate word matches will enable the use of this statistic in a more
precise manner for sequence comparison, database searches, and identification
of transcription factor binding sites.


  We consider the problem of estimating a sparse multi-response regression
function, with an application to expression quantitative trait locus (eQTL)
mapping, where the goal is to discover genetic variations that influence
gene-expression levels. In particular, we investigate a shrinkage technique
capable of capturing a given hierarchical structure over the responses, such as
a hierarchical clustering tree with leaf nodes for responses and internal nodes
for clusters of related responses at multiple granularity, and we seek to
leverage this structure to recover covariates relevant to each
hierarchically-defined cluster of responses. We propose a tree-guided group
lasso, or tree lasso, for estimating such structured sparsity under
multi-response regression by employing a novel penalty function constructed
from the tree. We describe a systematic weighting scheme for the overlapping
groups in the tree-penalty such that each regression coefficient is penalized
in a balanced manner despite the inhomogeneous multiplicity of group
memberships of the regression coefficients due to overlaps among groups. For
efficient optimization, we employ a smoothing proximal gradient method that was
originally developed for a general class of structured-sparsity-inducing
penalties. Using simulated and yeast data sets, we demonstrate that our method
shows a superior performance in terms of both prediction errors and recovery of
true sparsity patterns, compared to other methods for learning a
multivariate-response regression.


  We study the noise characteristics of stochastic oscillations in protein
number dynamics of simple genetic oscillatory systems. Using the
three-component negative feedback transcription regulatory system called the
repressilator as a prototypical example, we quantify the degree of fluctuations
in oscillation periods and amplitudes, as well as the noise propagation along
the regulatory cascade in the stable oscillation regime via dynamic Monte Carlo
simulations. For the single protein-species level, the fluctuation in the
oscillation amplitudes is found to be larger than that of the oscillation
periods, the distributions of which are reasonably described by the Weibull
distribution and the Gaussian tail, respectively. Correlations between
successive periods and between successive amplitudes, respectively, are
measured to assess the noise propagation properties, which are found to decay
faster for the amplitude than for the period. The local fluctuation property is
also studied.


  The genetic code has been shown to be very error robust compared to randomly
selected codes, but to be significantly less error robust than a certain code
found by a heuristic algorithm. We formulate this optimisation problem as a
Quadratic Assignment Problem and thus verify that the code found by the
heuristic is the global optimum. We also argue that it is strongly misleading
to compare the genetic code only with codes sampled from the fixed block model,
because the real code space is orders of magnitude larger. We thus enlarge the
space from which random codes can be sampled from approximately 2.433 x 10^18
codes to approximately 5.908 x 10^45 codes. We do this by leaving the fixed
block model, and using the wobble rules to formulate the characteristics
acceptable for a genetic code. By relaxing more constraints three larger spaces
are also constructed. Using a modified error function, the genetic code is
found to be more error robust compared to a background of randomly generated
codes with increasing space size. We point out that these results do not
necessarily imply that the code was optimized during evolution for error
minimization, but that other mechanisms could explain this error robustness.


  We study the survival of a prey that is hunted by N predators. The predators
perform independent random walks on a square lattice with V sites and start a
direct chase whenever the prey appears within their sighting range. The prey is
caught when a predator jumps to the site occupied by the prey. We analyze the
efficacy of a lazy, minimal-effort evasion strategy according to which the prey
tries to avoid encounters with the predators by making a hop only when any of
the predators appears within its sighting range; otherwise the prey stays
still. We show that if the sighting range of such a lazy prey is equal to 1
lattice spacing, at least 3 predators are needed in order to catch the prey on
a square lattice. In this situation, we establish a simple asymptotic relation
ln(Pev)(t) \sim (N/V)2ln(Pimm(t)) between the survival probabilities of an
evasive and an immobile prey. Hence, when the density of the predators is low
N/V<<1, the lazy evasion strategy leads to the spectacular increase of the
survival probability. We also argue that a short-sighting prey (its sighting
range is smaller than the sighting range of the predators) undergoes an
effective superdiffusive motion, as a result of its encounters with the
predators, whereas a far-sighting prey performs a diffusive-type motion.


  In a recent paper List, Elsholtz and Seeley [Phil. Trans. Roy. Soc. B. 364
(2009) 755] have devised an agent-based model of the the nest-choice dynamics
in swarms of honeybees, and have concluded that both interdependence and
independence are needed for the bees to reach a consensus on the best nest
site. We here present a simplified version of the model which can be treated
analytically with the tools of statistical physics and which largely has the
same features as the original dynamics. Based on our analytical approaches it
is possible to characterize the co-ordination outcome exactly on the
deterministic level, and to a good approximation if stochastic effects are
taken into account, reducing the need for computer simulations on the
agent-based level. In the second part of the paper we present a spatial
extension, and show that transient non-trivial patterns emerge, before
consensus is reached. Approaches in terms of Langevin equations for continuous
field variables are discussed.


  Today's sequencing technology allows sequencing an individual genome within a
few weeks for a fraction of the costs of the original Human Genome project.
Genomics labs are faced with dozens of TB of data per week that have to be
automatically processed and made available to scientists for further analysis.
This paper explores the potential and the limitations of using relational
database systems as the data processing platform for high-throughput genomics.
In particular, we are interested in the storage management for high-throughput
sequence data and in leveraging SQL and user-defined functions for data
analysis inside a database system. We give an overview of a database design for
high-throughput genomics, how we used a SQL Server database in some
unconventional ways to prototype this scenario, and we will discuss some
initial findings about the scalability and performance of such a more
database-centric approach.


  Phylogenetic data arising on two possibly different tree topologies might be
mixed through several biological mechanisms, including incomplete lineage
sorting or horizontal gene transfer in the case of different topologies, or
simply different substitution processes on characters in the case of the same
topology. Recent work on a 2-state symmetric model of character change showed
such a mixture model has non-identifiable parameters, and thus it is
theoretically impossible to determine the two tree topologies from any amount
of data under such circumstances. Here the question of identifiability is
investigated for 2-tree mixtures of the 4-state group-based models, which are
more relevant to DNA sequence data. Using algebraic techniques, we show that
the tree parameters are identifiable for the JC and K2P models. We also prove
that generic substitution parameters for the JC mixture models are
identifiable, and for the K2P and K3P models obtain generic identifiability
results for mixtures on the same tree. This indicates that the full
phylogenetic signal remains in such mixtures, and that the 2-state symmetric
result is thus a misleading guide to the behavior of other models.


  We study the dynamics of intracellular calcium oscillations in the presence
of proteins that bind calcium on multiple sites and that are generally believed
to act as passive calcium buffers in cells. We find that multisite
calcium-binding proteins set a sharp threshold for calcium oscillations. Even
with high concentrations of calcium-binding proteins, internal noise, which
shows up spontaneously in cells in the process of calcium wave formation, can
lead to self-oscillations. This produces oscillatory behaviors strikingly
similar to those observed in real cells. In addition, for given intracellular
concentrations of both calcium and calcium-binding proteins the regularity of
these oscillations changes and reaches a maximum as a function noise variance,
and the overall system dynamics displays stochastic coherence. We conclude that
calcium-binding proteins may have an important and active role in cellular
communication.


  Ring topologies of repressing genes have qualitatively different long-term
dynamics if the number of genes is odd (they oscillate) or even (they exhibit
bistability). However, these attractors may not fully explain the observed
behavior in transient and stochastic environments such as the cell. We show
here that even repressilators possess quasi-stable, travelling-wave periodic
solutions that are reachable, long-lived and robust to parameter changes. These
solutions underlie the sustained oscillations observed in even rings in the
stochastic regime, even if these circuits are expected to behave as switches.
The existence of such solutions can also be exploited for control purposes:
operation of the system around the quasi-stable orbit allows us to turn on and
off the oscillations reliably and on demand. We illustrate these ideas with a
simple protocol based on optical interference that can induce oscillations
robustly both in the stochastic and deterministic regimes.


  In this study, an efficient method to immunize modular networks (i.e.,
networks with community structure) is proposed. The immunization of networks
aims at fragmenting networks into small parts with a small number of removed
nodes. Its applications include prevention of epidemic spreading, intentional
attacks on networks, and conservation of ecosystems. Although preferential
immunization of hubs is efficient, good immunization strategies for modular
networks have not been established. On the basis of an immunization strategy
based on the eigenvector centrality, we develop an analytical framework for
immunizing modular networks. To this end, we quantify the contribution of each
node to the connectivity in a coarse-grained network among modules. We verify
the effectiveness of the proposed method by applying it to model and real
networks with modular structure.


  We propose a stochastic model for evolution. Births and deaths of species
occur with constant probabilities. Each new species is associated with a
fitness sampled from the uniform distribution on [0,1]. Every time there is a
death event then the type that is killed is the one with the smallest fitness.
We show that there is a sharp phase transition when the birth probability is
larger than the death probability. The set of species with fitness higher than
a certain critical value approach an uniform distribution. On the other hand
all the species with fitness less than the critical disappear after a finite
(random) time.


  The explanation is that in aortic tissue (both diseased and nondiseased) a
BAK1 pseudogene is expressed; while in the matching blood samples the actual
BAK1 gene is expressed. This explanation was reached after we realized that
BAK1 has two edited copies in human genome. These copies are probably BAK1
pseudogenes. One copy belongs to chromosome 11 (NG_005599.3) and the other to
chromosome 20 (NC_000850.5). The first copy has frameshifts which means that
probably it does not express any functional protein; by other hand, the
chromosome 20 copy has no frameshifts and what is more important contains all
the reported polymorphisms.


  On 11 June the World Health Organization officially raised the phase of
pandemic alert (with regard to the new H1N1 influenza strain) to level 6. We
use a global structured metapopulation model integrating mobility and
transportation data worldwide in order to estimate the transmission potential
and the relevant model parameters we used the data on the chronology of the
2009 novel influenza A(H1N1). The method is based on the maximum likelihood
analysis of the arrival time distribution generated by the model in 12
countries seeded by Mexico by using 1M computationally simulated epidemics. An
extended chronology including 93 countries worldwide seeded before 18 June was
used to ascertain the seasonality effects. We found the best estimate R0 = 1.75
(95% CI 1.64 to 1.88) for the basic reproductive number. Correlation analysis
allows the selection of the most probable seasonal behavior based on the
observed pattern, leading to the identification of plausible scenarios for the
future unfolding of the pandemic and the estimate of pandemic activity peaks in
the different hemispheres. We provide estimates for the number of
hospitalizations and the attack rate for the next wave as well as an extensive
sensitivity analysis on the disease parameter values. We also studied the
effect of systematic therapeutic use of antiviral drugs on the epidemic
timeline. The analysis shows the potential for an early epidemic peak occurring
in October/November in the Northern hemisphere, likely before large-scale
vaccination campaigns could be carried out. We suggest that the planning of
additional mitigation policies such as systematic antiviral treatments might be
the key to delay the activity peak inorder to restore the effectiveness of the
vaccination programs.


  A model of mutation rate evolution for multiple loci under arbitrary
selection is analyzed. Results are obtained using techniques from Karlin (1982)
that overcome the weak selection constraints needed for tractability in prior
studies of multilocus event models. A multivariate form of the reduction
principle is found: reduction results at individual loci combine topologically
to produce a surface of mutation rate alterations that are neutral for a new
modifier allele. New mutation rates survive if and only if they fall below this
surface - a generalization of the hyperplane found by Zhivotovsky et al. (1994)
for a multilocus recombination modifier. Increases in mutation rates at some
loci may evolve if compensated for by decreases at other loci. The strength of
selection on the modifier scales in proportion to the number of germline cell
divisions, and increases with the number of loci affected. Loci that do not
make a difference to marginal fitnesses at equilibrium are not subject to the
reduction principle, and under fine tuning of mutation rates would be expected
to have higher mutation rates than loci in mutation-selection balance. Other
results include the nonexistence of 'viability analogous, Hardy-Weinberg'
modifier polymorphisms under multiplicative mutation, and the sufficiency of
average transmission rates to encapsulate the effect of modifier polymorphisms
on the transmission of loci under selection. A conjecture is offered regarding
situations, like recombination in the presence of mutation, that exhibit
departures from the reduction principle. Constraints for tractability are:
tight linkage of all loci, initial fixation at the modifier locus, and mutation
distributions comprising transition probabilities of reversible Markov chains.


  The mean time for a diffusing ligand to activate a target protein located on
the surface of a microdomain can regulate cellular signaling. When the ligand
switches between various states induced by chemical interactions or
conformational changes, while target activation occurs in only one state, this
activation time is affected. We investigate this dynamics using new equations
for the sojourn times spent in each state. For two states, we obtain exact
solutions in dimension one, and asymptotic ones confirmed by Brownian
simulations in dimension 3. We find that the activation time is quite sensitive
to changes of the switching rates, which can be used to modulate signaling.
Interestingly, our analysis reveals that activation can be fast although the
ligand spends most of the time 'hidden' in the non-activating state. Finally,
we obtain a new formula for the narrow escape time in the presence of
switching.


  Filopodia are bundles of actin filaments that extend out ahead of the leading
edge of a crawling cell to probe its upcoming environment. {\it In vitro}
experiments [D. Vignjevic {\it et al.}, J. Cell Biol. {\bf 160}, 951 (2003)]
have determined the minimal ingredients required for the formation of filopodia
from the dendritic-like morphology of the leading edge. We model these
experiments using kinetic aggregation equations for the density of growing
bundle tips. In mean field, we determine the bundle size distribution to be
broad for bundle sizes smaller than a characteristic bundle size above which
the distribution decays exponentially. Two-dimensional simulations
incorporating both bundling and cross-linking measure a bundle size
distribution that agrees qualitatively with mean field. The simulations also
demonstrate a nonmonotonicity in the radial extent of the dendritic region as a
function of capping protein concentration, as was observed in experiments, due
to the interplay between percolation and the ratcheting of growing filaments
off a spherical obstacle.


  In order to produce specific complex structures from a large set of similar
biochemical building blocks, many biochemical systems require high sensitivity
to small molecular differences. The first and most common model used to explain
this high specificity is kinetic proofreading, which has been extended to a
variety of systems from detection of DNA mismatch to cell signaling processes.
While the specification properties of the kinetic proofreading model are well
known and were studied in various contexts, very little is known about its
temporal behavior. In this work, we study the dynamical properties of discrete
stochastic two branch kinetic proofreading schemes. Using the Laplace transform
of the corresponding chemical master equation, we obtain an analytical solution
for the completion time distribution. In particular we provide expressions for
the specificity and the mean and the variance of the process completion times.
We also show that, for a wide range of parameters a process distinguishing
between two different products can be reduced to a much simpler three point
process. Our results allow for the systematic study of the interplay between
specificity and completion times as well as testing the validity of the kinetic
proofreading model in biological systems.


  The possible mechanisms of communications between distant bio-systems by
means of optical and UV photons are studied. It is argued that their main
production mechanism is owed to the biochemical reactions, occurring during the
cell division.. In the proposed model the bio-systems perform such
communications, radiating the photons in form of short periodic bursts, which
were observed experimentally for fish and frog eggs1. For experimentally
measured photon rates the communication algorithm is supposedly similar to the
exchange of binary encoded data in computer net via optical channels


  The genome-wide recombination rate ($RR$) of a species is often described by
one parameter, the ratio between total genetic map length ($G$) and physical
map length ($P$), measured in centiMorgans per Megabase (cM/Mb). The value of
this parameter varies greatly between species, but the cause for these
differences is not entirely clear. A constraining factor of overall $RR$ in a
species, which may cause increased $RR$ for smaller chromosomes, is the
requirement of at least one chiasma per chromosome (or chromosome-arm) per
meiosis. In the present study, we quantify the relative excess of recombination
events on smaller chromosomes by a linear regression model, which relates the
genetic length of chromosomes to their physical length. We find for several
species that the two-parameter regression, $G= G_0 + k \cdot P$ provides a
better characterization of the relationship between genetic and physical map
length than the one-parameter regression that runs through the origin. A
non-zero intercept ($G_0$) indicates a relative excess of recombination on
smaller chromosomes in a genome. Given $G_0$, the parameter $k$ predicts the
increase of genetic map length over the increase of physical map length. The
observed values of $G_0$ have a similar magnitude for diverse species, whereas
$k$ varies by two orders of magnitude. The implications of this strategy for
the genetic maps of human, mouse, rat, chicken, honeybee, worm and yeast are
discussed.


  Background: By a new concept called "phase diagram", we compare two commonly
used genotype-based tests for case-control genetic analysis, one is a
Cochran-Armitage trend test (CAT test at $x=0.5$, or CAT0.5) and another
(called MAX2) is the maximization of two chi-square test results: one from the
two-by-two genotype count table that combines the baseline homozygotes and
heterozygotes, and another from the table that combines heterozygotes with risk
homozygotes. CAT0.5 is more suitable for multiplicative disease models and MAX2
is better for dominant/recessive models.
  Methods: We define the CAT0.5-MAX2 phase diagram on the disease model space
such that regions where MAX2 is more powerful than CAT0.5 are separated from
regions where the CAT0.5 is more powerful, and the task is to choose the
appropriate parameterization to make the separation possible.
  Results: We find that using the difference of allele frequencies ($\delta_p$)
and the difference of Hardy-Weinberg disequilibrium coefficients
($\delta_\epsilon$) can separate the two phases well, and the phase boundaries
are determined by the angle $tan^{-1}(\delta_p/\delta_\epsilon)$, which is an
improvement over the disease model selection using $\delta_\epsilon$ only.
  Conclusions: We argue that phase diagrams similar to the one for CAT0.5-MAX2
have graphical appeals in understanding power performance of various tests,
clarifying simulation schemes, summarizing case-control datasets, and guessing
the possible mode of inheritance.


  Two non-integer parameters are defined for MAX statistics, which are maxima
of $d$ simpler test statistics. The first parameter, $d_{MAX}$, is the
fractional number of tests, representing the equivalent numbers of independent
tests in MAX. If the $d$ tests are dependent, $d_{MAX} < d$. The second
parameter is the fractional degrees of freedom $k$ of the chi-square
distribution $\chi^2_k$ that fits the MAX null distribution. These two
parameters, $d_{MAX}$ and $k$, can be independently defined, and $k$ can be
non-integer even if $d_{MAX}$ is an integer. We illustrate these two parameters
using the example of MAX2 and MAX3 statistics in genetic case-control studies.
We speculate that $k$ is related to the amount of ambiguity of the model
inferred by the test. In the case-control genetic association, tests with low
$k$ (e.g. $k=1$) are able to provide definitive information about the disease
model, as versus tests with high $k$ (e.g. $k=2$) that are completely uncertain
about the disease model. Similar to Heisenberg's uncertain principle, the
ability to infer disease model and the ability to detect significant
association may not be simultaneously optimized, and $k$ seems to measure the
level of their balance.


  The issues we attempt to tackle here are what the first peptides did look
like when they emerged on the primitive earth, and what simple catalytic
activities they fulfilled. We conjecture that the early functional peptides
were short (3 to 8 amino acids long), were made of those amino acids, Gly, Ala,
Val and Asp, that are abundantly produced in many prebiotic synthesis
experiments and observed in meteorites, and that the neutralization of Asp's
negative charge is achieved by metal ions. We further assume that some traces
of these prebiotic peptides still exist, in the form of active sites in
present-day proteins. Searching these proteins for prebiotic peptide candidates
led us to identify three main classes of motifs, bound mainly to Mg^{2+} ions:
D(F/Y)DGD corresponding to the active site in RNA polymerases, DGD(G/A)D
present in some kinds of mutases, and DAKVGDGD in dihydroxyacetone kinase. All
three motifs contain a DGD submotif, which is suggested to be the common
ancestor of all active peptides. Moreover, all three manipulate phosphate
groups, which was probably a very important biological function in the very
first stages of life. The statistical significance of our results is supported
by the frequency of these motifs in today's proteins, which is three times
higher than expected by chance, with a P-value of 3 10^{-2}. The implications
of our findings in the context of the appearance of life and the possibility of
an experimental validation are discussed.


  A fundamental question for evolutionary biology is why rates of evolution
vary dramatically between proteins. Perhaps surprisingly, it is controversial
how much a protein's functional importance affects its rate of evolution. In
most studies, functional importance has been measured on the coarse scale of
protein knock-outs, while evolutionary rate has been measured on the fine scale
of amino acid substitutions. Here we introduce dynamical influence, a finer
measure of protein functional importance. To measure dynamical influence, we
first use detailed biochemical models of particular reaction networks to
measure the influence of each reaction rate constant on network dynamics. We
then define the dynamical influence of a protein to be the average influence of
the rate constants for all reactions it is involved in.
  Using models of a dozen biochemical systems and sequence data from
vertebrates, we show that dynamical influence and evolutionary rate are
negatively correlated; proteins with greater dynamical influence evolve more
slowly. We also show that proteins with greater dynamical influence are not
more likely to be essential. This suggests that there are many cellular
reactions whose presence is essential for life, but whose quantitative rate is
relatively unimportant to fitness. We also provide evidence that the effect of
dynamical influence on evolutionary rate is independent of protein expression
level, expression specificity, gene compactness, and reaction degree. Dynamical
influence offers a finer view of functional importance, and our results suggest
that focusing on essentiality may have previously led to an underestimation of
the role function plays in protein evolution.


  Background: Regions with copy number variations (in germline cells) or copy
number alteration (in somatic cells) are of great interest for human disease
gene mapping and cancer studies. They represent a new type of mutation and are
larger-scaled than the single nucleotide polymorphisms. Using genotyping
microarray for copy number variation detection has become standard, and there
is a need for improving analysis methods. Results: We apply the cumulative plot
to the detection of regions with copy number variation/alteration, on samples
taken from a chronic lymphocytic leukemia patient. Two sets of whole-genome
genotyping of 317k single nucleotide polymorphisms, one from the normal cell
and another from the cancer cell, are analyzed. We demonstrate the utility of
cumulative plot in detecting a 9Mb (9 x 10^6 bases) hemizygous deletion and 1Mb
homozygous deletion on chromosome 13. We also show the possibility to detect
smaller copy number variation/alteration regions below the 100kb range.
Conclusions: As a graphic tool, the cumulative plot is an intuitive and a
scale-free (window-less) way for detecting copy number variation/alteration
regions, especially when such regions are small.


  {\bf Background}: Several features are known to correlate with the GC-content
in the human genome, including recombination rate, gene density and distance to
telomere. However, by testing for pairwise correlation only, it is impossible
to distinguish direct associations from indirect ones and to distinguish
between causes and effects. {\bf Results}: We use partial correlations to
construct partially directed graphs for the following four variables:
GC-content, recombination rate, exon density and distance-to-telomere.
Recombination rate and exon density are unconditionally uncorrelated, but
become inversely correlated by conditioning on GC-content. This pattern
indicates a model where recombination rate and exon density are two independent
causes of GC-content variation. {\bf Conclusions}: Causal inference and
graphical models are useful methods to understand genome evolution and the
mechanisms of isochore evolution in the human genome.


  The bifurcation diagram of a model stochastic differential equation with
delayed feedback is presented. We are motivated by recent research on
stochastic effects in models of transcriptional gene regulation. We start from
the normal form for a pitchfork bifurcation, and add multiplicative or
parametric noise and linear delayed feedback. The latter is sufficient to
originate a Hopf bifurcation in that region of parameters in which there is a
sufficiently strong negative feedback. We find a sharp bifurcation in parameter
space, and define the threshold as the point in which the stationary
distribution function p(x) changes from a delta function at the trivial state
x=0 to p(x) ~ x^alpha at small x (with alpha = -1 exactly at threshold). We
find that the bifurcation threshold is shifted by fluctuations relative to the
deterministic limit by an amount that scales linearly with the noise intensity.
Analytic calculations of the bifurcation threshold are also presented in the
limit of small delay tau -> 0 that compare quite favorably with the numerical
solutions even for tau = 1.


  Some important biomolecules (for instance, bacterial FtsZ and eukaryotic DNA)
are known to posses spontaneous (intrinsic) curvature. Using a simple extension
of the wormlike chain model, we study the response of a weakly bending filament
in two dimensions to a pulling force applied at its ends (a configuration
common in classical in-vitro experiments and relevant to several in-vivo cell
cases). The spontaneous curvature of such a chain or filament can in general be
arc-length dependent and we study a case of sinusoidal variation, from which an
arbitrary case can be reconstructed via Fourier transformation. We obtain
analytic results for the force-extension relationship and the width of
transverse fluctuations. We show that spontaneous-curvature undulations can
affect the force-extension behavior even in relatively flexible filaments with
a persistence length smaller than the contour length.


  Evolutionary adaptation is often likened to climbing a hill or peak. While
this process is simple for fitness landscapes where mutations are independent,
the interaction between mutations (epistasis) as well as mutations at loci that
affect more than one trait (pleiotropy) are crucial in complex and realistic
fitness landscapes. We investigate the impact of epistasis and pleiotropy on
adaptive evolution by studying the evolution of a population of asexual haploid
organisms (haplotypes) in a model of N interacting loci, where each locus
interacts with K other loci. We use a quantitative measure of the magnitude of
epistatic interactions between substitutions, and find that it is an increasing
function of K. When haplotypes adapt at high mutation rates, more epistatic
pairs of substitutions are observed on the line of descent than expected. The
highest fitness is attained in landscapes with an intermediate amount of
ruggedness that balance the higher fitness potential of interacting genes with
their concomitant decreased evolvability. Our findings imply that the synergism
between loci that interact epistatically is crucial for evolving genetic
modules with high fitness, while too much ruggedness stalls the adaptive
process.


  Significant volumes of knowledge have been accumulated in recent years
linking subtle genetic variations to a wide variety of medical disorders from
Cystic Fibrosis to mental retardation. Nevertheless, there are still great
challenges in applying this knowledge routinely in the clinic, largely due to
the relatively tedious and expensive process of DNA sequencing. Since the
genetic polymorphisms that underlie these disorders are relatively rare in the
human population, the presence or absence of a disease-linked polymorphism can
be thought of as a sparse signal. Using methods and ideas from compressed
sensing and group testing, we have developed a cost-effective genotyping
protocol. In particular, we have adapted our scheme to a recently developed
class of high throughput DNA sequencing technologies, and assembled a
mathematical framework that has some important distinctions from 'traditional'
compressed sensing ideas in order to address different biological and technical
constraints.


  We study the distribution of the maximum of a set of random fitnesses with
fixed number of mutations in a model of biological evolution. The fitness
variables are not independent and the correlations can be varied via a
parameter $\ell=1,...,L$. We present analytical calculations for the following
three solvable cases: (i) one-step mutants with arbitrary $\ell$ (ii) weakly
correlated fitnesses with $\ell=L/2$ (iii) strongly correlated fitnesses with
$\ell=2$. In all these cases, we find that the limit distribution for the
maximum fitness is not of the standard Gumbel form.


  We study the evolution of cooperation in the prisoner's dilemma game, whereby
a coevolutionary rule is introduced that molds the random topology of the
interaction network in two ways. First, existing links are deleted whenever a
player adopts a new strategy or its degree exceeds a threshold value, and
second, new links are added randomly after a given number of game iterations.
These coevolutionary processes correspond to the generic formation of new and
deletion of existing links that, especially in human societies, appear
frequently as a consequence of ongoing socialization, change of lifestyle or
death. Due to the counteraction of deletions and additions of links the initial
heterogeneity of the interaction network is qualitatively preserved, and thus
cannot be held responsible for the observed promotion of cooperation. Indeed,
the coevolutionary rule evokes the spontaneous emergence of a powerful
multilevel selection mechanism, which despite of the sustained random topology
of the evolving network, maintains cooperation across the whole span of
defection temptation values.


  Confidence is an essential ingredient of success in a wide range of domains
ranging from job performance and mental health, to sports, business, and
combat. Some authors have suggested that not just confidence but
overconfidence-believing you are better than you are in reality-is advantageous
because it serves to increase ambition, morale, resolve, persistence, or the
credibility of bluffing, generating a self-fulfilling prophecy in which
exaggerated confidence actually increases the probability of success. However,
overconfidence also leads to faulty assessments, unrealistic expectations, and
hazardous decisions, so it remains a puzzle how such a false belief could
evolve or remain stable in a population of competing strategies that include
accurate, unbiased beliefs. Here, we present an evolutionary model showing
that, counter-intuitively, overconfidence maximizes individual fitness and
populations will tend to become overconfident, as long as benefits from
contested resources are sufficiently large compared to the cost of competition.
In contrast, "rational" unbiased strategies are only stable under limited
conditions. The fact that overconfident populations are evolutionarily stable
in a wide range of environments may help to explain why overconfidence remains
prevalent today, even if it contributes to hubris, market bubbles, financial
collapses, policy failures, disasters, and costly wars.


  Proteomics can be defined as the large-scale analysis of proteins. Due to the
complexity of biological systems, it is required to concatenate various
separation techniques prior to mass spectrometry. These techniques, dealing
with proteins or peptides, can rely on chromatography or electrophoresis. In
this review, the electrophoretic techniques are under scrutiny. Their
principles are recalled, and their applications for peptide and protein
separations are presented and critically discussed. In addition, the features
that are specific to gel electrophoresis and that interplay with mass
spectrometry (i.e., protein detection after electrophoresis, and the process
leading from a gel piece to a solution of peptides) are also discussed.


  We propose a new model describing the production and the establishment of the
stable gradient of the Bicoid protein along the antero-posterior axis of the
embryo of \textit{Drosophila}. In this model, we consider that \textit{bicoid}
mRNA diffuses along the antero-posterior axis of the embryo and the protein is
produced in the ribosomes localized near the syncytial nuclei. Bicoid protein
stays localized near the syncytial nuclei as observed in experiments. We
calibrate the parameters of the mathematical model with experimental data taken
during the cleavage stages 11 to 14 of the developing embryo of
\textit{Drosophila}. We obtain good agreement between the experimental and the
model gradients, with relative errors in the range 5-8%. The inferred diffusion
coefficient of \textit{bicoid} mRNA is in the range 4.6\times
10^{-12}-1.5\times 10^{-11}m^2s^{-1}, in agreement with the theoretical
predictions and experimental measurements for the diffusion of macromolecules
in the cytoplasm. We show that the model based on the mRNA diffusion hypothesis
is consistent with the known observational data, supporting the recent
experimental findings of the gradient of \textit{bicoid} mRNA in
\textit{Drosophila} [Spirov et al. (2009) Development 136:605-614].


  We present a general methodology in order to build mathematical models of
genetic regulatory networks. This approach is based on the mass action law and
on the Jacob and Monod operon model. The mathematical models are built
symbolically by the \emph{Mathematica} software package \emph{GeneticNetworks}.
This package accepts as input the interaction graphs of the transcriptional
activators and repressors and, as output, gives the mathematical model in the
form of a system of ordinary differential equations. All the relevant
biological parameters are chosen automatically by the software. Within this
framework, we show that threshold effects in biology emerge from the catalytic
properties of genes and its associated conservation laws. We apply this
methodology to the segment patterning in \emph{Drosophila} early development
and we calibrate and validate the genetic transcriptional network responsible
for the patterning of the gap proteins Hunchback and Knirps, along the
antero-posterior axis of the \emph{Drosophila} embryo. This shows that
patterning at the gap genes stage is a consequence of the relations between the
transcriptional regulators and their initial conditions along the embryo.


  The standard genetic code is known to be robust to translation errors and
point mutations. We studied how small modifications of the standard code affect
its robustness. The robustness was assessed in terms of a proper stability
function, the negative variations of which correspond to a more robust code.
The fraction of more robust codes obtained under small modifications appeared
to be unexpectedly high, ranged from 0.1 to 0.4 depending on the choice of
stability function and code modifications, yet significantly lower than the
corresponding fraction in the random codes (about a half). In this sense the
standard code ought to be considered distinctly non-random in accordance with
previous observations. The distribution of the negative variations of stability
function revealed very abrupt drop beyond one standard deviation, much sharper
than for Gaussian distribution or for the random codes with the same number of
codons in the sets coding for amino acids or stop-codons. This behavior holds
for both the standard code as a whole and its binary NRN-NYN, NWN-NSN, and
NMN-NKN blocks. Previously, it has been proved that such binary block structure
is necessary for the robustness of a code and is inherent to the standard
genetic code. The modifications of the standard code corresponding to more
robust coding may be related to the different variants of the code. These
effects may also contribute to the rates of replacements of amino acids. The
observed features demonstrate the joint impact of random factors and natural
selection during evolution of the genetic code.


  We discuss a realistic scenario, accounting for the existence of
sub-micrometric protein domains in cell membranes. At the biological level,
such membrane domains have been shown to be specialized, in order to perform a
determined biological task, in the sense that they gather one or a few protein
species out of the hundreds of different ones that a cell membrane may contain.
By analyzing the balance between mixing entropy and protein affinities, we
propose that such protein sorting in distinct domains can be explained without
appealing to pre-existing lipidic micro-phase separations, as in the lipid raft
scenario. We show that the proposed scenario is compatible with known physical
interactions between membrane proteins, even if thousands of different species
coexist.


  Recent work by Erik Volz [arXiv:0705.2092] has shown how to calculate the
growth and eventual decay of an SIR epidemic on a static random network,
assuming infection and recovery each happen at constant rates. This calculation
allows us to account for effects due to heterogeneity in degree that are
neglected in the standard mass-action SIR equations. In this note we offer an
alternate derivation which arrives at a simpler -- though equivalent -- system
of governing equations to that of Volz. This new derivation is more closely
connected to the underlying physical processes, and the resulting equations are
of comparable complexity to the mass-action SIR equations. We further show that
earlier derivations of the final size of epidemics on networks can be
reproduced using the same approach, thereby providing a common framework for
calculating both the dynamics and the final size of an epidemic spreading on a
random network.


  Metamorphic proteins like Lymphotactin are a notable exception of the
empirical principle that structured natural proteins possess a unique three
dimensional structure. In particular, the human chemokine lymphotactin protein
(Ltn) exists in two distinct conformations (one monomeric and one dimeric)
under physiological conditions. In this work we use a Ca Go-model to show how
this very peculiar behavior can be reproduced. From the study of the
thermodynamics and of the kinetics we characterize the interconversion
mechanism. In particular, this takes place through the docking of the two
chains living in a third monomeric, partially unfolded, state which shows a
residual structure involving a set of local contacts common to the two native
conformations. The main feature of two-fold proteins appears to be the sharing
of a common set of local contacts between the two distinct folds as confirmed
by the study of two designed two-fold proteins. Metamorphic proteins may be
more common than expected.


  Biological cells are able to accurately sense chemicals with receptors at
their surfaces, allowing cells to move towards sources of attractant and away
from sources of repellent. The accuracy of sensing chemical concentration is
ultimately limited by the random arrival of particles at the receptors by
diffusion. This fundamental physical limit is generally considered to be the
Berg & Purcell limit [H.C. Berg and E.M. Purcell, Biophys. J. {\bf 20}, 193
(1977)]. Here we derive a lower limit by applying maximum likelihood to the
time series of receptor occupancy. The increased accuracy stems from solely
considering the unoccupied time intervals - disregarding the occupied time
intervals as these do not contain any information about the external particle
concentration, and only decrease the accuracy of the concentration estimate.
Receptors which minimize the bound time intervals achieve the highest possible
accuracy. We discuss how a cell could implement such an optimal sensing
strategy by absorbing or degrading bound particles.


  Bundles of filamentous polymers are primary structural components of a broad
range of cytoskeletal structures, and their mechanical properties play key
roles in cellular functions ranging from locomotion to mechanotransduction and
fertilization. We give a detailed derivation of a wormlike bundle model as a
generic description for the statics and dynamics of polymer bundles consisting
of semiflexible polymers interconnected by crosslinking agents. The elastic
degrees of freedom include bending as well as twist deformations of the
filaments and shear deformation of the crosslinks. We show that a competition
between the elastic properties of the filaments and those of the crosslinks
leads to renormalized effective bend and twist rigidities that become
mode-number dependent. The strength and character of this dependence is found
to vary with bundle architecture, such as the arrangement of filaments in the
cross section and pretwist. We discuss two paradigmatic cases of bundle
architecture, a uniform arrangement of filaments as found in F-actin bundles
and a shell-like architecture as characteristic for microtubules. Each
architecture is found to have its own universal ratio of maximal to minimal
bending rigidity, independent of the specific type of crosslink induced
filament coupling; our predictions are in reasonable agreement with available
experimental data for microtubules. Moreover, we analyze the predictions of the
wormlike bundle model for experimental observables such as the tangent-tangent
correlation function and dynamic response and correlation functions. Finally,
we analyze the effect of pretwist (helicity) on the mechanical properties of
bundles. We predict that microtubules with different number of protofilaments
should have distinct variations in their effective bending rigidity.


  A bacterial colony may develop a small number of cells genetically identical
to, but phenotypically different from other normally growing bacteria. These
so-called persister cells keep themselves in a dormant state and thus are
insensitive to antibiotic treatment, resulting in serious problems of drug
resistance. In this paper, we proposed a novel strategy to "kill" persister
cells by triggering them to switch, in a fast and synchronized way, into
normally growing cells that are susceptible to antibiotics. The strategy is
based on resonant activation (RA), a well-studied phenomenon in physics where
the internal noise of a system can constructively facilitate fast and
synchronized barrier crossings. Through stochastic Gilliespie simulation with a
generic toggle switch model, we demonstrated that RA exists in the phenotypic
switching of a single bacterium. Further, by coupling single cell level and
population level simulations, we showed that with RA, one can greatly reduce
the time and total amount of antibiotics needed to sterilize a bacterial
population. We suggest that resonant activation is a general phenomenon in
phenotypic transition, and can find other applications such as cancer therapy.


  We introduce a framework to investigate ab-initio the dynamics of rare
thermally activated reactions. The electronic degrees of freedom are described
at the quantum-mechanical level in the Born-Oppenheimer approximation, while
the nuclear degrees of freedom are coupled to a thermal bath, through a
Langevin equation. This method is based on the path integral representation for
the stochastic dynamics and yields the time evolution of both nuclear and
electronic degrees of freedom, along the most probable reaction pathways,
without spending computational time to explore metastable states. This approach
is very efficient and allows to study thermally activated reactions which
cannot be simulated using ab-initio molecular dynamics techniques. As a first
illustrative application, we characterize the dominant pathway in the
cyclobutene to butadiene reaction.


  We analyze large systems of interacting proteins, using techniques from the
non-equilibrium statistical mechanics of disordered many-particle systems.
Apart from protein production and removal, the most relevant microscopic
processes in the proteome are complex formation and dissociation, and the
microscopic degrees of freedom are the evolving concentrations of unbound
proteins (in multiple post-translational states) and of protein complexes. Here
we only include dimer-complexes, for mathematical simplicity, and we draw the
network that describes which proteins are reaction partners from an ensemble of
random graphs with an arbitrary degree distribution. We show how generating
functional analysis methods can be used successfully to derive closed equations
for dynamical order parameters, representing an exact macroscopic description
of the complex formation and dissociation dynamics in the infinite system
limit. We end this paper with a discussion of the possible routes towards
solving the nontrivial order parameter equations, either exactly (in specific
limits) or approximately.


  Although widely used in practice, the behavior and accuracy of the popular
module identification technique called modularity maximization is not well
understood in practical contexts. Here, we present a broad characterization of
its performance in such situations. First, we revisit and clarify the
resolution limit phenomenon for modularity maximization. Second, we show that
the modularity function Q exhibits extreme degeneracies: it typically admits an
exponential number of distinct high-scoring solutions and typically lacks a
clear global maximum. Third, we derive the limiting behavior of the maximum
modularity Q_max for one model of infinitely modular networks, showing that it
depends strongly both on the size of the network and on the number of modules
it contains. Finally, using three real-world metabolic networks as examples, we
show that the degenerate solutions can fundamentally disagree on many, but not
all, partition properties such as the composition of the largest modules and
the distribution of module sizes. These results imply that the output of any
modularity maximization procedure should be interpreted cautiously in
scientific contexts. They also explain why many heuristics are often successful
at finding high-scoring partitions in practice and why different heuristics can
disagree on the modular structure of the same network. We conclude by
discussing avenues for mitigating some of these behaviors, such as combining
information from many degenerate solutions or using generative models.


  We consider an asexual biological population of constant size $N$ evolving in
discrete time under the influence of selection and mutation. Beneficial
mutations appear at rate $U$ and their selective effects $s$ are drawn from a
distribution $g(s)$. After introducing the required models and concepts of
mathematical population genetics, we review different approaches to computing
the speed of logarithmic fitness increase as a function of $N$, $U$ and $g(s)$.
We present an exact solution of the infinite population size limit and provide
an estimate of the population size beyond which it is valid. We then discuss
approximate approaches to the finite population problem, distinguishing between
the case of a single selection coefficient, $g(s) = \delta(s - s_b)$, and a
continuous distribution of selection coefficients. Analytic estimates for the
speed are compared to numerical simulations up to population sizes of order
$10^{300}$.


  We review enzymatic systems which involve biocatalytic reactions utilized for
information processing (biocomputing). Extensive ongoing research in
biocomputing, mimicking Boolean logic gates has been motivated by potential
applications in biotechnology and medicine. Furthermore, novel sensor concepts
have been contemplated with multiple inputs processed biochemically before the
final output is coupled to transducing "smart-material" electrodes and other
systems. These applications have warranted recent emphasis on networking of
biocomputing gates. First few-gate networks have been experimentally realized,
including coupling, for instance, to signal-responsive electrodes for signal
readout. In order to achieve scalable, stable network design and functioning,
considerations of noise propagation and control have been initiated as a new
research direction. Optimization of single enzyme-based gates for avoiding
analog noise amplification has been explored, as were certain
network-optimization concepts. We review and exemplify these developments, as
well as offer an outlook for possible future research foci. The latter include
design and uses of non-Boolean network elements, e.g., filters, as well as
other developments motivated by potential novel sensor and biotechnology
applications.


  We investigate how perturbations propagate up and down a food chain with and
without self-interaction and omnivory. A source of perturbation is a shift in
death rate of a trophic level, and the measure of perturbation is the
difference between the perturbed and unperturbed steady state populations. For
Lotka-Volterra food chains with linear functional response, we show
analytically that both intraspecific competition and intraguild predation can
either dampen or enhance the propagation of perturbations, thus stabilizing or
destabilizing the food web. The direction of the effect depend on the position
of the source of perturbation, as well as on the position of the additional
competitive and predatory links . These conclusions are confirmed numerically
for a food chain with more realistic Type-II functional response. Our results
support the positions of both sides in the long-standing debate on the effect
of intraspecific competition and omnivory on the stability of trophic systems.


  Concerns about biodiversity and the long-term stability of forest ecosystems
have lead to changing attitudes with respect to plantations. These artificial
communities are ubiquitous, yet provide reduced habitat value in comparison to
their naturally established counterparts, key factors being high density,
homogeneous spatial structure, and their even-sized/aged nature.
  Transformation (manipulation of plantations to produce stands more
reminiscent of natural ones) represents a major challenge for forest managers,
and the shift from even- to uneven-aged stand management is far from simple.
Tree species' attributes, too, vary dramatically. This study aims to identify
generic aspects of forest population dynamics, in order to understand the
temporal evolution of single species forest stands from an initial plantation
lattice structure, into an ``old-growth'' state.
  This is achieved using a size-structured model, which implements simple rules
to simulate the growth and fate of individuals in a spatial arena. We
parameterise and discuss model behaviour in the context of available plantation
and semi-natural data from Caledonian Scots Pine stands. We move on to
illustrate the use of the model in predicting the outcome of silviculture, and
discuss results in the context of stand transformation.


  The presence of clusters of rare codons is known to negatively impact the
efficiency and accuracy of protein production. In this paper, we demonstrate a
statistical method of identifying such clusters in the coding sequence of a
gene. Using E. coli as our model organism, we show that genes having denser
clusters tend to have lower protein yields.