  Boolean networks are discrete dynamical systems in which the state (zero or
one) of each node is updated at each time t to a state determined by the states
at time t-1 of those nodes that have links to it. When these systems are used
to model genetic control, the case of 'canalizing' update rules is of
particular interest. A canalizing rule is one for which a node state at time
$t$ is determined by the state at time t-1 of a single one of its inputs when
that inputting node is in its canalizing state. Previous work on the
order/disorder transition in Boolean networks considered complex, non-random
network topology. In the current paper we extend this previous work to account
for canalizing behavior.


  Prevalence of cooperation within groups of selfish individuals is puzzling in
that it contradicts with the basic premise of natural selection. Favoring
players with higher fitness, the latter is key for understanding the challenges
faced by cooperators when competing with defectors. Evolutionary game theory
provides a competent theoretical framework for addressing the subtleties of
cooperation in such situations, which are known as social dilemmas. Recent
advances point towards the fact that the evolution of strategies alone may be
insufficient to fully exploit the benefits offered by cooperative behavior.
Indeed, while spatial structure and heterogeneity, for example, have been
recognized as potent promoters of cooperation, coevolutionary rules can extend
the potentials of such entities further, and even more importantly, lead to the
understanding of their emergence. The introduction of coevolutionary rules to
evolutionary games implies, that besides the evolution of strategies, another
property may simultaneously be subject to evolution as well. Coevolutionary
rules may affect the interaction network, the reproduction capability of
players, their reputation, mobility or age. Here we review recent works on
evolutionary games incorporating coevolutionary rules, as well as give a
didactic description of potential pitfalls and misconceptions associated with
the subject. In addition, we briefly outline directions for future research
that we feel are promising, thereby particularly focusing on dynamical effects
of coevolutionary rules on the evolution of cooperation, which are still widely
open to research and thus hold promise of exciting new discoveries.


  We present a minimal physical model for the flagellar motor that enables
bacteria to swim. Our model explains the experimentally measured torque-speed
relationship of the proton-driven E. coli motor at various pH and temperature
conditions. In particular, the dramatic drop of torque at high rotation speeds
(the "knee") is shown to arise from saturation of the proton flux. Moreover, we
show that shot noise in the proton current dominates the diffusion of motor
rotation at low loads. This suggests a new way to probe the discreteness of the
energy source, analogous to measurements of charge quantization in
superconducting tunnel junctions.


  This is a preliminary article stating and proving a new maximum entropy
theorem. The entropies that we consider can be used as measures of
biodiversity. In that context, the question is: for a given collection of
species, which frequency distribution(s) maximize the diversity? The theorem
provides the answer. The chief surprise is that although we are dealing with
not just a single entropy, but a one-parameter family of entropies, there is a
single distribution maximizing all of them simultaneously.


  Computing the similarity between two protein structures is a crucial task in
molecular biology, and has been extensively investigated. Many protein
structure comparison methods can be modeled as maximum clique problems in
specific k-partite graphs, referred here as alignment graphs. In this paper, we
propose a new protein structure comparison method based on internal distances
(DAST) which is posed as a maximum clique problem in an alignment graph. We
also design an algorithm (ACF) for solving such maximum clique problems. ACF is
first applied in the context of VAST, a software largely used in the National
Center for Biotechnology Information, and then in the context of DAST. The
obtained results on real protein alignment instances show that our algorithm is
more than 37000 times faster than the original VAST clique solver which is
based on Bron & Kerbosch algorithm. We furthermore compare ACF with one of the
fastest clique finder, recently conceived by Ostergard. On a popular benchmark
(the Skolnick set) we observe that ACF is about 20 times faster in average than
the Ostergard's algorithm.


  Gradual changes in exploitation, nutrient loading, etc. produce shifts
between alternative stable states (ASS) in ecosystems which, quite often, are
not smooth but abrupt or catastrophic. Early warnings of such catastrophic
regime shifts are fundamental for designing management protocols for
ecosystems. Here we study the spatial version of a popular ecological model,
involving a logistically growing single species subject to exploitation, which
is known to exhibit ASS. Spatial heterogeneity is introduced by a carrying
capacity parameter varying from cell to cell in a regular lattice. Transport of
biomass among cells is included in the form of diffusion. We investigate
whether different quantities from statistical mechanics -like the variance, the
two-point correlation function and the patchiness- may serve as early warnings
of catastrophic phase transitions between the ASS. In particular, we find that
the patch-size distribution follows a power law when the system is close to the
catastrophic transition. We also provide links between spatial and temporal
indicators and analyze how the interplay between diffusion and spatial
heterogeneity may affect the earliness of each of the observables. We find that
possible remedial procedures, which can be followed after these early signals,
are more effective as the diffusion becomes lower. Finally, we comment on
similarities and differences between these catastrophic shifts and paradigmatic
thermodynamic phase transitions like the liquid-vapour change of state for a
fluid like water.


  We study the enzymatic degradation of an elastic fiber under tension using an
an isotropic random-walk model, coupled with binding-unbinding reactions that
weaken the fiber. The fiber is represented by a chain of elastic springs in
series, surrounded by two layers of sites along which enzyme molecules can
diffuse. Through numerical simulations we show that the fiber stiffness
decreases exponentially with two distinct regimes. The time constant associated
with the first regime decreases with increasing applied force, which is in
agreement with published experimental data. In addition, a simple mean field
calculation allows us to partition the time constant into geometrical, chemical
and externally controllable factors, which is corroborated by the simulations.


  The role of the topology and its relation with the geometry of biopolymers
under different physical conditions is a nontrivial and interesting problem.
Aiming at understanding this issue for a related simpler system, we use Monte
Carlo methods to investigate the interplay between writhe and knotting of ring
polymers in good and poor solvents. The model that we consider is interacting
self-avoiding polygons on the simple cubic lattice. For polygons with fixed
knot type we find a writhe distribution whose average depends on the knot type
but is insensitive to the length $N$ of the polygon and to solvent conditions.
This "topological contribution" to the writhe distribution has a value that is
consistent with that of ideal knots. The standard deviation of the writhe
increases approximately as $\sqrt{N}$ in both regimes and this constitutes a
geometrical contribution to the writhe. If the sum over all knot types is
considered, the scaling of the standard deviation changes, for compact
polygons, to $\sim N^{0.6}$. We argue that this difference between the two
regimes can be ascribed to the topological contribution to the writhe that, for
compact chains, overwhelms the geometrical one thanks to the presence of a
large population of complex knots at relatively small values of $N$. For
polygons with fixed writhe we find that the knot distribution depends on the
chosen writhe, with the occurrence of achiral knots being considerably
suppressed for large writhe. In general, the occurrence of a given knot thus
depends on a nontrivial interplay between writhe, chain length, and solvent
conditions.


  A biological regulatory network can be modeled as a discrete function that
contains all available information on network component interactions. From this
function we can derive a graph representation of the network structure as well
as of the dynamics of the system. In this paper we introduce a method to
identify modules of the network that allow us to construct the behavior of the
given function from the dynamics of the modules. Here, it proves useful to
distinguish between dynamical and structural modules, and to define network
modules combining aspects of both. As a key concept we establish the notion of
symbolic steady state, which basically represents a set of states where the
behavior of the given function is in some sense predictable, and which gives
rise to suitable network modules. We apply the method to a regulatory network
involved in T helper cell differentiation.


  The chemotactic pathway allows bacteria to respond and adapt to environmental
changes, by tuning the tumbling and running motions that are due to clockwise
and counterclockwise rotations of their flagella. The pathway is tightly
regulated by feedback mechanisms governed by the phosphorylation and
methylation of several proteins. In this paper, we present a detailed
mechanistic model for chemotaxis, that considers all of its transmembrane and
cytoplasmic components, and their mutual interactions. Stochastic simulations
of the dynamics of a pivotal protein, CheYp, are performed by means of tau
leaping algorithm. This approach is then used to investigate the interplay
between the stochastic fluctuations of CheYp amount and the number of cellular
flagella. Our results suggest that the combination of these factors might
represent a relevant component for chemotaxis. Moreover, we study the pathway
under various conditions, such as different methylation levels and ligand
amounts, in order to test its adaptation response. Some issues for future work
are finally discussed.


  Appropriate regulation of gene expression is essential to ensure that protein
synthesis occurs in a selective manner. The control of transcription is the
most dominant type of regulation mediated by a complex of molecules such as
transcription factors. In general, regulatory molecules are of two types:
activator and repressor. Activators promote the initiation of transcription
whereas repressors inhibit transcription. In many cases, they regulate the gene
transcription on binding the promoter mutually exclusively and the observed
gene expression response is either graded or binary. In experiments, the gene
expression response is quantified by the amount of proteins produced on varying
the concentration of an external inducer molecules in the cell. In this paper,
we study a gene regulatory network where activators and repressors both bind
the same promoter mutually exclusively. The network is modeled by assuming that
the gene can be in three possible states: repressed, unregulated and active. An
exact analytical expression for the steady-state probability distribution of
protein levels is then derived. The exact result helps to explain the
experimental observations that in the presence of activator molecules the
response is graded at all inducer levels whereas in the presence of both
activator and repressor molecules, the response is graded at low and high
inducer levels and binary at an intermediate inducer level.


  Regulation of mRNA decay is a critical component of global cellular
