  Motivation: Millions of genes in the modern species belong to only thousands
of `gene families'. A gene family includes instances of the same gene in
different species (orthologs) and duplicate genes in the same species
(paralogs). Genes are gained and lost during evolution. With advances in
sequencing technology, researchers are able to investigate the important roles
of gene duplications and losses in adaptive evolution. Because of gene complex
evolution, ortholog identification is a basic but difficult task in comparative
genomics. A key method for the task is to use an explicit model of the
evolutionary history of the genes being studied, called the gene (family) tree.
It compares the gene tree with the evolutionary history of the species in which
the genes reside, called the species tree, using the procedure known as tree
reconciliation. Reconciling binary gene and specific trees is simple. However,
both gene and species trees may be non-binary in practice and thus tree
reconciliation presents challenging problems. Here, non-binary gene and species
tree reconciliation is studied in a binary refinement model.
  Results: The problem of reconciling arbitrary gene and species trees is
proved NP-hard even for the duplication cost. We then present the first
efficient method for reconciling a non-binary gene tree and a non-binary
species tree. It attempts to find binary refinements of the given gene and
species trees that minimize reconciliation cost. Our algorithms have been
implemented into a software to support quick automated analysis of large data
sets.
  Availability: The program, together with the source code, is available at its
online server http://phylotoo.appspot.com.


  The rock-paper-scissors game is a model example of the on-going cyclic
turnover typical of many ecosystems, ranging from the terrestrial and aquatic
to the microbial. Here we explore the evolution of a rock-paper-scissors system
where three species compete for space. The species are allowed to mutate and
change the speed by which they invade one another. In the case when all species
have similar mutation rates, we observe a perpetual arms race where no single
species prevails. When only two species mutate, their aggressions increase
indefinitely until the ecosystem collapses and only the non-mutating species
survives. Finally we show that when only one species mutates, group selection
removes individual predators with the fastest growth rates, causing the growth
rate of the species to stabilize. We explain this group selection
quantitatively.


  Bundles of filaments and motors are central to contractility in cells. The
classic example is striated muscle, where actomyosin contractility is mediated
by highly organized sarcomeres which act as fundamental contractile units.
However, many contractile bundles in vivo and in vitro lack sarcomeric
organization. Here we propose a model for how contractility can arise in
actomyosin bundles without sarcomeric organization and validate its predictions
with experiments on a reconstituted system. In the model, internal stresses in
frustrated arrangements of motors with diverse velocities cause filaments to
buckle, leading to overall shortening. We describe the onset of buckling in the
presence of stochastic actin-myosin detachment and predict that
buckling-induced contraction occurs in an intermediate range of motor
densities. We then calculate the size of the "contractile units" associated
with this process. Consistent with these results, our reconstituted actomyosin
bundles contract at relatively high motor density, and we observe buckling at
the predicted length scale.


  Familial Searching is the process of searching in a DNA database for
relatives of a certain individual. It is well known that in order to evaluate
the genetic evidence in favour of a certain given form of relatedness between
two individuals, one needs to calculate the appropriate likelihood ratio, which
is in this context called a Kinship Index. Suppose that the database contains,
for a given type of relative, at most one related individual. Given prior
probabilities for being the relative for all persons in the database, we derive
the likelihood ratio for each database member in favour of being that relative.
This likelihood ratio takes all the Kinship Indices between the target
individual and the members of the database into account. We also compute the
corresponding posterior probabilities. We then discuss two methods to select a
subset from the database that contains the relative with a known probability,
or at least a useful lower bound thereof. One method needs prior probabilities
and yields posterior probabilities, the other does not. We discuss the relation
between the approaches, and illustrate the methods with familial searching
carried out in the Dutch National DNA Database.


  We examine the properties of a recently proposed model for antigenic
variation in malaria which incorporates multiple epitopes and both long-lasting
and transient immune responses. We show that in the case of a vanishing decay
rate for the long-lasting immune response, the system exhibits the so-called
"bifurcations without parameters" due to the existence of a hypersurface of
equilibria in the phase space. When the decay rate of the long-lasting immune
response is different from zero, the hypersurface of equilibria degenerates,
and a multitude of other steady states are born, many of which are related by a
permutation symmetry of the system. The robustness of the fully symmetric state
of the system was investigated by means of numerical computation of transverse
Lyapunov exponents. The results of this exercise indicate that for a vanishing
decay of long-lasting immune response, the fully symmetric state is not robust
in the substantial part of the parameter space, and instead all variants
develop their own temporal dynamics contributing to the overall time evolution.
At the same time, if the decay rate of the long-lasting immune response is
increased, the fully symmetric state can become robust provided the growth rate
of the long-lasting immune response is rapid.


  We introduce a formalism for computing bond percolation properties of a class
of correlated and clustered random graphs. This class of graphs is a
generalization of the Configuration Model where nodes of different types are
connected via different types of hyperedges, edges that can link more than 2
nodes. We argue that the multitype approach coupled with the use of clustered
hyperedges can reproduce a wide spectrum of complex patterns, and thus enhances
our capability to model real complex networks. As an illustration of this
claim, we use our formalism to highlight unusual behaviors of the size and
composition of the components (small and giant) in a synthetic, albeit
realistic, social network.


  In this note we provide the necessary and sufficient conditions to uniquely
reconstruct an oncogenetic tree.


  Microarrays are made it possible to simultaneously monitor the expression
profiles of thousands of genes under various experimental conditions.
Identification of co-expressed genes and coherent patterns is the central goal
in microarray or gene expression data analysis and is an important task in
Bioinformatics research. In this paper, K-Means algorithm hybridised with
Cluster Centre Initialization Algorithm (CCIA) is proposed Gene Expression
Data. The proposed algorithm overcomes the drawbacks of specifying the number
of clusters in the K-Means methods. Experimental analysis shows that the
proposed method performs well on gene Expression Data when compare with the
traditional K- Means clustering and Silhouette Coefficients cluster measure.


  Assessing the statistical power to detect susceptibility variants plays a
critical role in GWA studies both from the prospective and retrospective points
of view. Power is empirically estimated by simulating phenotypes under a
disease model H1. For this purpose, the "gold" standard consists in simulating
genotypes given the phenotypes (e.g. Hapgen). We introduce here an alternative
approach for simulating phenotypes under H1 that does not require generating
new genotypes for each simulation. In order to simulate phenotypes with a fixed
total number of cases and under a given disease model, we suggest three
algorithms: i) a simple rejection algorithm; ii) a numerical Markov Chain
Monte-Carlo (MCMC) approach; iii) and an exact and efficient backward sampling
algorithm. In our study, we validated the three algorithms both on a
toy-dataset and by comparing them with Hapgen on a more realistic dataset. As
an application, we then conducted a simulation study on a 1000 Genomes Project
dataset consisting of 629 individuals (314 cases) and 8,048 SNPs from
Chromosome X. We arbitrarily defined an additive disease model with two
susceptibility SNPs and an epistatic effect. The three algorithms are
consistent, but backward sampling is dramatically faster than the other two.
Our approach also gives consistent results with Hapgen. Using our application
data, we showed that our limited design requires a biological a priori to limit
the investigated region. We also proved that epistatic effects can play a
significant role even when simple marker statistics (e.g. trend) are used. We
finally showed that the overall performance of a GWA study strongly depends on
the prevalence of the disease: the larger the prevalence, the better the power.


  Large multidimensionality of high-throughput datasets pertaining to cell
signaling and gene regulation renders it difficult to extract mechanisms
underlying the complex kinetics involving various biochemical compounds (e.g.,
proteins, lipids). Data-driven models often circumvent this difficulty by using
pair correlations of the protein expression levels to produce a small numbers
(<10) of principal components, each a linear combination of the concentrations,
to successfully model how cells respond to different stimuli. However, it is
not understood if this reduction is specific to a particular biological system
or to nature of the stimuli used in these experiments. We study temporal
changes in pair correlations described by the covariance matrix between
different molecular species that evolve following deterministic mass action
kinetics in large biologically relevant reaction networks and show that this
dramatic reduction of dimensions (from hundreds to <5) arises from the strong
correlations between different species at any time and is in sensitive of the
form of the nonlinear interactions, network architecture and values of rate
constants and concentrations over a wide range. We relate temporal changes in
the eigenvalue spectrum of the covariance matrix to low-dimensional, local
changes in directions of the trajectory embedded in much larger dimensions
using elementary differential geometry. We illustrate how to extract
biologically relevant insights such as identifying significant time scales and
groups of correlated chemical species from our analysis. Our work provides for
the first time a theoretical underpinning for the successful experimental
analysis and points to way to extract mechanisms from large- scale high
throughput data sets.


  Many populations in nature are fragmented: they consist of local populations
occupying separate patches. A local population is prone to extinction due to
the shot noise of birth and death processes. A migrating population from
another patch can dramatically delay the extinction. What is the optimal
migration rate that minimizes the extinction risk of the whole population? Here
we answer this question for a connected network of model habitat patches with
different carrying capacities.


  Resource competition in heterogeneous environments is still an unresolved
problem of theoretical ecology. In this article I analyze competition between
two phytoplankton species in a deep water column, where the distributions of
main resources (light and a limiting nutrient) have opposing gradients and
co-limitation by both resources causes a deep biomass maximum. Assuming that
the species have a trade-off in resource requirements and the water column is
weakly mixed, I apply the invasion threshold analysis (Ryabov and Blasius 2011)
to determine relations between environmental conditions and phytoplankton
composition. Although species deplete resources in the interior of the water
column, the resource levels at the bottom and surface remain high. As a result,
the slope of resources gradients becomes a new crucial factor which, rather
than the local resource values, determines the outcome of competition. The
value of resource gradients nonlinearly depend on the density of consumers.
This leads to complex relationships between environmental parameters and
species composition. In particular, it is shown that an increase of both the
incident light intensity or bottom nutrient concentrations favors the best
light competitors, while an increase of the turbulent mixing or background
turbidity favors the best nutrient competitors. These results might be
important for prediction of species composition in deep ocean.


  Both genomic stability and sustenance of day-to-day life rely on efficient
and accurate readout of the genetic code. Single-molecule experiments show that
transcription and replication are highly stochastic and irregular processes,
with the polymerases frequently pausing and even reversing direction. While
such behavior is recognized as stemming from a sophisticated proofreading
mechanism during replication, the origin and functional significance of
irregular transcription dynamics remain controversial. Here, we theoretically
examine the implications of RNA polymerase backtracking and transcript cleavage
on transcription rates and fidelity. We illustrate how an extended state space
for backtracking provides entropic fidelity enhancements that, together with
additional fidelity checkpoints, can account for physiological error rates. To
explore the competing demands of transcription fidelity, nucleotide
triphosphate (NTP) consumption and transcription speed in a physiologically
relevant setting, we establish an analytically framework for evaluating
transcriptional performance at the level of extended sequences. Using this
framework, we reveal a mechanism by which moderately irregular transcription
results in astronomical gains in the rate at which extended high-fidelity
transcripts can be produced under physiological conditions.


  A simple way to model phenotypic evolution is to assume that after splitting,
the trait values of the sister species diverge as independent Brownian motions.
Relying only on a prior distribution for the underlying species tree
(conditioned on the number, n, of extant species) we study the random vector
(X_1,...,X_n) of the observed trait values. In this paper we derive compact
formulae for the variance of the sample mean and the mean of the sample
variance for the vector (X_1,...,X_n).
  The key ingredient of these formulae is the correlation coefficient between
two trait values randomly chosen from (X_1,...,X_n). This interspecies
correlation coefficient takes into account not only variation due to the random
sampling of two species out of n and the stochastic nature of Brownian motion
but also the uncertainty in the phylogenetic tree. The latter is modeled by a
(supercritical or critical) conditioned branching process. In the critical case
we modify the Aldous-Popovic model by assuming a proper prior for the time of
origin.


  Protein function frequently involves conformational changes with large
amplitude on timescales which are difficult and computationally expensive to
access using molecular dynamics. In this paper, we report on the combination of
three computationally inexpensive simulation methods-normal mode analysis using
the elastic network model, rigidity analysis using the pebble game algorithm,
and geometric simulation of protein motion-to explore conformational change
along normal mode eigenvectors. Using a combination of ELNEMO and FIRST/FRODA
software, large-amplitude motions in proteins with hundreds or thousands of
residues can be rapidly explored within minutes using desktop computing
resources. We apply the method to a representative set of six proteins covering
a range of sizes and structural characteristics and show that the method
identifies specific types of motion in each case and determines their amplitude
limits.


  Known results on the moments of the distribution generated by the two-locus
Wright-Fisher diffusion model and a duality between the diffusion process and
the ancestral process with recombination are briefly summarized. A numerical
methods for computing moments by a Markov chain Monte Carlo and a method to
compute closed-form expressions of the moments are presented. By using the
duality argument properties of the ancestral recombination graph are studied in
terms of the moments.


  Modeling stochasticity in gene regulatory networks is an important and
complex problem in molecular systems biology. To elucidate intrinsic noise,
several modeling strategies such as the Gillespie algorithm have been used
successfully. This paper contributes an approach as an alternative to these
classical settings. Within the discrete paradigm, where genes, proteins, and
other molecular components of gene regulatory networks are modeled as discrete
variables and are assigned as logical rules describing their regulation through
interactions with other components. Stochasticity is modeled at the biological
function level under the assumption that even if the expression levels of the
input nodes of an update rule guarantee activation or degradation there is a
probability that the process will not occur due to stochastic effects. This
approach allows a finer analysis of discrete models and provides a natural
setup for cell population simulations to study cell-to-cell variability. We
applied our methods to two of the most studied regulatory networks, the outcome
of lambda phage infection of bacteria and the p53-mdm2 complex.


  The fact that individuals will most likely behave differently in different
situations begets the introduction of conditional strategies. Inspired by this,
we study the evolution of cooperation in the spatial public goods game, where
besides unconditional cooperators and defectors, also different types of
conditional cooperators compete for space. Conditional cooperators will
contribute to the public good only if other players within the group are likely
to cooperate as well, but will withhold their contribution otherwise. Depending
on the number of other cooperators that are required to elicit cooperation of a
conditional cooperator, the latter can be classified in as many types as there
are players within each group. We find that the most cautious cooperators, such
that require all other players within a group to be conditional cooperators,
are the undisputed victors of the evolutionary process, even at very low
synergy factors. We show that the remarkable promotion of cooperation is due
primarily to the spontaneous emergence of quarantining of defectors, which
become surrounded by conditional cooperators and are forced into isolated
convex "bubbles" from where they are unable to exploit the public good. This
phenomenon can be observed only in structured populations, thus adding to the
relevance of pattern formation for the successful evolution of cooperation.


  We study a model of competition for resource through a chemostat-type model
where species consume the common resource that is constantly supplied. We
assume that the species and resources are characterized by a continuous trait.
As already proved, this model, although more complicated than the usual
Lotka-Volterra direct competition model, describes competitive interactions
leading to concentrated distributions of species in continuous trait space.
Here we assume a very fast dynamics for the supply of the resource and a fast
dynamics for death and uptake rates. In this regime we show that factors that
are independent of the resource competition become as important as the
competition efficiency and that the direct competition model is a good
approximation of the chemostat. Assuming these two timescales allows us to
establish a mathematically rigorous proof showing that our resource-competition
model with continuous traits converges to a direct competition model. We also
show that the two timescales assumption is required to mathematically justify
the corresponding classic result on a model consisting of only finite number of
species and resources (MacArthur, R. Theor. Popul. Biol. 1970:1, 1-11). This is
performed through asymptotic analysis, introducing different scales for the
resource renewal rate and the uptake rate. The mathematical difficulty relies
in a possible initial layer for the resource dynamics. The chemostat model
comes with a global convex Lyapunov functional. We show that the particular
form of the competition kernel derived from the uptake kernel, satisfies a
positivity property which is known to be necessary for the direct competition
model to enjoy the related Lyapunov functional.


  A mechanism of double strand breaking (DSB) in DNA due to the action of two
electrons is considered. These are the electrons produced in the vicinity of
DNA molecules due to ionization of water molecules with a consecutive emission
of two electrons, making such a mechanism possible. This effect qualitatively
solves a puzzle of large yields of DSBs following irradiation of DNA molecules.
The transport of secondary electrons, including the additional electrons, is
studied in relation to the assessment of radiation damage due to incident ions.
This work is a stage in the inclusion of Auger mechanism and like effects into
the multiscale approach to ion-beam cancer therapy.


  This article reformulates a common illness-death model in terms of a new
system of stochastical differential equations (SDEs). The SDEs are used to
estimate epidemiological characteristics and burden of systemic lupus
erythematosus in England and Wales in 1995.


  Some years ago a cellular automata model was proposed to describe the
evolution of the immune repertoire of B cells and antibodies based on Jerne's
immune network theory and shape-space formalism. Here we investigate if the
networks generated by this model in the different regimes can be classified as
complex networks. We have found that in the chaotic regime the network has
random characteristics with large, constant values of clustering coefficients,
while in the ordered phase, the degree distribution of the network is
exponential and the clustering coefficient exhibits power law behavior. In the
transition region we observed a mixed behavior (random-like and exponential) of
the degree distribution as opposed to the scale-free behavior reported for
other biological networks. Randomness and low connectivity in the active sites
allow for rapid changes in the connectivity distribution of the immune network
in order to include and/or discard information and generate a dynamic memory.
However it is the availability of the low concentration nodes to change rapidly
without driving the system to pathological states that allow the generation of
dynamic memory and consequently a reproduction of immune system behavior in
mice. Although the overall behavior of degree correlation is positive, there is
an interplay between assortative and disassortative mixing in the stable and
transition regions regulated by a threshold value of the node degree, which
achieves a maximum value on the transition region and becomes totally
assortative in the chaotic regime.


  The enzyme FoF1-ATP synthase provides the 'chemical energy currency'
adenosine triphosphate (ATP) for living cells. Catalysis is driven by
mechanochemical coupling of subunit rotation within the enzyme with
conformational changes in the three ATP binding sites. Proton translocation
through the membrane-bound Fo part of ATP synthase powers a 10-step rotary
motion of the ring of c subunits. This rotation is transmitted to the gamma and
epsilon subunits of the F1 part. Because gamma and epsilon subunits rotate in
120 deg steps, we aim to unravel this symmetry mismatch by real time monitoring
subunit rotation using single-molecule Forster resonance energy transfer
(FRET). One fluorophore is attached specifically to the F1 motor, another one
to the Fo motor of the liposome-reconstituted enzyme. Photophysical artifacts
due to spectral fluctuations of the single fluorophores are minimized by a
previously developed duty cycle-optimized alternating laser excitation scheme
(DCO-ALEX). We report the detection of reversible elastic deformations between
the rotor parts of Fo and F1 and estimate the maximum angular displacement
during the load-free rotation using Monte Carlo simulations


  FoF1-ATP synthases in Escherichia coli (E. coli) bacteria are membrane-bound
enzymes which use an internal proton-driven rotary double motor to catalyze the
synthesis of adenosine triphosphate (ATP). According to the 'chemiosmotic
hypothesis', a series of proton pumps generate the necessary pH difference plus
an electric potential across the bacterial plasma membrane. These proton pumps
are redox-coupled membrane enzymes which are possibly organized in
supercomplexes, as shown for the related enzymes in the mitochondrial inner
membrane. We report diffusion measurements of single fluorescent FoF1-ATP
synthases in living E. coli by localization microscopy and single enzyme
tracking to distinguish a monomeric enzyme from a supercomplex-associated form
in the bacterial membrane. For quantitative mean square displacement (MSD)
analysis, the limited size of the observation area in the membrane with a
significant membrane curvature had to be considered. The E. coli cells had a
diameter of about 500 nm and a length of about 2 to 3 \mum. Because the surface
coordinate system yielded different localization precision, we applied a
sliding observation window approach to obtain the diffusion coefficient D =
0.072 \mum2/s of FoF1-ATP synthase in living E. coli cells.


  The relationship of base pair openings to DNA flexibility is examined.
Published experimental data on the temperature dependence of the persistence
length by two different groups are well described in terms of an inhomogeneous
Kratky-Porot model with soft and hard joints, corresponding to open and closed
base pairs, and sequence-dependent statistical information about the state of
each pair provided by a Peyrard-Bishop-Dauxois (PBD) model calculation with no
freely adjustable parameters.


  A strain of Halomonas bacteria, GFAJ-1, has been reported to be able to use
arsenate as a nutrient when phosphate is limiting, and to specifically
incorporate arsenic into its DNA in place of phosphorus. However, we have found
that arsenate does not contribute to growth of GFAJ-1 when phosphate is
limiting and that DNA purified from cells grown with limiting phosphate and
abundant arsenate does not exhibit the spontaneous hydrolysis expected of
arsenate ester bonds. Furthermore, mass spectrometry showed that this DNA
contains only trace amounts of free arsenate and no detectable covalently bound
arsenate.


  How self-incompatibility systems are maintained in plant populations is still
a debated issue. Theoretical models predict that self-incompatibility systems
break down according to the intensity of inbreeding depression and number of
S-alleles. Other studies have explored the function of asexual reproduction in
the maintenance of self-incompatibility. However, the population genetics of
partially asexual, self-incompatible populations are poorly understood and
previous studies have failed to consider all possible effects of asexual
reproduction or could only speculate on those effects. In this study, we
investigated how partial asexuality may affect genetic diversity at the S-locus
and fitness in small self-incompatible populations. A genetic model including
an S-locus and a viability locus was developed to perform forward simulations
of the evolution of populations of various sizes. Drift combined with partial
asexuality produced a decrease in the number of alleles at the S-locus. In
addition, an excess of heterozygotes was present in the population, causing an
increase in mutation load. This heterozygote excess was enhanced by the
self-incompatibility system in small populations. In addition, in highly
asexual populations, individuals produced asexually had some fitness advantages
over individuals produced sexually, because sexual reproduction produces
homozygotes of the deleterious allele, contrary to asexual reproduction. Our
results suggest that future research on the function of asexuality for the
maintenance of self-incompatibility will need to (1) account for whole-genome
fitness (mutation load generated by asexuality, self-incompatibility and drift)
and (2) acknowledge that the maintenance of self-incompatibility may not be
independent of the maintenance of sex itself.


  Knowing which individuals can be more efficient in spreading a pathogen
throughout a determinate environment is a fundamental question in disease
control. Indeed, over the last years the spread of epidemic diseases and its
relationship with the topology of the involved system have been a recurrent
topic in complex network theory, taking into account both network models and
real-world data. In this paper we explore possible correlations between the
heterogeneous spread of an epidemic disease governed by the
susceptible-infected-recovered (SIR) model, and several attributes of the
originating vertices, considering Erd\"os-R\'enyi (ER), Barab\'asi-Albert (BA)
and random geometric graphs (RGG), as well as a real case of study, the US Air
Transportation Network that comprises the US 500 busiest airports along with
inter-connections. Initially, the heterogeneity of the spreading is achieved
considering the RGG networks, in which we analytically derive an expression for
the distribution of the spreading rates among the established contacts, by
assuming that such rates decay exponentially with the distance that separates
the individuals. Such distribution is also considered for the ER and BA models,
where we observe topological effects on the correlations. In the case of the
airport network, the spreading rates are empirically defined, assumed to be
directly proportional to the seat availability. Among both the theoretical and
the real networks considered, we observe a high correlation between the total
epidemic prevalence and the degree, as well as the strength and the
accessibility of the epidemic sources. For attributes such as the betweenness
centrality and the $k$-shell index, however, the correlation depends on the
topology considered.


  We study the relative importance of "top-speed" (long-term growth rate) and
"acceleration" (how quickly the long-term growth rate can be reached) in the
evolutionary race to increase population size. We observe that fitness alone
does not capture growth rate: robustness, a property of neutral network shape,
combines with fitness to include the effect of deleterious mutations, giving
growth rate. Similarly, we show that growth rate alone does not capture
population size: regularity, a different property of neutral network shape,
combines with growth rate to include the effect of higher depletion rates early
on, giving size. Whereas robustness is a function of the principal eigenvalue
of the neutral network adjacency matrix, regularity is a function of the
principal eigenvector. We show that robustness is not correlated with
regularity, and observe in silico the selection for regularity by evolving RNA
ribozymes. Despite having smaller growth rates, the more regular ribozymes have
the biggest populations.


  Dengue viral infections show unique infection patterns arising from its four
serot- ypes, (DENV-1,2,3,4). Its effects range from simple fever in primary
infections to potentially fatal secondary infections. We analytically and
numerically analyse virus dynamics and humoral response in a host during
primary and secondary dengue infection for long periods using micro-epidemic
models. The models presented here incorporate time delays, antibody dependent
enhancement (ADE), a dynamic switch and a correlation factor between different
DENV serotypes. We find that the viral load goes down to undetectable levels
within 7-14 days as is observed for dengue infection, in both cases. For
primary infection, the stability analysis of steady states shows interesting
dependence on the time delay involved in the production of antibodies from
plasma cells. We demonstrate the existence of a critical value for the immune
response parameter, beyond which the infection gets completely cured. For
secondary infections with a different serotype, the homologous antibody
production is enhanced due to the influence of heterologous antibodies. The
antibody production is also controlled by the correlation factor, which is a
measure of similarities between the different DENV serotypes involved. Our
results agree with clinically observed humoral responses for primary and
secondary infections.


  Recent experimental evidence [Gruji\'c et al., PLoS ONE 5, e13749 (2010)] on
the spatial Prisoner's Dilemma suggests that players choosing to cooperate or
not on the basis of their previous action and the actions of their neighbors
coexist with steady defectors and cooperators. We here study the coexistence of
these three strategies in the multiplayer iterated Prisoner's Dilemma by means
of the replicator dynamics. We consider groups with n = 2, 3, 4 and 5 players
and compute the payoffs to every type of player as the limit of a Markov chain
where the transition probabilities between actions are found from the
corresponding strategies. We show that for group sizes up to n = 4 there exists
an interior point in which the three strategies coexist, the corresponding
basin of attraction decreasing with increasing number of players, whereas we
have not been able to locate such a point for n = 5. We analytically show that
in the infinite n limit no interior points can arise. We conclude by discussing
the implications of this theoretical approach on the behavior observed in
experiments.


  We investigate the problem of signal transduction via a descriptive analysis
of the spatial organization of the complement of proteins exerting a certain
function within a cellular compartment. We propose a scheme to assign a
numerical value to individual proteins in a protein interaction network by
means of a simple optimization algorithm. We test our procedure against
datasets focusing on the proteomes in the neurite and soma compartments.


  DNA replication is an essential process in biology and its timing must be
robust so that cells can divide properly. Random fluctuations in the formation
of replication starting points, called origins, and the subsequent activation
of proteins lead to variations in the replication time. We analyse these
stochastic properties of DNA and derive the positions of origins corresponding
to the minimum replication time. We show that under some conditions the
minimization of replication time leads to the grouping of origins, and relate
this to experimental data in a number of species showing origin grouping.


  The Moran process models the spread of genetic mutations through a
population. A mutant with relative fitness $r$ is introduced into a population
and the system evolves, either reaching fixation (in which every individual is
a mutant) or extinction (in which none is). In a widely cited paper (Nature,
2005), Lieberman, Hauert and Nowak generalize the model to populations on the
vertices of graphs. They describe a class of graphs (called "superstars"), with
a parameter $k$. Superstars are designed to have an increasing fixation
probability as $k$ increases. They state that the probability of fixation tends
to $1-r^{-k}$ as graphs get larger but we show that this claim is untrue as
stated. Specifically, for $k=5$, we show that the true fixation probability (in
the limit, as graphs get larger) is at most $1-1/j(r)$ where
$j(r)=\Theta(r^4)$, contrary to the claimed result. We do believe that the
qualitative claim of Lieberman et al.\ --- that the fixation probability of
superstars tends to 1 as $k$ increases --- is correct, and that it can probably
be proved along the lines of their sketch. We were able to run larger computer
simulations than the ones presented in their paper. However, simulations on
graphs of around 40,000 vertices do not support their claim. Perhaps these
graphs are too small to exhibit the limiting behaviour.


  GU "wobble" pair may contain either the enol form of uracil and the keto form
of guanine or the enol form of guanine and the keto form of uracil. A direct
proof of uracil existence in enol form can be obtained by using NMR
spectroscopy. In HNN-COSY experiments the GU pair imitates the GC pair due to
the enolization of uracil.


  Motivation: Cell-biological processes are regulated through a complex network
of interactions between genes and their products. The processes, their
activating conditions, and the associated transcriptional responses are often
unknown. Organism-wide modeling of network activation can reveal unique and
shared mechanisms between physiological conditions, and potentially as yet
unknown processes. We introduce a novel approach for organism-wide discovery
and analysis of transcriptional responses in interaction networks. The method
searches for local, connected regions in a network that exhibit coordinated
transcriptional response in a subset of conditions. Known interactions between
genes are used to limit the search space and to guide the analysis. Validation
on a human pathway network reveals physiologically coherent responses,
functional relatedness between physiological conditions, and coordinated,
context-specific regulation of the genes. Availability: Implementation is
freely available in R and Matlab at http://netpro.r-forge.r-project.org


  Systemic properties of living cells are the result of molecular dynamics
governed by so-called genetic regulatory networks (GRN). These networks capture
all possible features of cells and are responsible for the immense levels of
adaptation characteristic to living systems. At any point in time only small
subsets of these networks are active. Any active subset of the GRN leads to the
expression of particular sets of molecules (expression modes). The subsets of
active networks change over time, leading to the observed complex dynamics of
expression patterns. Understanding of this dynamics becomes increasingly
important in systems biology and medicine. While the importance of
transcription rates and catalytic interactions has been widely recognized in
modeling genetic regulatory systems, the understanding of the role of
degradation of biochemical agents (mRNA, protein) in regulatory dynamics
remains limited. Recent experimental data suggests that there exists a
functional relation between mRNA and protein decay rates and expression modes.
In this paper we propose a model for the dynamics of successions of sequences
of active subnetworks of the GRN. The model is able to reproduce key
characteristics of molecular dynamics, including homeostasis, multi-stability,
periodic dynamics, alternating activity, differentiability, and self-organized
critical dynamics. Moreover the model allows to naturally understand the
mechanism behind the relation between decay rates and expression modes. The
model explains recent experimental observations that decay-rates (or turnovers)
vary between differentiated tissue-classes at a general systemic level and
highlights the role of intracellular decay rate control mechanisms in cell
differentiation.


  Theoretical studies of protein folding on lattice models relie on the
assumption that water close to amino-acids is always in thermal equilibrium all
along the folding pathway. Within this framework, it has always been considered
that out-of-equilibrium properties, such as folding time, could be evaluated
equivalently from an averaging over a collection of trajectories of the protein
with water described eitherexplicitly or through a mean-field approach. To
critically assess this hypothesis, we built a two-dimensional lattice model of
a protein in interaction with water molecules that can adopt a wide range of
conformations. This microscopic description of the solvent has been used
further to derive an effective model by averaging over all the degrees of
freedom of the solvent. At thermal equilbrium, the two descriptions are
rigourously equivalent, predicting the same folded conformation of the protein.
The model allows exact calculations of some relaxation properties using the
master equations associated to both solvent descriptions. The kinetic patterns
associated to the folding pathways are remarkably different. In this work we
demonstrate, that an effective description of the solvent can not described
properly the folding pathway of a protein. The microscopic solvent model, that
describes correctly the microscopic routes, appears to be the only candidate to
study folding kinetics.


  Thermophilic enzymes can operate at higher temperatures but show reduced
activities at room temperature. They are in general more stable during
preparation and, accordingly, are considered to be more rigid in structure.
Crystallization is often easier compared to proteins from bacteria growing at
ambient temperatures, especially for membrane proteins. The ATP-producing
enzyme FoF1-ATP synthase from thermoalkaliphilic Caldalkalibacillus thermarum
strain TA2.A1 is driven by a Fo motor consisting of a ring of 13 c-subunits. We
applied a single-molecule F\"orster resonance energy transfer (FRET) approach
using duty cycle-optimized alternating laser excitation (DCO-ALEX) to monitor
the expected 13-stepped rotary Fo motor at work. New FRET transition histograms
were developed to identify the smaller step sizes compared to the 10-stepped Fo
motor of the Escherichia coli enzyme. Dwell time analysis revealed the
temperature and the LDAO dependence of the Fo motor activity on the single
molecule level. Back-and-forth stepping of the Fo motor occurs fast indicating
a high flexibility in the membrane part of this thermophilic enzyme.


  Massively parallel DNA sequencing techniques are rapidly changing the
dynamics of phylogenetic study design by exponentially increasing the discovery
of phylogenetically useful loci. This increase in the number of phylogenetic
markers potentially provides researchers the opportunity to select subsets of
loci best-addressing particular phylogenetic hypotheses based on objective
measures of performance over different time scales. Investigators may also want
to determine the power of particular phylogenetic markers relative to each
other. However, currently available tools are designed to evaluate a small
number of markers and are not well-suited to screening hundreds or thousands of
candidate loci across the genome. TAPIR is an alternative implementation of
Townsend's estimate of phylogenetic informativeness (PI) that enables rapid
estimation and summary of PI when applied to data sets containing hundreds to
thousands of candidate, phylogenetically informative loci.


  Several indices that measure the degree of balance of a rooted phylogenetic
tree have been proposed so far in the literature. In this work we define and
study a new index of this kind, which we call the total cophenetic index: the
sum, over all pairs of different leaves, of the depth of their least common
ancestor. This index makes sense for arbitrary trees, can be computed in linear
time and it has a larger range of values and a greater resolution power than
other indices like Colless' or Sackin's. We compute its maximum and minimum
values for arbitrary and binary trees, as well as exact formulas for its
expected value for binary trees under the Yule and the uniform models of
evolution. As a byproduct of this study, we obtain an exact formula for the
expected value of the Sackin index under the uniform model, a result that seems
to be new in the literature.


  MicroRNAs can affect the protein translation using nine mechanistically
different mechanisms, including repression of initiation and degradation of the
transcript. There is a hot debate in the current literature about which
mechanism and in which situations has a dominant role in living cells. The
worst, same experimental systems dealing with the same pairs of mRNA and miRNA
can provide ambiguous evidences about which is the actual mechanism of
translation repression observed in the experiment. We start with reviewing the
current knowledge of various mechanisms of miRNA action and suggest that
mathematical modeling can help resolving some of the controversial
interpretations. We describe three simple mathematical models of miRNA
translation that can be used as tools in interpreting the experimental data on
the dynamics of protein synthesis. The most complex model developed by us
includes all known mechanisms of miRNA action. It allowed us to study possible
dynamical patterns corresponding to different miRNA-mediated mechanisms of
translation repression and to suggest concrete recipes on determining the
dominant mechanism of miRNA action in the form of kinetic signatures. Using
computational experiments and systematizing existing evidences from the
literature, we justify a hypothesis about co-existence of distinct
miRNA-mediated mechanisms of translation repression. The actually observed
mechanism will be that acting on or changing the limiting "place" of the
translation process. The limiting place can vary from one experimental setting
to another. This model explains the majority of existing controversies
reported.


  In an experimental study of single enzyme reactions, it has been proposed
that the rate constants of the enzymatic reactions fluctuate randomly,
according to a given distribution. To quantify the uncertainty arising from
random rate constants, it is necessary to investigate how one can simulate such
a biochemical system. To do this, we will take the Gillespie's stochastic
simulation algorithm for simulating the evolution of the state of a chemical
system, and study a modification of the algorithm that incorporates the random
rate constants. In addition to simulating the waiting time of each reaction
step, the modified algorithm also involves simulating the random fluctuation of
the rate constant at each reaction time. We consider the modified algorithm in
a general framework, then specialize it to two contrasting physical models, one
in which the fluctuations occur on a much faster time scale than the reaction
step, and the other in which the fluctuations occur much more slowly. The
latter case was applied to the single enzyme reaction system, using in part the
Metropolis-Hastings algorithm to enact the given distribution on the random
rate constants. The modified algorithm is shown to produce simulation outputs
that are corroborated by the experimental results. It is hoped that this
modified algorithm can subsequently be used as a tool for the estimation or
calibration of parameters in the system using experimental data.


  We developed mathematical models that simulate community dynamics derived
from a series of perturbation experiments. These experiments were performed in
an Atlantic Forest stream. The three trophic level community was submitted to
two combinations of press perturbation experiment. In the first, the top
trophic level was removed by electrical exclusion. In the second configuration,
the top trophic level was still excluded plus a group of species from the
second trophic level was inhibited, also by electrical pulses. Experiments were
repeated at different times to increase confidence in the observed mechanisms.
Community responses indicated the existence of cascading interactions similar
to a classic trophic cascade. The community was composed of Macrobrachium
shrimps at the top trophic level, Ephemeroptera and Chironomidae larval insects
at the second trophic level and periphyton (= epilithon = biofilm) at the first
trophic level. The shrimp exclusion caused an increase in ephemeropterans and
chironomids and led the periphyton to decrease in the exclusion zone. When
ephemeropterans were inhibited together with the exclusion of shrimps, the
chironomids and the periphyton increased. We modelled the community
interactions by means of differential equation systems, simulating the three
configurations: 1. natural condition, 2. shrimp exclusion condition and 3.
shrimp-exclusion-plus-ephemeropteran-inhibition condition. All parameters were
calculated via experimental results and the stability of the models was
determined by their matrix eigenvalues. The models successfully reproduced the
qualitative responses of the community. They proved to have attraction points,
which enables us to predict that the system has some ability to return to its
previous condition after local perturbation.


  The fundamental law for protein folding is the Thermodynamic Principle: the
amino acid sequence of a protein determines its native structure and the native
structure has the minimum Gibbs free energy. If all chemical problems can be
answered by quantum mechanics, there should be a quantum mechanics derivation
of Gibbs free energy formula G(X) for every possible conformation X of the
protein. We apply quantum statistics to derive such a formula. For simplicity,
only monomeric self folding globular proteins are covered. We point out some
immediate applications of the formula. We show that the formula explains the
observed phenomena very well. It gives a unified explanation to both folding
and denaturation; it explains why hydrophobic effect is the driving force of
protein folding and clarifies the role played by hydrogen bonding; it explains
the successes and deficients of various surface area models. The formula also
gives a clear kinetic force of the folding: Fi(X) = - \nablaxi G(X). This also
gives a natural way to perform the ab initio prediction of protein structure,
minimizing G(X) by Newton's fastest desciending method.


  In biochemical signaling, information is often encoded in oscillatory
signals. However, the advantages of such a coding strategy over an amplitude
encoding scheme of constant signals remain unclear. Here we study the dynamics
of a simple model gene promoter in response to oscillating and constant
transcription factor signals. We find that in biologically-relevant parameter
regimes an oscillating input can produce a more constant protein level than a
constant input. Our results suggest that oscillating signals may be used to
minimize noise in gene regulation.


  Electropermeabilization is a method by which molecules can be introduced in
living cells through the plasma membrane. Here, we show that
electropermeabilization affects the lateral mobility of Rae-1, a GPi anchored
protein. Our results shows that 10 to 20 % of the membrane surface is occupied
by defaults or pores. Gene electrotransfert also affects the lateral mobility
of Rae-1. Furthermore, we clearly show that once inserted into the plasma
membrane, DNA totaly excludes Rae-1, indicating that DNA molecules are tightly
packed together to forms huge aggregates.


  In this paper we summarize the main parts of the first exit time theory
developed in connection to the life table data and the resulting theoretical
and applied issues. Several new tools arise from the development of this theory
and especially the Health State Function and some important characteristics of
this function. Special attention has being done in the presentation of the
health state function along with the well established theory for the Force of
Mortality and the related applications as are the life tables and the
estimation of life expectancies. A main part of this work is the formulation of
the appropriate non-linear analysis program including a model which provides an
almost perfect fit to life table data. This model, proposed in 1995 is now
expanded as to include the mortality excess for the age group from 15-30 years.
A version of the program is given in Excel and provided at the website:
http://www.cmsim.net


  We study the classic Susceptible-Infected-Recovered (SIR) model for the
spread of an infectious disease. In this stochastic process, there are two
competing mechanism: infection and recovery. Susceptible individuals may
contract the disease from infected individuals, while infected ones recover
from the disease at a constant rate and are never infected again. Our focus is
the behavior at the epidemic threshold where the rates of the infection and
recovery processes balance. In the infinite population limit, we establish
analytically scaling rules for the time-dependent distribution functions that
characterize the sizes of the infected and the recovered sub-populations. Using
heuristic arguments, we also obtain scaling laws for the size and duration of
the epidemic outbreaks as a function of the total population. We perform
numerical simulations to verify the scaling predictions and discuss the
consequences of these scaling laws for near-threshold epidemic outbreaks.


  Tree rearrangement operations typically induce a metric on the space of
phylogenetic trees. One important property of these metrics is the size of the
neighbourhood, that is, the number of trees exactly one operation from a given
tree. We present an expression for the size of the TBR (tree bisection and
reconnection) neighbourhood, thus answering a question first posed in [Annals
of Combinatorics, 5, 2001 1-15].


  Starfish typically show pentameral symmetry, and they are typically similar
in shape to a pentagram. Although starfish can evolve and live with other
numbers of arms, the dominant species always show pentameral symmetry. We used
mathematical and physical methods to analyze the superiority of starfish with
five arms in comparison with those with a different number of arms with respect
to detection, turning over, autotomy and adherence. In this study, we
determined that starfish with five arms, although slightly inferior to others
in one or two aspects, exhibit the best performance when the four
aforementioned factors are considered together. In addition, five-armed
starfish perform best on autotomy, which is crucially important for starfish
survival. This superiority contributes to the dominance of five-armed starfish
in evolution, which is consistent with the practical situation. Nevertheless,
we can see some flexibility in the number and conformation of arms. The
analyses performed in our research will be of great help in unraveling the
mysteries of dominant shapes and structures.


  We present a differential rate equation model of chiral polymerization based
on a simple copolymerization scheme in which the enantiomers are added to, or
removed from, the homochiral or heterochiral chains (reversible stepwise
isodesmic growth or dissociation). The model is set up for closed systems and
takes into account the corresponding thermodynamic constraints implied by the
reversible monomer attachments, while obeying a constant mass constraint. In
its simplest form, the model depends on a single variable rate constant, the
maximum chain length N, and the initial concentrations. We have fit the model
to the experimental data from the Rehovot group on lattice-controlled chiral
amplification of oligopeptides. We find in all the chemical systems employed
except for one, that the model fits the measured relative abundances of the
oligopetides with higher degrees of correlation than from a purely random
polymerization process.


  Phylogenetic mixture models, in which the sites in sequences undergo
different substitution processes along the same or different trees, allow the
description of heterogeneous evolutionary processes. As data sets consisting of
longer sequences become available, it is important to understand such models,
for both theoretical insights and use in statistical analyses. Some recent
articles have highlighted disturbing "mimicking" behavior in which a
distribution from a mixture model is identical to one arising on a different
tree or trees. Other works have indicated such problems are unlikely to occur
in practice, as they require very special parameter choices.
  After surveying some of these works on mixture models, we give several new
results. In general, if the number of components in a generating mixture is not
too large and we disallow zero or infinite branch lengths, then it cannot mimic
the behavior of a non-mixture on a different tree. On the other hand, if the
mixture model is locally over-parameterized, it is possible for a phylogenetic
mixture model to mimic distributions of another tree model. Though theoretical
questions remain, these sorts of results can serve as a guide to when the use
of mixture models in either ML or Bayesian frameworks is likely to lead to
statistically consistent inference, and when mimicking due to heterogeneity
should be considered a realistic possibility.


  We present a unified mathematical approach to epidemiological models with
parametric heterogeneity, i.e., to the models that describe individuals in the
population as having specific parameter (trait) values that vary from one
individuals to another. This is a natural framework to model, e.g.,
heterogeneity in susceptibility or infectivity of individuals. We review, along
with the necessary theory, the results obtained using the discussed approach.
In particular, we formulate and analyze an SIR model with distributed
susceptibility and infectivity, showing that the epidemiological models for
closed populations are well suited to the suggested framework. A number of
known results from the literature is derived, including the final epidemic size
equation for an SIR model with distributed susceptibility. It is proved that
the bottom up approach of the theory of heterogeneous populations with
parametric heterogeneity allows to infer the population level description,
which was previously used without a firm mechanistic basis; in particular, the
power law transmission function is shown to be a consequence of the initial
gamma distributed susceptibility and infectivity. We discuss how the general
theory can be applied to the modeling goals to include the heterogeneous
contact population structure and provide analysis of an SI model with
heterogeneous contacts. We conclude with a number of open questions and
promising directions, where the theory of heterogeneous populations can lead to
important simplifications and generalizations.


  Motivated by questions in biological classification, we discuss some
elementary combinatorial and computational properties of certain set systems
that generalize hierarchies, namely, 'patchworks', 'weak patchworks', 'ample
patchworks' and 'saturated patchworks' and also outline how these concepts
relate to an apparently new 'duality theory' for cluster systems that is based
on the fundamental concept of 'compatibility' of clusters.


  The innate immune system, acting as the first line of host defense, senses
and adapts to foreign challenges through complex intracellular and
intercellular signaling networks. Endotoxin tolerance and priming elicited by
macrophages are classic examples of the complex adaptation of innate immune
cells. Upon repetitive exposures to different doses of bacterial endotoxin
(lipopolysaccharide) or other stimulants, macrophages show either suppressed or
augmented inflammatory responses compared to a single exposure to the
stimulant. Endotoxin tolerance and priming are critically involved in both
immune homeostasis and the pathogenesis of diverse inflammatory diseases.
However, the underlying molecular mechanisms are not well understood. By means
of a computational search through the parameter space of a coarse-grained
three-node network with a two-stage Metropolis sampling approach, we enumerated
all the network topologies that can generate priming or tolerance. We
discovered three major mechanisms for priming (pathway synergy, suppressor
deactivation, activator induction) and one for tolerance (inhibitor
persistence). These results not only explain existing experimental
observations, but also reveal intriguing test scenarios for future experimental
studies to clarify mechanisms of endotoxin priming and tolerance.


  This paper presents a novel method to segment/decode DNA sequences based on
n-grams statistical language model. Firstly, we find the length of most DNA
'words' is 12 to 15 bps by analyzing the genomes of 12 model species. Then we
design an unsupervised probability based approach to segment the DNA sequences.
The benchmark of segmenting method is also proposed.


  To elucidate the possible role of structural waters in stabilizing bent DNA,
various conformations of AT-containing decamers, (A5T5)2 and A10:T10, were
studied by Monte Carlo simulations. The duplexes were constrained to reproduce
the NMR inter-proton distances for the A-tracts at two temperatures: 5 and 35C.
Analysis of the water shell structures revealed a strong correlation between
the groove widths on the one hand, and the types of hydration patterns and
their probabilities on the other hand. Depending on the minor groove width, the
following patterns were observed in this groove: either an interstrand
"hydration spine", or interstrand two-water bridges, or a double ribbon of
intrastrand sugar-base "water strings". Hydration shell in the major groove is
less regular than in the minor groove, which agrees with crystallographic data.
As in previous studies, energetically advantageous hydration is found for the
A-tract conformations with narrow minor groove and wide major groove (B'-like
DNA), known to be correlated with DNA bending and curvature. In addition, our
calculations indicate that the advantage of such DNA conformations is coupled
with the increase in adenine N7 hydration in the major groove. As the major
groove is opened wider, its hydration shell is enriched by energetically
favorable "trident waters" hydrogen bonded to three hydrophilic centers
belonging to two adjacent adenines: two N7 atoms and one amino-proton H(N6).
Based on these results, we suggest that formation of the novel hydration
pattern in the major groove is one of the factors responsible for stabilization
of the B'-like conformations of A-tracts at low temperature in the absense of
dehydrating agents, leading in turn to the strong intrinsic curvature of DNA
containing alternating A-tracts and "mixed" GC-rich sequences.


  Rock-paper-scissors (RPS) game is a nice model to study the biodiversity in
ecosystem. However, the previous studies only consider the nearest- neighbor-
interaction among the species. In this paper, taking the long range migration
into account, the effects of the interplay between nearest-neighbor-interaction
and long-range-interaction of Levy flight obey the power law distance
distribution with the exponent h (-0.3<h<-0.1) in spatial RPS game is
investigated. Taking the probability of long range Levy flight and the power
exponent as parameters, the coexistence conditions of three species are found.
The critical curves for stable coexistence of three species in the parameters
space are presented. It is also found that long-range-interaction with Levy
flight has interesting effects on the final spatiotemporal pattern of the
system. The results reveal that the long-range-interaction of Levy flight
exhibit pronounced effects on biodiversity of ecosystem.


  The cellular phenotype is described by a complex network of molecular
interactions. Elucidating network properties that distinguish disease from the
healthy cellular state is therefore of critical importance for gaining
systems-level insights into disease mechanisms and ultimately for developing
improved therapies. By integrating gene expression data with a protein
interaction network to induce a stochastic dynamics on the network, we here
demonstrate that cancer cells are characterised by an increase in the dynamic
network entropy, compared to cells of normal physiology. Using a fundamental
relation between the macroscopic resilience of a dynamical system and the
uncertainty (entropy) in the underlying microscopic processes, we argue that
cancer cells will be more robust to random gene perturbations. In addition, we
formally demonstrate that gene expression differences between normal and cancer
tissue are anticorrelated with local dynamic entropy changes, thus providing a
systemic link between gene expression changes at the nodes and their local
network dynamics. In particular, we also find that genes which drive
cell-proliferation in cancer cells and which often encode oncogenes are
associated with reductions in the dynamic network entropy. In summary, our
results support the view that the observed increased robustness of cancer cells
to perturbation and therapy may be due to an increase in the dynamic network
entropy that allows cells to adapt to the new cellular stresses. Conversely,
genes that exhibit local flux entropy decreases in cancer may render cancer
cells more susceptible to targeted intervention and may therefore represent
promising drug targets.


  We study the dynamics of quantum excitations inside macromolecules which can
undergo conformational transitions. In the first part of the paper, we use the
path integral formalism to rigorously derive a set of coupled equations of
motion which simultaneously describe the molecular and quantum transport
dynamics, and obey the fluctuation/dissipation relationship. We also introduce
an algorithm which yields the most probable molecular and quantum transport
pathways in rare, thermally-activated reactions. In the second part of the
paper, we apply this formalism to simulate the propagation of a charge during
the collapse of a polymer from an initial stretched conformation to a final
globular state. We find that the charge dynamics is quenched when the chain
reaches a molten globule state. Using random matrix theory we show that this
transition is due to an increase of quantum localization driven by dynamical
disorder.


  Forces in the actin cytoskeleton are generated by small groups of
non-processive myosin II motors for which stochastic effects are highly
relevant. Using a crossbridge model with the assumptions of fast powerstroke
kinetics and equal load sharing between equivalent states, we derive a one-step
master equation for the activity of a finite-sized ensemble of mechanically
coupled myosin II motors. For constant external load, this approach yields
analytical results for duty ratio and force-velocity relation as a function of
ensemble size. We find that stochastic effects cannot be neglected for ensemble
sizes below 15. The one-step master equation can be used also for efficient
computer simulations with linear elastic external load and reveals the sequence
of build-up of force and ensemble rupture that is characteristic for
reconstituted actomyosin contractility.


  We study the influence of complex graphs on the metastability and fixation
properties of a set of evolutionary processes. In the framework of evolutionary
game theory, where the fitness and selection are frequency-dependent and vary
with the population composition, we analyze the dynamics of snowdrift games
(characterized by a metastable coexistence state) on scale-free networks. Using
an effective diffusion theory in the weak selection limit, we demonstrate how
the scale-free structure affects the system's metastable state and leads to
anomalous fixation. In particular, we analytically and numerically show that
the probability and mean time of fixation are characterized by stretched
exponential behaviors with exponents depending on the network's degree
distribution.


  While microscopic organisms can use gradient-based search to locate
resources, this strategy can be poorly suited to the sensory signals available
to macroscopic organisms. We propose a framework that models search-decision
making in cases where sensory signals are infrequent, subject to large
fluctuations, and contain little directional information. Our approach
simultaneously models an organism's intrinsic movement behavior (e.g. Levy
walk) while allowing this behavior to be adjusted based on sensory data. We
find that including even a simple model for signal response can dominate other
features of random search and greatly improve search performance. In
particular, we show that a lack of signal is not a lack of information.
Searchers that receive no signal can quickly abandon target-poor regions. Such
phenomena naturally give rise to the area-restricted search behavior exhibited
by many searching organisms.


  The topological analysis of biological networks has been a prolific topic in
network science during the last decade. A persistent problem with this approach
is the inherent uncertainty and noisy nature of the data. One of the cases in
which this situation is more marked is that of transcriptional regulatory
networks (TRNs) in bacteria. The datasets are incomplete because regulatory
pathways associated to a relevant fraction of bacterial genes remain unknown.
Furthermore, direction, strengths and signs of the links are sometimes unknown
or simply overlooked. Finally, the experimental approaches to infer the
regulations are highly heterogeneous, in a way that induces the appearance of
systematic experimental-topological correlations. And yet, the quality of the
available data increases constantly. In this work we capitalize on these
advances to point out the influence of data (in)completeness and quality on
some classical results on topological analysis of TRNs, specially regarding
modularity at different levels. In doing so, we identify the most relevant
factors affecting the validity of previous findings, highlighting important
caveats to future prokaryotic TRNs topological analysis.


  We present determinant criteria for the preclusion of non-degenerate multiple
steady states in networks of interacting species. A network is modeled as a
system of ordinary differential equations in which the form of the species
formation rate function is restricted by the reactions of the network and how
the species influence each reaction. We characterize families of so-called
power-law kinetics for which the associated species formation rate function is
injective within each stoichiometric class and thus the network cannot exhibit
multistationarity. The criterion for power-law kinetics is derived from the
determinant of the Jacobian of the species formation rate function. Using this
characterization we further derive similar determinant criteria applicable to
general sets of kinetics. The criteria are conceptually simple, computationally
tractable and easily implemented. Our approach embraces and extends previous
work on multistationarity, such as work in relation to chemical reaction
networks with dynamics defined by mass-action or non-catalytic kinetics, and
also work based on graphical analysis of the interaction graph associated to
the system. Further, we interpret the criteria in terms of circuits in the
so-called DSR-graph


  The ability of genetic isolation to block gene flow plays a key role in the
speciation of sexually reproducing organisms. This paper analyses the hybrid
zone dynamics affected by "weak" Haldane's rule, namely the incomplete hybrids
inferiority (sterility/inviability) against the heterogametic (XY or ZW) sex
caused by a Dobzhansky-Muller incompatibility. Different strengths of
incompatibility, dispersal and density-dependent regulation are considered; and
the gene flow and clinal structures of allele frequencies in the presence of
short-range dispersal (the stepping-stone model) are examined. I show that a
weak heterogametic hybrid incompatibility could constitute a substantial
barrier that could reduce gene flow and result in non-coincident and discordant
clines of alleles. It is found that the differential gene flow is more
pronounced under a stronger density-dependent regulation. This study provides a
mechanistic explanation for how an adaptive mutation, which may only have a
marginal fitness effect, could set a gene up as an evolutionary hot-spot.


  It is now well accepted that cellular responses to materials in a biological
medium reflect greatly the adsorbed biomolecular layer, rather than the
material itself. Here, we study by molecular dynamic simulations the
competitive protein adsorption on a surface (Vroman effect), i.e. the
non-monotonic behavior of the amount of protein adsorbed on a surface in
contact with plasma as a function of contact time and plasma concentration. We
find a complex behavior, with regimes during which small and large proteins are
not necessarily competing between them, but are both competing with others in
solution. We show how the effect can be understood, controlled and inverted.


  We consider a nonlinear structured population model with a distributed
recruitment term. The question of the existence of non-trivial steady states
can be treated (at least!) in three different ways. One approach is to study
spectral properties of a parametrized family of unbounded operators. The
alternative approach, on which we focus here, is based on the reformulation of
the problem as an integral equation. In this context we introduce a density
dependent net reproduction rate and discuss its relationship to a biologically
meaningful quantity. Finally, we briefly discuss a third approach, which is
based on the finite rank approximation of the recruitment operator.


  Modeling sequence evolution on phylogenetic trees is a useful technique in
computational biology. Especially powerful are models which take account of the
heterogeneous nature of sequence evolution according to the "grammar" of the
encoded gene features. However, beyond a modest level of model complexity,
manual coding of models becomes prohibitively labor-intensive. We demonstrate,
via a set of case studies, the new built-in model-prototyping capabilities of
XRate (macros and Scheme extensions). These features allow rapid implementation
of phylogenetic models which would have previously been far more
labor-intensive. XRate's new capabilities for lineage-specific models,
ancestral sequence reconstruction, and improved annotation output are also
discussed. XRate's flexible model-specification capabilities and computational
efficiency make it well-suited to developing and prototyping phylogenetic
grammar models. XRate is available as part of the DART software package:
http://biowiki.org/DART .


  Motivated by the lack of a suitable constructive framework for analyzing
popular stochastic models of Systems Biology, we devise conditions for
existence and uniqueness of solutions to certain jump stochastic differential
equations (SDEs). Working from simple examples we find reasonable and explicit
assumptions on the driving coefficients for the SDE representation to make
sense. By `reasonable' we mean that stronger assumptions generally do not hold
for systems of practical interest. In particular, we argue against the
traditional use of global Lipschitz conditions and certain common growth
restrictions. By `explicit', finally, we like to highlight the fact that the
various constants occurring among our assumptions all can be determined once
the model is fixed.
  We show how basic long time estimates and some limit results for
perturbations can be derived in this setting such that these can be contrasted
with the corresponding estimates from deterministic dynamics. The main
complication is that the natural path-wise representation is generated by a
counting measure with an intensity that depends nonlinearly on the state.


  The Multiple Sequence Alignment (MSA) is a computational abstraction that
represents a partial summary either of indel history, or of structural
similarity. Taking the former view (indel history), it is possible to use
formal automata theory to generalize the phylogenetic likelihood framework for
finite substitution models (Dayhoff's probability matrices and Felsenstein's
pruning algorithm) to arbitrary-length sequences. In this paper, we report
results of a simulation-based benchmark of several methods for reconstruction
of indel history. The methods tested include a relatively new algorithm for
statistical marginalization of MSAs that sums over a stochastically-sampled
ensemble of the most probable evolutionary histories. For mammalian
evolutionary parameters on several different trees, the single most likely
history sampled by our algorithm appears less biased than histories
reconstructed by other MSA methods. The algorithm can also be used for
alignment-free inference, where the MSA is explicitly summed out of the
analysis. As an illustration of our method, we discuss reconstruction of the
evolutionary histories of human protein-coding genes.


  Echinoderms take many forms of symmetry. Pentameral symmetry is the major
form and the other forms are derived from it. However, the ancestors of
echinoderms, which originated from Cambrian period, were believed to be
bilaterians. Echinoderm larvae are bilateral during their early development.
During embryonic development of starfish and sea urchins, the position and the
developmental sequence of each arm are fixed, implying an auxological
anterior/posterior axis. Starfish also possess the Hox gene cluster, which
controls symmetrical development. Overall, echinoderms are thought to have a
bilateral developmental mechanism and process. In this article, we focused on
adult starfish behaviors to corroborate its bilateral tendency. We weighed
their central disk and each arm to measure the position of the center of
gravity. We then studied their turning-over behavior, crawling behavior and
fleeing behavior statistically to obtain the center of frequency of each
behavior. By joining the center of gravity and each center of frequency, we
obtained three behavioral symmetric planes. These behavioral bilateral
tendencies might be related to the A/P axis during the embryonic development of
the starfish. It is very likely that the adult starfish is, to some extent,
bilaterian because it displays some bilateral propensity and has a definite
behavioral symmetric plane. The remainder of bilateral symmetry may have
benefited echinoderms during their evolution from the Cambrian period to the
present.


  We monitored twelve acoustically-tagged small pelagic fish (Selar
crumenophthalmus) around a floating object in shallow water, playing the role
of a coastal fish aggregating device (FAD). We characterized the response of
the tagged-fish aggregation to varying current strengths and daylight. We found
that the current induced a displacement of the aggregation upstream of the FAD,
at distances that were increasing with the current strength. We gave evidence,
of an expansion and a higher coordination in the aggregation at dusk, with
increasing swimming speed, distances among congeners and alignment. We
discussed possible scenarios where fish polarization increases at dusk and
proposed complementary measurements in future experiments that could confirm
our findings.


  The concept of a fitness landscape is a powerful metaphor that offers insight
into various aspects of evolutionary processes and guidance for the study of
evolution. Until recently, empirical evidence on the ruggedness of these
landscapes was lacking, but since it became feasible to construct all possible
genotypes containing combinations of a limited set of mutations, the number of
studies has grown to a point where a classification of landscapes becomes
possible. The aim of this review is to identify measures of epistasis that
allow a meaningful comparison of fitness landscapes and then apply them to the
empirical landscapes to discern factors that affect ruggedness. The various
measures of epistasis that have been proposed in the literature appear to be
equivalent. Our comparison shows that the ruggedness of the empirical landscape
is affected by whether the included mutations are beneficial or deleterious and
by whether intra- or intergenic epistasis is involved. Finally, the empirical
landscapes are compared to landscapes generated with the Rough Mt.\ Fuji model.
Despite the simplicity of this model, it captures the features of the
experimental landscapes remarkably well.


  Biochemical reaction systems in mesoscopic volume, under sustained
environmental chemical gradient(s), can have multiple stochastic attractors.
Two distinct mechanisms are known for their origins: ($a$) Stochastic
single-molecule events, such as gene expression, with slow gene on-off
dynamics; and ($b$) nonlinear networks with feedbacks. These two mechanisms
yield different volume dependence for the sojourn time of an attractor. As in
the classic Arrhenius theory for temperature dependent transition rates, a
landscape perspective provides a natural framework for the system's behavior.
However, due to the nonequilibrium nature of the open chemical systems, the
landscape, and the attractors it represents, are all themselves {\em emergent
properties} of complex, mesoscopic dynamics. In terms of the landscape, we show
a generalization of Kramers' approach is possible to provide a rate theory. The
emergence of attractors is a form of self-organization in the mesoscopic
system; stochastic attractors in biochemical systems such as gene regulation
and cellular signaling are naturally inheritable via cell division.
Delbr\"{u}ck-Gillespie's mesoscopic reaction system theory, therefore, provides
a biochemical basis for spontaneous isogenetic switching and canalization.


  Identifying spatial and temporal patterns can reveal the driving factors that
govern the behavior of fish in their environment. In this study, we
characterized the spatial and temporal occupation of 37 acoustically tagged
bigeye scads (Selar Crumenophthalmus) in an array of shallow Fish Aggregating
Devices (FADs) to clarify the mechanism that leads fish to associate with FADs.
A comparison of the number of visits and residence times exhibited by the fish
at different FADs revealed a strong variability over the array of FADS, with
the emergence of a leading FAD, which recorded the majority of visits and
retained the fish for a longer period of time. We found diel variability in the
residence times, with fish associated at daytime and exploring the array of
FADs at nighttime. We demonstrated that this diel temporal pattern was
amplified in the leading FAD. We identified a 24-hour periodicity for a subset
of individuals aggregated to the leading FAD, thus suggesting that those fish
were able to find this FAD after night excursions. The modeling of fish
movements based on a Monte Carlo sampling of inter-FAD transitions revealed
that the observed spatial heterogeneity in the number of visits could not be
explained through simple array-connectivity arguments. Similarly, we
demonstrated that the high residence times recorded at the leading FAD were not
due to the spatial arrangement of individual fish having different associative
characters. We discussed the relationships between these patterns of
association with the FADs, the exploration of the FAD array and the possible
effects of social interactions and environmental factors.


  We examine the budding of a nanoscale particle through a lipid bilayer using
molecular dynamics simulations, free energy calculations, and an elastic
theory, with the aim of determining the extent to which equilibrium elasticity
theory can describe the factors that control the mechanism and efficiency of
budding. The particle is a smooth sphere which experiences attractive
interactions to the lipid head groups. Depending on the parameters, we observe
four classes of dynamical trajectories: particle adhesion to the membrane,
stalled partially wrapped states, budding followed by scission, and membrane
rupture. In most regions of parameter space we find that the elastic theory
agrees nearly quantitatively with the simulated phase behavior as a function of
adhesion strength, membrane bending rigidity, and particle radius. However, at
parameter values near the transition between particle adhesion and budding, we
observe long-lived partially wrapped states which are not captured by existing
elastic theories. These states could constrain the accessible system parameters
for those enveloped viruses or drug delivery vehicles which rely on exo- or
endocytosis for membrane transport.


  Chemical purity of RNA samples is critical for high-precision studies of RNA
folding and catalytic behavior, but such purity may be compromised by
photodamage accrued during ultraviolet (UV) visualization of gel-purified
samples. Here, we quantitatively assess the breadth and extent of such damage
by using reverse transcription followed by single-nucleotide-resolution
capillary electrophoresis. We detected UV-induced lesions across a dozen
natural and artificial RNAs including riboswitch domains, other non-coding
RNAs, and artificial sequences; across multiple sequence contexts, dominantly
at but not limited to pyrimidine doublets; and from multiple lamps that are
recommended for UV shadowing in the literature. Most strikingly, irradiation
time-courses reveal detectable damage within a few seconds of exposure, and
these data can be quantitatively fit to a 'skin effect' model that accounts for
the increased exposure of molecules near the top of irradiated gel slices. The
results indicate that 200-nucleotide RNAs subjected to 20 seconds or less of UV
shadowing can incur damage to 20% of molecules, and the molecule-by-molecule
distribution of these lesions is more heterogeneous than a Poisson
distribution. Photodamage from UV shadowing is thus likely a widespread but
unappreciated cause of artifactual heterogeneity in quantitative and
single-molecule-resolution RNA biophysical measurements.


  Migrations have played an important role in shaping the genetic diversity of
human populations. Understanding genomic data thus requires careful modeling of
historical gene flow. Here we consider the effect of relatively recent
population structure and gene flow, and interpret genomes of individuals that
have ancestry from multiple source populations as mosaics of segments
originating from each population. We propose general and tractable models for
describing the evolution of these patterns of local ancestry and their impact
on genetic diversity. We focus on the length distribution of continuous
ancestry tracts, and the variance in total ancestry proportions among
individuals. The proposed models offer improved agreement with Wright-Fisher
simulation data when compared to state-of-the art models, and can be used to
infer various demographic parameters in gene flow models. Considering HapMap
African-American (ASW) data, we find that a model with two distinct phases of
`European' gene flow significantly improves the modeling of both tract lengths
and ancestry variances.


  Introduction : Mass spectrometry approaches are very attractive to detect
