  Quantifying interactions in DNA microarrays is of central importance for a
better understanding of their functioning. Hybridization thermodynamics for
nucleic acid strands in aqueous solution can be described by the so-called
nearest-neighbor model, which estimates the hybridization free energy of a
given sequence as a sum of dinucleotide terms. Compared with its solution
counterparts, hybridization in DNA microarrays may be hindered due to the
presence of a solid surface and of a high density of DNA strands. We present
here a study aimed at the determination of hybridization free energies in DNA
microarrays. Experiments are performed on custom Agilent slides. The solution
contains a single oligonucleotide. The microarray contains spots with a perfect
matching complementary sequence and other spots with one or two mismatches: in
total 1006 different probe spots, each replicated 15 times per microarray. The
free energy parameters are directly fitted from microarray data. The
experiments demonstrate a clear correlation between hybridization free energies
in the microarray and in solution. The experiments are fully consistent with
the Langmuir model at low intensities, but show a clear deviation at
intermediate (non-saturating) intensities. These results provide new
interesting insights for the quantification of molecular interactions in DNA
microarrays.


  The Electron Microscopy Data Bank (EMDB) is a rapidly growing repository for
the dissemination of structural data from single-particle reconstructions of
supramolecular protein assemblies including motors, chaperones, cytoskeletal
assemblies, and viral capsids. While the static structure of these assemblies
provides essential insight into their biological function, their conformational
dynamics and mechanics provide additional important information regarding the
mechanism of their biological function. Here, we present an unsupervised
computational framework to analyze and store for public access the
conformational dynamics of supramolecular protein assemblies deposited in the
EMDB. Conformational dynamics are analyzed using normal mode analysis in the
finite element framework, which is used to compute equilibrium thermal
fluctuations, cross-correlations in molecular motions, and strain energy
distributions for 452 of the 681 entries stored in the EMDB at present. Results
for the viral capsid of hepatitis B, ribosome-bound termination factor RF2, and
GroEL are presented in detail and validated with all-atom based models. The
conformational dynamics of protein assemblies in the EMDB may be useful in the
interpretation of their biological function, as well as in the classification
and refinement of EM-based structures.


  Most bacteriophages are known to inject their double-stranded DNA into
bacteria upon receptor binding in an essentially spontaneous way. This downhill
thermodynamic process from the intact virion toward the empty viral capsid plus
released DNA is made possible by the energy stored during active packaging of
the genome into the capsid. Only indirect measurements of this energy have been
available until now using either single-molecule or osmotic suppression
techniques. In this paper, we describe for the first time the use of isothermal
titration calorimetry to directly measure the heat released (or equivalently
the enthalpy) during DNA ejection from phage lambda, triggered in solution by a
solubilized receptor. Quantitative analyses of the results lead to the
identification of thermodynamic determinants associated with DNA ejection. The
values obtained were found to be consistent with those previously predicted by
analytical models and numerical simulations. Moreover, the results confirm the
role of DNA hydration in the energetics of genome confinement in viral capsids.


  The Gene Ontology (GO) provides a knowledge base to effectively describe
proteins. However, measuring similarity between proteins based on GO remains a
challenge. In this paper, we propose a new similarity measure, information
coefficient similarity measure (SimIC), to effectively integrate both the
information content (IC) of GO terms and the structural information of GO
hierarchy to determine the similarity between proteins. Testing on yeast
proteins, our results show that SimIC efficiently addresses the shallow
annotation issue in GO, thus improves the correlations between GO similarities
of yeast proteins and their expression similarities as well as between GO
similarities of yeast proteins and their sequence similarities. Furthermore, we
demonstrate that the proposed SimIC is superior in predicting yeast protein
interactions. We predict 20484 yeast protein-protein interactions (PPIs)
between 2462 proteins based on the high SimIC values of biological process (BP)
and cellular component (CC). Examining the 214 MIPS complexes in our predicted
PPIs shows that all members of 159 MIPS complexes can be found in our PPI
predictions, which is more than those (120/214) found in PPIs predicted by
relative specificity similarity (RSS). Integrating IC and structural
information of GO hierarchy can improve the effectiveness of the semantic
similarity measure of GO terms. The new SimIC can effectively correct the
effect of shallow annotation, and then provide an effective way to measure
similarity between proteins based on Gene Ontology.


  Quorum sensing is the mechanism by which bacteria communicate and synchronize
group behaviors. Quantitative information on parameters such as the copy number
of particular quorum-sensing proteins should contribute strongly to
understanding how the quorum-sensing network functions. Here we show that the
copy number of the master regulator protein LuxR in Vibrio harveyi, can be
determined in vivo by exploiting small-number fluctuations of the protein
distribution when cells undergo division. When a cell divides, both its volume
and LuxR protein copy number N are partitioned with slight asymmetries. We have
measured the distribution functions describing the partitioning of the protein
fluorescence and the cell volume. The fluorescence distribution is found to
narrow systematically as the LuxR population increases while the volume
partitioning is unchanged. Analyzing these changes statistically, we have
determined that N = 80-135 dimers at low cell density and 575 dimers at high
cell density. In addition, we have measured the static distribution of LuxR
over a large (3,000) clonal population. Combining the static and time-lapse
experiments, we determine the magnitude of the Fano factor of the distribution.
This technique has broad applicability as a general, in vivo technique for
measuring protein copy number and burst size.


  The significant role of space in maintaining species coexistence and
determining community structure and function is well established. However,
community ecology studies have mainly focused on simple competition and
predation systems, and the relative impact of positive interspecific
interactions in shaping communities in a spatial context is not well
understood. Here we employ a spatially explicit metacommunity model to
investigate the effect of local dispersal on the structure and function of
communities in which species are linked through an interaction web comprising
mutualism, competition and exploitation. Our results show that function,
diversity and interspecific interactions of locally linked communities undergo
a phase transition with changes in the rate of species dispersal. We find that
low spatial interconnectedness favors the spontaneous emergence of strongly
mutualistic communities which are more stable but less productive and diverse.
On the other hand, high spatial interconnectedness promotes local biodiversity
at the expense of local stability and supports communities with a wide range of
interspecific interactions. We argue that investigations of the relationship
between spatial processes and the self-organization of complex interaction webs
are critical to understanding the geographic structure of interactions in real
landscapes.


  An approach for simulating bionanosystems, such as viruses and ribosomes, is
presented. This calibration-free approach is based on an all-atom description
for bionanosystems, a universal interatomic force field, and a multiscale
perspective. The supramillion-atom nature of these bionanosystems prohibits the
use of a direct molecular dynamics approach for phenomena like viral structural
transitions or self-assembly that develop over milliseconds or longer. A key
element of these multiscale systems is the cross-talk between, and consequent
strong coupling of, processes over many scales in space and time. We elucidate
the role of interscale cross-talk and overcome bionanosystem simulation
difficulties with automated construction of order parameters (OPs) describing
supra-nanometer scale structural features, construction of OP dependent
ensembles describing the statistical properties of atomistic variables that
ultimately contribute to the entropies driving the dynamics of the OPs, and the
derivation of a rigorous equation for the stochastic dynamics of the OPs. Since
the atomic scale features of the system are treated statistically, several
ensembles are constructed that reflect various experimental conditions. The
theory provides a basis for a practical, quantitative bionanosystem modeling
approach that preserves the cross-talk between the atomic and nanoscale
features. A method for integrating information from nanotechnical experimental
data in the derivation of equations of stochastic OP dynamics is also
introduced.


  We consider DNA codes based on the nearest-neighbor (stem) similarity model
which adequately reflects the "hybridization potential" of two DNA sequences.
Our aim is to present a survey of bounds on the rate of DNA codes with respect
to a thermodynamically motivated similarity measure called an additive stem
similarity. These results yield a method to analyze and compare known samples
of the nearest neighbor "thermodynamic weights" associated to stacked pairs
that occurred in DNA secondary structures.


  The evolutionary effect of recombination depends crucially on the epistatic
interactions between linked loci. A paradigmatic case where recombination is
known to be strongly disadvantageous is a two-locus fitness landscape dis-
playing reciprocal sign epistasis with two fitness peaks of unequal height.
Focusing on the occurrence of bistability in the equilibrium solutions, we
consider here the deterministic, haploid two-locus model with reversible mu-
tations, selection and recombination. We find analytic formulae for the criti-
cal recombination probability rc above which two stable stationary solutions
appear which are localized on each of the two fitness peaks. We also derive the
stationary genotype frequencies in various parameter regimes. When the
recombination rate is close to rc and the fitness difference between the two
peaks is small, we obtain a compact description in terms of a cubic polyno-
mial which is analogous to the Landau theory of physical phase transitions.


  Here we study the effects of adopting different strategies against different
opponent instead of adopting the same strategy against all of them in the
prisoner dilemma structured in well-mixed populations. We consider an
evolutionary process in which strategies that provide reproductive success are
imitated and players replace one of their worst interactions by the new one. We
set individuals in a well-mixed population so that network reciprocity effect
is excluded and we analyze both synchronous and asynchronous updates. As a
consequence of the replacement rule, we show that mutual cooperation is never
destroyed and the initial fraction of mutual cooperation is a lower bound for
the level of cooperation. We show by simulation and mean-field analysis that
for synchronous update cooperation dominates while for asynchronous update only
cooperations associated to the initial mutual cooperations are maintained. As a
side effect of the replacement rule, an "implicit punishment" mechanism comes
up in a way that exploitations are always neutralized providing evolutionary
stability for cooperation.


  It is proposed that the ability of humans to flourish in diverse environments
