  It is well known that the general transcription factors (GTF) specifically
recognize correct TATA boxes, distinguishing them from many others. Employing
the principles of determinacy analysis (mathematical theory of rules) we
analyzed a fragment of human chromosome 1 DNA sequence and identified specific
genetic markers (IG-markers = Irreducible Genetic markers) in the nearest
proximity to TATAAA motifs. The IG-markers enable determining the exact
location of any TATAAA motif within the investigated DNA fragment. Based on our
data we hypothesize that the GTF recognize the
{\guillemotleft}true{\guillemotright} transcriptional start TATA box by means
of IG-markers. The math method described here is universal and can be used to
find IG-markers that will provide, like a global navigation satellite system,
for the specific location of any distinct sequence motif within larger DNA
sequence content.


  In a recent paper it was shown that, for chemical reaction networks
possessing a subtle structural property called concordance, dynamical behavior
of a very circumscribed (and largely stable) kind is enforced, so long as the
kinetics lies within the very broad and natural weakly monotonic class. In
particular, multiple equilibria are precluded, as are degenerate positive
equilibria. Moreover, under certain circumstances, also related to concordance,
all real eigenvalues associated with a positive equilibrium are negative.
Although concordance of a reaction network can be decided by readily available
computational means, we show here that, when a nondegenerate network's
Species-Reaction Graph satisfies certain mild conditions, concordance and its
dynamical consequences are ensured. These conditions are weaker than earlier
ones invoked to establish kinetic system injectivity, which, in turn, is just
one ramification of network concordance. Because the Species-Reaction Graph
resembles pathway depictions often drawn by biochemists, results here expand
the possibility of inferring signicant dynamical information directly from
standard biochemical reaction diagrams.


  We introduce a family of Rock-Paper-Scissors type models with $Z_N$ symmetry
($N$ is the number of species) and we show that it has a very rich structure
with many completely different phases. We study realizations which lead to the
formation of domains, where individuals of one or more species coexist,
separated by interfaces whose (average) dynamics is curvature driven. This type
of behavior, which might be relevant for the development of biological
complexity, leads to an interface network evolution and pattern formation
similar to the ones of several other nonlinear systems in condensed matter and
cosmology.


  In this work we study a modified Susceptible-Infected-Susceptible (SIS) model
in which the infection rate $\lambda$ decays exponentially with the number of
reinfections $n$, saturating after $n=l$. We find a critical decaying rate
$\epsilon_{c}(l)$ above which a finite fraction of the population becomes
permanently infected. From the mean-field solution and computer simulations on
hypercubic lattices we find evidences that the upper critical dimension is 6
like in the SIR model, which can be mapped in ordinary percolation.


  Ever since the disorder of proteins is the main cause for many diseases. As
compared with other disorders, the major reason that causes disease is of
structural inability of many proteins. The potentially imminent availability of
recent datasets helps one to discover the protein disorders, however in
majority of cases, the stability of proteins depend on the carbon content.
Addressing this distinct feature, it is possible to hit upon the carbon
distribution along the sequence and can easily recognize the stable nature of
protein. There are certain reported mental disorders which fall in to this
category. Regardless, such kind of disorder prone protein FMR1p (Fragile X
mental retardation 1 protein) is identified as the main cause for the disease
Fragile X syndrome. This paper deals with the identification of defects in the
FMR1 protein sequence considering the carbon contents along the sequence. This
attempt is to evaluate the stability of proteins, accordingly the protein
disorders in order to improvise the certain Biological functions of proteins to
prevent disease. The transition of the disorder to order protein involves
careful considerations and can be achieved by detecting the unstable region
that lacks hydrophobicity. This work focuses the low carbon content in the FMR1
protein so as to attain the stable status in future to reduce the morbidity
rate caused by Fragile X syndrome for the society.


  The cytosol state in living cell is treated as homogeneous phase equilibrium
with a special feature: the pressure of one phase is positive and the pressure
of the other is negative. From this point of view the cytosol is neither
solution nor gel (or sol as a whole) regardless its components (water and
dissolved substances). This is its unique capability for selecting, sorting and
transporting reagents to the proper place of the living cell without a
so-called "pipeline". To base this statement the theoretical investigation of
the conditions of equilibrium and stability of the medium with alternative-sign
pressure is carried out under using the thermodynamic laws and the Gibbs'
equilibrium criterium.


  Cooperation is crucial for the remarkable evolutionary success of the human
species. Not surprisingly, some individuals are willing to bare additional
costs in order to punish defectors. Current models assume that, once set, the
fine and cost of punishment do not change over time. Here we show that relaxing
this assumption by allowing players to adapt their sanctioning efforts in
dependence on the success of cooperation can explain both, the spontaneous
emergence of punishment, as well as its ability to deter defectors and those
unwilling to punish them with globally negligible investments. By means of
phase diagrams and the analysis of emerging spatial patterns, we demonstrate
that adaptive punishment promotes public cooperation either through the
invigoration of spatial reciprocity, the prevention of the emergence of cyclic
dominance, or through the provision of competitive advantages to those that
sanction antisocial behavior. Presented results indicate that the process of
self-organization significantly elevates the effectiveness of punishment, and
they reveal new mechanisms by means of which this fascinating and widespread
social behavior could have evolved.


  The statistical mechanics of DNA denaturation under fixed linking number is
qualitatively different from that of the unconstrained DNA. Quantitatively
different melting scenarios are reached from two alternative assumptions,
namely, that the denatured loops are formed in expense of 1) overtwist, 2)
supercoils. Recent work has shown that the supercoiling mechanism results in a
BEC-like picture where a macroscopic loop appears at Tc and grows steadily with
temperature, while the nature of the denatured phase for the overtwisting case
has not been studied. By extending an earlier result, we show here that a
macroscopic loop appears in the overtwisting scenario as well. We calculate its
size as a function of temperature and show that the fraction of the total sum
of microscopic loops decreases above Tc, with a cusp at the critical point.


  In this article the Gordan theorem is applied to the thermodynamics of a
chemical reaction network at steady state. From a theoretical viewpoint it is
equivalent to the Clausius formulation of the second law for the out of
equilibrium steady states of chemical networks, i.e. it states that the
exclusion (presence) of closed reactions loops makes possible (impossible) the
definition of a thermodynamic potential and vice versa. On the computational
side, it reveals that calculating reactions free energy and searching
infeasible loops in flux states are dual problems whose solutions are
alternatively inconsistent. The relevance of this result for applications is
discussed with an example in the field of constraints-based modeling of
cellular metabolism where it leads to efficient and scalable methods to afford
the energy balance analysis.


  We propose a general method for predicting potentially good folders from a
given number of amino acid sequences. Our approach is based on the calculation
of the rate of convergence of each amino acid chain towards the native
structure using only the very initial parts of the dynamical trajectories. It
does not require any preliminary knowledge of the native state and can be
applied to different kinds of models, including atomistic descriptions. We
tested the method within both the lattice and off-lattice model frameworks and
obtained several so far unknown good folders. The unbiased algorithm also
allows to determine the optimal folding temperature and takes at least 3--4
orders of magnitude less time steps than those needed to compute folding times.


  Certain regulatory proteins influence the polymerization dynamics of
microtubules by inducing catastrophe with a rate that depends on the
microtubule length. Using a discrete formulation, here we show that, for a
catastrophe rate proportional to the microtubule length, the steady-state
probability distributions of length decay much faster with length than an
exponential decay as seen in the absence of these proteins.


  In 2009, for the first time in Cape Verde, an outbreak of dengue was reported
and over twenty thousand people were infected. Only a few prophylactic measures
were taken. The effects of vector control on disease spreading, such as
insecticide (larvicide and adulticide) and mechanical control, as well as an
hypothetical vaccine, are estimated through simulations with the Cape Verde
data.


  Without mutation and migration, evolutionary dynamics ultimately leads to the
extinction of all but one species. Such fixation processes are well understood
and can be characterized analytically with methods from statistical physics.
However, many biological arguments focus on stationary distributions in a
mutation-selection equilibrium. Here, we address the equilibration time
required to reach stationarity in the presence of mutation, this is known as
the mixing time in the theory of Markov processes. We show that mixing times in
evolutionary games have the opposite behaviour from fixation times when the
intensity of selection increases: In coordination games with bistabilities, the
fixation time decreases, but the mixing time increases. In coexistence games
with metastable states, the fixation time increases, but the mixing time
decreases. Our results are based on simulations and the WKB approximation of
the master equation.


  We study DNA self-assembly and DNA computation using a coarse-grained DNA
model within the directional dynamic bonding framework {[}C. Svaneborg, Comp.
Phys. Comm. 183, 1793 (2012){]}. In our model, a single nucleotide or domain is
represented by a single interaction site. Complementary sites can reversibly
hybridize and dehybridize during a simulation. This bond dynamics induces a
dynamics of the angular and dihedral bonds, that model the collective effects
of chemical structure on the hybridization dynamics. We use the DNA model to
perform simulations of the self-assembly kinetics of DNA tetrahedra, an
icosahedron, as well as strand displacement operations used in DNA computation.


  Adaptive dynamics is a widely used framework for modeling long-term evolution
of continuous phenotypes. It is based on invasion fitness functions, which
determine selection gradients and the canonical equation of adaptive dynamics.
Even though the derivation of the adaptive dynamics from a given invasion
fitness function is general and model-independent, the derivation of the
invasion fitness function itself requires specification of an underlying
ecological model. Therefore, evolutionary insights gained from adaptive
dynamics models are generally model-dependent. Logistic models for symmetric,
frequency-dependent competition are widely used in this context. Such models
have the property that the selection gradients derived from them are gradients
of scalar functions, which reflects a certain gradient property of the
corresponding invasion fitness function. We show that any adaptive dynamics
model that is based on an invasion fitness functions with this gradient
property can be transformed into a generalized symmetric competition model.
This provides a precise delineation of the generality of results derived from
competition models. Roughly speaking, to understand the adaptive dynamics of
the class of models satisfying a certain gradient condition, one only needs a
complete understanding of the adaptive dynamics of symmetric,
frequency-dependent competition. We show how this result can be applied to
number of basic issues in evolutionary theory.


  We present McGenus, an algorithm to predict RNA secondary structures with