protein panels in a sensitive and high speed way. MS can be coupled to many
proteomic separation techniques. However, controlling technological variability
on these analytical chains is a critical point. Adequate information processing
is mandatory for data analysis to take into account the complexity of the
analysed mixture, to improve the measurement reliability and to make the
technology user friendly. Therefore we develop a hierarchical parametric
probabilistic model of the LC-MS analytical chain including the technological
variability. We introduce a Bayesian reconstruction methodology to recover the
protein biomarkers content in a robust way. We will focus on the digestion step
since it brings a major contribution to technological variability. Method : In
this communication, we introduce a hierarchical model of the LC-MS analytical
chain. Such a chain is a cascade of molecular events depicted by a graph
structure, each node being associated to a molecular state such as protein,
peptide and ion and each branch to a molecular processing such as digestion,
ionisation and LC-MS separation. This molecular graph defines a hierarchical
mixture model. We extend the Bayesian statistical framework we have introduced
previously [1] to this hierarchical description. As an example, we will
consider the digestion step. We describe the digestion process on a pair of
peptides within the targeted protein as a Bernoulli random process associated
with a cleavage probability controlled by the digestion kinetic law.


  In most biological studies and processes, cell proliferation and population
dynamics play an essential role. Due to this ubiquity, a multitude of
mathematical models has been developed to describe these processes. While the
simplest models only consider the size of the overall populations, others take
division numbers and labeling of the cells into account. In this work, we
present a modeling and computational framework for proliferating cell
population undergoing symmetric cell division. In contrast to existing models,
the proposed model incorporates both, the discrete age structure and continuous
label dynamics. Thus, it allows for the consideration of division number
dependent parameters as well as the direct comparison of the model prediction
with labeling experiments, e.g., performed with Carboxyfluorescein succinimidyl
ester (CFSE). We prove that under mild assumptions the resulting system of
coupled partial differential equations (PDEs) can be decomposed into a system
of ordinary differential equations (ODEs) and a set of decoupled PDEs, which
reduces the computational effort drastically. Furthermore, the PDEs are solved
analytically and the ODE system is truncated, which allows for the prediction
of the label distribution of complex systems using a low-dimensional system of
ODEs. In addition to modeling of labeling dynamics, we link the label-induced
fluorescence to the measure fluorescence which includes autofluorescence. For
the resulting numerically challenging convolution integral, we provide an
analytical approximation. This is illustrated by modeling and simulating a
proliferating population with division number dependent proliferation rate.


  The binding of small metal ions to complex macromolecular structures is