adaptation to changing environments. The corresponding changes in mRNA
lifetimes can be coordinated with changes in mRNA transcription rates to
fine-tune gene expression. Current approaches for measuring mRNA lifetimes can
give rise to secondary effects due to transcription inhibition and require
separate experiments to estimate changes in mRNA transcription rates. Here, we
propose an approach for simultaneous determination of changes in mRNA
transcription rate and lifetime using regulatory small RNAs to control mRNA
decay. We analyze a stochastic model for coupled degradation of mRNAs and sRNAs
and derive exact results connecting RNA lifetimes and transcription rates to
mean abundances. The results obtained show how steady-state measurements of RNA
levels can be used to analyze factors and processes regulating changes in mRNA
transcription and decay.


  Recent molecular dynamics simulations (Herce and Garcia, PNAS, 104: 20805
(2007)) have suggested that the arginine-rich HIV Tat peptides might be able to
translocate by destabilizing and inducing transient pores in phospholipid
bilayers. In this pathway for peptide translocation, arginine residues play a
fundamental role not only in the binding of the peptide to the surface of the
membrane but also in the destabilization and nucleation of transient pores
across the bilayer, despite being charged and highly hydrophilic. Here we
present a molecular dynamics simulation of a peptide composed of nine arginines
(Arg-9) that shows that this peptide follows the same translocation pathway
previously found for the Tat peptide. We test this hypothesis experimentally by
measuring ionic currents across phospholipid bilayers and cell membranes
through the pores induced by Arg-9 peptides. We find that Arg-9 peptides, in
the presence of an electrostatic potential gradient, induce ionic currents
across planar phospholipid bilayers, as well as in cultured osteosarcoma cells
and human smooth muscle cells freshly isolated from the umbilical artery. Our
results suggest that the mechanism of action of Arg-9 peptide involves the
creation of transient pores in lipid bilayers and cell membranes.


  Accurate reconstruction of phylogenies remains a key challenge in
evolutionary biology. Most biologically plausible formulations of the problem
are formally NP-hard, with no known efficient solution. The standard in
practice are fast heuristic methods that are empirically known to work very
well in general, but can yield results arbitrarily far from optimal. Practical
exact methods, which yield exponential worst-case running times but generally
much better times in practice, provide an important alternative. We report
progress in this direction by introducing a provably optimal method for the
weighted multi-state maximum parsimony phylogeny problem. The method is based
on generalizing the notion of the Buneman graph, a construction key to
efficient exact methods for binary sequences, so as to apply to sequences with
arbitrary finite numbers of states with arbitrary state transition weights. We
implement an integer linear programming (ILP) method for the multi-state
problem using this generalized Buneman graph and demonstrate that the resulting
method is able to solve data sets that are intractable by prior exact methods
in run times comparable with popular heuristics. Our work provides the first
method for provably optimal maximum parsimony phylogeny inference that is
practical for multi-state data sets of more than a few characters.


  Feldman and Karlin conjectured that the number of isolated fixed points for