pseudoknots. The method is based on a classification of RNA structures
according to their topological genus. McGenus can treat sequences of up to 1000
bases and performs an advanced stochastic search of their minimum free energy
structure allowing for non trivial pseudoknot topologies. Specifically, McGenus
employs a multiple Markov chain scheme for minimizing a general scoring
function which includes not only free energy contributions for pair stacking,
loop penalties, etc. but also a phenomenological penalty for the genus of the
pairing graph. The good performance of the stochastic search strategy was
successfully validated against TT2NE which uses the same free energy
parametrization and performs exhaustive or partially exhaustive structure
search, albeit for much shorter sequences (up to 200 bases). Next, the method
was applied to other RNA sets, including an extensive tmRNA database, yielding
results that are competitive with existing algorithms. Finally, it is shown
that McGenus highlights possible limitations in the free energy scoring
function. The algorithm is available as a web-server at
http://ipht.cea.fr/rna/mcgenus.php .


  This article presents an algorithm that allows modeling of biological
networks in a qualitative framework with continuous time. Mathematical modeling
is used as a systems biology tool to answer biological questions, and more
precisely, to validate a network that describes biological observations and to
predict the effect of perturbations.
  We propose a modeling approach that is intrinsically continuous in time. The
algorithm presented here fills the gap between qualitative and quantitative
modeling. It is based on continuous time Markov process applied on a Boolean
state space. In order to describe the temporal evolution, we explicitly specify
the transition rates for each node. For that purpose, we built a language that
can be seen as a generalization of Boolean equations. The values of transition
rates have a natural interpretation: it is the inverse of the time for the
transition to occur. Mathematically, this approach can be translated in a set
of ordinary differential equations on probability distributions; therefore, it
can be seen as an approach in between quantitative and qualitative.
  We developed a C++ software, MaBoSS, that is able to simulate such a system
by applying Kinetic Monte-Carlo (or Gillespie algorithm) in the Boolean state
space. This software, parallelized and optimized, computes temporal evolution
of probability distributions and can also estimate stationary distributions.
Applications of Boolean Kinetic Monte-Carlo have been demonstrated for two
qualitative models: a toy model and a published p53/Mdm2 model. Our approach
allows to describe kinetic phenomena which were difficult to handle in the
original models. In particular, transient effects are represented by time
dependent probability distributions, interpretable in terms of cell
populations.


  We use the worm-like chain model to study supercoiling of DNA under tension
and torque. The model reproduces experimental data for a much broader range of
forces, salt concentrations and contour lengths than previous approaches. Our
theory shows, for the first time, how the behavior of the system is controlled
by a multi-plectoneme phase in a wide range of parameters. This phase does not
only affect turn-extension curves but also leads to a non-constant torque in
the plectonemic phase. Shortcomings from previous models and inconsistencies
between experimental data are resolved in our theory without the need of
adjustable parameters.


  Metabotropic glutamate receptor 5 (mGluR5) null mutant (-/-) mice have been
reported to totally lack the rein- forcing or locomotor stimulating effects of
cocaine. We tested mGluR5 -/- and +/+ mice for their locomotor and conditioned
place- preference response to cocaine. Unlike the previous finding, here we
show that compared to mGluR5 +/+ mice, -/- mice exhibit no difference in the
locomotor response to low to moderate doses of cocaine (10 or 20 mg/kg). A high
dose of cocaine (40 mg/kg) resulted in a blunted rather than absent locomo- tor
response. We tested mGluR5 -/- and +/+ mice for conditioned place-preference to
cocaine and found no group differences at a conditioning dose of 10 mg/kg,
suggesting normal conditioned rewarding properties of cocaine. These results
differ substantially from Chiamulera et al. (2001) and replicates Olsen et al.,
(2010), who found normal cocaine place-preference in mGluR5 -/- mice at 5
mg/kg. Our results indicate mGluR5 receptors exert a modulatory rather than
necessary role in cocaine-induced locomotor stimulation and exert no effect on
the conditioned rewarding effects of cocaine.


  We simulate the nonlocal Stokesian hydrodynamics of an elastic filament which
is active due a permanent distribution of stresslets along its contour. A
bending instability of an initially straight filament spontaneously breaks flow
symmetry and leads to autonomous filament motion which, depending on
conformational symmetry, can be translational or rotational. At high ratios of
activity to elasticity, the linear instability develops into nonlinear
fluctuating states with large amplitude deformations. The dynamics of these
states can be qualitatively understood as a superposition of translational and
rotational motion associated with filament conformational modes of opposite
symmetry. Our results can be tested in molecular-motor filament mixtures,
synthetic chains of autocatalytic particles, or other linearly connected
systems where chemical energy is converted to mechanical energy in a fluid
environment.


  Curvatures in the most probable rupture force ($f^*$) versus log-loading rate
($\log{r_f}$) observed in dynamic force spectroscopy (DFS) on biomolecular
complexes are interpreted using a one-dimensional free energy profile with
multiple barriers or a single barrier with force-dependent transition state.
Here, we provide a criterion to select one scenario over another. If the
rupture dynamics occurs by crossing a single barrier in a physical free energy
profile describing unbinding, the exponent $\nu$, from $(1-
f^*/f_c)^{1/\nu}\sim(\log r_f)$ with $f_c$ being a critical force in the
absence of force, is restricted to $0.5 \leq \nu \leq 1$. For biotin-ligand
complexes and leukocyte-associated antigen-1 bound to intercellular adhesion
molecules, which display large curvature in the DFS data, fits to experimental
data yield $\nu<0.5$, suggesting that ligand unbinding is associated with
multiple-barrier crossing.


  The Price equation partitions total evolutionary change into two components.
The first component provides an abstract expression of natural selection. The
second component subsumes all other evolutionary processes, including changes
during transmission. The natural selection component is often used in
applications. Those applications attract widespread interest for their
simplicity of expression and ease of interpretation. Those same applications
attract widespread criticism by dropping the second component of evolutionary
change and by leaving unspecified the detailed assumptions needed for a
complete study of dynamics. Controversies over approximation and dynamics have
nothing to do with the Price equation itself, which is simply a mathematical
equivalence relation for total evolutionary change expressed in an alternative
form. Disagreements about approach have to do with the tension between the
relative valuation of abstract versus concrete analyses. The Price equation's
greatest value has been on the abstract side, particularly the invariance
relations that illuminate the understanding of natural selection. Those
abstract insights lay the foundation for applications in terms of kin
selection, information theory interpretations of natural selection, and
partitions of causes by path analysis. I discuss recent critiques of the Price
equation by Nowak and van Veelen.


  The tryptophan metabolite, kynurenic acid (KYNA), is classically known to be
an antagonist of ionotropic glutamate receptors. Within the last decade several
reports have been published suggesting that KYNA also blocks nicotinic
acetylcholine receptors (nAChRs) containing the \alpha7 subunit (\alpha7*).
Most of these reports involve either indirect measurements of KYNA effects on
\alpha7 nAChR function, or are reports of KYNA effects in complicated in vivo
systems. However, a recent report investigating KYNA interactions with \alpha7
nAChRs failed to detect an interaction using direct measurements of \alpha7
nAChRs function. Further, it showed that a KYNA blockade of \alpha7 nAChR
stimulated GABA release (an indirect measure of \alpha7 nAChR function) was not
due to KYNA blockade of the \alpha7 nAChRs. The current study measured the
direct effects of KYNA on \alpha7-containing nAChRs expressed on interneurons
in the hilar and CA1 stratum radiatum regions of the mouse hippocampus and on
interneurons in the CA1 region of the rat hippocampus. Here we show that KYNA
does not block \alpha7* nACHRs using direct patch-clamprecording of \alpha7
currents in adult brain slices.


  X-ray crystallography, NMR (Nuclear Magnetic Resonance) spectroscopy, and
dual polarization interferometry, etc are indeed very powerful tools to
determine the 3D structures of proteins (including the membrane proteins),
though they are time-consuming and costly. However, for some proteins, due to
their unstable, noncrystalline and insoluble nature, these tools cannot work.
Under this condition, mathematical and physical theoretical methods and
computational approaches allow us to obtain a description of the protein 3D
structure at a submicroscopic level. This Chapter presents some practical and
useful mathematical optimization computational approaches to produce 3D
structures of the Prion AGAAAAGA Amyloid Fibrils, from a potential energy
minimization point of view.
  X-ray crystallography finds the X-ray final structure of a protein, which
usually need refinements in order to produce a better structure. The
computational methods presented in this Chapter can be also acted as a tool for
the refinements.


  Many forms of cell motility rely on Brownian ratchet mechanisms that involve
multiple stochastic processes. We present a computational and theoretical study
of the nonequilibrium statistical dynamics of such a many-body ratchet, in the
specific form of a growing polymer gel that pushes a diffusing obstacle. We
find that oft-neglected correlations among constituent filaments impact
steady-state kinetics and significantly deplete the gel's density within
molecular distances of its leading edge. These behaviors are captured
quantitatively by a self-consistent theory for extreme fluctuations in
filaments' spatial distribution.


  In this paper, we investigated a density-dependent reaction-diffusion