typically dominated by strong local interactions of the ion with its nearest
ligands. For this reason, it is often possible to understand the microscopic
origin of ion binding selectivity by considering simplified reduced models
comprised of only the nearest ion-coordinating ligands. Although the main
ingredients underlying simplified reduced models are intuitively clear, a
formal statistical mechanical treatment is nonetheless necessary in order to
draw meaningful conclusions about complex macromolecular systems. By
construction, reduced models only treat the ion and the nearest coordinating
ligands explicitly. The influence of the missing atoms from the protein or the
solvent is incorporated indirectly. Quasi-chemical theory offers one example of
how to carry out such a separation in the case of ion solvation in bulk
liquids, and in several ways, a statistical mechanical formulation of reduced
binding site models for macromolecules is expected to follow a similar route.
Here, some critical issues with recent theories of reduced binding site models
are examined.


  The history of life on Earth and in other potential life-bearing planetary
platforms is deeply linked to the history of the universe. Since life as we
know it relies on chemical elements forged in dying heavy stars, the universe
needs to be old enough for stars to form and evolve. Current cosmological
theory indicates that the universe is 13.7$\pm 0.13$ billion years old and that
the first stars formed hundreds of millions of years after the big bang. At
least some stars formed with stable planetary systems wherein a set of
biochemical reactions leading to life could have taken place. In this lecture,
I argue that we can divide cosmological history into four ages, from the big
bang to intelligent life. The Physical Age describes the origin of the
universe, of matter, of cosmic nucleosynthesis, as well as the formation of the
first stars and galaxies. The Chemical Age begun when heavy stars provided the
raw ingredients for life through stellar nucleosynthesis and describes how
heavier chemical elements collected in nascent planets and moons to give rise
to prebiotic biomolecules. The Biological Age describes the origin of early
life, its evolution through Darwinian natural selection, and the emergence of
complex multicellular life forms. Finally, the Cognitive Age describes how
complex life evolved into intelligent life capable of self-awareness and of
developing technology through the directed manipulation of energy and
materials. We conclude discussing whether we are the rule or the exception.


  A key open question in the study of life is the origin of biomolecular