and evolve complex cultures reflects the following two underlying cognitive
transitions. The transition from the coarse-grained associative memory of Homo
habilis to the fine-grained memory of Homo erectus enabled limited
representational redescription of perceptually similar episodes, abstraction,
and analytic thought, the last of which is modeled as the formation of states
and of lattices of properties and contexts for concepts. The transition to the
modern mind of Homo sapiens is proposed to have resulted from onset of the
capacity to spontaneously and temporarily shift to an associative mode of
thought conducive to interaction amongst seemingly disparate concepts, modeled
as the forging of conjunctions resulting in states of entanglement. The fruits
of associative thought became ingredients for analytic thought, and vice versa.
The ratio of associative pathways to concepts surpassed a percolation threshold
resulting in the emergence of a self-modifying, integrated internal model of
the world, or worldview.


  In the context of the current financial crisis, when more companies are
facing bankruptcy or insolvency, the paper aims to find methods to identify
distressed firms by using financial ratios. The study will focus on identifying
a group of Romanian listed companies, for which financial data for the year
2008 were available. For each company a set of 14 financial indicators was
calculated and then used in a principal component analysis, followed by a
cluster analysis, a logit model, and a CHAID classification tree.


  A critical task in systems biology is the identification of genes that
interact to control cellular processes by transcriptional activation of a set
of target genes. Many methods have been developed to use statistical
correlations in high-throughput datasets to infer such interactions. However,
cellular pathways are highly cooperative, often requiring the joint effect of
many molecules, and few methods have been proposed to explicitly identify such
higher-order interactions, partially due to the fact that the notion of
multivariate statistical dependency itself remains imprecisely defined. We
define the concept of dependence among multiple variables using maximum entropy
techniques and introduce computational tests for their identification.
Synthetic network results reveal that this procedure uncovers dependencies even
in undersampled regimes, when the joint probability distribution cannot be
reliably estimated. Analysis of microarray data from human B cells reveals that
third-order statistics, but not second-order ones, uncover relationships
between genes that interact in a pathway to cooperatively regulate a common set
of targets.


  We investigate the stability properties of two different classes of metabolic
cycles using a combination of analytical and computational methods. Using
principles from structural kinetic modeling (SKM), we show that the stability
of metabolic networks with certain structural regularities can be studied using
a combination of analytical and computational techniques. We then apply these
techniques to a class of single input, single output metabolic cycles, and find
that the cycles are stable under all conditions tested. Next, we extend our
analysis to a small autocatalytic cycle, and determine parameter regimes within
which the cycle is very likely to be stable. We demonstrate that analytical
methods can be used to understand the relationship between kinetic parameters
and stability, and that results from these analytical methods can be confirmed
with computational experiments. In addition, our results suggest that elevated
metabolite concentrations and certain crucial saturation parameters can
strongly affect the stability of the entire metabolic cycle. We discuss our
results in light of the possibility that evolutionary forces may select for
metabolic network topologies with a high intrinsic probability of being stable.
Furthermore, our conclusions support the hypothesis that certain types of
metabolic cycles may have played a role in the development of primitive
metabolism despite the absence of regulatory mechanisms.


  In computational molecular biology, gene regulatory binding sites prediction
in whole genome remains a challenge for the researchers. Now a days, the genome
wide regulatory binding site prediction tools required either direct pattern
sequence or weight matrix. Although there are known transcription factor
binding sites databases available for genome wide prediction but no tool is
available which can construct different weight matrices as per need of user or
tools available for large data set scanning by first aligning the input
upstream or promoter sequences and than construct the matrices in different
level and file format. Considering this, we developed a DNA MATRIX tool for
searching putative regulatory binding sites in gene upstream sequences. This
tool uses the simple biological rule based heuristic algorithm for weight
matrix construction, which can be transformed into different formats after
motif alignment and therefore provides the possibility to identify the most
potential conserved binding sites in the regulated genes. The user may
construct and save specific weight or frequency matrices in different form and
file formats based on user based selection of conserved aligned block of short
sequences ranges from 6 to 20 base pairs and prior nucleotide frequency before
weight scoring.


  Over the years, data mining has attracted most of the attention from the
research community. The researchers attempt to develop faster, more scalable
algorithms to navigate over the ever increasing volumes of spatial gene
expression data in search of meaningful patterns. Association rules are a data
mining technique that tries to identify intrinsic patterns in spatial gene
expression data. It has been widely used in different applications, a lot of
algorithms introduced to discover these rules. However Priori like algorithms
has been used to find positive association rules. In contrast to positive
rules, negative rules encapsulate relationship between the occurrences of one
set of items with absence of the other set of items. In this paper, an
algorithm for mining negative association rules from spatial gene expression
data is introduced. The algorithm intends to discover the negative association
rules which are complementary to the association rules often generated by
Priori like algorithm. Our study shows that negative association rules can be
discovered efficiently from spatial gene expression data.


  Test experiments of hybridization in DNA microarrays show systematic
deviations from the equilibrium isotherms. We argue that these deviations are
due to the presence of a partially hybridized long-lived state, which we
include in a kinetic model. Experiments confirm the model predictions for the
intensity vs. free energy behavior. The existence of slow relaxation phenomena
has important consequences for the specificity of microarrays as devices for
the detection of a target sequence from a complex mixture of nucleic acids.


  Bayesian inference is now a leading technique for reconstructing phylogenetic
trees from aligned sequence data. In this short note, we formally show that the
maximum posterior tree topology provides a statistically consistent estimate of
a fully-resolved evolutionary tree under a wide variety of conditions. This
includes the inference of gene trees from aligned sequence data across the
entire parameter range of branch lengths, and under general conditions on
priors in models where the usual `identifiability' conditions hold. We extend
this to the inference of species trees from sequence data, where the gene trees
constitute `nuisance parameters', as in the program *BEAST. This note also
addresses earlier concerns raised in the literature questioning the extent to
which statistical consistency for Bayesian methods might hold in general.


  Motivation: Protein sequence world is discrete as 20 amino acids (AA) while
its structure world is continuous, though can be discretized into structural
alphabets (SA). In order to reveal the relationship between sequence and
structure, it is interesting to consider both AA and SA in a joint space.
However, such space has too many parameters, so the reduction of AA is
necessary to bring down the parameter numbers.
  Result: We've developed a simple but effective approach called entropic
clustering based on selecting the best mutual information between a given
reduction of AAs and SAs. The optimized reduction of AA into two groups leads
to hydrophobic and hydrophilic. Combined with our SA, namely conformational
letter (CL) of 17 alphabets, we get a joint alphabet called hydropathy
conformational letter (hp-CL). A joint substitution matrix with (17*2)*(17*2)
indices is derived from FSSP. Moreover, we check the three coding systems, say
AA, CL and hp-CL against a large database consisting proteins from family to
fold, with their performance on the TopK accuracy of both similar fragment pair
(SFP) and the neighbor of aligned fragment pair (AFP). The TopK selection is
according to the score calculated by the coding system's substitution matrix.
Finally, embedding hp-CL in a pairwise alignment algorithm, say CLeFAPS, to
replace the original CL, will get an improvement on the HOMSTRAD benchmark.


  RNA viruses exist in large intra-host populations which display great
genotypic and phenotypic diversity. We analyze a model of viral competition
between two different viral strains infecting a constantly replenished cell
pool, in which we assume a trade-off between the virus' colonization skills
(cell killing ability or virulence) and its local competition skills
(replication performance within coinfected cells). We characterize the
conditions that allow for viral spread by means of the basic reproductive
number and show that a local coexistence equilibrium exists, which is
asymptotically stable. At this equilibrium, the less virulent competitor has a
reproductive advantage over the more virulent colonizer. The equilibria at
which one strain outcompetes the other one are unstable, i.e., a second viral
strain is always able to permanently invade. One generalization of the model is
to consider multiple viral strains, each one displaying a different virulence.
However, to account for the large phenotypic diversity in viral populations, we
consider a continuous spectrum of virulences and present a continuum limit of
this multiple viral strains model that describes the time evolution of an
initial continuous distribution of virulence. We provide a proof of the
existence of solutions of the model's equations and present numerical
approximations of solutions for different initial distributions. Our
simulations suggest that initial continuous distributions of virulence evolve
towards a stationary distribution that is extremely skewed in favor of
competitors. Consequently, collective virulence attenuation takes place. This
finding may contribute to understanding the phenomenon of virulence
attenuation, which has been reported in previous experimental studies.


  Motivation : Molecular signatures for diagnosis or prognosis estimated from
large-scale gene expression data often lack robustness and stability, rendering
their biological interpretation challenging. Increasing the signature's
interpretability and stability across perturbations of a given dataset and, if
possible, across datasets, is urgently needed to ease the discovery of
important biological processes and, eventually, new drug targets. Results : We
propose a new method to construct signatures with increased stability and
easier interpretability. The method uses a gene network as side interpretation
and enforces a large connectivity among the genes in the signature, leading to
signatures typically made of genes clustered in a few subnetworks. It combines
the recently proposed graph Lasso procedure with a stability selection
procedure. We evaluate its relevance for the estimation of a prognostic
signature in breast cancer, and highlight in particular the increase in
interpretability and stability of the signature.


  The incidence of Dengue epidemiologic disease has grown in recent decades. In
this paper an application of optimal control in Dengue epidemics is presented.
The mathematical model includes the dynamic of Dengue mosquito, the affected
persons, the people's motivation to combat the mosquito and the inherent social
cost of the disease, such as cost with ill individuals, educations and sanitary
campaigns. The dynamic model presents a set of nonlinear ordinary differential
equations. The problem was discretized through Euler and Runge Kutta schemes,
and solved using nonlinear optimization packages. The computational results as
well as the main conclusions are shown.


  A 62 Myr periodicity is superimposed on other longer-term trends in fossil