equation, $u_t = (u^{m})_{xx} + u - u^{m}$. This equation is known as the
extension of the Fisher or Kolmogoroff-Petrovsky-Piscounoff equation which is
widely used in the population dynamics, combustion theory and plasma physics.
By employing the suitable transformation, this equation was mapped to the
anomalous diffusion equation where the nonlinear reaction term was eliminated.
Due to its simpler form, some exact self-similar solutions with the compact
support have been obtained. The solutions, evolving from an initial state,
converge to the usual traveling wave at a certain transition time. Hence, it is
quite clear the connection between the self-similar solution and the traveling
wave solution from these results. Moreover, the solutions were found in the
manner that either propagates to the right or propagates to the left.
Furthermore, the two solutions form a symmetric solution, expanding in both
directions. The application on the spatiotemporal pattern formation in
biological population has been mainly focused.


  The chemotaxis pathway of the bacterium Rhodobacter sphaeroides has many
similarities to the well-studied pathway in Escherichia coli. It exhibits
robust adaptation and has several homologues of the latter's chemotaxis
proteins. Recent theoretical results have been able to correctly predict that
the chemotactic response of Escherichia coli exhibits the same output behavior
in response to scaled ligand inputs, a dynamic property known as fold-change
detection (FCD), or input-scale invariance. In this paper, we present
theoretical assumptions on the R. sphaeroides chemotaxis sensing dynamics that
can be analytically shown to yield FCD behavior in a specific ligand
concentration range. Based on these assumptions, we construct two models of the
full chemotaxis pathway that are able to reproduce experimental time-series
data from earlier studies. To test the validity of our assumptions, we propose
a series of experiments in which our models predict robust FCD behavior where
earlier models do not. In this way, we illustrate how a dynamic phenotype such
as FCD can be used for the purposes of discriminating between two models that
reproduce the same experimental time-series data.


  The phage $\lambda$ infection of an \textit{E. coli} cell has become a
paradigm for understanding the molecular processes involved in gene expression
and cell signaling. This system provides an example of a genetic switch, as
cells with identical DNA choose either of two cell cycles: a lysogenic cycle,
in which the phage genome is incorporated into the host and copied by the host;
or a lytic cycle, resulting in the death of the cell and a burst of viruses.
The robustness of this switch is remarkable; although the first stages of the
lysogenic and lytic cycles are identical, a lysogen rarely spontaneously flips,
and external stressors or instantaneous cell conditions are required to induce
flipping. In particular, the cell fate decision can depend on the populations
of two proteins, cI and Cro, as well as their oligomerization and subsequent
binding affinities to three DNA sites. These processes in turn govern the rates
at which RNAp transcribes the cI and Cro genes to produce more of their
respective proteins.
  In this work, a dynamical model of the non-equilibrium statistical mechanics
is revisited and generalized. The low number of proteins and other sources of
noise are non-negligble and corrections to the kinetics are essential to
understanding the stability. To this end, general integral forms for
advection-diffusion equations have been developed and numerically solved for a
variety of mutants and assumptions about the state of the cells. These
solutions quantify the probabilistic and flux landscapes of the ensembles'
evolution in concentration space and are used to predict the populations of the
cell states, entropy production, passage times, and potential barriers of wild
type and mutant bacteria to illuminate some structure of the configuration
space from which Nature naturally selects.


  Advanced Monte Carlo simulations are used to study the effect of nano-slit
confinement on metric and topological properties of model DNA chains. We
consider both linear and circularised chains with contour lengths in the
1.2--4.8 $\mu$m range and slits widths spanning continuously the 50--1250nm
range. The metric scaling predicted by de Gennes' blob model is shown to hold
for both linear and circularised DNA up to the strongest levels of confinement.
More notably, the topological properties of the circularised DNA molecules have
two major differences compared to three-dimensional confinement. First, the
overall knotting probability is non-monotonic for increasing confinement and
can be largely enhanced or suppressed compared to the bulk case by simply
varying the slit width. Secondly, the knot population consists of knots that
are far simpler than for three-dimensional confinement. The results suggest
that nano-slits could be used in nano-fluidic setups to produce DNA rings
having simple topologies (including the unknot) or to separate heterogeneous
ensembles of DNA rings by knot type.


  This thesis aims at the logical analysis of discrete processes, in particular
of such generated by gene regulatory networks. States, transitions and
operators from temporal logics are expressed in the language of Formal Concept
Analysis. By the attribute exploration algorithm, an expert or a computer
program is enabled to validate a minimal and complete set of implications, e.g.
by comparison of predictions derived from literature with observed data. Here,
these rules represent temporal dependencies within gene regulatory networks
including coexpression of genes, reachability of states, invariants or possible
causal relationships. This new approach is embedded into the theory of
universal coalgebras, particularly automata, Kripke structures and Labelled
Transition Systems. A comparison with the temporal expressivity of Description
Logics is made. The main theoretical results concern the integration of
background knowledge into the successive exploration of the defined data
structures (formal contexts). Applying the method a Boolean network from
literature modelling sporulation of Bacillus subtilis is examined. Finally, we
developed an asynchronous Boolean network for extracellular matrix formation
and destruction in the context of rheumatoid arthritis.


  Mutations create the genetic diversity on which selective pressures can act,
yet also create structural instability in proteins. How, then, is it possible
for organisms to ameliorate mutation-induced perturbations of protein stability
while maintaining biological fitness and gaining a selective advantage? Here we
used a new technique of site-specific chromosomal mutagenesis to introduce a
selected set of mostly destabilizing mutations into folA - an essential
chromosomal gene of E. coli encoding dihydrofolate reductase (DHFR) - to
determine how changes in protein stability, activity and abundance affect
fitness. In total, 27 E.coli strains carrying mutant DHFR were created. We
found no significant correlation between protein stability and its catalytic
activity nor between catalytic activity and fitness in a limited range of
variation of catalytic activity observed in mutants. The stability of these
mutants is strongly correlated with their intracellular abundance; suggesting
that protein homeostatic machinery plays an active role in maintaining
intracellular concentrations of proteins. Fitness also shows a significant
correlation with intracellular abundance of soluble DHFR in cells growing at
30oC. At 42oC, on the other hand, the picture was mixed, yet remarkable: a few
strains carrying mutant DHFR proteins aggregated rendering them nonviable, but,
intriguingly, the majority exhibited fitness higher than wild type. We found
that mutational destabilization of DHFR proteins in E. coli is counterbalanced
at 42oC by their soluble oligomerization, thereby restoring structural
stability and protecting against aggregation.


  Inference of new biological knowledge, e.g., prediction of protein function,
from protein-protein interaction (PPI) networks has received attention in the
post-genomic era. A popular strategy has been to cluster the network into
functionally coherent groups of proteins and predict protein function from the
clusters. Traditionally, network research has focused on clustering of nodes.
However, why favor nodes over edges, when clustering of edges may be preferred?
For example, nodes belong to multiple functional groups, but clustering of
nodes typically cannot capture the group overlap, while clustering of edges
can. Clustering of adjacent edges that share many neighbors was proposed
recently, outperforming different node clustering methods. However, since some
biological processes can have characteristic "signatures" throughout the
network, not just locally, it may be of interest to consider edges that are not
necessarily adjacent. Hence, we design a sensitive measure of the "topological
similarity" of edges that can deal with edges that are not necessarily
adjacent. We cluster edges that are similar according to our measure in
different baker's yeast PPI networks, outperforming existing node and edge
clustering approaches.


  Recent analysis of empirical data [F. Radicchi, A. Baronchelli & L.A.N.
Amaral. PloS ONE 7, e029910 (2012)] showed that humans adopt L\'evy flight
strategies when exploring the bid space in on-line auctions. A game theoretical
model proved that the observed L\'evy exponents are nearly optimal, being close
to the exponent value that guarantees the maximal economical return to players.
Here, we rationalize these findings by adopting an evolutionary perspective. We
show that a simple evolutionary process is able to account for the empirical
measurements with the only assumption that the reproductive fitness of a player
is proportional to her search ability. Contrarily to previous modeling, our
approach describes the emergence of the observed exponent without resorting to
any strong assumptions on the initial searching strategies. Our results
generalize earlier research, and open novel questions in cognitive, behavioral
and evolutionary sciences.


  The bioinformatical methods to detect lateral gene transfer events are mainly
based on functional coding DNA characteristics. In this paper, we propose the
use of DNA traits not depending on protein coding requirements. We introduce
several semilocal variables that depend on DNA primary sequence and that
reflect thermodynamic as well as physico-chemical magnitudes that are able to
tell apart the genome of different organisms. After combining these variables
in a neural classificator, we obtain results whose power of resolution go as
far as to detect the exchange of genomic material between bacteria that are
phylogenetically close.


  This is an extended abstract of lectures delivered at the 4th Warsaw Summer
School on Statistical Physics in Kazimierz Dolny, June 25-July 2, 2011.


  We argue that protein native state structures reside in a novel "phase" of
matter which confers on proteins their many amazing characteristics. This phase
arises from the common features of all globular proteins and is characterized
by a sequence-independent free energy landscape with relatively few low energy
minima with funnel-like character. The choice of a sequence that fits well into
one of these predetermined structures facilitates rapid and cooperative
folding. Our model calculations show that this novel phase facilitates the
formation of an efficient route for sequence design starting from random
peptides.


  Recently the A/H1N1-2009 virus pandemic appeared in Mexico and in other
nations. We present a study of this pandemic in the Mexican case using the SIR
model to describe epidemics. This model is one of the simplest models but it
has been a successful description of some epidemics of closed populations. We
consider the data for the Mexican case and use the SIR model to make some
predictions. Then, we generalize the SIR model in order to describe the spatial
dynamics of the disease. We make a study of the spatial and temporal spread of
the infected population with model parameters that are consistent with temporal
SIR model parameters obtained by fitting to the Mexican case.


  Food webs with intraguild predation (IGP) are widespread in natural habitats.