deterministic models of viability selection and recombination among n possible
haplotypes has an upper bound of 2^n - 1. Here a proof is provided. The upper
bound of 3^{n-1} obtained by Lyubich et al. (2001) using Bezout's Theorem
(1779) is reduced here to 2^n through a change of representation that reduces
the third-order polynomials to second order. A further reduction to 2^n - 1 is
obtained using the homogeneous representation of the system, which yields
always one solution `at infinity'. While the original conjecture was made for
systems of viability selection and recombination, the results here generalize
to viability selection with any arbitrary system of bi-parental transmission,
which includes recombination and mutation as special cases. An example is
constructed of a mutation-selection system that has 2^n - 1 fixed points given
any n, which shows that 2^n - 1 is the sharpest possible upper bound that can
be found for the general space of selection and transmission coefficients.


  We show that strategy independent adaptations of random interaction networks
can induce powerful mechanisms, ranging from the Red Queen to group selection,
that promote cooperation in evolutionary social dilemmas. These two mechanisms
emerge spontaneously as dynamical processes due to deletions and additions of
links, which are performed whenever players adopt new strategies and after a
certain number of game iterations, respectively. The potency of cooperation
promotion, as well as the mechanism responsible for it, can thereby be tuned
via a single parameter determining the frequency of link additions. We thus
demonstrate that coevolving random networks may evoke an appropriate mechanism
for each social dilemma, such that cooperation prevails even by highly
unfavorable conditions.


  We use computer simulations to study a model, first proposed by Wales [1],
for the reversible and monodisperse self-assembly of simple icosahedral virus
capsid structures. The success and efficiency of assembly as a function of
thermodynamic and geometric factors can be qualitatively related to the
potential energy landscape structure of the assembling system. Even though the
model is strongly coarse-grained, it exhibits a number of features also
observed in experiments, such as sigmoidal assembly dynamics, hysteresis in
capsid formation and numerous kinetic traps. We also investigate the effect of
macromolecular crowding on the assembly dynamics. Crowding agents generally
reduce capsid yields at optimal conditions for non-crowded assembly, but may
increase yields for parameter regimes away from the optimum. Finally, we
generalize the model to a larger triangulation number T = 3, and observe more
complex assembly dynamics than that seen for the original T = 1 model.


  Since its discovery of the living fossil in 1938, the coelacanth (Latimeria
chalumnae) has generally been considered to be the closest living relative of
the land vertebrates, and this is still the prevailing opinion in most general
biology textbooks. However, the origin of tetrapods has been the subject of
intense debate for decades. The three principal hypothesis (lungfish-tetrapod,
coelacanth-tetrapod, or lungfish-coelacanth sister group) have been proposed.
We used the maximum gene-support tree approach to analyze 43 nuclear genes
encoding amino acid residues, and compared the results of concatenation and
majority-rule tree approaches. The results inferred with three common
phylogenetic methods and three genome-scale approaches consistently rejected
the hypothesis that the coelacanth is the closest living relative of tetrapods.


  This paper has been withdrawn by the author(s), due t a modification in Eqn.
1.


  Summary: GeneSupport implements a genome-scale algorithm: Maximum
Gene-Support Tree to estimate species tree from gene trees based on multilocus
sequences. It provides a new option for multiple genes to infer species tree.
It is incorporated into popular phylogentic program: PHYLIP package with the
same usage and user interface. It is suitable for phylogenetic methods such as
maximum parsimony, maximum likelihood, Baysian and neighbour-joining, which is
used to reconstruct single gene trees firstly with a variety of phylogenetic
inference programs.


  Many pathogens spread primarily via direct contact between infected and
susceptible hosts. Thus, the patterns of contacts or contact network of a
population fundamentally shapes the course of epidemics. While there is a
robust and growing theory for the dynamics of single epidemics in networks, we
know little about the impacts of network structure on long term epidemic or
endemic transmission. For seasonal diseases like influenza, pathogens
repeatedly return to populations with complex and changing patterns of
susceptibility and immunity acquired through prior infection. Here, we develop
two mathematical approaches for modeling consecutive seasonal outbreaks of a
partially-immunizing infection in a population with contact heterogeneity.
Using methods from percolation theory we consider both leaky immunity, where
all previously infected individuals gain partial immunity, and perfect
immunity, where a fraction of previously infected individuals are fully immune.
By restructuring the epidemiologically active portion of their host population,
such diseases limit the potential of future outbreaks. We speculate that these
dynamics can result in evolutionary pressure to increase infectiousness.


  Macromolecules occupy between 34 and 44% of the cell cytoplasm, about half
the maximum pack- ing density of spheres in three dimension. Yet, there is no
clear understanding of what is special about this value. To address this
fundamental question we investigate the effect of macromolecular crowding on
cell metabolism. We develop a cell scale flux balance model capturing the main
features of cell metabolism at different nutrient uptakes and macromolecular
densities. Using this model we show there are two metabolic regimes at low and
high nutrient uptakes. The latter regime is charac- terized by an optimal
cytoplasmatic density where the increase of reaction rates by confinement and
the decrease by diffusion slow-down balance. More important, the predicted
optimal density is in the range of the experimentally determined density of E.
coli. We conclude that cells have evolved to a cytoplasmatic density resulting
in the maximum metabolic rate given the nutrient availability and
macromolecular crowding effects and report a flux balance model accounting for
its effect.


  We present a framework for analyzing luminescence regulation during quorum
sensing in the bioluminescent bacterium \emph{Vibrio harveyi}. Using a
simplified model for signal transduction in the quorum sensing pathway, we
identify key dimensionless parameters that control the system's response. These
parameters are estimated using experimental data on luminescence phenotypes for
different mutant strains. The corresponding model predictions are consistent
with results from other experiments which did not serve as inputs for
determining model parameters. Furthermore, the proposed framework leads to
novel testable predictions for luminescence phenotypes and for responses of the
network to different perturbations.


  The processes, resulting in the transcription of RNA, are intrinsically
noisy. It was observed experimentally that the synthesis of mRNA molecules is
driven by short, burst-like, events. An accurate prediction of the protein
level often requires one to take these fluctuations into account. Here, we
consider the stochastic model of gene expression regulated by small RNAs. Small
RNA post-transcriptional regulation is achieved by base-pairing with mRNA. We
show that in a strong pairing limit the mRNA steady state distribution can be
determined exactly. The obtained analytical results reveal large discrepancies
in the deterministic description, due to transcriptional bursting. The proposed
approach is expected to work for a large variety of multi-species models.


  At the onset of X Chromosomes Inactivation, the vital process whereby female
mammal cells equalize X products with respect to males, the X chromosomes are
colocalized along their Xic (X-Inactivation Center) regions. The mechanism
inducing recognition and pairing of the X's remains, though, elusive. Starting
from recent discoveries on the molecular factors and on the DNA sequences (the
so-called ``pairing sites'') involved, we dissect the mechanical basis of Xic
colocalization by using a Statistical Physics model. We show that soluble DNA
specific binding molecules, as those experimentally identified, can be indeed
sufficient to induce the spontaneous colocalization of the homologous
chromosomes, but only when their concentration, or chemical affinity, rises
above a threshold value, as a consequence of a thermodynamic phase transition.
We derive the likelihood of pairing and its probability distribution.
Chromosome dynamics has two stages: an initial independent Brownian diffusion
followed, after a characteristic time scale, by recognition and pairing.
Finally, we investigate the effects of DNA deletion/insertions in the region of
pairing sites and compare model predictions to available experimental data.


  This paper explores coevolution and governance of common goods using models
of coevolving biospheres, in which adapting populations must collectively
regulate their planet's climate or face extinction. The results support the
Gaia hypothesis against challenges based on the tragedy of the commons: model
creatures are often able to work together to maintain the common good (a
suitable climate) without being undermined by "free riders." A long-term
dynamics appears in which communities that cannot sustain Gaian cooperation
give way to communities that can. This result provides an argument why a Gaia
scenario should generally be observed, rather than a tragedy of the commons
scenario. Second, a close look at how communities fail reveals failures that do
not fit the tragedy of the commons framework and are better described in terms
of conflict between differently positioned parties, with power over different
aspects of the system. In the context of Norgaard's work, all these
observations can be read as narratives of coevolution relevant to social
communities as well as ecological ones, contrasting with pessimistic scenarios
about common governance and supporting respect for traditional arrangements and
restraint in intervention.


  Many physical and chemical processes, such as folding of biopolymers, are
best described as dynamics on large combinatorial energy landscapes. A concise
approximate description of dynamics is obtained by partitioning the
micro-states of the landscape into macro-states. Since most landscapes of
interest are not tractable analytically, the probabilities of transitions
between macro-states need to be extracted numerically from the microscopic
ones, typically by full enumeration of the state space. Here we propose to
approximate transition probabilities by a Markov chain Monte-Carlo method. For
landscapes of the number partitioning problem and an RNA switch molecule we
show that the method allows for accurate probability estimates with
significantly reduced computational cost.


  Many mechanisms, functions and structures of life have been unraveled.
However, the fundamental driving force that propelled chemical evolution and
led to life has remained obscure. The 2nd law of thermodynamics, written as an
equation of motion, reveals that elemental abiotic matter evolves from the
equilibrium via chemical reactions that couple to external energy toward
complex biotic non-equilibrium systems. Each time a new mechanism of energy
transduction emerges, e.g., by random variation in syntheses, evolution prompts
by punctuation and settles to a stasis when the accessed free energy has been
consumed. The evolutionary course toward an increasingly larger energy
transduction system accumulates a diversity of energy transduction mechanisms,
i.e., species. The rate of entropy increase is identified as the fitness
criterion among the diverse mechanisms which places the theory of evolution by
natural selection on the fundamental thermodynamic principle with no
demarcation line between inanimate and animate.


  The architecture of the eukaryotic genome is characterized by a high degree
of spatial organization. Chromosomes occupy preferred territories correlated to
their state of activity and, yet, displace their genes to interact with remote
sites in complex patterns requiring the orchestration of a huge number of DNA
loci and molecular regulators. Far from random, this organization serves
crucial functional purposes, but its governing principles remain elusive. By
computer simulations of a Statistical Mechanics model, we show how
architectural patterns spontaneously arise from the physical interaction
between soluble binding molecules and chromosomes via collective thermodynamics
mechanisms. Chromosomes colocalize, loops and territories form and find their
relative positions as stable thermodynamic states. These are selected by
"thermodynamic switches" which are regulated by concentrations/affinity of
soluble mediators and by number/location of their attachment sites along
chromosomes. Our "thermodynamic switch model" of nuclear architecture, thus,
explains on quantitative grounds how well known cell strategies of upregulation
of DNA binding proteins or modification of chromatin structure can dynamically
shape the organization of the nucleus.


  An equation, strongly reminiscent of Fisher's equation, is used to model the
response of tsetse populations to proposed control measures in the vicinity of
a game reserve. The model assumes movement is by diffusion and that growth is
logistic. This logistic growth is dependent on an historical population, in
contrast to Fisher's equation which bases it on the present population. The
model therefore takes into account the fact that new additions to the adult fly
population are, in actual fact, the descendents of a population which existed
one puparial duration ago, furthermore, that this puparial duration is
temperature dependent. Artificially imposed mortality is modelled as a
proportion at a constant rate. Fisher's equation is also solved as a formality.
  The temporary imposition of a 2 % $\mathrm{day}^{-1}$ mortality everywhere
outside the reserve for a period of 2 years will have no lasting effect on the
influence of the reserve on either the Glossina austeni or the G. brevipalpis
populations, although it certainly will eradicate tsetse from poor habitat,
outside the reserve. A 5 $\mathrm{km}$-wide barrier with a minimum mortality of
4 % $\mathrm{day}^{-1}$, throughout, will succeed in isolating a worst-case, G.
austeni population and its associated trypanosomiasis from the surrounding
areas. A more optimistic estimate of its mobility suggests a mortality of 2 %
$\mathrm{day}^{-1}$ will suffice. For a given target-related mortality, more
mobile species are found to be more vulnerable to eradication than more
sedentary species, while the opposite is true for containment.


  Natively unfolded proteins lack a well defined three dimensional structure
but have important biological functions, suggesting a re-assignment of the
structure-function paradigm. Many proteins have amino acidic compositions
compatible both with the folded and unfolded status, and belong to a twilight
zone between order and disorder. This makes difficult a dichotomic
classification of protein sequences into folded and natively unfolded ones. In
this methodological paper dichotomic folding indexes are considered:
hydrophobicity-charge, mean packing, mean pairwise energy, Poodle-W and a new
global index, that is called here gVSL2, based on the local disorder predictor
VSL2. The performance of these indexes is evaluated on different datasets.
Poodle-W, gVSL2 and mean pairwise energy have good performance and stability in
all the datasets considered and are combined into a strictly unanimous
combination score SSU, that leaves proteins unclassified when the consensus of
all combined indexes is not reached. The unclassified proteins: i) belong to an
overlap region in the vector space of amino acidic compositions occupied by
both folded and unfolded proteins; ii) are composed by approximately the same
number of order-promoting and disorder-promoting amino acids; iii) have a mean
flexibility intermediate between that of folded and that of unfolded proteins.
These proteins reasonably have physical properties intermediate between those
of folded and those of natively unfolded proteins and their structural
properties and evolutionary history are worth to be investigated.


  Modern high-throughput gene perturbation screens are key technologies at the
forefront of genetic research. Combined with rich phenotypic descriptors they
enable researchers to observe detailed cellular reactions to experimental
perturbations on a genome-wide scale. This review surveys the current
state-of-the-art in analyzing perturbation screens from a network point of
view. We describe approaches to make the step from the parts list to the wiring
diagram by using phenotypes for network inference and integrating them with
complementary data sources. The first part of the review describes methods to
analyze one- or low-dimensional phenotypes like viability or reporter activity;
the second part concentrates on high-dimensional phenotypes showing global
changes in cell morphology, transcriptome or proteome.


  Phylogenetic trees are widely used to display estimates of how groups of
species evolved. Each phylogenetic tree can be seen as a collection of
clusters, subgroups of the species that evolved from a common ancestor. When
phylogenetic trees are obtained for several data sets (e.g. for different
genes), then their clusters are often contradicting. Consequently, the set of
all clusters of such a data set cannot be combined into a single phylogenetic
tree. Phylogenetic networks are a generalization of phylogenetic trees that can
be used to display more complex evolutionary histories, including reticulate
events such as hybridizations, recombinations and horizontal gene transfers.
Here we present the new CASS algorithm that can combine any set of clusters
into a phylogenetic network. We show that the networks constructed by CASS are
usually simpler than networks constructed by other available methods. Moreover,
we show that CASS is guaranteed to produce a network with at most two
reticulations per biconnected component, whenever such a network exists. We
have implemented CASS and integrated it in the freely available Dendroscope
software.


  Networks are useful for describing systems of interacting objects, where the
nodes represent the objects and the edges represent the interactions between
them. The applications include chemical and metabolic systems, food webs as
well as social networks. Lately, it was found that many of these networks
display some common topological features, such as high clustering, small
average path length (small world networks) and a power-law degree distribution
(scale free networks). The topological features of a network are commonly
related to the network's functionality. However, the topology alone does not
account for the nature of the interactions in the network and their strength.
Here we introduce a method for evaluating the correlations between pairs of
nodes in the network. These correlations depend both on the topology and on the
functionality of the network. A network with high connectivity displays strong
correlations between its interacting nodes and thus features small-world
functionality. We quantify the correlations between all pairs of nodes in the
network, and express them as matrix elements in the correlation matrix. From
this information one can plot the correlation function for the network and to
extract the correlation length. The connectivity of a network is then defined
as the ratio between this correlation length and the average path length of the
network. Using this method we distinguish between a topological small world and
a functional small world, where the latter is characterized by long range
correlations and high connectivity. Clearly, networks which share the same
topology, may have different connectivities, based on the nature and strength
of their interactions. The method is demonstrated on metabolic networks, but
can be readily generalized to other types of networks.


  Starting with the results of Li et al. in 1992 there is valuable interest in
finding long range correlations in dna sequences since it raises questions
about the role of introns and intron-containing genes. In the present paper we
studied two sequences. We applied the method of the recurrence quantification
analysis (rqa) that was introduced by Zbilut and Webber in 1994. The
significant result that we have here is that both Lmax and Laminarity exhibit
very large values in non coding respect to coding sequences. Therefore we
suggest that there the claimed higher long range correlations of introns
respect to exons from many authors may be explained here in reason of such
found higher values of Lmax and of Laminarity.


  Chromatin immunoprecipitation combined with DNA microarrays (ChIP-chip) is an
assay for DNA-protein-binding or post-translational chromatin/histone
modifications. As with all high-throughput technologies, it requires a thorough
bioinformatic processing of the data for which there is no standard yet. The
primary goal is the reliable identification and localization of genomic regions
that bind a specific protein. The second step comprises comparison of binding
profiles of functionally related proteins, or of binding profiles of the same
protein in different genetic backgrounds or environmental conditions.
Ultimately, one would like to gain a mechanistic understanding of the effects
of DNA binding events on gene expression. We present a free, open-source R
package Starr that, in combination with the package Ringo, facilitates the
comparative analysis of ChIP-chip data across experiments and across different
microarray platforms. Core features are data import, quality assessment,
normalization and visualization of the data, and the detection of ChIP-enriched
genomic regions. The use of common Bioconductor classes ensures the
compatibility with other R packages. Most importantly, Starr provides methods
for integration of complementary genomics data, e.g., it enables systematic
investigation of the relation between gene expression and dna binding.


  Lattice protein models, as the Hydrophobic-Polar (HP) model, are a common
abstraction to enable exhaustive studies on structure, function, or evolution
of proteins. A main issue is the high number of optimal structures, resulting
from the hydrophobicity-based energy function applied. We introduce an
equivalence relation on protein structures that correlates to the energy
function. We discuss the efficient enumeration of optimal representatives of
the corresponding equivalence classes and the application of the results.


  The simulation of a protein's folding process is often done via stochastic
local search, which requires a procedure to apply structural changes onto a
given conformation. Here, we introduce a constraint-based approach to enumerate
lattice protein structures according to k-local moves in arbitrary lattices.
Our declarative description is much more flexible for extensions than standard
operational formulations. It enables a generic calculation of k-local neighbors
in backbone-only and side chain models. We exemplify the procedure using a
simple hierarchical folding scheme.


  Why reaction rate constants for enzymatic reactions are typically inversely
proportional to fractional power exponents of solvent viscosity remains to be
already a thirty years old puzzle. Available interpretations of the phenomenon
invoke to either a modification of 1. the conventional Kramers' theory or that
of 2. the Stokes law. We show that there is an alternative interpretation of
the phenomenon at which neither of these modifications is in fact
indispensable. We reconcile 1. and 2. with the experimentally observable
dependence. We assume that an enzyme solution in solvent with or without
cosolvent molecules is an ensemble of samples with different values of the
viscosity for the movement of the system along the reaction coordinate. We
assume that this viscosity consists of the contribution with the weight $q$
from cosolvent molecules and that with the weight $1-q$ from protein matrix and
solvent molecules. We introduce heterogeneity in our system with the help of a
distribution over the weight $q$. We verify the obtained solution of the
integral equation for the unknown function of the distribution by direct
substitution. All parameters of the model are related to experimentally
observable values. General formalism is exemplified by the analysis of
literature experimental data for oxygen escape from hemerythin.


  Many human diseases are associated with protein aggregation and fibrillation.
We present experiments on in vitro glucagon fibrillation using total internal
reflection fluorescence microscopy, providing real-time measurements of
single-fibril growth. We find that amyloid fibrils grow in an intermittent
fashion, with periods of growth followed by long pauses. The observed
exponential distributions of stop and growth times support a Markovian model,
in which fibrils shift between the two states with specific rates. Remarkably,
the probability of being in the growing (stopping) state is very close to 1/4
(3/4) in all experiments, even if the rates vary considerably. This finding
suggests the presence of 4 independent conformations of the fibril tip; we
discuss this possibility in terms of the existing structural knowledge.


  Demographic noise has profound effects on evolutionary and population
dynamics, as well as on chemical reaction systems and models of epidemiology.
Such noise is intrinsic and due to the discreteness of the dynamics in finite
populations. We here show that similar noise-sustained trajectories arise in
game dynamical learning, where the stochasticity has a different origin: agents
sample a finite number of moves of their opponents in-between adaptation
events. The limit of infinite batches results in deterministic modified
replicator equations, whereas finite sampling leads to a stochastic dynamics.
The characteristics of these fluctuations can be computed analytically using
methods from statistical physics, and such noise can affect the attractors
significantly, leading to noise-sustained cycling or removing periodic orbits
of the standard replicator dynamics.


  Recently much attention has been devoted to the construction of phylogenetic
networks which generalize phylogenetic trees in order to accommodate complex
evolutionary processes. Here we present an efficient, practical algorithm for
reconstructing level-1 phylogenetic networks - a type of network slightly more
general than a phylogenetic tree - from triplets. Our algorithm has been made
publicly available as the program LEV1ATHAN. It combines ideas from several
known theoretical algorithms for phylogenetic tree and network reconstruction
with two novel subroutines. Namely, an exponential-time exact and a greedy
algorithm both of which are of independent theoretical interest. Most
importantly, LEV1ATHAN runs in polynomial time and always constructs a level-1
network. If the data is consistent with a phylogenetic tree, then the algorithm
constructs such a tree. Moreover, if the input triplet set is dense and, in
addition, is fully consistent with some level-1 network, it will find such a
network. The potential of LEV1ATHAN is explored by means of an extensive
simulation study and a biological data set. One of our conclusions is that
LEV1ATHAN is able to construct networks consistent with a high percentage of
input triplets, even when these input triplets are affected by a low to
moderate level of noise.


  Gene regulatory networks typically have low in-degrees, whereby any given
gene is regulated by few of the genes in the network. What mechanisms might be
responsible for these low in-degrees? Starting with an accepted framework of
the binding of transcription factors to DNA, we consider a simple model of gene
regulatory dynamics. In this model, we show that the constraint of having a
given function leads to the emergence of minimum connectivities compatible with
function. We exhibit mathematically this behavior within a limit of our model
and show that it also arises in the full model. As a consequence, functionality
in these gene networks is parsimonious, i.e., is concentrated on a sparse
number of interactions as measured for instance by their essentiality. Our
model thus provides a simple mechanism for the emergence of sparse regulatory
networks, and leads to very heterogeneous effects of mutations.


  Evolutionary prisoner's dilemma games are studied with players located on
square lattice and random regular graphs defining four neighbors for each one.
The players follow one of the three strategies: tit-for-tat, unconditional
cooperation, and defection. The simplified payoff matrix is characterized by
two parameters: the temptation $b$ to choose defection, and the cost $c$ of
inspection reducing the income of tit-for-tat. The strategy imitation from one
of the neighbors is controlled by pairwise comparison at a fixed level of
noise. Using Monte Carlo simulations and the extended versions of pair
approximation we have evaluated the $b-c$ phase diagrams indicating a rich
plethora of phase transitions between stationary coexistence, absorbing and
oscillatory states, including continuous and discontinuous phase transitions.
By reasonable costs the tit-for-tat strategy prevents extinction of cooperators
across the whole span of $b$ values determining the prisoner's dilemma game,
irrespective of the interaction graph structure. We also demonstrate that the
system can exhibit a repetitive succession of oscillatory and stationary states
upon changing a single payoff value, which highlights the remarkable
sensitivity of cyclical interactions on parameters that define the strength of
dominance.


  Design Theory, a branch of mathematics, was born out of the experimental
statistics research of the population geneticist R. A. Fisher and of Indian
mathematical statisticians in the 1930s. The field combines elements of
combinatorics, finite projective geometries, Latin squares, and a variety of
further mathematical structures, brought together in surprising ways. This
essay will present these structures and ideas as well as how the field came
together, in itself an interesting story.


  The multi-species coalescent provides an elegant theoretical framework for
estimating species trees and species demographics from genetic markers.
Practical applications of the multi-species coalescent model are, however,
limited by the need to integrate or sample over all gene trees possible for
each genetic marker. Here we describe a polynomial-time algorithm that computes
the likelihood of a species tree directly from the markers under a finite-sites
model of mutation, effectively integrating over all possible gene trees. The
method applies to independent (unlinked) biallelic markers such as well-spaced
single nucleotide polymorphisms (SNPs), and we have implemented it in SNAPP, a
Markov chain Monte-Carlo sampler for inferring species trees, divergence dates,
and population sizes. We report results from simulation experiments and from an
analysis of 1997 amplified fragment length polymorphism (AFLP) loci in 69
individuals sampled from six species of {\em Ourisia} (New Zealand native
foxglove).


  We present a combined mean-field and simulation approach to different models
describing the dynamics of classes formed by elements that can appear,
disappear or copy themselves. These models, related to a paradigm
duplication-innovation model known as Chinese Restaurant Process, are devised
to reproduce the scaling behavior observed in the genome-wide repertoire of
protein domains of all known species. In view of these data, we discuss the
qualitative and quantitative differences of the alternative model formulations,
focusing in particular on the roles of element loss and of the specificity of
empirical domain classes.


  The study of fluctuations in gene regulatory networks is extended to the case
of Gaussian colored noise. Firstly, the solution of the corresponding Langevin
equation with colored noise is expressed in terms of an Ito integral. Then, two
important lemmas concerning the variance of an Ito integral and the covariance
of two Ito integrals are shown. Based on the lemmas, we give the general
formulae for the variances and covariance of molecular concentrations for a
regulatory network near a stable equilibrium explicitly. Two examples, the gene
auto-regulatory network and the toggle switch, are presented in details. In
general, it is found that the finite correlation time of noise reduces the
fluctuations and enhances the correlation between the fluctuations of the
molecular components.


  Classical descriptions of enzyme kinetics ignore the physical nature of the
intracellular environment. Main implicit assumptions behind such approaches are
that reactions occur in compartment volumes which are large enough so that
molecular discreteness can be ignored and that molecular transport occurs via
diffusion. Starting from a master equation description of enzyme reaction
kinetics and assuming metabolic steady-state conditions, we derive novel
mesoscopic rate equations which take into account (i) the intrinsic molecular
noise due to the low copy number of molecules in intracellular compartments
(ii) the physical nature of the substrate transport process, i.e. diffusion or
vesicle-mediated transport. These equations replace the conventional
macroscopic and deterministic equations in the context of intracellular
kinetics. The latter are recovered in the limit of infinite compartment
volumes. We find that deviations from the predictions of classical kinetics are
pronounced (hundreds of percent in the estimate for the reaction velocity) for
enzyme reactions occurring in compartments which are smaller than approximately
200nm, for the case of substrate transport to the compartment being mediated
principally by vesicle or granule transport and in the presence of competitive
enzyme inhibitors. This has implications for the common approach of modelling
large intracellular reaction networks using ordinary differential equations and
also for the calculation of the effective dosage of competitive inhibitor
drugs.


  Recent studies attracted the attention on the inherent structure landscape
(ISL) approach as a reduced description of proteins allowing to map their full
thermodynamic properties. However, the analysis has been so far limited to a
single topology of a two-state folding protein, and the simplifying assumptions
of the method have not been examined. In this work, we construct the
thermodynamics of four two-state folding proteins of different sizes and
secondary structure by MD simulations using the ISL method, and critically
examine possible limitations of the method. Our results show that the ISL
approach correctly describes the thermodynamics function, such as the specific
heat, on a qualitative level. Using both analytical and numerical methods, we
show that some quantitative limitations cannot be overcome with enhanced
sampling or the inclusion of harmonic corrections.


  We explore the complex dynamical behavior of two simple predator-prey models
of biological coevolution that on the ecological level account for
interspecific and intraspecific competition, as well as adaptive foraging
behavior. The underlying individual-based population dynamics are based on a
ratio-dependent functional response [W.M. Getz, J. Theor. Biol. 108, 623
(1984)]. Analytical results for fixed-point population sizes in some simple
communities are derived and discussed. In long kinetic Monte Carlo simulations
we find quite robust, approximate 1/f noise in species diversity and population
sizes, as well as power-law distributions for the lifetimes of individual
species and the durations of periods of relative evolutionary stasis. Adaptive
foraging enhances coexistence of species and produces a metastable
low-diversity phase and a stable high-diversity phase.


  Metabolic networks perform some of the most fundamental functions in living
cells, including energy transduction and building block biosynthesis. While
these are the best characterized networks in living systems, understanding
their evolutionary history and complex wiring constitutes one of the most
fascinating open questions in biology, intimately related to the enigma of
life's origin itself. Is the evolution of metabolism subject to general
principles, beyond the unpredictable accumulation of multiple historical
accidents? Here we search for such principles by applying to an artificial
chemical universe some of the methodologies developed for the study of genome
scale models of cellular metabolism. In particular, we use metabolic flux
constraint-based models to exhaustively search for artificial chemistry
pathways that can optimally perform an array of elementary metabolic functions.
Despite the simplicity of the model employed, we find that the ensuing pathways
display a surprisingly rich set of properties, including the existence of
autocatalytic cycles and hierarchical modules, the appearance of universally
preferable metabolites and reactions, and a logarithmic trend of pathway length
as a function of input/output molecule size. Some of these properties can be
derived analytically, borrowing methods previously used in cryptography. In
addition, by mapping biochemical networks onto a simplified carbon atom
reaction backbone, we find that several of the properties predicted by the
artificial chemistry model hold for real metabolic networks. These findings
suggest that optimality principles and arithmetic simplicity might lie beneath
some aspects of biochemical complexity.


  Context. So far, theoretical explanations for body-size patterns in periodic
population dynamics have received little attention. In particular, tuning and
testing of allometric models on empirical data and regressions has not been
carried out yet. Here, oscillations expected from a one-species (delayed
logistic) and a two-species (Rosenzweig-MacArthur) model were compared to
cycles observed in laboratory experiments and field surveys for a wide range of
invertebrates and vertebrates. The parameters in the equations were linked to
body mass, using a consistent set of allometric relationships that was
calibrated on 230 regressions. Oscillation period and amplitude predicted by
the models were validated with data taken from literature. Results. The
collected data showed that cycle times of herbivores scaled to species body
mass with a slope up to 1/4 as expected from the models. With exception of
aquatic herbi-detritivores, intercepts were observed at the level calculated by
the two-species model. Remarkably, oscillation periods were size-independent
for predatory invertebrates, fishes, birds and mammals. Average cycles were of
4 to 5 years, similar to those predicted by the one-species model with a
size-independent delay of 1 year. The consistent difference between herbivores
and the carnivores could be explained by the models from the small parameter
space for consumer-resource cycles in generalist predators. As expected,
amplitudes recorded in the field did not scale to size. Observed oscillation
periods were generally within a factor of about 2 from the values expected from
the models. This demonstrates that a set of slopes and intercepts for age and
density parameters applicable to a wide range of species allows a reasonable
estimate of independently measured cycle times.


  The mechanical unfolding of a simple RNA hairpin and of a 236--bases portion
of the Tetrahymena thermophila ribozyme is studied by means of an Ising--like
model. Phase diagrams and free energy landscapes are computed exactly and
suggest a simple two--state behaviour for the hairpin and the presence of
intermediate states for the ribozyme. Nonequilibrium simulations give the
possible unfolding pathways for the ribozyme, and the dominant pathway
corresponds to the experimentally observed one.


  The Yule (pure-birth) model is the simplest null model of speciation; each
lineage gives rise to a new lineage independently with the same rate $\lambda$.
We investigate the expected length of an edge chosen at random from the
resulting evolutionary tree. In particular, we compare the expected length of a
randomly selected edge with the expected length of a randomly selected pendant
edge. We provide some exact formulae, and show how our results depend slightly
on whether the depth of the tree or the number of leaves is conditioned on, and
whether $\lambda$ is known or is estimated using maximum likelihood.


  A simple explanation for the symmetry and degeneracy of the genetic code has
been suggested. An alternative to the wobble hypothesis has been proposed. This
hypothesis offers explanations for: i) the difference between thymine and
uracil, ii) encoding of tryptophan by only one codon, iii) why E. coli have no
inosine in isoleucine tRNA, but isoleucine is encoded by three codons. The
facts revealed in this study offer a new insight into physical mechanisms of
the functioning of the genetic code.


  Precise temporal coordination of gene expression is crucial for many
developmental processes. One central question in developmental biology is how
such coordinated expression patterns are robustly controlled. During embryonic
development of the Drosophila central nervous system, neural stem cells called
neuroblasts sequentially express a group of genes in a definite order, which
generates the diversity of cell types. By producing all possible regulatory
networks of these genes and examining their expression dynamics numerically, we
identify requisite regulations and predict an unknown factor to reproduce known
expression profiles caused by loss-of-function or overexpression of the genes
in vivo, as well as in the wild type. We then evaluate the stability of the
actual Drosophila network for sequential expression. This network shows the
highest robustness against parameter variations and gene expression
fluctuations among the possible networks that reproduce the expression
profiles. We propose a regulatory module composed of three kinds of regulations
which is responsible for precise sequential expression. The present study
suggests an underlying principle on how biological systems are robustly
designed under functional constraint.


  We study a modified version of the stochastic
susceptible-infected-refractory-susceptible (SIRS) model by employing a
nonlinear (exponential) reinforcement in the contagion rate and no diffusion.
We run simulations for complete and random graphs as well as d-dimensional
hypercubic lattices (for d=3,2,1). For weak nonlinearity, a continuous
nonequilibrium phase transition between an absorbing and an active phase is
obtained, such as in the usual stochastic SIRS model [Joo and Lebowitz, Phys.
Rev. E 70, 036114 (2004)]. However, for strong nonlinearity, the nonequilibrium
transition between the two phases can be discontinuous for d>=2, which is
confirmed by well-characterized hysteresis cycles and bistability. Analytical
mean-field results correctly predict the overall structure of the phase
diagram. Furthermore, contrary to what was observed in a model of phase-coupled
stochastic oscillators with a similar nonlinearity in the coupling [Wood et
al., Phys. Rev. Lett. 96, 145701 (2006)], we did not find a transition to a
stable (partially) synchronized state in our nonlinearly pulse-coupled
excitable elements. For long enough refractory times and high enough
nonlinearity, however, the system can exhibit collective excitability and
unstable stochastic oscillations.


  The Kinesin family of motor proteins are involved in a variety of cellular
processes that transport materials and generate force. With recent advances in
experimental techniques, such as optical tweezers which can probe individual
molecules, there has been an increasing interest in understanding the
mechanisms by which motor proteins convert chemical energy into mechanical
work. Here we present a mathematical model for the chemistry and three
dimensional mechanics of the Kinesin motor protein which captures many of the
force dependent features of the motor. For the elasticity of the tether that
attaches cargo to the motor we develop a method for deriving the non-linear
force-extension relationship from optical trap data. For the Kinesin heads,
cargo, and microscope stage we formulate a three dimensional Brownian Dynamics
model that takes into account excluded volume interactions. To efficiently
compute statistics from the model an algorithm is proposed that uses a two step
protocol that separates the simulation of the mechanical features of the model
from the chemical kinetics of the model. Using this approach for a bead
transported by the motor, the force dependent average velocity and randomness
parameter are computed and compared with the experimental data.


  Multiple genome alignment remains a challenging problem. Effects of
recombination including rearrangement, segmental duplication, gain, and loss
can create a mosaic pattern of homology even among closely related organisms.
We describe a method to align two or more genomes that have undergone
large-scale recombination, particularly genomes that have undergone substantial
amounts of gene gain and loss (gene flux). The method utilizes a novel
alignment objective score, referred to as a sum-of-pairs breakpoint score. We
also apply a probabilistic alignment filtering method to remove erroneous
alignments of unrelated sequences, which are commonly observed in other genome
alignment methods. We describe new metrics for quantifying genome alignment
accuracy which measure the quality of rearrangement breakpoint predictions and
indel predictions. The progressive genome alignment algorithm demonstrates
markedly improved accuracy over previous approaches in situations where genomes
have undergone realistic amounts of genome rearrangement, gene gain, loss, and
duplication. We apply the progressive genome alignment algorithm to a set of 23
completely sequenced genomes from the genera Escherichia, Shigella, and
Salmonella. The 23 enterobacteria have an estimated 2.46Mbp of genomic content
conserved among all taxa and total unique content of 15.2Mbp. We document
substantial population-level variability among these organisms driven by
homologous recombination, gene gain, and gene loss. Free, open-source software
implementing the described genome alignment approach is available from
http://gel.ahabs.wisc.edu/mauve .


  We describe a modification of the TAP method for purification and analysis of
multiprotein complexes, termed here DEF-TAP (for Differential Elution
Fractionation after Tandem Affinity Purification). Its essential new feature is
the use for last purification step of 6XHis-Ni++ interaction, which is
resistant to a variety of harsh washing conditions, including high ionic
strength and presence of organic solvents. This allows us to use various
fractionation schemes before the protease digestion, which is expected to
improve the coverage of the analysed protein mixture and also to provide an
additional insight into the structure of the purified macromolecular complex
and the nature of protein-protein interactions involved. We illustrate our new
approach by analysis of soluble nuclear complexes containing histone H4
purified from HeLa cells. In particular, we observed different fractionation
patterns of HAT1 and RbAp46 proteins as compared to RbAp48 protein, all
identified as interaction partners of H4 histone. In addition, we report all
components of the licensing MCM2-7 complex and the apoptosis-related DAXX
protein among the interaction partners of the soluble H4. Finally, we show that
HAT1 requires N-terminal tail of H4 for its stable association with this
histone.


  We consider the concept of mutual information in ecological networks, and use
this idea to analyse the Tangled Nature model of co-evolution. We show that
this measure of correlation has two distinct behaviours depending on how we
define the network in question: if we consider only the network of viable
species this measure increases, whereas for the whole system it decreases. It
is suggested that these are complimentary behaviours that show how ecosystems
can become both more stable and better adapted.


  Individual expression profiles from EBV transformed cell lines are an
emerging resource for genomic investigation. In this study we characterize the
effects of age, sex, and genetic variation on gene expression by surveying
public datasets of such profiles. We establish that the expression space of
cell lines maintains genetic as well as non-germline information, in an
individual-specific and cross-tissue manner. Age of donor is associated with
the expression of 949 genes in the derived cell line. Age-associated genes
include over-representation of immune-related genes, specifically MHC Class I
genes, a phenomenon that replicates across tissues and organisms. Sex
associated genes in these cell lines include likely candidates, such as genes
that escape X-inactivation,testis specific expressed genes, androgen and
estrogen specific genes, but also gene families previously unknown to be sex
associated such as common microRNA targets (MIR-490, V_ARP1_01, MIR-489).
Finally, we report 494 transcripts whose expression levels are associated with
a genetic variant in cis, overlapping and validating previous reports.
Incorporating age in analysis of association facilitates additional discovery
of trans-acting regulatory genetic variants. Our findings promote expression
profiling of transformed cell lines as a vehicle for understanding cellular
systems beyond the specific lines.


  Evolutionary game theory has been successfully used to investigate the
dynamics of systems, in which many entities have competitive interactions. From
a physics point of view, it is interesting to study conditions under which a
coordination or cooperation of interacting entities will occur, be it spins,
particles, bacteria, animals, or humans. Here, we analyze the case, where the
entities are heterogeneous, particularly the case of two populations with
conflicting interactions and two possible states. For such systems, explicit
mathematical formulas will be determined for the stationary solutions and the
associated eigenvalues, which determine their stability. In this way, four
different types of system dynamics can be classified, and the various kinds of
phase transitions between them will be discussed. While these results are
interesting from a physics point of view, they are also relevant for social,
economic, and biological systems, as they allow one to understand conditions
for (1) the breakdown of cooperation, (2) the coexistence of different
behaviors ("subcultures"), (2) the evolution of commonly shared behaviors
("norms"), and (4) the occurrence of polarization or conflict. We point out
that norms have a similar function in social systems that forces have in
physics.


  Mucin glycoprotein consist of tandem repeating glycosylated regions flanked
by non-repetitive protein domains with little glycosylation. These
non-repetitive domains are involved in the pH dependent gelation of gastric
mucin, which is essential to protecting the stomach from autodigestion. We have
examined the folding of the non-repetitive sequence of von Willebrand factor
vWF-C1 domain (67 amino acids) and PGM 2X (242 amino acids) at neutral and low
pH using Discrete Molecular Dynamics. A four-bead protein model with hydrogen
bonding and amino acid-specific hydrophobic/hydrophilic and electrostatic
interactions of side chains) was used. The simulations reveal that the distant
N- and C-terminal regions form salt-bridges at neutral pH giving a relatively
compact folded structure. At low pH, the salt bridges break giving a more open
and extended structure. The calculated average value of the beta-strand
increases from 0.23 at neutral pH to 0.36 at low pH in very good agreement with
CD data. Simulations of vWF C1 show 4-6 beta strands separated by turns/loops
and we found that pH did not affect significantly the folded structure. The
average beta-strand structure of 0.32 was again in very good agreement with the
CD results.


  We explore aspects of the community structures generated by a simple
predator-prey model of biological coevolution, using large-scale kinetic Monte
Carlo simulations. The model accounts for interspecies and intraspecies
competition for resources, as well as adaptive foraging behavior. It produces a
metastable low-diversity phase and a stable high-diversity phase. The
structures and joint indegree-outdegree distributions of the food webs
generated in the latter phase are discussed.


  The idea of 'date' and 'party' hubs has been influential in the study of
protein-protein interaction networks. Date hubs display low co-expression with
their partners, whilst party hubs have high co-expression. It was proposed that
party hubs are local coordinators whereas date hubs are global connectors. Here
we show that the reported importance of date hubs to network connectivity can
in fact be attributed to a tiny subset of them. Crucially, these few, extremely
central, hubs do not display particularly low expression correlation,
undermining the idea of a link between this quantity and hub function. The
date/party distinction was originally motivated by an approximately bimodal
distribution of hub co-expression; we show that this feature is not always
robust to methodological changes. Additionally, topological properties of hubs
do not in general correlate with co-expression. Thus, we suggest that a
date/party dichotomy is not meaningful and it might be more useful to conceive
of roles for protein-protein interactions rather than individual proteins. We
find significant correlations between interaction centrality and the functional
similarity of the interacting proteins.


  A new theoretical survey of proteins' resistance to constant speed stretching
is performed for a set of 17 134 proteins as described by a structure-based
model. The proteins selected have no gaps in their structure determination and
consist of no more than 250 amino acids. Our previous studies have dealt with
7510 proteins of no more than 150 amino acids. The proteins are ranked
according to the strength of the resistance. Most of the predicted top-strength
proteins have not yet been studied experimentally. Architectures and folds
which are likely to yield large forces are identified. New types of potent
force clamps are discovered. They involve disulphide bridges and, in
particular, cysteine slipknots. An effective energy parameter of the model is
estimated by comparing the theoretical data on characteristic forces to the
corresponding experimental values combined with an extrapolation of the
theoretical data to the experimental pulling speeds. These studies provide
guidance for future experiments on single molecule manipulation and should lead
to selection of proteins for applications. A new class of proteins, involving
cystein slipknots, is identified as one that is expected to lead to the
strongest force clamps known. This class is characterized through molecular
dynamics simulations.


  The nuclear-cytoplasmic shuttling of NF-kB is characterized by damped
oscillations of the nuclear concentration with a time period of around 1-2
hours. The NF-kB network contains several feedback loops modulating the overall
response of NF-kB activity. While IkBa is known to drive and IkBe is known to
dampen the oscillations, the precise role of A20 negative feedback remains to
be elucidated. Here we propose a model of the NF-kB system focusing on three
negative feedback loops (IkBa, IkBe and A20) which capture the experimentally
observed responses in wild-type and knockout cells. We find that A20, like
IkBe, efficiently dampens the oscillations albeit through a distinct mechanism.
In addition, however, we have discovered a new functional role of A20 by which
it controls the oscillation period of nuclear NF-kB. The design based on three
nested feedback loops allows independent control of period and amplitude decay
in the oscillatory response. Based on these results we predict that adjusting
the expression level of A20, e.g. by siRNA, the period can be changed by up to
a factor 2.


  We introduce a coarse-grained rigid nucleotide model of DNA that reproduces
the basic thermodynamics of short strands: duplex hybridization,
single-stranded stacking and hairpin formation, and also captures the essential
structural properties of DNA: the helical pitch, persistence length and
torsional stiffness of double-stranded molecules, as well as the comparative
flexibility of unstacked single strands. We apply the model to calculate the
detailed free-energy landscape of one full cycle of DNA 'tweezers', a simple
machine driven by hybridization and strand displacement.


  We have simulated the evolution of sexually reproducing populations composed
of individuals represented by diploid genomes. A series of eight bits formed an
allele occupying one of 128 loci of one haploid genome (chromosome). The
environment required a specific activity of each locus, this being the sum of
the activities of both alleles located at the corresponding loci on two
chromosomes. This activity is represented by the number of bits set to zero. In
a constant environment the best fitted individuals were homozygous with
alleles' activities corresponding to half of the environment requirement for a
locus (in diploid genome two alleles at corresponding loci produced a proper
activity). Changing the environment under a relatively low recombination rate
promotes generation of more polymorphic alleles. In the heterozygous loci,
alleles of different activities complement each other fulfilling the
environment requirements. Nevertheless, the genetic pool of populations evolves
in the direction of a very restricted number of complementing haplotypes and a
fast changing environment kills the population. If simulations start with all
loci heterozygous, they stay heterozygous for a long time.


  This paper is a review of the VBL project, where we develop and test a
numerical simulator of tumor spheroids.


  Tree reconstruction methods are often judged by their accuracy, measured by
how close they get to the true tree. Yet most reconstruction methods like ML do
not explicitly maximize this accuracy. To address this problem, we propose a
Bayesian solution. Given tree samples, we propose finding the tree estimate
which is closest on average to the samples. This ``median'' tree is known as
the Bayes estimator (BE). The BE literally maximizes posterior expected
accuracy, measured in terms of closeness (distance) to the true tree. We
discuss a unified framework of BE trees, focusing especially on tree distances
which are expressible as squared euclidean distances. Notable examples include
Robinson--Foulds distance, quartet distance, and squared path difference. Using
simulated data, we show Bayes estimators can be efficiently computed in
practice by hill climbing. We also show that Bayes estimators achieve higher
accuracy, compared to maximum likelihood and neighbor joining.


  We theorize that phylogenetic profiles provide a quantitative method that can
relate the structural and functional properties of proteins, as well as their
evolutionary relationships. A key feature of phylogenetic profiles is the
interoperable data format (e.g. alignment information, physiochemical
information, genomic information, etc). Indeed, we have previously demonstrated
Position Specific Scoring Matrices (PSSMs) are an informative M-dimension which
can be scored from quantitative measure of embedded or unmodified sequence
alignments. Moreover, the information obtained from these alignments is
informative, even in the twilight zone of sequence similarity (<25%
identity)(1-5). Although powerful, our previous embedding strategy suffered
from contaminating alignments(embedded AND unmodified) and computational
expense. Herein, we describe the logic and algorithmic process for a heuristic
embedding strategy (Adaptive GDDA-BLAST, Ada-BLAST). Ada-BLAST on average up to
~19-fold faster and has similar sensitivity to our previous method. Further, we
provide data demonstrating the benefits of embedded alignment measurements for
isolating secondary structural elements and the classifying
transmembrane-domain structure/function. We theorize that sequence-embedding is
one of multiple ways that low-identity alignments can be measured and
incorporated into high-performance PSSM-based phylogenetic profiles.


  Just as physicists strive to develop a TOE (theory of everything), which
explains and unifies the physical laws of the universe, the life-scientist
wishes to uncover the TOE as it relates to cellular systems. This can only be
achieved with a quantitative platform that can comprehensively deduce and
relate protein structure, functional, and evolution of genomes and proteomes in
a comparative fashion. Were this perfected, proper analyses would start to
uncover the underlying physical laws governing the emergent behavior of
biological systems and the evolutionary pressures responsible for functional
innovation. In the near term, such methodology would allow the vast quantities
of uncharacterized (e.g. metagenomic samples) primary amino acid sequences to
be rapidly decoded. Analogous to natural products found in the Amazon, genomes
of living organisms contain large numbers of proteins that would prove useful
as new therapeutics for human health, energy sources, and/or waste management
solutions if they could be identified and characterized. We previously
theorized that phylogenetic profiles could provide a quantitative platform for
obtaining unified measures of structure, function, and evolution (SF&E)(1). In
the present manuscript, we present data that support this theory and
demonstrates how refinements of our analysis algorithms improve the performance
of phylogenetic profiles for deriving structural/functional relationships.


  In this study, we explored the utility of a descriptive and predictive
bionetwork model for phospholipase C-coupled calcium signaling pathways, built
with non-kinetic experimental information. Boolean models generated from these
data yield oscillatory activity patterns for both the endoplasmic reticulum
resident inositol-1,4,5-trisphosphate receptor (IP3R) and the plasma-membrane
resident canonical transient receptor potential channel 3 (TRPC3). These
results are specific as randomization of the Boolean operators ablates
oscillatory pattern formation. Furthermore, knock-out simulations of the IP3R,
TRPC3, and multiple other proteins recapitulate experimentally derived results.
The potential of this approach can be observed by its ability to predict
previously undescribed cellular phenotypes using in vitro experimental data.
Indeed our cellular analysis of the developmental and calcium-regulatory
protein, DANGER1a, confirms the counter-intuitive predictions from our Boolean
models in two highly relevant cellular models. Based on these results, we
theorize that with sufficient legacy knowledge and/or computational biology
predictions, Boolean networks provide a robust method for predictive-modeling
of any biological system.


  We have simulated the evolution of age structured populations whose
individuals represented by their diploid genomes were distributed on a square
lattice. The environmental conditions on the whole territory changed
simultaneously in the same way by switching on or off some requirements.
Mutations accumulated in the genes dispensable during a given period of time
were neutral, but they could cause a genetic death of individuals if the
environment required their functions again. Populations survived due to
retaining some surplus of genetic information in the individual genomes. The
changes of the environment caused the fluctuations of the population size.
Since the simulations were performed with individuals spatially distributed on
the lattice and the maximal distance between mating partners was set as a
parameter of the model, the inbreeding coefficient in populations changed
unevenly, following the fluctuation of population size and enhancing the
speciation phenomena.


  We construct a global bifurcation diagram of the plane differential system $$
{l} \dot x = x(1-x)-x y/(a+x^2), \dot y = y(\delta-\beta y/x), x(t)>0, y(t)>0,
a>0, \delta>0, \beta>0, $$ which describes the predator-prey interaction.


  We examine the mechanism of bundling of cytoskeletal actin filaments by two
representative bundling proteins, fascin and espin. Small-angle X-ray studies
show that increased binding from linkers drives a systematic \textit{overtwist}
of actin filaments from their native state, which occurs in a linker-dependent
fashion. Fascin bundles actin into a continuous spectrum of intermediate twist
states, while espin only allows for untwisted actin filaments and
fully-overtwisted bundles. Based on a coarse-grained, statistical model of
protein binding, we show that the interplay between binding geometry and the
intrinsic \textit{flexibility} of linkers mediates cooperative binding in the
bundle. We attribute the respective continuous/discontinous bundling mechanisms
of fascin/espin to differences in the stiffness of linker bonds themselves.


  The complex pattern of presence and absence of many genes across different
species provides tantalising clues as to how genes evolved through the
processes of gene genesis, gene loss and lateral gene transfer (LGT). The
extent of LGT, particularly in prokaryotes, and its implications for creating a
`network of life' rather than a `tree of life' is controversial. In this paper,
we formally model the problem of quantifying LGT, and provide exact
mathematical bounds, and new computational results. In particular, we
investigate the computational complexity of quantifying the extent of LGT under
the simple models of gene genesis, loss and transfer on which a recent
heuristic analysis of biological data relied. Our approach takes advantage of a
relationship between LGT optimization and graph-theoretical concepts such as
tree width and network flow.


  It is a long-standing question in origin-of-life research whether the