biodiversity. This cycle can be discerned in marine data based on the Sepkoski
compendium, the Paleobiology Database, and the Fossil Record 2. The signal also
exists in changes in sea level/sediment, but is much weaker than in
biodiversity itself. A significant excess of 19 previously identified
Phanerozoic mass extinctions occur on the declining phase of the 62 Myr cycle.
appearance of the signal in sampling-standardized biodiversity data, it is
likely not to be a sampling artifact, but either a consequence of sea-level
changes or an additional effect of some common cause for them both. In either
case, it is intriguing why both changes would have a regular pattern.


  Recent work has highlighted deep connections between sequence-length
requirements for high-probability phylogeny reconstruction and the related
problem of the estimation of ancestral sequences. In [Daskalakis et al.'09],
building on the work of [Mossel'04], a tight sequence-length requirement was
obtained for the CFN model. In particular the required sequence length for
high-probability reconstruction was shown to undergo a sharp transition (from
$O(\log n)$ to $\hbox{poly}(n)$, where $n$ is the number of leaves) at the
"critical" branch length $\critmlq$ (if it exists) of the ancestral
reconstruction problem.
  Here we consider the GTR model. For this model, recent results of [Roch'09]
show that the tree can be accurately reconstructed with sequences of length
$O(\log(n))$ when the branch lengths are below $\critksq$, known as the
Kesten-Stigum (KS) bound. Although for the CFN model $\critmlq = \critksq$, it
is known that for the more general GTR models one has $\critmlq \geq \critksq$
with a strict inequality in many cases. Here, we show that this phenomenon also
holds for phylogenetic reconstruction by exhibiting a family of symmetric
models $Q$ and a phylogenetic reconstruction algorithm which recovers the tree
from $O(\log n)$-length sequences for some branch lengths in the range
$(\critksq,\critmlq)$. Second we prove that phylogenetic reconstruction under
GTR models requires a polynomial sequence-length for branch lengths above
$\critmlq$.


  The search for interesting association rules is an important topic in
knowledge discovery in spatial gene expression databases. The set of admissible
rules for the selected support and confidence thresholds can easily be
extracted by algorithms based on support and confidence, such as Apriori.
However, they may produce a large number of rules, many of them are
uninteresting. The challenge in association rule mining (ARM) essentially
becomes one of determining which rules are the most interesting. Association
rule interestingness measures are used to help select and rank association rule
patterns. Besides support and confidence, there are other interestingness
measures, which include generality reliability, peculiarity, novelty,
surprisingness, utility, and applicability. In this paper, the application of
the interesting measures entropy and variance for association pattern discovery
from spatial gene expression data has been studied. In this study the fast
mining algorithm has been used which produce candidate itemsets and it spends
less time for calculating k-supports of the itemsets with the Boolean matrix
pruned, and it scans the database only once and needs less memory space.
Experimental results show that using entropy as the measure of interest for the
spatial gene expression data has more diverse and interesting rules.


  The dynamical characterization of proteins is crucial to understand protein
function. From a microscopic point of view, protein dynamics is governed by the
local atomic interactions that, in turn, trigger the functional conformational
changes. Unfortunately, the relationship between local atomic fluctuations and
global protein rearrangements is still elusive. Here, atomistic molecular
dynamics simulations in conjunction with complex network analysis show that
fast peptide relaxations effectively build the backbone of the global
free-energy landscape, providing a connection between local and global atomic
rearrangements. A minimum-spanning-tree representation, built on the base of
transition gradients networks, results in a high resolution mapping of the
system dynamics and thermodynamics without requiring any a priori knowledge of
the relevant degrees of freedom. These results suggest the presence of a local
mechanism for the high communication efficiency generally observed in complex
systems.


  There has been concern in the literature about the methodology of using
secondary calibration timepoints when estimating evolutionary divergence dates.
Such timepoints are divergence time estimates that have been derived from one
molecular data set on the basis of a primary external calibration timepoint,
and which are then used independently on a second data set. Logically, the
primary and secondary calibration points must be mutually consistent, in the
sense that it must be possible to predict each time point from the other.
However, the attempt by Shaul and Graur (2002, Gene 300: 59-61) to assess the
reliability of secondary timepoints is flawed because they presented time
estimates without presenting confidence intervals on those estimates, and so it
was not possible to make any explicit hypothesis tests of divergence times.
Also, they inappropriately excluded some of the data, which leads to a very
biased estimate of one of the divergence times. Here, I present a re-analysis
of the same data set, with more appropriate methodology, and come to the
conclusion that no inconsistencies are involved. However, it is clear from the
analysis that molecular data often have such large confidence intervals that
they are uninformative, and thus cannot be used for reliable hypothesis tests.


  Huelsenbeck and Rannala (2004, Systematic Biology 53, 904-913) presented a
series of simulations in order to assess the extent to which the bayesian
posterior probabilities associated with phylogenetic trees represent the
standard frequentist statistical interpretation. They concluded that when the
analysis model matches the generating model then the bayesian posterior
probabilities are correct, but that the probabilities are much too large when
the model is under-specified and slightly too small when the model is
over-specified. Here, I take issue with the first conclusion, and instead
contend that their simulation data show that the posterior probabilities are
still slightly too large even when the models match. Furthermore, I suggest
that the data show that the degree of this over-estimation increases as the
sequence length increases, and that it might increase as model complexity
increases. I also provide some comments on the authors' conclusions concerning
whether bootstrap proportions over- or under-estimate the true probabilities.


  Evolutionary game dynamics describes not only frequency dependent genetical
evolution, but also cultural evolution in humans. In this context, successful
strategies spread by imitation. It has been shown that the details of strategy
update rules can have a crucial impact on evolutionary dynamics in theoretical
models and e.g. significantly alter the level of cooperation in social
dilemmas. But what kind of strategy update rules can describe imitation
dynamics in humans? Here, we present a way to measure such strategy update
rules in a behavioral experiment. We use a setting in which individuals are
virtually arranged on a spatial lattice. This produces a large number of
different strategic situations from which we can assess strategy updating. Most
importantly, spontaneous strategy changes corresponding to mutations or
exploration behavior are more frequent than assumed in many models. Our
experimental approach to measure properties of the update mechanisms used in
theoretical models will be useful for mathematical models of cultural
evolution.


  The ability to locate and annotate mitochondrial genes is an important
practical issue, given the rapidly increasing number of mitogenomes appearing
in the public databases. Unfortunately, tRNA genes in Metazoan mitochondria
have proved to be problematic because they often vary in number (genes missing
or duplicated) and also in the secondary structure of the transcribed tRNAs (T
or D arms missing). I have performed a series of comparative analyses of the
tRNA genes of a broad range of Metazoan mitogenomes in order to address this
issue. I conclude that no single computer program is necessarily capable of
finding all of the tRNA genes in any given mitogenome, and that use of both the
ARWEN and DOGMA programs is sometimes necessary because they produce
complementary false negatives. There are apparently a very large number of
erroneous annotations in the databased mitogenome sequences, including missed
genes, wrongly annotated locations, false complements, and inconsistent
criteria for assigning the 5' and 3' boundaries; and I have listed many of
these. The extent of overlap between genes is often greatly exaggerated due to
inconsistent annotations, although notable overlaps involving tRNAs are
apparently real. Finally, three novel hypotheses were examined and found to
have support from the comparative analyses: (1) some organisms have mitogenomic
locations that simultaneously code for multiple tRNAs; (2) some organisms have
mitogenomic locations that simultaneously code for tRNAs and proteins (but not
rRNAs); and (3) one group of nematodes has several genes that code for tRNAs
lacking both the D and T arms.


  This paper investigates the stability of the power-law steady state often
observed in marine ecosystems. Three dynamical systems are considered,
describing the abundance of organisms as a function of body mass and time: a
"jump-growth" equation, a first order approximation which is the widely used
McKendrick-von Foerster equation, and a second order approximation which is the
McKendrick-von Foerster equation with a diffusion term. All of these yield a
power-law steady state. We derive, for the first time, the eigenvalue spectrum
for the linearised evolution operator, under certain constraints on the
parameters. This provides new knowledge of the stability properties of the
power-law steady state. It is shown analytically that the steady state of the
McKendrick-von Foerster equation without the diffusion term is always unstable.
Furthermore, numerical plots show that eigenvalue spectra of the McKendrick-von
Foerster equation with diffusion give a good approximation to those of the
jump-growth equation. The steady state is more likely to be stable with a low
preferred predator : prey mass ratio, a large diet breadth and a high feeding
efficiency.


  Many types of cells can sense external ligand concentrations with
cell-surface receptors at extremely high accuracy. Interestingly, ligand-bound
receptors are often internalized, a process also known as receptor-mediated
endocytosis. While internalization is involved in a vast number of important
functions for the life of a cell, it was recently also suggested to increase
the accuracy of sensing ligand as the overcounting of the same ligand molecules
is reduced. Here we show, by extending simple ligand-receptor models to
out-of-equilibrium thermodynamics, that internalization increases the accuracy
with which cells can measure ligand concentrations in the external environment.
Comparison with experimental rates of real receptors demonstrates that our
model has indeed biological significance.


  A mechanism for creating amino acid enantiomerism that always selects the
same global chirality is identified, and subsequent chemical replication and
galactic mixing that would populate the galaxy with the predominant species is
described. This involves: (1) the spin of the 14N in the amino acids, or in
precursor molecules from which amino acids might be formed, coupling to the
chirality of the molecules; 2) the neutrinos emitted from the supernova,
together with magnetic field from the nascent neutron star or black hole formed
from the supernova selectively destroying one orientation of the 14N, and thus
selecting the chirality associated with the other 14N orientation; (3) chemical
evolution, by which the molecules replicate and evolve to more complex forms of
a single chirality on a relatively short timescale; and (4) galactic mixing on
a longer timescale mixing the selected molecules throughout the galaxy.


  This study elaborates some examples of a simple evolutionary stochastic rate
process where the population rate of change depends on the distribution of
properties--so different cohorts change at different rates. We investigate the
effect on the evolution arising from parametrized perturbations of uniformity
for the initial inhomogeneity. The information geometric neighbourhood system
yields also solutions for a wide range of other initial inhomogeneity
distributions, including approximations to truncated Gaussians of arbitrarily
small variance and distributions with pronounced extreme values. It is found
that, under quite considerable alterations in the shape and variance of the
initial distribution of inhomogeneity in unfitness, the decline of the mean
does change markedly with the variation in starting conditions, but the net
population evolution seems surprisingly stable.


  We generalize the concept of the geometric phase in stochastic kinetics to a
noncyclic evolution. Its application is demonstrated on kinetics of the
Michaelis-Menten reaction. It is shown that the nonperiodic geometric phase is
responsible for the correction to the Michaelis-Menten law when parameters,
such as a substrate concentration, are changing with time. We apply these ideas
to a model of chemical reactions in a bacterial culture of a growing size,
where the geometric correction qualitatively changes the outcome of the
reaction kinetics.


  Community assembly is studied using individual-based multispecies models. The
models have stochastic population dynamics with mutation, migration, and
extinction of species. Mutants appear as a result of mutation of the resident
species, while migrants have no correlation with the resident species. It is
found that the dynamics of community assembly with mutations are quite
different from the case with migrations. In contrast to mutation models, which
show intermittent dynamics of quasi-steady states interrupted by sudden
reorganizations of the community, migration models show smooth and gradual
renewal of the community. As a consequence, instead of the $1/f$ diversity
fluctuations found for the mutation models, $1/f^2$, random-walk like
fluctuations are observed for the migration models. In addition, a
characteristic species-lifetime distribution is found: a power law that is cut
off by a "skewed" distribution in the long-lifetime regime. The latter has a
longer tail than a simple exponential function, which indicates an
age-dependent species-mortality function. Since this characteristic profile has
been observed, both in fossil data and in several other mathematical models, we
conclude that it is a universal feature of macroevolution.


  Chromatin looping is a major epigenetic regulatory mechanism in higher
eukaryotes. Besides its role in transcriptional regulation, chromatin loops
have been proposed to play a pivotal role in the segregation of entire
chromosomes. The detailed topological and entropic forces between loops still
remain elusive. Here, we quantitatively determine the potential of mean force
between the centers of mass of two ring polymers, i.e. loops. We find that the
transition from a linear to a ring polymer induces a strong increase in the
entropic repulsion between these two polymers. On top, topological interactions
such as the non-catenation constraint further reduce the number of accessible
conformations of close-by ring polymers by about 50%, resulting in an
additional effective repulsion. Furthermore, the transition from linear to ring
polymers displays changes in the conformational and structural properties of
the system. In fact, ring polymers adopt a markedly more ordered and aligned
state than linear ones. The forces and accompanying changes in shape and
alignment between ring polymers suggest an important regulatory function of
such a topology in biopolymers. We conjecture that dynamic loop formation in
chromatin might act as a versatile control mechanism regulating and maintaining
different local states of compaction and order.


  RNAnet provides a bridge between two widely used Human gene databases.
Ensembl describes DNA sequences and transcripts but not experimental gene
expression. Whilst NCBI's GEO contains actual expression levels from Human
samples. RNAnet provides immediate access to thousands of Affymetrix HG-U133
2plus GeneChip measurements Homo sapiens genes in most medically interesting
tissues.
  Without RNAnet comparison across experiments in GEO is very labour intensive
requiring man months of effort to down load and clean data. With RNAnet anyone
can access cleaned quantile normalised data in seconds. It can be used to data
mine patterns of co-expression. The network of strongly correlated genes is
huge but sparse. Thousands of genes interact strongly with thousands of others.
Conversely there are tens of thousands of genes which interact strongly with
less than 100 others. I.e. RNAnet gives new views for RNA Systems Biology.


  We examine a spatially discrete reaction diffusion model based on the
interactions that create a periodic pattern in the Drosophila eye imaginal
disc. This model is capable of generating a regular hexagonal pattern of gene
expression behind a moving front, as observed in the fly system. In order to
better understand the novel switch and template mechanism behind this pattern
formation, we present here a detailed study of the model's behavior in one
dimension, using a combination of analytic methods and numerical searches of
parameter space. We find that patterns are created robustly provided that there
is an appropriate separation of timescales and that self-activation is
sufficiently strong, and we derive expressions in this limit for the front
speed and the pattern wavelength. Moving fronts in pattern-forming systems near
an initial linear instability generically select a unique pattern, but our
model operates in a strongly nonlinear regime where the final pattern depends
on the initial conditions as well as on parameter values. Our work highlights
the important role that cellularization and cell-autonomous feedback can play
in biological pattern formation.


  Hidden Markov models are traditionally decoded by the Viterbi algorithm which
finds the highest probability state path in the model. In recent years, several
limitations of the Viterbi decoding have been demonstrated, and new algorithms
have been developed to address them
\citep{Kall2005,Brejova2007,Gross2007,Brown2010}.
  In this paper, we propose a new efficient highest expected reward decoding
algorithm (HERD) that allows for uncertainty in boundaries of individual
sequence features. We demonstrate usefulness of our approach on jumping HMMs
for recombination detection in viral genomes.


  We study the effects of demographic stochasticity on the long-term dynamics
of biological coevolution models of community assembly. The noise is induced in
order to check the validity of deterministic population dynamics. While
mutualistic communities show little dependence on the stochastic population
fluctuations, predator-prey models show strong dependence on the stochasticity,
indicating the relevance of the finiteness of the populations. For a
predator-prey model, the noise causes drastic decreases in diversity and total
population size. The communities that emerge under influence of the noise
consist of species strongly coupled with each other and have stronger linear
stability around the fixed-point populations than the corresponding noiseless
model. The dynamics on evolutionary time scales for the predator-prey model are
also altered by the noise. Approximate $1/f$ fluctuations are observed with
noise, while $1/f^{2}$ fluctuations are found for the model without demographic
noise.


  With the completion of human genome mapping, the focus of scientists seeking
to explain the biological complexity of living systems is shifting from
analyzing the individual components (such as a particular gene or biochemical
reaction) to understanding the set of interactions amongst the large number of
components that results in the different functions of the organism. To this
end, the area of systems biology attempts to achieve a "systems-level"
description of biology by focusing on the network of interactions instead of
the characteristics of its isolated parts. In this article, we briefly describe
some of the emerging themes of research in "network" biology, looking at
dynamical processes occurring at the two different length scales of within the
cell and between cells, viz., the intra-cellular signaling network and the
nervous system. We show that focusing on the systems-level aspects of these
problems allows one to observe surprising and illuminating common themes
amongst them.


  We analyze gene co-expression network under the random matrix theory
framework. The nearest neighbor spacing distribution of the adjacency matrix of
this network follows Gaussian orthogonal statistics of random matrix theory
(RMT). Spectral rigidity test follows random matrix prediction for a certain
range, and deviates after wards. Eigenvector analysis of the network using
inverse participation ratio (IPR) suggests that the statistics of bulk of the
eigenvalues of network is consistent with those of the real symmetric random
matrix, whereas few eigenvalues are localized. Based on these IPR calculations,
we can divide eigenvalues in three sets; (A) The non-degenerate part that
follows RMT. (B) The non-degenerate part, at both ends and at intermediate
eigenvalues, which deviate from RMT and expected to contain information about
{\it important nodes} in the network. (C) The degenerate part with $zero$
eigenvalue, which fluctuates around RMT predicted value. We identify nodes
corresponding to the dominant modes of the corresponding eigenvectors and
analyze their structural properties.


  The seed bank dynamics of the three co-occurring obligate-seeder (i.e.
fire-sensitive) Proteaceae species, Banksia ericifolia, Banksia marginata and
Petrophile pulchella, were examined at sites of varying time since the most
recent fire (i.e. plant age) in the Sydney region. Significant variation among
species was found in the number of cones produced, the position of the cones
within the canopy, the percentage of barren cones produced (Banksia species
only), the number of follicles/bracts produced per cone, and the number of
seeds lost/released due to spontaneous fruit rupture. Thus, three different
regeneration strategies were observed, highlighting the variation in
reproductive strategies of co-occurring Proteaceae species. Ultimately, B.
marginata potentially accumulated a seed bank of c. 3000 seeds per plant after
20 years, with c. 1500 seeds per plant for P. pulchella and c. 500 for B.
ericifolia. Based on these data, B. marginata and B. ericifolia require a
minimum fire-free period of 8-10 years, with 7-8 years for P. pulchella, to
allow for an adequate seed bank to accumulate and thus ensure local persistence
of these species in fire-prone habitats.


  The regulation of the cell state is a complex process involving several
components. These complex dynamics can be modeled using Boolean networks,
allowing us to explain the existence of different cell states and the
transition between them. Boolean models have been introduced both as specific
examples and as ensemble or distribution network models. However, current
ensemble Boolean network models do not make a systematic distinction between
different cell components such as epigenetic factors, gene and transcription
factors. Consequently, we still do not understand their relative contributions
in controlling the cell fate. In this work we introduce and study higher order
Boolean networks, which feature an explicit distinction between the different
cell components and the types of interactions between them. We show that the
stability of the cell state dynamics can be determined solving the eigenvalue
problem of a matrix representing the regulatory interactions and their
strengths. The qualitative analysis of this problem indicates that, in addition
to the classification into stable and chaotic regimes, the cell state can be
simple or complex depending on whether it can be deduced from the independent
study of its components or not. Finally, we illustrate how the model can be
expanded considering higher levels and higher order dynamics.


  Gene-gene interactions have long been recognized to be fundamentally
important to understand genetic causes of complex disease traits. At present,
identifying gene-gene interactions from genome-wide case-control studies is
computationally and methodologically challenging. In this paper, we introduce a
simple but powerful method, named `BOolean Operation based Screening and
Testing'(BOOST). To discover unknown gene-gene interactions that underlie
complex diseases, BOOST allows examining all pairwise interactions in
genome-wide case-control studies in a remarkably fast manner. We have carried
out interaction analyses on seven data sets from the Wellcome Trust Case
Control Consortium (WTCCC). Each analysis took less than 60 hours on a standard
3.0 GHz desktop with 4G memory running Windows XP system. The interaction
patterns identified from the type 1 diabetes data set display significant
difference from those identified from the rheumatoid arthritis data set, while
both data sets share a very similar hit region in the WTCCC report. BOOST has
also identified many undiscovered interactions between genes in the major
histocompatibility complex (MHC) region in the type 1 diabetes data set. In the
coming era of large-scale interaction mapping in genome-wide case-control
studies, our method can serve as a computationally and statistically useful
tool.


  Populations of competing biological species exhibit a fascinating interplay
between the nonlinear dynamics of evolutionary selection forces and random
fluctuations arising from the stochastic nature of the interactions. The
processes leading to extinction of species, whose understanding is a key
component in the study of evolution and biodiversity, are influenced by both of
these factors.
  In this paper, we investigate a class of stochastic population dynamics
models based on generalized Lotka-Volterra systems. In the case of neutral
stability of the underlying deterministic model, the impact of intrinsic noise
on the survival of species is dramatic: it destroys coexistence of interacting
species on a time scale proportional to the population size. We introduce a new
method based on stochastic averaging which allows one to understand this
extinction process quantitatively by reduction to a lower-dimensional effective
dynamics. This is performed analytically for two highly symmetrical models and
can be generalized numerically to more complex situations. The extinction
probability distributions and other quantities of interest we obtain show
excellent agreement with simulations.


  It is well known among phylogeneticists that adding an extra taxon (e.g.
species) to a data set can alter the structure of the optimal phylogenetic tree
in surprising ways. However, little is known about this "rogue taxon" effect.
In this paper we characterize the behavior of balanced minimum evolution (BME)
phylogenetics on data sets of this type using tools from polyhedral geometry.
First we show that for any distance matrix there exist distances to a "rogue
taxon" such that the BME-optimal tree for the data set with the new taxon does
not contain any nontrivial splits (bipartitions) of the optimal tree for the
original data. Second, we prove a theorem which restricts the topology of
BME-optimal trees for data sets of this type, thus showing that a rogue taxon
cannot have an arbitrary effect on the optimal tree. Third, we construct
polyhedral cones computationally which give complete answers for BME rogue
taxon behavior when our original data fits a tree on four, five, and six taxa.
We use these cones to derive sufficient conditions for rogue taxon behavior for
four taxa, and to understand the frequency of the rogue taxon effect via
simulation.


  The development of systemic approaches in biology has put emphasis on
identifying genetic modules whose behavior can be modeled accurately so as to
gain insight into their structure and function. However most gene circuits in a
cell are under control of external signals and thus quantitative agreement
between experimental data and a mathematical model is difficult. Circadian
biology has been one notable exception: quantitative models of the internal
clock that orchestrates biological processes over the 24-hour diurnal cycle
have been constructed for a few organisms, from cyanobacteria to plants and
mammals. In most cases, a complex architecture with interlocked feedback loops
has been evidenced. Here we present first modeling results for the circadian
clock of the green unicellular alga Ostreococcus tauri. Two plant-like clock
genes have been shown to play a central role in Ostreococcus clock. We find
that their expression time profiles can be accurately reproduced by a minimal
model of a two-gene transcriptional feedback loop. Remarkably, best adjustment
of data recorded under light/dark alternation is obtained when assuming that
the oscillator is not coupled to the diurnal cycle. This suggests that coupling
to light is confined to specific time intervals and has no dynamical effect
when the oscillator is entrained by the diurnal cycle. This intringuing
property may reflect a strategy to minimize the impact of fluctuations in
daylight intensity on the core circadian oscillator, a type of perturbation
that has been rarely considered when assessing the robustness of circadian
clocks.


  Cell differentiation in multicellular organisms is a complex process whose
mechanism can be understood by a reductionist approach, in which the individual
processes that control the generation of different cell types are identified.
Alternatively, a large scale approach in search of different organizational
features of the growth stages promises to reveal its modular global structure
with the goal of discovering previously unknown relations between cell types.
Here we sort and analyze a large set of scattered data to construct the network
of human cell differentiation (NHCD) based on cell types (nodes) and
differentiation steps (links) from the fertilized egg to a crying baby. We
discover a dynamical law of critical branching, which reveals a fractal
regularity in the modular organization of the network, and allows us to observe
the network at different scales. The emerging picture clearly identifies
clusters of cell types following a hierarchical organization, ranging from
sub-modules to super-modules of specialized tissues and organs on varying
scales. This discovery will allow one to treat the development of a particular
cell function in the context of the complex network of human development as a
whole. Our results point to an integrated large-scale view of the network of
cell types systematically revealing ties between previously unrelated domains
in organ functions.


  In the human genomes, recombination frequency between homologous chromosomes
during meiosis is highly correlated with their physical length while it differs
significantly when their coding density is considered. Furthermore, it has been
observed that the recombination events are distributed unevenly along the
chromosomes. We have found that many of such recombination properties can be
predicted by computer simulations of population evolution based on the Monte
Carlo methods. For example, these simulations have shown that the probability
of acceptance of the recombination events by selection is higher at the ends of
chromosomes and lower in their middle parts. The regions of high coding density
are more prone to enter the strategy of haplotype complementation and to form
clusters of genes which are "recombination deserts". The phenomenon of
switching in-between the purifying selection and haplotype complementation has
a phase transition character, and many relations between the effective
population size, coding density, chromosome size and recombination frequency
are those of the power law type.


  The functional correlation of missense mutations which cause disease remains
a challenge to understanding the basis of genetic diseases. This is
particularly true for proteins related to diseases for which there are no
available three dimensional structures. One such disease is Shwachman Diamond
syndrome SDS OMIM 260400, a multi system disease arising from loss of
functional mutations. The Homo sapiens Shwachman Bodian Diamond Syndrome gene
hSBDS is responsible for SDS. hSBDS is expressed in all tissues and encodes a
protein of 250 amino acids SwissProt accession code Q9Y3A5. Sequence analysis
of disease associated alleles has identified more than 20 different mutations
in affected individuals. While a number of these mutations have been described
as leading to the loss of protein function due to truncation, translation or
surface epitope association, the structural basis for these mutations has yet
to be determined due to the lack of a three-dimensional structure for SBDS.


  Ensembl's human non-coding and protein coding genes are used to automatically
find DNA pattern motifs. The Backus-Naur form (BNF) grammar for regular
expressions (RE) is used by genetic programming to ensure the generated strings
are legal. The evolved motif suggests the presence of Thymine followed by one
or more Adenines etc. early in transcripts indicate a non-protein coding gene.
  Keywords: pseudogene, short and microRNAs, non-coding transcripts, systems
biology, machine learning, Bioinformatics, motif, regular expression, strongly
typed genetic programming, context-free grammar.


  In this paper we consider a mathematical model for the evolution and collapse
of the Easter Island society, starting from the fifth century until the last
period of the society collapse (fifteen century). Based on historical reports,
the available primary sources consisted almost exclusively on the trees. We
describe the inhabitants and the resources as an isolated system and both
considered as dynamic variables. A mathematical analysis about why the
structure of the Easter Island community collapse is performed. In particular,
we analyze the critical values of the fundamental parameters driving the
interaction humans-environment and consequently leading to the collapse. The
technological parameter, quantifying the exploitation of the resources, is
calculated and applied to the case of other extinguished civilization (Cop\'an
Maya) confirming, with a sufficiently precise estimation, the consistency of
the adopted model.


  A concentration difference of particles across a membrane perforated by pores
will induce a diffusive flux. If the diffusing objects are of the same length
scale as the the pores, diffusion may not be simple, objects can move into the
pore in a configuration that requires them to back up in order to continue
forward. A configuration that blocks flow through the pore may be statistically
preferred, an attracting metastable state of the system. This effect is purely
kinetic, and not dependent on potentials, friction or dissipation. We discuss
several geometries which generate this effect, and introduce a heuristic model
which captures the qualitative features.


  We investigate the reversible disassembly of VOV1-ATPase in life yeast cells
by time resolved confocal FRET imaging. VOV1-ATPase in the vacuolar membrane
pumps protons from the cytosol into the vacuole. VOV1-ATPase is a rotary
biological nanomotor driven by ATP hydrolysis. The emerging proton gradient is
used for transport processes as well as for pH and Ca2+ homoeostasis in the
cell. Activity of the VOV1-ATPase is regulated through assembly / disassembly
processes. During starvation the two parts of VOV1-ATPase start to disassemble.
This process is reversed after addition of glucose. The exact mechanisms are
unknown. To follow the disassembly / reassembly in vivo we tagged two subunits
C and E with different fluorescent proteins. Cellular distributions of C and E
were monitored using a duty cycle-optimized alternating laser excitation scheme
(DCO-ALEX) for time resolved confocal FRET-FLIM measurements.


  Flux analysis is a class of constraint-based approaches to the study of
biochemical reaction networks: they are based on determining the reaction flux
configurations compatible with given stoichiometric and thermodynamic
constraints. One of its main areas of application is the study of cellular
metabolic networks. We briefly and selectively review the main approaches to
this problem and then, building on recent work, we provide a characterization
of the productive capabilities of the metabolic network of the bacterium E.coli
in a specified growth medium in terms of the producible biochemical species.
While a robust and physiologically meaningful production profile clearly
emerges (including biomass components, biomass products, waste etc.), the
underlying constraints still allow for significant fluctuations even in key
metabolites like ATP and, as a consequence, apparently lay the ground for very
different growth scenarios.


  In this work we analyse the topological and dynamical properties of a simple
model of complex food webs, namely the niche model. In order to underline
competition among species, we introduce "prey" and "predators" weighted overlap
graphs derived from the niche model and compare synthetic food webs with real
data. Doing so, we find new tests for the goodness of synthetic food web models
and indicate a possible direction of improvement for existing ones. We then
exploit the weighted overlap graphs to define a competition kernel for
Lotka-Volterra population dynamics and find that for such a model the stability
of food webs decreases with its ecological complexity.


  As the behavior of a system composed of cyclically competing species is
strongly influenced by the presence of fluctuations, it is of interest to study
cyclic dominance in low dimensions where these effects are the most prominent.
We here discuss rock-paper-scissors games on a one-dimensional lattice where
the interaction rates and the mobility can be species dependent. Allowing only
single site occupation, we realize mobility by exchanging individuals of
different species. When the interaction and swapping rates are symmetric, a
strongly enhanced swapping rate yields an increased mixing of the species,
leading to a mean-field like coexistence even in one-dimensional systems. This
coexistence is transient when the rates are asymmetric, and eventually only one
species will survive. Interestingly, in our spatial games the dominating
species can differ from the species that would dominate in the corresponding
nonspatial model. We identify different regimes in the parameter space and
construct the corresponding dynamical phase diagram.


  The inability to resolve deep node relationships of highly divergent/rapidly
evolving protein families is a major factor that stymies evolutionary studies.
In this manuscript, we propose a Multiple Sequence Alignment (MSA) independent
method to infer evolutionary relationships. We previously demonstrated that
phylogenetic profiles built using position specific scoring matrices (PSSMs)
are capable of constructing informative evolutionary histories(1;2). In this
manuscript, we theorize that PSSMs derived specifically from the query
sequences used to construct the phylogenetic tree will improve this method for
the study of rapidly evolving proteins. To test this theory, we performed
phylogenetic analyses of a benchmark protein superfamily (reverse
transcriptases (RT)) as well as simulated datasets. When we compare the results
obtained from our method, PHYlogenetic ReconstructioN (PHYRN), with other MSA
dependent methods, we observe that PHYRN provides a 4- to 100-fold increase in
accurate measurements at deep nodes. As phylogenetic profiles are used as the
information source, rather than MSA, we propose PHYRN as a paradigm shift in
studying evolution when MSA approaches fail. Perhaps most importantly, due to
the improvements in our computational approach and the availability of vast
amount of sequencing data, PHYRN is scalable to thousands of sequences. Taken
together with PHYRNs adaptability to any protein family, this method can serve
as a tool for resolving ambiguities in evolutionary studies of rapidly
evolving/highly divergent protein families.


  We introduce a method to convert an ensemble of sequences of symbols into a
weighted directed network whose nodes are motifs, while the directed links and
their weights are defined from statistically significant co-occurences of two
motifs in the same sequence. The analysis of communities of networks of motifs
is shown to be able to correlate sequences with functions in the human proteome
database, to detect hot topics from online social dialogs, to characterize
trajectories of dynamical systems, and might find other useful applications to
process large amount of data in various fields.


  The past decade has witnessed the development and success of coarse-grained
network models of proteins for predicting many equilibrium properties related
to collective modes of motion. Curiously, the results are usually robust
towards the different cutoff distances used for constructing the residue
networks from the knowledge of the experimental coordinates. In this study, we
present a systematical study of network construction, and their effect on the
predicted properties. Probing bond orientational order around each residue, we
propose a natural partitioning of the interactions into an essential and a
residual set. In this picture the robustness originates from the way with which
new contacts are added so that an unusual local orientational order builds up.
These residual interactions have a vanishingly small effect on the force
vectors on each residue. The stability of the overall force balance then
translates into the Hessian as small shifts in the slow modes of motion and an
invariance of the corresponding eigenvectors. We introduce a rescaled version
of the Hessian matrix and point out a link between the matrix Frobenius norm
based on spectral stability arguments. A recipe for the optimal choice of
partitioning the interactions into essential and residual components is
prescribed. Implications for the study of biologically relevant properties of
proteins are discussed with specific examples.


  The unfolding of pandemic influenza A(H1N1) for Fall 2009 in the Northern
Hemisphere is still uncertain. Plans for vaccination campaigns and vaccine
trials are underway, with the first batches expected to be available early
October. Several studies point to the possibility of an anticipated pandemic
peak that could undermine the effectiveness of vaccination strategies. Here we
use a structured global epidemic and mobility metapopulation model to assess
the effectiveness of massive vaccination campaigns for the Fall/Winter 2009.
Mitigation effects are explored depending on the interplay between the
predicted pandemic evolution and the expected delivery of vaccines. The model
is calibrated using recent estimates on the transmissibility of the new A(H1N1)
influenza. Results show that if additional intervention strategies were not
used to delay the time of pandemic peak, vaccination may not be able to
considerably reduce the cumulative number of cases, even when the mass
vaccination campaign is started as early as mid-October. Prioritized
vaccination would be crucial in slowing down the pandemic evolution and
reducing its burden.


  Mitochondria are the power house of living cells, where the synthesis of the
chemical "energy currency" adenosine triphosphate (ATP) occurs. Oxidative
phosphorylation by a series of membrane protein complexes I to IV, that is, the
electron transport chain, is the source of the electrochemical potential
difference or proton motive force (PMF) of protons across the inner
mitochondrial membrane. The PMF is required for ATP production by complex V of
the electron transport chain, i.e. by FoF1-ATP synthase. Destroying cytochrome
C oxidase (COX; complex IV) in Photodynamic Therapy (PDT) is achieved by the
cationic photosensitizer Pt(II)-TMPyP. Electron microscopy revealed the
disruption of the mitochondrial christae as a primary step of PDT. Time
resolved phosphorescence measurements identified COX as the binding site for
Pt(II)-TMPyP in living HeLa cells. As this photosensitizer competed with
cytochrome C in binding to COX, destruction of COX might not only disturb ATP
synthesis but could expedite the release of cytochrome C to the cytosol
inducing apoptosis.


  Many biochemical processes can successfully be described by dynamical systems
allowing some form of switching when, depending on their initial conditions,
solutions of the dynamical system end up in different regions of state space
(associated with different biochemical functions). Switching is often realized
by a bistable system (i.e. a dynamical system allowing two stable steady state
solutions) and, in the majority of cases, bistability is established
numerically. In our point of view this approach is too restrictive, as, one the
one hand, due to predominant parameter uncertainty numerical methods are
generally difficult to apply to realistic models originating in Systems
Biology. And on the other hand switching already arises with the occurrence of
a saddle type steady state (characterized by a Jacobian where exactly one
Eigenvalue is positive and the remaining eigenvalues have negative real part).
Consequently we derive conditions based on linear inequalities that allow the
analytic computation of states and parameters where the Jacobian derived from a
mass action network has a defective zero eigenvalue so that -- under certain
genericity conditions -- a saddle-node bifurcation occurs. Our conditions are
applicable to general mass action networks involving at least one conservation
relation, however, they are only sufficient (as infeasibility of linear
inequalities does not exclude defective zero eigenvalues).


  In many interesting cases the reconstruction of a correct phylogeny is
blurred by high mutation rates and/or horizontal transfer events. As a
consequence a divergence arises between the true evolutionary distances and the
differences between pairs of taxa as inferred from available data, making the
phylogenetic reconstruction a challenging problem. Mathematically this
divergence translates in a loss of additivity of the actual distances between
taxa. In distance-based reconstruction methods, two properties of additive
distances were extensively exploited as antagonist criteria to drive phylogeny
reconstruction: on the one hand a local property of quartets, i.e., sets of
four taxa in a tree, the four-points condition; on the other hand a recently
proposed formula that allows to write the tree length as a function of the
distances between taxa, the Pauplin's formula. Here we introduce a new
reconstruction scheme, that exploits in a unified framework both the
four-points condition and the Pauplin's formula. We propose, in particular, a
new general class of distance-based Stochastic Local Search algorithms, which
reduces in a limit case to the minimization of the Pauplin's length. When
tested on artificially generated phylogenies our Stochastic Big-Quartet
Swapping algorithmic scheme significantly outperforms state-of-art
distance-based algorithms in cases of deviation from additivity due to high
rate of back mutations. A significant improvement is also observed with respect
to the state-of-art algorithms in case of high rate of horizontal transfer.


  Culture evolves following a process that is akin to biological evolution,
although with some significant differences. At the same time culture has often
a collective good value for human groups. This paper studies culture in an
evolutionary perspective, with a focus on the implications of group definition
for the coexistence of different cultures. A model of cultural evolution is
presented where agents interacts in an artificial environment. The belonging to
a specific memetic group is a major factor allowing agents to exploit different
environmental niches with, as a result, the coexistence of different cultures
in the same environment.


  Investigations of natural variation among cells within a population are
essential for understanding the stochastic nature of tissue cell deformation
under applied load. In the existing literature, the population variation of
single-cell creep compliance has so far been modeled universally by using a
log-normal distribution. Here we use optical stretching, a non-contact and
relatively high-throughput technique for probing cell mechanics, to accumulate
a sufficient data set that demonstrates robustly that cell compliance varies
according to the similar but distinct gamma distribution. Additionally, we
re-examine existing simulations that were originally proposed to justify a
log-normal fit, and show that in fact these simulation data also correspond to
the gamma distribution. Finally, we propose a general stochastic differential
equation that analytically predicts a gamma distribution of creep compliance
during cell stretching, as well as the Gaussian distribution of cell recovery
that we observe experimentally upon removal of applied load. The population
variation is well characterized by just a single parameter in each of the creep
and recovery regimes. We expect our correction of a phenomenological
distribution fit, enabled by an expansive data set for mesenchymal stem cells,
to enable the development of more accurate constitutive laws to describe
cytoskeletal deformation. These findings thus serve to replace an empirical
distribution with a better-fitting model that rests on a more solid
experimental and theoretical foundation, and also provides a basis to predict
and understand the stochastic nature of the mechanical response of individual
cells within populations.


  Recent experiments revealing possible nanoscale electrostatic interactions in
force generation at kinetochores for chromosome motions have prompted
speculation regarding possible models for interactions between positively
charged molecules in kinetochores and negative charge on C-termini near the
plus ends of microtubules. A clear picture of how kinetochores establish and
maintain a dynamic coupling to microtubules for force generation during the
complex motions of mitosis remains elusive. The current paradigm of molecular
cell biology requires that specific molecules, or molecular geometries, for
force generation be identified. However, it is possible to account for mitotic
motions within a classical electrostatics approach in terms of experimentally
known cellular electric charge interacting over nanometer distances. These
charges are modeled as bound surface and volume continuum charge distributions.
Electrostatic consequences of intracellular pH changes during mitosis may
provide a master clock for the events of mitosis.


  We consider stochastic matrix models for population driven by random
environments which form a Markov chain. The top Lyapunov exponent $a$, which
describes the long-term growth rate, depends smoothly on the demographic
parameters (represented as matrix entries) and on the parameters that define
the stochastic matrix of the driving Markov chain. The derivatives of $a$ --
the "stochastic elasticities" -- with respect to changes in the demographic
parameters were derived by \cite{tuljapurkar1990pdv}. These results are here
extended to a formula for the derivatives with respect to changes in the Markov
chain driving the environments. We supplement these formulas with rigorous
bounds on computational estimation errors, and with rigorous derivations of
both the new and the old formulas.


  Microtubules are highly regulated dynamic elements of the cytoskeleton of
eukaryotic cells. One of the regulation mechanisms observed in living cells is
the severing by the proteins katanin and spastin. We introduce a model for the
dynamics of microtubules in the presence of randomly occurring severing events.
Under the biologically motivated assumption that the newly created plus end
undergoes a catastrophe, we investigate the steady state length distribution.
We show that the presence of severing does not affect the number of
microtubules, regardless of the distribution of severing events. In the special
case in which the microtubules cannot recover from the depolymerizing state (no
rescue events) we derive an analytical expression for the length distribution.
In the general case we transform the problem into a single ODE that is solved
numerically.


  The population dynamics theory of B cells in a typical germinal center could
play an important role in revealing how affinity maturation is achieved.
However, the existing models encountered some conflicts with experiments. To
resolve these conflicts, we present a coarse-grained model to calculate the B
cell population development in affinity maturation, which allows a
comprehensive analysis of its parameter space to look for optimal values of
mutation rate, selection strength, and initial antibody-antigen binding level
that maximize the affinity improvement. With these optimized parameters, the
model is compatible with the experimental observations such as the ~100-fold
affinity improvements, the number of mutations, the hypermutation rate, and the
"all or none" phenomenon. Moreover, we study the reasons behind the optimal
parameters. The optimal mutation rate, in agreement with the hypermutation rate
in vivo, results from a tradeoff between accumulating enough beneficial
mutations and avoiding too many deleterious or lethal mutations. The optimal
selection strength evolves as a balance between the need for affinity
improvement and the requirement to pass the population bottleneck. These
findings point to the conclusion that germinal centers have been optimized by
evolution to generate strong affinity antibodies effectively and rapidly. In
addition, we study the enhancement of affinity improvement due to B cell
migration between germinal centers. These results could enhance our
understandings to the functions of germinal centers.


  Cells can often choose among several stably heritable phenotypes. Examples
are the expression of genes in eukaryotic cells where long chromosomal regions
can adopt persistent and heritable silenced or active states, that may be
associated with positive feedback in dynamic modification of nucleosomes. We
generalize this mechanism in terms of bistability associated with valleys in an
epigenetic landscape. A transfer matrix method was used to rigorously follow
the system through the disruptive process of cell division. This combined
treatment of noisy dynamics both between and during cell division provides an
efficient way to calculate the stability of alternative states in a broad range
of epigenetic systems.


  New high-throughput sequencing technologies have made it possible to pursue
the advent of genome-wide transcriptomics. That progress combined with the
recent discovery of regulatory non-coding RNAs (ncRNAs) has necessitated fast
and accurate algorithms to predict RNA-RNA interaction probability and
structure. Although there are algorithms to predict minimum free energy
interaction secondary structure for two nucleic acids, little work has been
done to exploit the information invested in the base pair probabilities to
improve interaction structure prediction. In this paper, we present an
algorithm to predict the Hamming centroid of the Boltzmann ensemble of
interaction structures. We also present an efficient algorithm to sample
interaction structures from the ensemble. Our sampling algorithm uses a
balanced scheme for traversing indices which improves the running time of the
Ding-Lawrence sampling algorithm. The Ding-Lawrence sampling algorithm has
$O(n^2m^2)$ time complexity whereas our algorithm has $O((n+m)^2\log(n+m))$
time complexity, in which $n$ and $m$ are the lengths of input strands. We
implemented our algorithm in a new version of {\tt piRNA} and compared our
structure prediction results with competitors. Our centroid prediction
outperforms competitor minimum-free-energy prediction algorithms on average.


  The constraints imposed by nano- and microscale confinement on the
conformational degrees of freedom of thermally fluctuating biopolymers are
utilized in contemporary nano-devices to specifically elongate and manipulate
single chains. A thorough theoretical understanding and quantification of the
statistical conformations of confined polymer chains is thus a central concern
in polymer physics. We present an analytical calculation of the radial
distribution function of harmonically confined semiflexible polymers in the
weakly bending limit. Special emphasis has been put on a proper treatment of
global modes, i.e. the possibility of the chain to perform global movements
within the channel. We show that the effect of these global modes significantly
impacts the chain statistics in cases of weak and intermediate confinement.
Comparing our analytical model to numerical data from Monte Carlo simulations
we find excellent agreement over a broad range of parameters.


  Recent experiments [J. van Mameren et al. PNAS 106, 18231 (2009)] provide a
detailed spatial picture of overstretched DNA, showing that under certain
conditions the two strands of the double helix separate at about 65 pN. It was
proposed that this observation rules out the existence of an elongated,
hybridized form of DNA ('S-DNA'). Here we argue that the S-DNA picture is
consistent with the observation of unpeeling during overstretching. We
demonstrate that assuming the existence of S-DNA does not imply DNA
overstretching to consist of the complete or near-complete conversion of the
molecule from B- to S-form. Instead, this assumption implies in general a more
complex dynamic coexistence of hybridized and unhybridized forms of DNA. We
argue that such coexistence can rationalize several recent experimental
observations.


  Twisting and bending deformations are crucial to the biological functions of
microfilaments such as DNA molecules. Although continuum-rod models have
emerged as efficient tools to describe the nonlinear dynamics of these
deformations, a major roadblock in the continuum-mechanics-based description of
microfilaments is the accurate modeling of the constitutive law, which follows
from its atomistic structure and bond-stiffnesses. Since first-principle
derivation of the constitutive law from atomistic structure is impractical and
so are direct experimental measurements due to the small length-scales, a
natural alternative is to estimate the constitutive law from discrete-structure
simulations such as molecular-dynamics (MD) simulations. In this paper, we
present a two-step inverse method for estimating the constitutive law using rod
theory and data generated from discrete-structure simulations. We illustrate
the method on a filament with an artificial and simplistic discrete-structure.
We simulate its deformation in response to a prescribed loading using a
multi-body dynamics (MBD) solver. Using data generated from the MBD solver, we
first estimate the curvature of the filament and subsequently use it in the
two-step method to estimate the effective constitutive-law relationship between
the restoring moment and curvature. Finally, we also illustrate how the
estimated constitutive law can be tested under independent loading conditions.


  Aggregates of twisted protein fibers, such as sickle hemoglobin and actin,
are important examples of biopolymers in which elastic interactions play a
crucial role in determining the (metastable) bundle radii. Here, we present a
corrected version of analysis on the stability of sickle hemoglobin fibers
using the classic nucleation theory.


  Using the perturbation-response scanning (PRS) technique, we study a set of
23 proteins that display a variety of conformational motions upon ligand
binding (e.g. shear, hinge, allosteric). In most cases, PRS determines residues
that may be manipulated to achieve the resulting conformational change. PRS
reveals that for some proteins, binding induced conformational change may be
achieved through the perturbation of residues scattered throughout the protein,
whereas in others, perturbation of specific residues confined to a highly
specific region are necessary. Correlations between the experimental and
calculated atomic displacements are always better or equivalent to those
obtained from a modal analysis of elastic network models. Furthermore, best
correlations obtained by the latter approach do not always appear in the most
collective modes. We show that success of the modal analysis depends on the
lack of redundant paths that exist in the protein. PRS thus demonstrates that
several relevant modes may simultaneously be induced by perturbing a single
select residue on the protein. We also illustrate the biological relevance of
applying PRS on the GroEL and ADK structures in detail, where we show that the
residues whose perturbation lead to the precise conformational changes usually
correspond to those experimentally determined to be functionally important.


  Intracellular transmission of information via chemical and transcriptional
networks is thwarted by a physical limitation: the finite copy number of the
constituent chemical species introduces unavoidable intrinsic noise. Here we
provide a method for solving for the complete probabilistic description of
intrinsically noisy oscillatory driving. We derive and numerically verify a
number of simple scaling laws. Unlike in the case of measuring a static
quantity, response to an oscillatory driving can exhibit a resonant frequency
which maximizes information transmission. Further, we show that the optimal
regulatory design is dependent on the biophysical constraints (i.e., the
allowed copy number and response time). The resulting phase diagram illustrates
under what conditions threshold regulation outperforms linear regulation.


  Dr. Hatchwell [2010] has proposed that the BAK1 gene variants were likely due
to sequencing of a processed gene on chromosome 20. However, in response, Dr.
Gottlieb and co-authors [2010] have argued that "some but not all of the
sequence changes present in the BAK1 sequence of our abdominal aorta samples
are also present in the chromosome 20 BAK1 sequence. However, all the AAA and
AA cDNA samples are identical to each other and different from chromosome 20
BAK1 sequence at amino acids 2 and 145". I have been following this discussion
because I have independently reached almost the same conclusion as Dr.
Hatchwell did [Yamagishi, 2009], and, unfortunately, the response from Dr.
Gottlieb and his co-authors seems to me to be unsatisfactory for the reasons
listed below


  Living cells are continually exposed to environmental signals that vary in
time. These signals are detected and processed by biochemical networks, which
are often highly stochastic. To understand how cells cope with a fluctuating
environment, we therefore have to understand how reliably biochemical networks
can transmit time-varying signals. To this end, we must understand both the
noise characteristics and the amplification properties of networks. In this
manuscript, we use information theory to study how reliably signalling cascades
employing autoregulation and feedback can transmit time-varying signals. We
calculate the frequency-dependence of the gain-to-noise ratio, which reflects
how reliably a network transmits signals at different frequencies. We find that
the gain-to-noise ratio may differ qualitatively from the power spectrum of the
output, showing that the latter does not directly reflect signaling
performance. Moreover, we find that auto-activation and auto-repression
increase and decrease the gain-to-noise ratio for all of frequencies,
respectively. Positive feedback specifically enhances information transmission
at low frequencies, while negative feedback increases signal fidelity at high
frequencies. Our analysis not only elucidates the role of autoregulation and
feedback in naturally-occurring biological networks, but also reveals design
principles that can be used for the reliable transmission of time-varying
signals in synthetic gene circuits.


  The collective dynamics of coupled two-dimensional chaotic maps on complex
networks is known to exhibit a rich variety of emergent properties which
crucially depend on the underlying network topology. We investigate the
collective motion of Chirikov standard maps interacting with time delay through
directed links of Gene Regulatory Network of bacterium Escherichia Coli.
Departures from strongly chaotic behavior of the isolated maps are studied in
relation to different coupling forms and strengths. At smaller coupling
intensities the network induces stable and coherent emergent dynamics. The
unstable behavior appearing with increase of coupling strength remains confined
within a connected sub-network. For the appropriate coupling, network exhibits
statistically robust self-organized dynamics in a weakly chaotic regime.


  In non-viral gene delivery, the variance of transgenic expression stems from
the low number of plasmids successfully transferred. Here, we experimentally
determine Lipofectamine- and PEI-mediated exogenous gene expression
distributions from single cell time-lapse analysis. Broad Poisson-like
distributions of steady state expression are observed for both transfection
agents, when used with synchronized cell lines. At the same time,
co-transfection analysis with YFP- and CFP-coding plasmids shows that multiple
plasmids are simultaneously expressed, suggesting that plasmids are delivered
in correlated units (complexes). We present a mathematical model of
transfection, where a stochastic, two-step process is assumed, with the first
being the low-probability entry step of complexes into the nucleus, followed by
the subsequent release and activation of a small number of plasmids from a
delivered complex. This conceptually simple model consistently predicts the
observed fraction of transfected cells, the cotransfection ratio and the
expression level distribution. It yields the number of efficient plasmids per
complex and elucidates the origin of the associated noise, consequently
providing a platform for evaluating and improving non-viral vectors.


  Recently, the behavior of different epidemic models and their relation both
to different types of geometries and to some biological models has been
revisited . Path equations representing the behavior of epidemic models and
their corresponding deviation vectors are examined. A comparison between paths
and their deviation vectors in Riemannian and Finslerian Geometries is
presented.


  Our understanding of the evolutionary process has gone a long way since the
publication, 150 years ago, of "On the origin of species" by Charles R. Darwin.
The XXth Century witnessed great efforts to embrace replication, mutation, and
selection within the framework of a formal theory, able eventually to predict
the dynamics and fate of evolving populations. However, a large body of
empirical evidence collected over the last decades strongly suggests that some
of the assumptions of those classical models necessitate a deep revision. The
viability of organisms is not dependent on a unique and optimal genotype. The
discovery of huge sets of genotypes (or neutral networks) yielding the same
phenotype --in the last term the same organism--, reveals that, most likely,
very different functional solutions can be found, accessed and fixed in a
population through a low-cost exploration of the space of genomes. The
'evolution behind the curtain' may be the answer to some of the current puzzles
that evolutionary theory faces, like the fast speciation process that is
observed in the fossil record after very long stasis periods.


  When a ligand that is bound to an integral membrane receptor is pulled, the
membrane and the underlying cytoskeleton can deform before either the membrane
delaminates from the cytoskeleton or the ligand detaches from the receptor. If
the membrane delaminates from the cytoskeleton, it may be further extruded and
form a membrane tether. We develop a phenomenological model for this processes
by assuming that deformations obey Hooke's law up to a critical force at which
the cell membrane locally detaches from the cytoskeleton and a membrane tether
forms. We compute the probability of tether formation and show that they can be
extruded only within an intermediate range of force loading rates and pulling
velocities. The mean tether length that arises at the moment of ligand
detachment is computed as are the force loading rates and pulling velocities
that yield the longest tethers.


  Individual resource intake rates are known to depend on both individual body
size and resource availability. Here, we have developed a model to integrate
these two drivers, accounting explicitly for the scaling of perceived resource
availability with individual body size. The model merges a Kleiber-like scaling
law with Holling functional responses into a single mathematical framework,
involving both body-size the density of resources.
  When the availability of resources is held constant the model predicts a
relationship between resource intake rates and body sizes whose log-log graph
is a concave curve. The significant deviation from a power law accounts for the
body size dependency of resource limitations. The model results are consistent
with data from both a laboratory experiment on benthic macro-invertebrates and
the available literature.


  Cooperation in the evolutionary snowdrift game with a self-questioning
updating mechanism is studied on annealed and quenched small-world networks
with directed couplings. Around the payoff parameter value $r=0.5$, we find a
size-invariant symmetrical cooperation effect. While generally suppressing
cooperation for $r>0.5$ payoffs, rewired networks facilitated cooperative
behavior for $r<0.5$. Fair amounts of noise were found to break the observed
symmetry and further weaken cooperation at relatively large values of $r$.
However, in the absence of noise, the self-questioning mechanism recovers
symmetrical behavior and elevates altruism even under large-reward conditions.
Our results suggest that an updating mechanism of this type is necessary to
stabilize cooperation in a spatially structured environment which is otherwise
detrimental to cooperative behavior, especially at high cost-to-benefit ratios.
Additionally, we employ component and local stability analyses to better
understand the nature of the manifested dynamics.


  We report on the nature of the thermal denaturation transition of homogeneous
DNA as determined from a renormalisation group analysis of the
Peyrard-Bishop-Dauxois model. Our approach is based on an analogy with the
phenomenon of critical wetting that goes further than previous qualitative
comparisons, and shows that the transition is continuous for the average
base-pair separation. However, since the range of universal critical behaviour
appears to be very narrow, numerically observed denaturation transitions may
look first-order, as it has been reported in the literature.


  Living cells provide a fluctuating, out-of-equilibrium environment in which
genes must coordinate cellular function. DNA looping, which is a common means
of regulating transcription, is very much a stochastic process; the loops arise
from the thermal motion of the DNA and other fluctuations of the cellular
environment. We present single-molecule measurements of DNA loop formation and
breakdown when an artificial fluctuating force, applied to mimic a fluctuating
cellular environment, is imposed on the DNA. We show that loop formation is
greatly enhanced in the presence of noise of only a fraction of $k_B T$, yet
find that hypothetical regulatory schemes that employ mechanical tension in the
DNA--as a sensitive switch to control transcription--can be surprisingly robust
due to a fortuitous cancellation of noise effects.


  We present stochastic, finite-population formulations of the Crow-Kimura and
Eigen models of quasispecies theory, for fitness functions that depend in an
arbitrary way on the number of mutations from the wild type. We include back
mutations in our description. We show that the fluctuation of the population
numbers about the average values are exceedingly large in these physical models
of evolution. We further show that horizontal gene transfer reduces by orders
of magnitude the fluctuations in the population numbers and reduces the
accumulation of deleterious mutations in the finite population due to Muller's
ratchet. Indeed the population sizes needed to converge to the infinite
population limit are often larger than those found in nature for smooth fitness
functions in the absence of horizontal gene transfer. These analytical results
are derived for the steady-state by means of a field-theoretic representation.
Numerical results are presented that indicate horizontal gene transfer speeds
up the dynamics of evolution as well.