Their adaptation and resilience behaviour is principal for understanding
restructuring of ecological communities. In spite of the importance of IGP food
webs their behaviour even for the simplest 3-species systems has not been fully
explored. One fundamental question is how an increase of diversity of the
lowest trophic level impacts the persistence of higher trophic levels in IGP
relationships. We analyze a 3-species food web model with a heterogeneous
resources and IGP. The model consists of two predators directly coupled via IGP
relation and indirectly via competition for resource. The resource is
subdivided into distinct subpopulations. Individuals in the subpopulations are
grazed at different rates by the predators. We consider two models: an IGP
module with immobilization by the top predator and an IGP module with species
turnover. We examine the effect of increasing enrichment and varying
immobilization (resource transfer) rate on a stable coexistence of predators
and resources. We explore how the predictions from the basic 3-species model
are altered when the IGP module is extended to multiple resource
subpopulations. We investigate which parameters support a robust coexistence in
the IGP system. For the case of multiple subpopulations of the resource we
present a numerical comparison of the percentage of food webs with stable
coexistence for different dimensionalities of the resource community. At low
immobilization (transfer) rates our model predicts a stable 3-species
coexistence only for intermediate enrichment meanwhile at high rates a large
set of stable equilibrium configurations is found for high enrichment as well.


  We consider population models in which the individuals reproduce, die and
also migrate in space. The population size scales according to some parameter
$N$, which can have different interpretations depending on the context. Each
individual is assigned a mass of 1/N and the total mass in the system is called
\emph{population density}. The dynamics has an intrinsic density regulation
mechanism that drives the population density towards an equilibrium. We show
that under a timescale separation between the \emph{slow} migration mechanism
and the \emph{fast} density regulation mechanism, the population dynamics
converges to a Fleming-Viot process as the scaling parameter $N$ approaches
$\infty$. We first prove this result for a basic model in which the birth and
death rates can only depend on the population density. In this case we obtain a
\emph{neutral} Fleming-Viot process. We then extend this model by including
position-dependence in the birth and death rates, as well as, offspring
dispersal and immigration mechanisms. We show how these extensions add
\emph{mutation} and \emph{selection} to the limiting Fleming-Viot process. All
the results are proved in a multi-type setting, where there are $q$ types of
individuals interacting with each other. We illustrate the usefulness of our
convergence result by discussing applications in population genetics and cell
biology.


  Coupling is a widely used technique in the theoretical study of interacting
stochastic processes. In this paper I present an example demonstrating its
usefulness also in the efficient computer simulation of such processes. I first
describe a basic coupling technique, applicable to all kinds of processes,
which allows trading memory use for a limited speedup. Next, I describe a
specialized variant of it, which can be used to speed up the simulation certain
kinds of processes satisfying a monotonicity criterion. This special algorithm
increases the speed by several orders of magnitude with only a modest increase
in memory usage.


  Heterogeneities in environmental conditions often induce corresponding
heterogeneities in the distribution of species. In the extreme case of a
localized patch of increased growth rates, reproducing populations can become
strongly concentrated at the patch despite the entropic tendency for population
to distribute evenly. Several deterministic mathematical models have been used
to characterize the conditions under which localized states can form, and how
they break down due to convective driving forces. Here, we study the
delocalization of a finite population in the presence of number fluctuations.
We find that any finite population delocalizes on sufficiently long time
scales. Depending on parameters, however, populations may remain localized for
a very long time. The typical waiting time to delocalization increases
exponentially with both population size and distance to the critical wind speed
of the deterministic approximation. We augment these simulation results by a
mathematical analysis that treats the reproduction and migration of individuals
as branching random walks subject to global constraints. For a particular
constraint, different from a fixed population size constraint, this model
yields a solvable first moment equation. We find that this solvable model
approximates very well the fixed population size model for large populations,
but starts to deviate as population sizes are small. The analytical approach
allows us to map out a phase diagram of the order parameter as a function of
the two driving parameters, inverse population size and wind speed. Our results
may be used to extend the analysis of delocalization transitions to different
settings, such as the viral quasi-species scenario.


  Biochemical reaction networks are subjected to large fluctuations
attributable to small molecule numbers, yet underlie reliable biological
functions. Most theoretical approaches describe them as purely deterministic or
stochastic dynamical systems, depending on which point of view is favored.
Here, we investigate the dynamics of a self-repressing gene using an
intermediate approach based on a moment closure approximation of the master
equation, which allows us to take into account the binary character of gene
activity. We thereby obtain deterministic equations that describe how
nonlinearity feeds back fluctuations into the mean-field equations, providing
insight into the interplay of determinism and stochasticity. This allows us to
identify regions of parameter space where fluctuations induce relatively
regular oscillations.


  The biological requirements for an ecosystem to develop and maintain species
diversity are in general unknown. Here we consider a model ecosystem of sessile
and mutually excluding organisms competing for space [Mathiesen et al. Phys.
Rev. Lett. 107, 188101 (2011)]. The competition is controlled by an interaction
network with fixed links chosen by a Bernoulli process. New species are
introduced in the system at a predefined rate. In the limit of small
introduction rates, the system becomes bistable and can undergo a phase
transition from a state of low diversity to high diversity. We suggest that
patches of isolated meta-populations formed by the collapse of cyclic relations
are essential for the transition to the state of high diversity.


  We formulate a dynamical model for microtubules interacting with a
catastrophe-inducing boundary. In this model microtubules are either waiting to
be nucleated, actively growing or shrinking, or stalled at the boundary. We
first determine the steady-state occupation of these various states and the
resultant length distribution. Next, we formulate the problem of the Mean
First-Passage Time to reach the boundary in terms of an appropriate set of
splitting probabilities and conditional Mean First-Passage Times, and solve
explicitly for these quantities using a differential equation approach. As an
application, we revisit a recently proposed search-and-capture model for the
interaction between microtubules and target chromosomes [Gopalakrishnan &
Govindan, Bull. Math. Biol. 73:2483--506 (2011)]. We show how our approach
leads to a direct and compact solution of this problem.


  Recent experimental and theoretical approaches have attempted to quantify the
physical organization (compaction and geometry) of the bacterial chromosome
with its complement of proteins (the nucleoid). The genomic DNA exists in a
complex and dynamic protein-rich state, which is highly organised at various
length scales. This has implications on modulating (when not enabling) the core
biological processes of replication, transcription, segregation. We overview
the progress in this area, driven in the last few years by new scientific ideas
and new interdisciplinary experimental techniques, ranging from high space- and
time-resolution microscopy to high-throughput genomics employing sequencing to
map different aspects of the nucleoid-related interactome. The aim of this
review is to present the wide spectrum of experimental and theoretical findings
coherently, from a physics viewpoint. We also discuss some attempts of
interpretation that unify different results, highlighting the role that
statistical and soft condensed matter physics, and in particular classic and
more modern tools from the theory of polymers, plays in describing this system
of fundamental biological importance, and pointing to possible directions for
future investigation.


  We analyze the distribution of RNA secondary structures given by the
Knudsen-Hein stochastic context-free grammar used in the prediction program
Pfold. We prove that the distribution of base pairs, helices and various types
of loops in RNA secondary structures in this probabilistic model is
asymptotically Gaussian, for a generic choice of the grammar probabilities. Our
proofs are based on singularity analysis of probability generating functions.
Finally, we use our results to discuss how this model reflects the properties
of some known ribosomal secondary structures.


  Previous human foraging experiments have shown that human groups routinely
undermatch environmental resources much like other animal species. In this
experiment, we test whether humans also selectively rely on others as
information sources when the environmental state is uncertain, and we also test
whether overt signals of other foragers' success influences group matching
behavior and group adaptation to a changing environment. The results show
evidence of reliance on social information in specific conditions, but
participants were primarily influenced by their individual assessments of food
location rather than the success of other foragers.


  We consider population dynamics on a network of patches, each of which has a
the same local dynamics, with different population scales (carrying
capacities). It is reasonable to assume that if the patches are coupled by very
fast migration the whole system will look like an individual patch with a large
effective carrying capacity. This is called a "well-mixed" system. We show
that, in general, it is not true that the well-mixed system has the same
dynamics as each local patch. Different global dynamics can emerge from
coupling, and usually must be figured out for each individual case. We give a
general condition which must be satisfied for well-mixed systems to have the
same dynamics as the constituent patches.


  We analyze a class of network motifs in which a short, two-node positive
feed- back motif is inserted in a three-node negative feedback loop. We
demonstrate that such networks can undergo a bifurcation to a state where a
stable fixed point and a stable limit cycle coexist. At the bifurcation point
the period of the oscillations diverges. Further, intrinsic noise can make the
system switch between oscillatory state and the stationary state spontaneously.
We find that this switching also occurs in previous models of circadian clocks
that use this combination of positive and negative feedback. Our results
suggest that real- life circadian systems may need specific regulation to
prevent or minimize such switching events.


  We study the stochastic dynamics of evolutionary games, and focus on the