information content of replicating molecules can be maintained in the presence
of replication errors. Extending standard quasispecies models of non-enzymatic
replication, we analyze highly specific enzymatic self-replication mediated
through an otherwise neutral recognition region, which leads to
frequency-dependent replication rates. We find a significant reduction of the
maximally tolerable error rate, because the replication rate of the fittest
molecules decreases with the fraction of functional enzymes. Our analysis is
extended to hypercyclic couplings as an example for catalytic networks.


  Evolutionary game dynamics is one of the most fruitful frameworks for
studying evolution in different disciplines, from Biology to Economics. Within
this context, the approach of choice for many researchers is the so-called
replicator equation, that describes mathematically the idea that those
individuals performing better have more offspring and thus their frequency in
the population grows. While very many interesting results have been obtained
with this equation in the three decades elapsed since it was first proposed, it
is important to realize the limits of its applicability. One particularly
relevant issue in this respect is that of non-mean-field effects, that may
arise from temporal fluctuations or from spatial correlations, both neglected
in the replicator equation. This review discusses these temporal and spatial
effects focusing on the non-trivial modifications they induce when compared to
the outcome of replicator dynamics. Alongside this question, the hypothesis of
linearity and its relation to the choice of the rule for strategy update is
also analyzed. The discussion is presented in terms of the emergence of
cooperation, as one of the current key problems in Biology and in other
disciplines.


  Protein translation is a multistep process which can be represented as a