homochirality: almost every life-form on Earth has exclusively levorotary amino
acids and dextrorotary sugars. Will the same handedness be preferred if life is
found elsewhere? We review some of the pertinent literature and discuss recent
results suggesting that life's homochirality resulted from sequential chiral
symmetry breaking triggered by environmental events. In one scenario,
autocatalytic prebiotic reactions undergo stochastic fluctuations due to
environmental disturbances. In another, chiral-selective polymerization
reaction rates influenced by environmental effects lead to substantial chiral
excess even in the absence of autocatalysis. Applying these arguments to other
potentially life-bearing platforms has implications to the search for
extraterrestrial life: we predict that a statistically representative sampling
of extraterrestrial stereochemistry will be racemic (chirally neutral) on
average.


  We discuss chemical reaction networks and metabolic pathways based on
stoichiometric network analysis, and introduce deformed toric ideal constraints
by the algebraic geometrical approach. This paper concerns steady state flux of
chemical reaction networks and metabolic pathways. With the deformed toric
ideal constraints, the linear combination parameters of extreme pathways are
automatically constrained without introducing ad hoc constraints. To illustrate
the effectiveness of such constraints, we discuss two examples of chemical
reaction network and metabolic pathway; in the former the flux and the
concentrations are constrained completely by deformed toric ideal constraints,
and in the latter, it is shown the deformed toric ideal constrains the linear
combination parameters of flux at least partially. Even in the latter case, the
flux and the concentrations are constrained completely with the additional
constraint that the total amount of enzyme is constant.


  This paper applies I.M. Gelfand's distinction between adequate and