so-called `stochastic slowdown' effect, previously observed in (Altrock et. al,
2010) for simple evolutionary dynamics. Slowdown here refers to the fact that a
beneficial mutation may take longer to fixate than a neutral one. More
precisely, the fixation time conditioned on the mutant taking over can show a
maximum at intermediate selection strength. We show that this phenomenon is
present in the prisoner's dilemma, and also discuss counterintuitive slowdown
and speedup in coexistence games. In order to establish the microscopic origins
of these phenomena, we calculate the average sojourn times. This allows us to
identify the transient states which contribute most to the slowdown effect, and
enables us to provide an understanding of slowdown in the takeover of a small
group of cooperators by defectors: Defection spreads quickly initially, but the
final steps to takeover can be delayed substantially. The analysis of
coexistence games reveals even more intricate behavior. In small populations,
the conditional average fixation time can show multiple extrema as a function
of the selection strength, e.g., slowdown, speedup, and slowdown again. We
classify two-player games with respect to the possibility to observe
non-monotonic behavior of the conditional average fixation time as a function
of selection strength.


  Although empirical and theoretical studies affirm that punishment can elevate
collaborative efforts, its emergence and stability remain elusive. By
peer-punishment the sanctioning is something an individual elects to do
depending on the strategies in its neighborhood. The consequences of
unsustainable efforts are therefore local. By pool-punishment, on the other
hand, where resources for sanctioning are committed in advance and at large,
the notion of sustainability has greater significance. In a population with
free-riders, punishers must be strong in numbers to keep the "punishment pool"
from emptying. Failure to do so renders the concept of institutionalized
sanctioning futile. We show that pool-punishment in structured populations is
sustainable, but only if second-order free-riders are sanctioned as well, and
to a such degree that they cannot prevail. A discontinuous phase transition
leads to an outbreak of sustainability when punishers subvert second-order
free-riders in the competition against defectors.


  Evolutionary game dynamics in finite populations provides a new framework to
understand the selection of traits with frequency-dependent fitness. Recently,
a simple but fundamental law of evolutionary dynamics, which we call {\sigma}
law, describes how to determine the selection between two competing strategies:
in most evolutionary processes with two strategies, A and B, strategy A is
favored over B in weak selection if and only if {\sigma}R + S > T + {\sigma}P.
This relationship holds for a wide variety of structured populations with
mutation rate and weak selection under certain assumptions. In this paper, we
propose a model of games based on a community-structured population and revisit
this law under the Moran process. By calculating the average payoffs of A and B
individuals with the method of effective sojourn time, we find that {\sigma}
features not only the structured population characteristics but also the
reaction rate between individuals. That's to say, an interaction between two
individuals are not uniform, and we can take {\sigma} as a reaction rate
between any two individuals with the same strategy. We verify this viewpoint by
the modified replicator equation with non-uniform interaction rates in a
simplified version of the prisoner's dilemma game (PDG).


  Recent genome-wide measurements of binding preferences of ~200 transcription
regulators in the vicinity of transcription start sites in yeast, have provided
a unique insight into the cis- regulatory code of a eukaryotic genome (Venters
et al., Mol. Cell 41, 480 (2011)). Here, we show that nonspecific transcription
factor (TF)-DNA binding significantly influences binding preferences of the
majority of transcription regulators in promoter regions of the yeast genome.
We show that promoters of SAGA-dominated and TFIID-dominated genes can be
statistically distinguished based on the landscape of nonspecific protein-DNA
binding free energy. In particular, we predict that promoters of SAGA-dominated
genes possess wider regions of reduced free energy compared to promoters of
TFIID-dominated genes. We also show that specific and nonspecific TF-DNA
binding are functionally linked and cooperatively influence gene expression in
yeast. Our results suggest that nonspecific TF-DNA binding is intrinsically
encoded into the yeast genome, and it may play a more important role in
transcriptional regulation than previously thought.


  Multiple-type branching processes that model the spread of infectious
diseases are investigated. In these stochastic processes, the disease goes
through multiple stages before it eventually disappears. We mostly focus on the
critical multistage Susceptible-Infected-Recovered (SIR) infection process. In
the infinite population limit, we compute the outbreak size distributions and
show that asymptotic results apply to more general multiple-type critical
branching processes. Finally using heuristic arguments and simulations we
establish scaling laws for a multistage SIR model in a finite population.


  Intra-cellular biochemical reactions exhibit a rich dynamical phenomenology
which cannot be explained within the framework of mean-field rate equations and
additive noise. Here, we show that the presence of metastable states and
radically different timescales are general features of a broad class of
autocatalyic reaction networks, and moreover, that this fact may be exploited
to gain analytical results. The latter point is demonstrated by a treatment of
the paradigmatic Togashi-Kaneko reaction, which has resisted theoretical
analysis for the last decade.


  The multispecies coalescent model describes the generation of gene trees from
a rooted metric species tree, and thus provides a framework for the inference
of species trees from sampled gene trees. We prove that the STAR method of Liu
et al., and generalizations of it, are statistically consistent methods of
topological species tree inference under this model. We discuss the impact of
gene tree sampling schemes for species tree inference using generalized STAR
methods, and reinterpret the original STAR as a consensus method based on
clades.


  Background: How to extract useful information from complex biological
networks is a major goal in many fields, especially in genomics and proteomics.
We have shown in several works that iterative hierarchical clustering, as
implemented in the UVCluster program, is a powerful tool to analyze many of
those networks. However, the amount of computation time required to perform
UVCluster analyses imposed significant limitations to its use.
  Methodology/Principal Findings: We describe the suite Jerarca, designed to
efficiently convert networks of interacting units into dendrograms by means of
iterative hierarchical clustering. Jerarca is divided into three main sections.
First, weighted distances among units are computed using up to three different
approaches: a more efficient version of UVCluster and two new, related
algorithms called RCluster and SCluster. Second, Jerarca builds dendrograms
based on those distances, using well-known phylogenetic algorithms, such as
UPGMA or Neighbor-Joining. Finally, Jerarca provides optimal partitions of the
trees using statistical criteria based on the distribution of intra- and
intercluster connections. Outputs compatible with the phylogenetic software
MEGA and the Cytoscape package are generated, allowing the results to be easily
visualized.
  Conclusions/Significance: The four main advantages of Jerarca in respect to
UVCluster are: 1) Improved speed of a novel UVCluster algorithm; 2) Additional,
alternative strategies to perform iterative hierarchical clustering; 3)
Automatic evaluation of the hierarchical trees to obtain optimal partitions;
and, 4) Outputs compatible with popular software such as MEGA and Cytoscape.


  In their 2008 and 2009 papers, Sumner and colleagues introduced the
"squangles" - a small set of Markov invariants for phylogenetic quartets. The
squangles are consistent with the general Markov model (GM) and can be used to
infer quartets without the need to explicitly estimate all parameters. As GM is
inhomogeneous and hence non-stationary, the squangles are expected to perform
well compared to standard approaches when there are changes in base-composition
amongst species. However, GM includes the IID assumption, so the squangles
should be confounded by data generated with invariant sites or with
rate-variation across sites. Here we implement the squangles in a least-squares
setting that returns quartets weighted by either confidence or internal edge
lengths; and use these as input into a variety of quartet-based supertree
methods. For the first time, we quantitatively investigate the robustness of
the squangles to the breaking of IID assumptions on both simulated and real
data sets; and we suggest a modification that improves the performance of the
squangles in the presence of invariant sites. Our conclusion is that the
squangles provide a novel tool for phylogenetic estimation that is
complementary to methods that explicitly account for rate-variation across
sites, but rely on homogeneous - and hence stationary - models.


  Several ionic liquids (ILs) are known to revert aggregation processes and
improve the in vitro refolding of denatured proteins. In this paper the
capacity of a particular class of ammonium based ILs to act as refolding
enhancers was tested using lysozyme as a model protein. Raman spectra of ILs
treated fibrillar lysozyme as well as lysozyme in its native and fibrillar
conformations were collected and carefully analyzed to characterize the
refolding extent under the effect of the IL interaction. Results obtained
confirm and largely extend the earlier knowledge on this class of protic ILs
and indicate Ethyl Ammonium Nitrate (EAN) as the most promising additive for
protein refolding. The experiment provides also the demonstration of the high
potentiality of Raman spectroscopy as a comprehensive diagnostic tool in this
field.


  Evolutionary experiments with microbes are a powerful tool to study mutations
and natural selection. These experiments, however, are often limited to the
well-mixed environments of a test tube or a chemostat. Since spatial
organization can significantly affect evolutionary dynamics, the need is
growing for evolutionary experiments in spatially structured environments. The
surface of a Petri dish provides such an environment, but a more detailed
understanding of microbial growth on Petri dishes is necessary to interpret
such experiments. We formulate a simple deterministic reaction-diffusion model,
which successfully predicts the spatial patterns created by two competing
species during colony expansion. We also derive the shape of these patterns
analytically without relying on microscopic details of the model. In
particular, we find that the relative fitness of two microbial strains can be
estimated from the logarithmic spirals created by selective sweeps. The theory
is tested with strains of the budding yeast Saccharomyces cerevisiae, for
spatial competitions with different initial conditions and for a range of
relative fitnesses. The reaction-diffusion model also connects the microscopic
parameters like growth rates and diffusion constants with macroscopic spatial
patterns and predicts the relationship between fitness in liquid cultures and
on Petri dishes, which we confirmed experimentally. Spatial sector patterns
therefore provide an alternative fitness assay to the commonly used liquid
culture fitness assays.


  Sexually reproducing populations with small number of individuals may go
extinct by stochastic fluctuations in sex determination, causing all their
members to become male or female in a generation. In this work we calculate the
time to extinction of isolated populations with fixed number $N$ of individuals
that are updated according to the Moran birth and death process. At each time
step, one individual is randomly selected and replaced by its offspring
resulting from mating with another individual of opposite sex; the offspring
can be male or female with equal probability. A set of $N$ time steps is called
a generation, the average time it takes for the entire population to be
replaced. The number k of females fluctuates in time, similarly to a random
walk, and extinction, which is the only asymptotic possibility, occurs when k=0
or k=N. We show that it takes only one generation for an arbitrary initial
distribution of males and females to approach the binomial distribution. This
distribution, however, is unstable and the population eventually goes extinct
in 2^N/N generations. We also discuss the robustness of these results against
bias in the determination of the sex of the offspring, a characteristic
promoted by infection by the bacteria Wolbachia in some arthropod species or by
temperature in reptiles.


  Molecular clock (MC) is a central concept of molecular evolution according to
which each gene evolves at a characteristic, near constant rate. Numerous
evolutionary studies have demonstrated the validity of MC but also have shown
that MC is substantially overdispersed, i.e. lineage-specific deviations of the
evolutionary rate of the given gene from the clock greatly exceed the
expectation from the sampling error. A fundamental observation of comparative
genomics that appears to complement the MC is that the distribution of
evolution rates across orthologous genes in pairs of related genomes remains
virtually unchanged throughout the evolution of life, from bacteria to mammals.
The conservation of this distribution implies that the relative evolution rates
of all genes remain nearly constant, or in other words, that evolutionary rates
of different genes are strongly correlated within each evolving genome. We
hypothesized that this correlation is not a simple consequence of MC but could
be better explained by a model we dubbed Universal PaceMaker (UPM) of genome
evolution. The UPM model posits that the rate of evolution changes
synchronously across genome-wide sets of genes in all evolving lineages. We
sought to differentiate between the MC and UPM models by fitting thousands of
phylogenetic trees for bacterial and archaeal genes to supertrees that reflect
the dominant trend of vertical descent in the evolution of archaea and bacteria
and that were constrained according to the two models. The goodness of fit for
the UPM model was better than the fit for the MC model, with overwhelming
statistical significance. These results reveal a universal pacemaker of genome
evolution that could have been in operation throughout the history of life.


  Boolean network models of molecular regulatory networks have been used
successfully in computational systems biology. The Boolean functions that
appear in published models tend to have special properties, in particular the
property of being nested canalizing, a concept inspired by the concept of
canalization in evolutionary biology. It has been shown that networks comprised
of nested canalizing functions have dynamic properties that make them suitable
for modeling molecular regulatory networks, namely a small number of (large)
attractors, as well as relatively short limit cycles.
  This paper contains a detailed analysis of this class of functions, based on
a novel normal form as polynomial functions over the Boolean field. The concept
of layer is introduced that stratifies variables into different classes
depending on their level of dominance. Using this layer concept a closed form
formula is derived for the number of nested canalizing functions with a given
number of variables. Additional metrics considered include Hamming weight, the
activity number of any variable, and the average sensitivity of the function.
It is also shown that the average sensitivity of any nested canalizing function
is between 0 and 2. This provides a rationale for why nested canalizing
functions are stable, since a random Boolean function in n variables has
average sensitivity n/2. The paper also contains experimental evidence that the
layer number is an important factor in network stability.


  Contacts between individuals serve as pathways where infections may
propagate. These contact patterns can be represented by network structures.
Static structures have been the common modeling paradigm but recent results
suggest that temporal structures play different roles to regulate the spread of
infections or infection-like dynamics. On temporal networks a vertex is active
only at certain moments and inactive otherwise such that a contact is not
continuously available. In several empirical networks, the time between two
consecutive vertex-activation events typically follows heterogeneous activity
(e.g. bursts). In this chapter, we present a simple and intuitive stochastic
model of a temporal network and investigate how epidemics co-evolves with the
temporal structures, focusing on the growth dynamics of the epidemics. The
model assumes no underlying topological structure and is only constrained by
the time between two consecutive events of vertex activation. The main
observation is that the speed of the infection spread is different in case of
heterogeneous and homogeneous temporal patterns but the differences depend on
the stage of the epidemics. In comparison to the homogeneous scenario, the
power law case results in a faster growth in the beginning but turns out to be
slower after a certain time, taking several time steps to reach the whole
network.


  Migrating cells possess intracellular gradients of Rho GTPases, but it is
unknown whether these shallow gradients themselves can induce motility. Here we
describe a new method to present cells with induced linear gradients of active,
endogenous Rac without receptor activation. Gradients as low as 15% were
sufficient to not only trigger cell migration up the synthetic gradient, but
also to induce both cell polarization and repolarization. Response kinetics
were inversely proportional to Rac gradient values, in agreement with a new
mathematical model, suggesting a role for natural input gradient amplification
upstream of Rac. Increases in Rac levels beyond a well-defined threshold
dramatically augmented polarization and decreased sensitivity to the gradient
value. The threshold was governed by initial cell polarity and PI3K activity,
supporting a role for both in defining responsiveness to natural or synthetic
Rac activation. Our methodology suggests a general way to investigate processes
regulated by intracellular signaling gradients.


  The copy number of any protein fluctuates among cells in a population;
characterizing and understanding these fluctuations is a fundamental problem in
biophysics. We show here that protein distributions measured under a broad
range of biological realizations collapse to a single non-Gaussian curve under
scaling by the first two moments. Moreover in all experiments the variance is
found to depend quadratically on the mean, showing that a single degree of
freedom determines the entire distribution. Our results imply that protein
fluctuations do not reflect any specific molecular or cellular mechanism, and
suggest that some buffering process masks these details and induces
universality.


  We explore a class of hybrid (piecewise deterministic) systems characterized
by a large number of individuals inhabiting an environment whose state is
described by a set of continuous variables. We use analytical and numerical
methods from nonequilibrium statistical mechanics to study the influence that
intrinsic noise has on the qualitative behavior of the system. We discuss the
application of these concepts to the case of semiarid ecosystems. Using a
system-size expansion we calculate the power spectrum of the fluctuations in
the system. This predicts the existence of noise-induced oscillations.


  Length-regulation of microtubules (MTs) is essential for many cellular
processes. Molecular motors like kinesin 8, which move along MTs and also act
as depolymerases, are known as key players in MT dynamics. However, the
regulatory mechanisms of length control remain elusive. Here, we investigate a
stochastic model accounting for the interplay between polymerization kinetics
and motor-induced depolymerization. We determine the dependence of MT length
and variance on rate constants and motor concentration. Moreover, our analyses
reveal how collective phenomena lead to a well-defined MT length.


  The aim of this paper is to contribute to the basic understanding of neuronal
polarization mechanisms by developing and studying a reaction-diffusion model
for protein activation and inactivation. In particular we focus on a feedback
loop between PI3 kinase and certain GTPases, and study its behaviour in
dependence of neurite lengths.
  We find that if an ultrasensitive activation is included, the model can
produce polarization at a critical length as observed in experiments. Symmetry
breaking to polarization in the longer neurite is found only if active
transport of a substance, in our case active PI3 kinase, is included into the
model.


  Protein synthesis is one of the most fundamental biological processes, which
consumes a significant amount of cellular resources. Despite existence of
multiple mathematical models of translation, varying in the level of
mechanistical details, surprisingly, there is no basic and simple chemical
kinetic model of this process, derived directly from the detailed kinetic
model. One of the reasons for this is that the translation process is
characterized by indefinite number of states, thanks to existence of polysomes.
We bypass this difficulty by applying a trick consisting in lumping multiple
states of translated mRNA into few dynamical variables and by introducing a
variable describing the pool of translating ribosomes. The simplest model can
be solved analytically under some assumptions. The basic and simple model can
be extended, if necessary, to take into account various phenomena such as the
interaction between translating ribosomes, limited amount of ribosomal units or
regulation of translation by microRNA. The model can be used as a building
block (translation module) for more complex models of cellular processes. We
demonstrate the utility of the model in two examples. First, we determine the
critical parameters of the single protein synthesis for the case when the
ribosomal units are abundant. Second, we demonstrate intrinsic bi-stability in
the dynamics of the ribosomal protein turnover and predict that a minimal
number of ribosomes should pre-exists in a living cell to sustain its protein
synthesis machinery, even in the absence of proliferation.


  In the present paper we discuss the clustering procedure in the case where
instead of a single metric we have a family of metrics. In this case we can
obtain a partially ordered graph of clusters which is not necessarily a tree.
We discuss a structure of a hypergraph above this graph. We propose two
definitions of dimension for hyperedges of this hypergraph and show that for
the multidimensional p-adic case both dimensions are reduced to the number of
p-adic parameters.
  We discuss the application of the hypergraph clustering procedure to the
construction of phylogenetic graphs in biology. In this case the dimension of a
hyperedge will describe the number of sources of genetic diversity.


  In this review, we discuss modularity and hierarchy in biological systems. We
review examples from protein structure, genetics, and biological networks of
modular partitioning of the geometry of biological space. We review theories to
explain modular organization of biology, with a focus on explaining how biology
may spontaneously organize to a structured form. That is, we seek to explain
how biology nucleated from among the many possibilities in chemistry. The
emergence of modular organization of biological structure will be described as
a symmetry-breaking phase transition, with modularity as the order parameter.
Experimental support for this description will be reviewed. Examples will be
presented from pathogen structure, metabolic networks, gene networks, and
protein-protein interaction networks. Additional examples will be presented
from ecological food networks, developmental pathways, physiology, and social
networks.


  The zebrafish (\emph{Danio rerio}) is one of the model animals for study of
immunology because the dynamics in the adaptive immune system of zebrafish are
similar to that in higher animals. In this work, we built a multi-scale model
to simulate the dynamics of B cells in the primary and secondary immune
responses of zebrafish. We use this model to explain the reported correlation
between VDJ usage of B cell repertoires in individual zebrafish. We use a delay
ordinary differential equation (ODE) system to model the immune responses in
the 6-month lifespan of a zebrafish. This mean field theory gives the number of
high affinity B cells as a function of time during an infection. The sequences
of those B cells are then taken from a distribution calculated by a
"microscopic" random energy model. This generalized $NK$ model shows that
mature B cells specific to one antigen largely possess a single VDJ
recombination. The model allows first-principles calculation of the
probability, $p$, that two zebrafish responding to the same antigen will select
the same VDJ recombination. This probability $p$ increases with the B cell
population size and the B cell selection intensity. The probability $p$
decreases with the B cell hypermutation rate. The multi-scale model predicts
correlations in the immune system of the zebrafish that are highly similar to
that from experiment.


  Many viruses evolve rapidly. For example, hemagglutinin of the H3N2 influenza
A virus evolves to escape antibody binding. This evolution of the H3N2 virus
means that people who have previously been exposed to an influenza strain may
be infected by a newly emerged virus. In this paper, we use Shannon entropy and
relative entropy to measure the diversity and selection pressure by antibody in
each amino acid site of H3 hemagglutinin between the 1992-1993 season and the
2009-2010 season. Shannon entropy and relative entropy are two independent
state variables that we use to characterize H3N2 evolution. The entropy method
estimates future H3N2 evolution and migration using currently available H3
hemagglutinin sequences. First, we show that the rate of evolution increases
with the virus diversity in the current season. The Shannon entropy of the
sequence in the current season predicts relative entropy between sequences in
the current season and those in the next season. Second, a global migration
pattern of H3N2 is assembled by comparing the relative entropy flows of
sequences sampled in China, Japan, the USA, and Europe. We verify this entropy
method by describing two aspects of historical H3N2 evolution. First, we
identify 54 amino acid sites in hemagglutinin that have evolved in the past to
evade the immune system. Second, the entropy method shows that epitopes A and B
in the top of hemagglutinin evolve most vigorously to escape antibody binding.
Our work provides a novel entropy-based method to predict and quantify future
H3N2 evolution and to describe the evolutionary history of H3N2.


  We consider the evolution of large but finite populations on arbitrary
fitness landscapes. We describe the evolutionary process by a Markov, Moran
process. We show that to $\mathcal O(1/N)$, the time-averaged fitness is lower
for the finite population than it is for the infinite population. We also show
that fluctuations in the number of individuals for a given genotype can be
proportional to a power of the inverse of the mutation rate. Finally, we show
that the probability for the system to take a given path through the fitness
landscape can be non-monotonic in system size.


  Influenza virus evolves to escape from immune system antibodies that bind to
it. We used free energy calculations with Einstein crystals as reference states
to calculate the difference of antibody binding free energy ($\Delta\Delta G$)
induced by amino acid substitution at each position in epitope B of the H3N2
influenza hemagglutinin, the key target for antibody. A substitution with
positive $\Delta\Delta G$ value decreases the antibody binding constant. On
average an uncharged to charged amino acid substitution generates the highest
$\Delta\Delta G$ values. Also on average, substitutions between small amino
acids generate $\Delta\Delta G$ values near to zero. The 21 sites in epitope B
have varying expected free energy differences for a random substitution.
Historical amino acid substitutions in epitope B for the A/Aichi/2/1968 strain
of influenza A show that most fixed and temporarily circulating substitutions
generate positive $\Delta\Delta G$ values. We propose that the observed pattern
of H3N2 virus evolution is affected by the free energy landscape, the mapping
from the free energy landscape to virus fitness landscape, and random genetic
drift of the virus. Monte Carlo simulations of virus evolution are presented to
support this view.


  The extent of coupling between the folding of a protein and its binding to a
substrate varies from protein to protein. Some proteins have highly structured
native states in solution, while others are natively disordered and only fold
fully upon binding. In this Letter, we use Monte Carlo simulations to
investigate how disordered polymer chains grafted around a binding site affect
the folding and binding of three model proteins. The protein that approaches
the substrate fully folded is more hindered during the binding process than
those whose folding and binding are cooperative. The polymer chains act as
localized crowding agents and can select correctly folded and bound
configurations in favor of non-specifically adsorbed states. The free energy
change for forming all intra-protein and protein-substrate contacts can depend
non-monotonically on the polymer length.


  Force field and first principles molecular dynamics simulations on complexes
of pig liver esterase (pig liver isoenzymes and a mutant) and selected
substrates (1-phenyl-1-ethyl acetate, 1- phenyl-2-butylacetate,
proline-{\beta}-naphthylamide and methyl butyrate) are presented. By restrained
force field simulations the access of the substrate to the hidden active site
was probed. For a few substrates spontaneous access to the active site via a
well defined entrance channel was found. The structure of the tetrahedral
intermediate was simulated for several substrates and our previous assignment
of GLU 452 instead of GLU 336 was confirmed. It was shown that the active site
readily adapts to the embedded substrate involving a varying number of
hydrophobic residues in the neighborhood. This puts into question key-lock
models for enantioselectivity. Ab initio molecular dynamics showed that the
structures we found for the tetrahedral intermediate in force field simulations
are consistent with the presumed mechanism of ester cleavage. Product release
from the active site as final step of the enzymatic reaction revealed to be
very slow and took already more than 20ns for the smallest product, methanol.


  Evolutionary branching is analysed in a stochastic, individual-based
population model under mutation and selection. In such models, the common
assumption is that individual reproduction and life career are characterised by
values of a trait, and also by population sizes, and that mutations lead to
small changes in trait value. Then, traditionally, the evolutionary dynamics is
studied in the limit of vanishing mutational step sizes. In the present
approach, small but non-negligible mutational steps are considered. By means of
theoretical analysis in the limit of infinitely large populations, as well as
computer simulations, we demonstrate how discrete mutational steps affect the
patterns of evolutionary branching. We also argue that the average time to the
first branching depends in a sensitive way on both mutational step size and
population size.


  H1N1 influenza causes substantial seasonal illness and was the subtype of the
2009 influenza pandemic. Precise measures of antigenic distance between the
vaccine and circulating virus strains help researchers design influenza
vaccines with high vaccine effectiveness. We here introduce a sequence-based
method to predict vaccine effectiveness in humans. Historical epidemiological
data show that this sequence-based method is as predictive of vaccine
effectiveness as hemagglutination inhibition (HI) assay data from ferret animal
model studies. Interestingly, the expected vaccine effectiveness is greater
against H1N1 than H3N2, suggesting a stronger immune response against H1N1 than
H3N2. The evolution rate of hemagglutinin in H1N1 is also shown to be greater
than that in H3N2, presumably due to greater immune selection pressure.


  Catastrophic regime shifts in complex natural systems may be averted through
advanced detection. Recent work has provided a proof-of-principle that many
systems approaching a catastrophic transition may be identified through the
lens of early warning indicators such as rising variance or increased return
times. Despite widespread appreciation of the difficulties and uncertainty
involved in such forecasts, proposed methods hardly ever characterize their
expected error rates. Without the benefits of replicates, controls, or
hindsight, applications of these approaches must quantify how reliable
different indicators are in avoiding false alarms, and how sensitive they are
to missing subtle warning signs. We propose a model based approach in order to
quantify this trade-off between reliability and sensitivity and allow
comparisons between different indicators. We show these error rates can be
quite severe for common indicators even under favorable assumptions, and also
illustrate how a model-based indicator can improve this performance. We
demonstrate how the performance of an early warning indicator varies in
different data sets, and suggest that uncertainty quantification become a more
central part of early warning predictions.


  Influenza has been circulating in the human population and has caused three
pandemics in the last century (1918 H1N1, 1957 H2N2, 1968 H3N2). The 2009
A(H1N1) was classified by the World Health Organization (WHO) as the fourth
pandemic. Influenza has a high evolution rate, which makes vaccine design
challenging. We here consider an approach for early detection of new dominant
strains. By clustering the 2009 A(H1N1) sequence data, we found two main
clusters. We then define a metric to detect the emergence of dominant strains.
We show on historical H3N2 data that this method is able to identify a cluster
around an incipient dominant strain before it becomes dominant. For example,
for H3N2 as of March 30, 2009, the method detects the cluster for the new
A/British Columbia/RV1222/2009 strain. This strain detection tool would appear
to be useful for annual influenza vaccine selection.


  This is supplementary information for arXiv:1204.4896, which also appeared in
Physical Biology 9, 026008 (2012). Main text abstract: Evolutionary experiments
with microbes are a powerful tool to study mutations and natural selection.
These experiments, however, are often limited to the well-mixed environments of
a test tube or a chemostat. Since spatial organization can significantly affect
evolutionary dynamics, the need is growing for evolutionary experiments in
spatially structured environments. The surface of a Petri dish provides such an
environment, but a more detailed understanding of microbial growth on Petri
dishes is necessary to interpret such experiments. We formulate a simple
deterministic reaction-diffusion model, which successfully predicts the spatial
patterns created by two competing species during colony expansion. We also
derive the shape of these patterns analytically without relying on microscopic
details of the model. In particular, we find that the relative fitness of two
microbial strains can be estimated from the logarithmic spirals created by
selective sweeps. The theory is tested with strains of the budding yeast
Saccharomyces cerevisiae, for spatial competitions with different initial
conditions and for a range of relative fitnesses. The reaction-diffusion model
also connects the microscopic parameters like growth rates and diffusion
constants with macroscopic spatial patterns and predicts the relationship
between fitness in liquid cultures and on Petri dishes, which we confirmed
experimentally. Spatial sector patterns therefore provide an alternative
fitness assay to the commonly used liquid culture fitness assays.


  In this paper we revisit and adapt to viral evolution an approach based on
the theory of branching process advanced by Demetrius, Schuster and Sigmund
("Polynucleotide evolution and branching processes", Bull. Math. Biol. 46
(1985) 239-262), in their study of polynucleotide evolution. By taking into
account beneficial effects we obtain a non-trivial multivariate generalization
of their single-type branching process model. Perturbative techniques allows us
to obtain analytical asymptotic expressions for the main global parameters of
the model which lead to the following rigorous results: (i) a new criterion for
"no sure extinction", (ii) a generalization and proof, for this particular
class of models, of the lethal mutagenesis criterion proposed by Bull,
Sanju\'an and Wilke ("Theory of lethal mutagenesis for viruses", J. Virology 18
(2007) 2930-2939), (iii) a new proposal for the notion of relaxation time with
a quantitative prescription for its evaluation, (iv) the quantitative
description of the evolution of the expected values in in four distinct
"stages": extinction threshold, lethal mutagenesis, stationary "equilibrium"
and transient. Finally, based on these quantitative results we are able to draw
some qualitative conclusions.


  Topological landscape is introduced for networks with functions defined on
the nodes. By extending the notion of gradient flows to the network setting,
critical nodes of different indices are defined. This leads to a concise and
hierarchical representation of the network. Persistent homology from
computational topology is used to design efficient algorithms for performing
such analysis. Applications to some examples in social and biological networks
are demonstrated, which show that critical nodes carry important information
about structures and dynamics of such networks.


  The network approach became a widely used tool to understand the behaviour of
complex systems in the last decade. We start from a short description of
structural rigidity theory. A detailed account on the combinatorial rigidity
analysis of protein structures, as well as local flexibility measures of
proteins and their applications in explaining allostery and thermostability is
given. We also briefly discuss the network aspects of cytoskeletal tensegrity.
Finally, we show the importance of the balance between functional flexibility
and rigidity in protein-protein interaction, metabolic, gene regulatory and
neuronal networks. Our summary raises the possibility that the concepts of
flexibility and rigidity can be generalized to all networks.


  The position of a reaction front, propagating into a metastable state,
fluctuates because of the shot noise of reactions and diffusion. A recent
theory [B. Meerson, P.V. Sasorov, and Y. Kaplan, Phys. Rev. E 84, 011147
(2011)] gave a closed analytic expression for the front diffusion coefficient
in the weak noise limit. Here we test this theory in stochastic simulations
involving reacting and diffusing particles on a one-dimensional lattice. We
also investigate a small noise-induced systematic shift of the front velocity
compared to the prediction from the spatially continuous deterministic
reaction-diffusion equation.


  As a generalization of deterministic, nonlinear conservative dynamical
systems, a notion of {\em canonical conservative dynamics} with respect to a
positive, differentiable stationary density $\rho(x)$ is introduced:
$\dot{x}=j(x)$ in which $\nabla\cdot\big(\rho(x)j(x)\big)=0$. Such systems have
a conserved "generalized free energy function" $F[u]$ $=$ $\int
u(x,t)\ln\big(u(x,t)/\rho(x)\big)dx$ in phase space with a density flow
$u(x,t)$ satisfying $\partial u_t =-\nabla\cdot(ju)$. Any general stochastic
diffusion process without detailed balance, in terms of its Fokker-Planck
equation, can be decomposed into a reversible diffusion process with detailed
balance and a canonical conservative dynamics. This decomposition can be
rigorously established in a function space with inner product defined as
$\langle\phi,\psi\rangle=\int\rho^{-1}(x)\phi(x)\psi(x)dx$. Furthermore, a law
for balancing $F[u]$ can be obtained: The non-positive $dF[u(x,t)]/dt$ $=$
$E_{in}(t)-e_p(t)$ where the "source" $E_{in}(t)\ge 0$ and the "sink"
$e_p(t)\ge 0$ are known as house-keeping heat and entropy production,
respectively. A reversible diffusion has $E_{in}(t)=0$. For a linear
(Ornstein-Uhlenbeck) diffusion process, our decomposition is equivalent to the
previous approaches developed by R. Graham and P. Ao, as well as the theory of
large deviations. In terms of two different formulations of time reversal for a
same stochastic process, the meanings of {\em dissipative} and {\em
conservative} stationary dynamics are discussed.


  Mining gene expression profiles has proven valuable for identifying
metagenes, defined as linear combinations of individual genes, serving as
surrogates of biological phenotypes. Typically, such metagenes are jointly
generated as the result of an optimization process for dimensionality
reduction. Here we present an unconstrained method for individually generating
metagenes that can point to the core of the underlying biological mechanisms.
We use an iterative process that starts from any seed gene and converges to one
of several precise attractor metagenes representing biomolecular events, such
as cell transdifferentiation or the presence of an amplicon. By analyzing six
rich gene expression datasets from three different cancer types, we identified
many such biomolecular events, some of which are present in all tested cancer
types. We focus on several such events including a stage-associated mesenchymal
transition and a grade-associated mitotic chromosomal instability.


  Background: The accumulation of deleterious mutations of a population
directly contributes to the fate as to how long the population would exist, a
process often described as Muller's ratchet with the absorbing phenomenon. The
key to understand this absorbing phenomenon is to characterize the decaying
time of the fittest class of the population. Adaptive landscape introduced by
Wright, a re-emerging powerful concept in systems biology, is used as a tool to
describe biological processes. To our knowledge, the dynamical behaviors for
Muller's ratchet over the full parameter regimes are not studied from the point
of the adaptive landscape. And the characterization of the absorbing phenomenon
is not yet quantitatively obtained without extraneous assumptions as well.
  Results: We describe the dynamical behavior of the population exposed to
Muller's ratchet in all parameters regimes by adaptive landscape. The adaptive
landscape has rich structures such as finite and infinite potential, real and
imaginary fixed points. We give the formula about the single click time with
finite and infinite potential. And we find the single click time increases with
selection rates and population size increasing, decreases with mutation rates
increasing. These results provide a new understanding of infinite potential. We
analytically demonstrate the adaptive and unadaptive states for the whole
parameters regimes. Interesting issues about the parameters regions with the
imaginary fixed points is demonstrated. This can help understand the biological
meaning of the critical points such as the intersection points of regimes. Most
importantly, we find that the absorbing phenomenon is characterized by the
adaptive landscape and the single click time without any extraneous
assumptions. These results suggest a graphical and quantitative framework to
study the absorbing phenomenon.


  The recently proposed close-packed motif for collagen is investigated using
first principles semi-empirical wave function theory and Kohn-Sham density
functional theory. Under these refinements the close-packed motif is shown to
be stable. For the case of the 7/2 motif a similar stability exists. The
electronic circular dichroism of the close-packed model has a significant
negative bias and a large signal. An interesting feature of the close-packed
structure is the existence of a central channel. Simulations show that, if
hydrogen atoms are placed in the cavity, a chain of molecular hydrogens is
formed suggesting a possible biological function for molecular hydrogen.


  Motivation
  The Burrows-Wheeler transform (BWT) is the foundation of many algorithms for
compression and indexing of text data, but the cost of computing the BWT of
very large string collections has prevented these techniques from being widely
applied to the large sets of sequences often encountered as the outcome of DNA
sequencing experiments. In previous work, we presented a novel algorithm that
allows the BWT of human genome scale data to be computed on very moderate
hardware, thus enabling us to investigate the BWT as a tool for the compression
of such datasets.
  Results
  We first used simulated reads to explore the relationship between the level
of compression and the error rate, the length of the reads and the level of
sampling of the underlying genome and compare choices of second-stage
compression algorithm.
  We demonstrate that compression may be greatly improved by a particular
reordering of the sequences in the collection and give a novel `implicit
sorting' strategy that enables these benefits to be realised without the
overhead of sorting the reads. With these techniques, a 45x coverage of real
human genome sequence data compresses losslessly to under 0.5 bits per base,
allowing the 135.3Gbp of sequence to fit into only 8.2Gbytes of space (trimming
a small proportion of low-quality bases from the reads improves the compression
still further).
  This is more than 4 times smaller than the size achieved by a standard