cascade of biochemical reactions (initiation, ribosome assembly, elongation,
etc.), the rate of which can be regulated by small non-coding microRNAs through
multiple mechanisms. It remains unclear what mechanisms of microRNA action are
most dominant: moreover, many experimental reports deliver controversal
messages on what is the concrete mechanism actually observed in the experiment.
Parker and Nissan (Parker and Nissan, RNA, 2008) demonstrated that it is
impossible to distinguish alternative biological hypotheses using the steady
state data on the rate of protein synthesis. For their analysis they used two
simple kinetic models of protein translation. In contrary, we show that
dynamical data allow to discriminate some of the mechanisms of microRNA action.
We demonstrate this using the same models as in (Parker and Nissan, RNA, 2008)
for the sake of comparison but the methods developed (asymptotology of
biochemical networks) can be used for other models. As one of the results of
our analysis, we formulate a hypothesis that the effect of microRNA action is
measurable and observable only if it affects the dominant system
(generalization of the limiting step notion for complex networks) of the
protein translation machinery. The dominant system can vary in different
experimental conditions that can partially explain the existing controversy of
some of the experimental data.


  Human-made structures are increasingly found in marine coastal habitats. The
aim of the present study was to explore whether urban coastal structures can
affect the genetic variation of hard-bottom species. We conducted a population
genetic analysis on the limpet Patella caerulea sampled in both natural and
artificial habitats along the Adriatic coast. Five microsatellite loci were
used to test for differences in genetic diversity and structure among samples.
Three microsatellite loci showed strong Hardy-Weinberg disequilibrium likely
linked with the presence of null alleles. Genetic diversity was significantly
higher in natural habitat than in artificial habitat. A weak but significant
differentiation over all limpet samples was observed, but not related to the
type of habitat. While the exact causes of the differences in genetic diversity
deserve further investigation, these results clearly point that the expansion
of urban structures can lead to genetic diversity loss at regional scales.


  We describe pattern formation in ecological systems using a version of the