non-adequate use of mathematical language in different contexts to the newly
opened window of model-based measurements of intracellular dynamics. The
specifics of geometry and dynamics on the mesoscale of cell physiology are
elaborated - in contrast to the familiar Newtonian mechanics and the more
recent, but by now also rather well established quantum field theories.
Examples are given originating from the systems biology of insulin secreting
pancreatic beta-cells and the mathematical challenges of an envisioned
non-invasive control of magnetic nanoparticles.


  Phylogenetic analysis of DNA or other data commonly gives rise to a
collection or sample of inferred evolutionary trees. Principal Components
Analysis (PCA) cannot be applied directly to collections of trees since the
space of evolutionary trees on a fixed set of taxa is not a vector space. This
paper describes a novel geometrical approach to PCA in tree-space that
constructs the first principal path in an analogous way to standard linear
Euclidean PCA. Given a data set of phylogenetic trees, a geodesic principal
path is sought that maximizes the variance of the data under a form of
projection onto the path. Due to the high dimensionality of tree-space and the
nonlinear nature of this problem, the computational complexity is potentially
very high, so approximate optimization algorithms are used to search for the
optimal path. Principal paths identified in this way reveal and quantify the
main sources of variation in the original collection of trees in terms of both
topology and branch lengths. The approach is illustrated by application to
simulated sets of trees and to a set of gene trees from metazoan (animal)
species.


  We consider the translocation of a one-dimensional polymer through a pore