BWT-based compressor (bzip2) on the untrimmed reads, but an important further
advantage of our approach is that it facilitates the building of compressed
full text indexes such as the FM-index on large-scale DNA sequence collections.


  MicroRNA-mediated regulation of gene expression is characterised by some
distinctive features that set it apart from unregulated and transcription
factor-regulated gene expression. Recently, a mathematical model has been
proposed to describe the dynamics of post-transcriptional regulation by
microRNAs. The model explains the observations made in single cell experiments
quite well. In this paper, we introduce some additional features into the model
and consider two specific cases. In the first case, a non-cooperative positive
feedback loop is included in the transcriptional regulation of the target gene
expression. In the second case, a stochastic version of the original model is
considered in which there are random transitions between the inactive and
active expression states of the gene. In the first case we show that
bistability is possible in a parameter regime, due to the presence of a
non-linear protein decay term in the gene expression dynamics. In the second
case, we derive the conditions for obtaining stochastic binary gene expression.
We find that this type of gene expression is more favourable in the case of
regulation by microRNAs as compared to the case of unregulated gene expression.
The theoretical predictions relating to binary gene expression are
experimentally testable.


  Living cell signaling systems include multistep biochemical signaling
reaction cascades (BSCs) comprising modifications of molecular signaling
proteins. Substantial data on BSCs have been accumulated in the field of
molecular biology and the analysis of signaling systems requires qualitative
evaluation. However, quantification of the information and channel capacity of
BSCs has not been focused on from the perspective of information theory. In the
current study, we aimed to derive basic equations for describing the channel
capacity and information density of BSCs using the fluctuation theorem. From
the results, channel capacity and information density can be described using
the average entropy production rate when the signaling system is minimally
redundant. The channel capacity could actually be calculated for the
mitogen-activated protein kinase BSC when it was minimally redundant. This
quantitative method of examination is applicable to the quantitative analysis
of BSCs.


  While worldwide census data provide statistical evidence that firmly link the
population density with several indicators of social welfare, the precise
mechanisms underlying these observations are largely unknown. Here we study the
impact of population density on the evolution of public cooperation in
structured populations and find that the optimal density is uniquely related to
the percolation threshold of the host graph irrespective of its topological
details. We explain our observations by showing that spatial reciprocity peaks
in the vicinity of the percolation threshold, when the emergence of a giant
cooperative cluster is hindered neither by vacancy nor by invading defectors,
thus discovering an intuitive yet universal law that links the population
density with social prosperity.


  This paper presents new results from a detailed study of the structure of
autocatalytic sets. We show how autocatalytic sets can be decomposed into
smaller autocatalytic subsets, and how these subsets can be identified and
classified. We then argue how this has important consequences for the
evolvability, enablement, and emergence of autocatalytic sets. We end with some
speculation on how all this might lead to a generalized theory of autocatalytic
sets, which could possibly be applied to entire ecologies or even economies.