classical Lotka-Volterra model characterized by a spatial scale which controls
the predator-prey interaction range. Analytical and simulational results show
that patterns can emerge in some regions of the parameters space where the
instability is driven by the range of the interaction. The individual-based
implementation captures realistic ecological features. In fact, spatial
structures emerge in an erratic oscillatory regime which can contemplate
predators' extinction.


  Few of food web theory hypotheses/predictions can be readily tested using
empirical data. An exception is represented by simple probabilistic models for
food web structure, for which the likelihood has been derived. Here I test the
performance of a more complex model for food web structure that is grounded in
the allometric scaling of interactions with body size and the theory of optimal
foraging (Allometric Diet Breadth Model - ADBM). This deterministic model has
been evaluated measuring the fraction of trophic relations correctly predicted.
I contrast this value with that produced by simpler models based on body sizes
and find that the data does not favor the more complex model: the information
on allometric scaling and optimal foraging does not significantly increase the
fit to the data. Also, I take a different approach and compute the p-value for
the fraction of trophic interactions correctly predicted by ADBM with respect
to three probabilistic null models. I find that the ADBM is clearly better at
predicting links than random graphs, but other models can do even better.
Although optimal foraging and allometric scaling could improve our
understanding of food webs, the models need to be ameliorated to find support
in the data.


  We continue former work on the modeling of potential effects of Gamma Ray