channel helped by a motor driven by a dichotomous noise with time exponential
correlation. We are interested in the study of the translocation time, mean
velocity and stall force of the system as a function of the mean driving
frequency. We find a monotonous translocation time, in contrast with the mean
velocity which shows a pronounced maximum at a given frequency. Interestingly,
the stall force shows a nonmonotonic behavior with the presence of a minimum.
The influence of the spring elastic constant to the mean translocation times
and velocities is also presented.


  For nearly 50 years microbiologists have been determining prokaryotic genome
relatedness by means of nucleic acid reassociation kinetics. These methods,
however, are technically challenging, difficult to reproduce, and - given the
time and resources it takes to generate a single data-point - not cost
effective. In the post genomic era, with the cost of sequencing whole
prokaryotic genomes no longer a limiting factor, we believed that
computationally predicting the output value from a traditional DNA-DNA
hybridization experiment using pair-wise comparisons of whole genome sequences
to be of value. While other computational whole-genome classification methods
exist, they predict values on widely different scales than DNA-DNA
hybridization, introducing yet another metric into the polyphasic approach of
defining microbial species. Our goal was to develop an in silico BLAST based
pipeline that would predict with a high level of certainty the value of the wet
lab-based DNA-DNA hybridization values. Here we report on one such method that
produces estimates that are both accurate and precise with respect to the
DNA-DNA hybridization values they are designed to emulate.