Bursts on Phanerozoic Earth. We focus on global biospheric effects of ozone
depletion and show a first modeling of the spectral reduction of light by NO2
formed in the stratosphere. We also illustrate the current complexities
involved in the prediction of how terrestrial ecosystems would respond to this
kind of burst. We conclude that more biological field and laboratory data are
needed to reach even moderate accuracy in this modeling


  Third political parties are influential in shaping American politics. In this
work we study the spread of a third party ideology in a voting population where
we assume that party members/activists are more influential in recruiting new
third party voters than non-member third party voters. The study uses an
epidemiological metaphor to develop a theoretical model with nonlinear ordinary
differential equations as applied to a case study, the Green Party. Considering
long-term behavior, we identify three threshold parameters in our model that
describe the different possible scenarios for the political party and its
spread. We also apply the model to the study of the Green Party's growth using
voting and registration data in six states and the District of Columbia to
identify and explain trends over the past decade. Our system produces a
backward bifurcation that helps identify conditions under which a sufficiently
dedicated activist core can enable a third party to thrive, under conditions
which would not normally allow it to arise. Our results explain the critical
role activists play in sustaining grassroots movements under adverse
conditions.


  Investigations of molecular bonds between single molecules and molecular
complexes by the dynamic force spectroscopy are subject to large fluctuations
at nanoscale and possible other aspecific binding, which mask the experimental
output. Big efforts are devoted to develop methods for effective selection of
the relevant experimental data, before taking the quantitative analysis of bond
parameters. Here we present a methodology which is based on the application of
graph theory. The force-distance curves corresponding to repeated pulling
events are mapped onto their correlation network (mathematical graph). On these
graphs the groups of similar curves appear as topological modules, which are
identified using the spectral analysis of graphs. We demonstrate the approach
by analyzing a large ensemble of the force-distance curves measured on:
ssDNA-ssDNA, peptide-RNA (system from HIV1), and peptide-Au surface. Within our
data sets the methodology systematically separates subgroups of curves which
are related to different intermolecular interactions and to spatial
arrangements in which the molecules are brought together and/or pulling speeds.
This demonstrates the sensitivity of the method to the spatial degrees of
freedom, suggesting potential applications in the case of large molecular
complexes and situations with multiple binding sites.


  We present an Evolutionary Placement Algorithm (EPA) for the rapid assignment
of sequence fragments (short reads) to branches of a given phylogenetic tree
under the Maximum Likelihood (ML) model. The accuracy of the algorithm is
evaluated on several real-world data sets and compared to placement by
pair-wise sequence comparison, using edit distances and BLAST.
  We test two versions of the placement algorithm, one slow and more accurate
where branch length optimization is conducted for each short read insertion and
a faster version where the branch lengths are approximated at the insertion
position. For the slow version, additional heuristic techniques are explored
that almost yield the same run time as the fast version, with only a small loss
of accuracy. When those additional heuristics are employed the run time of the
more accurate algorithm is comparable to that of a simple BLAST search for data
sets with a high number of short query sequences. Moreover, the accuracy of the
Evolutionary Placement Algorithm is significantly higher, in particular when
the taxon sampling of the reference topology is sparse or inadequate. Our
algorithm, which has been integrated into RAxML, therefore provides an equally
fast but more accurate alternative to BLAST for phylogeny-aware analysis of
short-read sequence data.


  The effect of target molecule depletion from the supernatant solution is
incorporated into a physico-chemical model of hybridisation on oligonucleotide
microarrays. Two possible regimes are identified: local depletion, in which
depletion by a given probe feature only affects that particular probe, and
global depletion, in which all features responding to a given target species
are affected. Examples are given of two existing spike-in data sets
experiencing measurable effects of target depletion. The first of these, from
an experiment by Suzuki et al. using custom built arrays with a broad range of
probe lengths and mismatch positions, is verified to exhibit local and not
global depletion. The second dataset, the well known Affymetrix HGU133a latin
square experiment is shown to be very well explained by a global depletion
model. It is shown that microarray calibrations relying on Langmuir isotherm
models which ignore depletion effects will significantly underestimate specific
target concentrations. It is also shown that a combined analysis of perfect
match and mismatch probe signals in terms of a simple graphical summary, namely
the hook curve method, can discriminate between cases of local and global
depletion.


  Parts of DNA sequences known as exons and introns play very different role in
coding and storage of genetic information. Here we show that their conducting
properties are also very different. Taking into account long-range correlations
among four basic nucleotides that form double-stranded DNA sequence, we
calculate electron localization length for exon and intron regions. Analyzing
different DNA molecules, we obtain that the exons have narrow bands of extended
states, unlike the introns where all the states are well localized. The band of
extended states is due to a specific form of the binary correlation function of
the sequence of basic DNA nucleotides.


  I tackle the problem of partitioning a sequence into homogeneous segments,
where homogeneity is defined by a set of Markov models. The problem is to study
the likelihood that a sequence is divided into a given number of segments.
Here, the moments of this likelihood are computed through an efficient
algorithm. Unlike methods involving Hidden Markov Models, this algorithm does
not require probability transitions between the models. Among many possible
usages of the likelihood, I present a maximum \textit{a posteriori} probability
criterion to predict the number of homogeneous segments into which a sequence
can be divided, and an application of this method to find CpG islands.


  Languages evolve over time in a process in which reproduction, mutation and
extinction are all possible, similar to what happens to living organisms. Using
this similarity it is possible, in principle, to build family trees which show
the degree of relatedness between languages.
  The method used by modern glottochronology, developed by Swadesh in the
1950s, measures distances from the percentage of words with a common historical
origin. The weak point of this method is that subjective judgment plays a
relevant role.
  Recently we proposed an automated method that avoids the subjectivity, whose
results can be replicated by studies that use the same database and that
doesn't require a specific linguistic knowledge. Moreover, the method allows a
quick comparison of a large number of languages.
  We applied our method to the Indo-European and Austronesian families,
considering in both cases, fifty different languages. The resulting trees are
similar to those of previous studies, but with some important differences in
the position of few languages and subgroups. We believe that these differences
carry new information on the structure of the tree and on the phylogenetic
relationships within families.


  The idea of measuring distance between languages seems to have its roots in
the work of the French explorer Dumont D'Urville (D'Urville 1832). He collected
comparative words lists of various languages during his voyages aboard the
Astrolabe from 1826 to1829 and, in his work about the geographical division of
the Pacific, he proposed a method to measure the degree of relation among
languages. The method used by modern glottochronology, developed by Morris
Swadesh in the 1950s (Swadesh 1952), measures distances from the percentage of
shared cognates, which are words with a common historical origin. Recently, we
proposed a new automated method which uses normalized Levenshtein distance
among words with the same meaning and averages on the words contained in a
list. Another classical problem in glottochronology is the study of the
stability of words corresponding to different meanings. Words, in fact, evolve
because of lexical changes, borrowings and replacement at a rate which is not
the same for all of them. The speed of lexical evolution is different for
different meanings and it is probably related to the frequency of use of the
associated words (Pagel et al. 2007). This problem is tackled here by an
automated methodology only based on normalized Levenshtein distance.


  A theory is presented of the longitudinal dispersion of DNA under equilibrium
confined in a nanochannel. Orientational fluctuations of the DNA chain build up
to give rise to substantial fluctuations of the coil in the longitudinal
direction of the channel. The translational and orientational degrees of
freedom of the polymer are described by the Green function satisfying the usual
Fokker-Planck equation. It is argued that this is analogous to the transport
equation occurring in the theory of convective diffusion of particles in pipe
flow. Moreover, Taylor's method may be used to reduce the Fokker-Planck
equation to a diffusion equation for long DNA although subtleties arise
connected with the orientational distribution of segments within the channel.
The longitudinal "step length" turns out to be proportional to the typical
angle of a DNA segment to the sixth power. The dispersion is underestimated
compared to experiment, probably because the harmonic approximation is used to
describe the polymer confinement.