  Electrode-immobilized glucose-6-phosphate dehydrogenase is used to catalyze
an enzymatic reaction which carries out the AND logic gate. This logic function
is considered here in the context of biocatalytic processes utilized for the
biocomputing applications for "digital" (threshold) sensing/actuation. We
outline the response functions desirable for such applications and report the
first experimental realization of a sigmoid-shape response in one of the
inputs. A kinetic model is developed and utilized to evaluate the extent to
which the experimentally realized gate is close to optimal.


  FADD (Fas-associated death domain) is the main death receptor adaptor
molecule that transmits apoptotic signal. Recently, FADD protein was shown to
be expressed both in the cytoplasm and nucleus of in vitro cell lines. In
contrast to the cytoplasmic FADD, the nuclear FADD was shown to protect cells
from apoptosis. However, in vivo subcellular localization of FADD was still
unknown. Here, we demonstrated that FADD protein was expressed in both
cytoplasmic and nuclear compartment in ex vivo thyroid cells demonstrating that
nuclear sublocalization of FADD protein was a relevant phenomenon occurring in
vivo. Moreover, we showed that in the nucleus of untransformed thyroid cells
FADD localized mainly on euchromatin. We confirmed the nuclear localization of
FADD in ex vivo liver and showed that in this organ FADD and MBD4 interact
together. These results demonstrate that FADD is physiologically expressed in
the nucleus of cells in at least two mouse organs. This particular localization
opens new possible role of FADD in vivo either asan inhibitor of cell death, or
as a transcription factor, or as a molecular link between apoptosis and genome
surveillance.


  Two-dimensional semiflexible polymer rings are studied both by imaging
circular DNA adsorbed on a mica surface and by Monte Carlo simulations of
phantom polymers as well as of polymers with finite thickness. Comparison of
size and shape of the different models over the full range of flexibilities
shows that excluded volume caused by finite thickness induces an anisotropic
increase of the main axes of the conformations, a change of shape, accomplished
by an enhanced correlation along the contour.


  By assuming that not only counter-ions but DNA molecules as well are
thermally distributed according to a Boltzmann law, we propose a modified
Poisson-Boltzmann equation at the classical level as starting point to compute
the effects of quantum fluctuations of the electric field on the interaction
among DNA-cation complexes. The latter are modeled here as infinite
one-dimensional wires ($\delta$-functions). Our goal is to single out such
quantum-vacuum-driven interaction from the counterion-induced and water-related
interactions. We obtain a universal, frustration-free Casimir-like (codimension
2) interaction that extensive numerical analysis show to be a good candidate to
explain the formation and stability of DNA aggregates. Such Casimir energy is
computed for a variety of configurations of up to 19 DNA strands in a hexagonal
array. It is found to be strongly many-body.


  Application of quantum principles to living cells requires a new
approximation of the full quantum mechanical description of intracellular
dynamics. We discuss what principal elements any such good approximation should
contain. As one such element, the notion of "Catalytic force" Cf is introduced.
Cf is the effect of the molecular target of catalysis on the catalytic
microenvironment that adjusts the microenvironment towards a state that
facilitates the catalytic act. This phenomenon is experimentally testable and
has an intriguing implication for biological organization and evolution, as it
amounts to "optimization without natural selection of replicators". Unlike the
statistical-mechanical approaches to self-organization, the Cf principle does
not encounter the problem of "tradeoff between stability and complexity" at the
level of individual cell. Physically, the Cf is considered as a harmonic-like
force of reaction, which keeps the state of the cell close to the ground state,
defined here as a state where enzymatic acts work most efficiently. Ground
state is subject to unitary evolution, and serves as a starting point in a
general strategy of quantum description of intracellular processes, termed here
"Euclidean approach". The next step of this strategy is transition from the
description of ground state to that one of growing state, and we suggest how it
can be accomplished using arguments from the fluctuation-dissipation theorem.
Finally, given that the most reliable and informative observable of an
individual cell is the sequence of its genome, we propose that the
non-classical correlations between individual molecular events at the single
cell level could be easiest to detect using high throughput DNA sequencing.


  Driven by the need for better models that allow one to shed light into the
question how life's diversity has evolved, phylogenetic networks have now
joined phylogenetic trees in the center of phylogenetics research. Like
phylogenetic trees, such networks canonically induce collections of
phylogenetic trees, clusters, and triplets, respectively. Thus it is not
surprising that many network approaches aim to reconstruct a phylogenetic
network from such collections. Related to the well-studied perfect phylogeny
problem, the following question is of fundamental importance in this context:
When does one of the above collections encode (i.e. uniquely describe) the
network that induces it? In this note, we present a complete answer to this
question for the special case of a level-1 (phylogenetic) network by
characterizing those level-1 networks for which an encoding in terms of one (or
equivalently all) of the above collections exists. Given that this type of
network forms the first layer of the rich hierarchy of level-k networks, k a
non-negative integer, it is natural to wonder whether our arguments could be
extended to members of that hierarchy for higher values for k. By giving
examples, we show that this is not the case.


  A promising new method for measuring intramolecular distances in solution
uses small-angle X-ray scattering interference between gold nanocrystal labels
(Mathew-Fenn et al, Science, 322, 446 (2008)). When applied to double stranded
DNA, it revealed that the DNA length fluctuations are strikingly strong and
correlated over at least 80 base pair steps. In other words, the DNA behaves as
accordion bellows, with distant fragments stretching and shrinking concertedly.
This hypothesis, however, disagrees with earlier experimental and computational
observations. This Letter shows that the discrepancy can be rationalized by
taking into account the cluster exclusion volume and assuming a moderate
long-range repulsion between them. The long-range interaction can originate
from an ion exclusion effect and cluster polarization in close proximity to the
DNA surface.


  We present numerical results on the temperature dependence of the
distribution of bubble lengths in DNA segments of various guanine-cytosine (GC)
concentrations. Base-pair openings are described by the Peyrard-Bishop-Dauxois
model and the corresponding thermal equilibrium distributions of bubbles are
obtained through Monte Carlo calculations for bubble sizes up to the order of a
hundred base pairs. The dependence of the parameters of bubble length
distribution on temperature and the GC content is investigated. We provide
simple expressions which approximately describe these relations. The variation
of the average bubble length is also presented. We find a temperature
dependence of the exponent c that appears in the distribution of bubble
lengths. If an analogous dependence exists in the loop entropy exponent of real
DNA, it may be relevant to understand overstretching in force-extension
experiments.


  A transition rate model of cargo transport by $N$ molecular motors is
proposed. Under the assumption of steady state, the force-velocity curve of
multi-motor system can be derived from the force-velocity curve of single
motor. Our work shows, in the case of low load, the velocity of multi-motor
system can decrease or increase with increasing motor number, which is
dependent on the single motor force-velocity curve. And most commonly, the
velocity decreases. This gives a possible explanation to some recent


  Increasingly, biologists are constructing evolutionary trees on large numbers
of overlapping sets of taxa, and then combining them into a `supertree' that
classifies all the taxa. In this paper, we ask how much coverage of the total
set of taxa is required by these subsets in order to ensure we have enough
information to reconstruct the supertree uniquely. We describe two results - a
combinatorial characterization of the covering subsets to ensure that at most
one supertree can be constructed from the smaller trees (whatever trees these
may be) and a more liberal analysis that asks only that the supertree is highly
likely to be uniquely specified by the tree structure on the covering subsets.


  Proteins can sometimes be knotted, and for many reasons the study of knotted
proteins is rapidly becoming very important. For example, it has been proposed
that a knot increases the stability of a protein. Knots may also alter
enzymatic activities and enhance binding. Moreover, knotted proteins may even
have some substantial biomedical significance in relation to illnesses such as
Parkinson's disease. But to a large extent the biological role of knots remains
a conundrum. In particular, there is no explanation why knotted proteins are so
scarce. Here we argue that knots are relatively rare because they tend to cause
swelling in proteins that are too short, and presently short proteins are
over-represented in the Protein Data Bank (PDB). Using Monte Carlo simulations
we predict that the figure-8 knot leads to the most compact protein
configuration when the number of amino acids is in the range of 200-600. For
the existence of the simplest knot, the trefoil, we estimate a theoretical
upper bound of 300-400 amino acids, in line with the available PDB data.


  Selection systems and the corresponding replicator equations model the
evolution of replicators with a high level of abstraction. In this paper we
apply novel methods of analysis of selection systems to the replicator
equations. To be suitable for the suggested algorithm the interaction matrix of
the replicator equation should be transformed; in particular the standard
singular value decomposition allows us to rewrite the replicator equation in a
convenient form. The original $n$-dimensional problem is reduced to the
analysis of asymptotic behavior of the solutions to the so-called escort
system, which in some important cases can be of significantly smaller dimension
than the original system. The Newton diagram methods are applied to study the
asymptotic behavior of the solutions to the escort system, when interaction
matrix has rank 1 or 2. A general replicator equation with the interaction
matrix of rank 1 is fully analyzed; the conditions are provided when the
asymptotic state is a polymorphic equilibrium. As an example of the system with
the interaction matrix of rank 2 we consider the problem from [Adams, M.R. and
Sornborger, A.T., J Math Biol, 54:357-384, 2007], for which we show, for
arbitrary dimension of the system and under some suitable conditions, that
generically one globally stable equilibrium exits on the 1-skeleton of the
simplex.


  A simple pillbox model with two adjustable parameters accounts for the
selectivity of both DEEA Ca channels and DEKA Na channels in many ionic
solutions of different composition and concentration. Only the side chains are
different in the model of the Ca and Na channels. Parameters are the same for
both channels in all solutions. 'Pauling' radii are used for ions. No
information from crystal structures is used in the model. Side chains are
grossly approximated as spheres. The predicted properties of the Na and Ca
channels are very different. How can such a simple model give such powerful
results when chemical intuition says that selectivity depends on the precise
relation of ions and side chains? We use Monte Carlo simulations of this model
that determine the most stable-lowest free energy-structure of the ions and
side chains. Structure is the computed consequence of the forces in this model.
The relationship of ions and side chains vary with ionic solution and are very
different in simulations of the Na and Ca channels. Selectivity is a
consequence of the 'induced fit' of side chains to ions and depends on the
flexibility (entropy) of the side chains as well as their location. The model
captures the relation of side chains and ions well enough to account for
selectivity of both Na channels and Ca channels in the wide range of conditions
measured in experiments. Evidently, the structures in the real Na and Ca
channels responsible for selectivity are self-organized, at their free energy
minimum. Oversimplified models are enough to account for selectivity if the
models calculate the 'most stable' structure as it changes from solution to
solution, and mutation to mutation.


  It has been recently shown that radical-ion pairs and their reactions are a
paradigm biological system manifesting non-trivial quantum effects, so far
invisible due to the phenomenological description of radical-ion-pair reactions
used until now. We here use the quantum-mechanically consistent master equation
describing magnetic-sensitive radical-ion-pair reactions to explain
experimental data [C. R. Timmel and K. B. Henbest, Phil. Trans. R. Soc. Lond. A
{\bf 362}, 2573 (2004); C. T. Rodgers, S. A. Norman, K. B. Henbest, C. R.
Timmel and P. J. Hore, J. Am. Chem. Soc. {\bf 129} 6746 (2007)] on the effect
of deuteration on the reaction yields. Anomalous behavior of radical-ion-pair
reactions after deuteration, i.e. data inconsistent with the predictions of the
phenomenological theory used so far, has been observed since the 70's and has
remained unexplained until now.


  Positive feedback regulation is ubiquitous in cell signaling networks, often
leading to binary outcomes in response to graded stimuli. However, the role of
such feedbacks in clustering, and in spatial spreading of activated molecules,
has come to be appreciated only recently. We focus on the latter, using a
simple model developed in the context of Ras activation with competing negative
and positive feedback mechanisms. We find that positive feedback, in the
presence of slow diffusion, results in clustering of activated molecules on the
plasma membrane, and rapid spatial spreading as the front of the cluster
propagates with a constant velocity (dependent on the feedback strength). The
advancing fronts of the clusters of the activated species are rough, with
scaling consistent with the Kardar-Parisi-Zhang (KPZ) equation in one
dimension. Our minimal model is general enough to describe signal transduction
in a wide variety of biological networks where activity in the
membrane-proximal region is subject to feedback regulation.


  We demonstrate that demographic noise can induce persistent spatial pattern
formation and temporal oscillations in the Levin-Segel predator-prey model for
plankton-herbivore population dynamics. Although the model exhibits a Turing
instability in mean field theory, demographic noise greatly enlarges the region
of parameter space where pattern formation occurs. To distinguish between
patterns generated by fluctuations and those present at the mean field level in
real ecosystems, we calculate the power spectrum in the noise-driven case and
predict the presence of fat tails not present in the mean field case. These
results may account for the prevalence of large-scale ecological patterns,
beyond that expected from traditional non-stochastic approaches.


  Stochastic effects may cause fade-out of an infectious disease in a
population immediately after an epidemic outbreak. We develop WKB theory to
determine the most probable path of the system toward epidemic fade-out and to
evaluate the fade-out probability. The most probable path is an instanton-like
orbit in the phase space of the underlying Hamiltonian flow.


  This paper addresses the problem of giving conditions for transcriptional
systems to be globally entrained to external periodic inputs. By using
contraction theory, a powerful tool from dynamical systems theory, it is shown
that certain systems driven by external periodic signals have the property that
all solutions converge to a fixed limit cycle. General results are proved, and
the properties are verified in the specific case of some models of
transcriptional systems. The basic mathematical results needed from contraction
theory are proved in the paper, making it self-contained.


  We investigate extinction of a long-lived self-regulating stochastic
population, caused by intrinsic (demographic) noise. Extinction typically
occurs via one of two scenarios depending on whether the absorbing state n=0 is
a repelling (scenario A) or attracting (scenario B) point of the deterministic
rate equation. In scenario A the metastable stochastic population resides in
the vicinity of an attracting fixed point next to the repelling point n=0. In
scenario B there is an intermediate repelling point n=n_1 between the
attracting point n=0 and another attracting point n=n_2 in the vicinity of
which the metastable population resides. The crux of the theory is WKB method
which assumes that the typical population size in the metastable state is
large. Starting from the master equation, we calculate the quasi-stationary
probability distribution of the population sizes and the (exponentially long)
mean time to extinction for each of the two scenarios. When necessary, the WKB
approximation is complemented (i) by a recursive solution of the
quasi-stationary master equation at small n and (ii) by the van Kampen
system-size expansion, valid near the fixed points of the deterministic rate
equation. The theory yields both entropic barriers to extinction and
pre-exponential factors, and holds for a general set of multi-step processes
when detailed balance is broken. The results simplify considerably for
single-step processes and near the characteristic bifurcations of scenarios A
and B.


  Proteins are regularly described with some general indices (mass fractal
dimension, surface fractal dimension, entropy, enthalpy, free energies,
hydrophobicity, denaturation temperature etc..), which are inherently
statistical in nature. These general indices emerge from innumerable (innately
context-dependent and time-dependent) interactions between various atoms of a
protein. Many a studies have been performed on the nature of these inter-atomic
interactions and the change of profile of atomic fluctuations that they cause.
However, we still do not know, under a given context, for a given duration of
time, how does a macroscopic biophysical property emerge from the cumulative
inter-atomic interactions. An exact answer to that question will involve
bridging the gap between nano-scale distinguishable atomic description and
macroscopic indistinguishable (statistical) measures, along the mesoscopic
scale of observation. In this work we propose a computationally implementable
mathematical model that derives expressions for observability of emergence of a
macroscopic biophysical property from a set of interacting (fluctuating) atoms.
Since most of the aforementioned interactions are non-linear in nature;
observability criteria are derived for both linear and the non-linear
descriptions of protein interior. The study assumes paramount importance in
21st-century biology, from both the theoretical and practical utilitarian point
of view.


  Boolean networks have been successfully used in modelling gene regulatory
networks. In this paper we propose a reduction method that reduces the
complexity of a Boolean network but keeps dynamical properties and topological
features and hence it makes the analysis easier; as a result, it allows for a
better understanding of the role of network topology on the dynamics. In
particular, we use the reduction method to study steady states of Boolean
models.


  We study the phase diagram of the standard pair approximation equations for
two different models in population dynamics, the
susceptible-infective-recovered-susceptible model of infection spread and a
predator-prey interaction model, on a network of homogeneous degree $k$. These
models have similar phase diagrams and represent two classes of systems for
which noisy oscillations, still largely unexplained, are observed in nature. We
show that for a certain range of the parameter $k$ both models exhibit an
oscillatory phase in a region of parameter space that corresponds to weak
driving. This oscillatory phase, however, disappears when $k$ is large. For
$k=3, 4$, we compare the phase diagram of the standard pair approximation
equations of both models with the results of simulations on regular random
graphs of the same degree. We show that for parameter values in the oscillatory
phase, and even for large system sizes, the simulations either die out or
exhibit damped oscillations, depending on the initial conditions. We discuss
this failure of the standard pair approximation model to capture even the
qualitative behavior of the simulations on large regular random graphs and the
relevance of the oscillatory phase in the pair approximation diagrams to
explain the cycling behavior found in real populations.


  Multisite covalent modification of proteins is omnipresent in eukaryotic
cells. A well-known example is the mitogen-activated protein kinase (MAPK)
cascade, where in each layer of the cascade a protein is phosphorylated at two
sites. It has long been known that the response of a MAPK pathway strongly
depends on whether the enzymes that modify the protein act processively or
distributively: distributive mechanism, in which the enzyme molecules have to
release the substrate molecules in between the modification of the two sites,
can generate an ultrasensitive response and lead to hysteresis and bistability.
We study by Green's Function Reaction Dynamics, a stochastic scheme that makes
it possible to simulate biochemical networks at the particle level and in time
and space, a dual phosphorylation cycle in which the enzymes act according to a
distributive mechanism. We find that the response of this network can differ
dramatically from that predicted by a mean-field analysis based on the chemical
rate equations. In particular, rapid rebindings of the enzyme molecules to the
substrate molecules after modification of the first site can markedly speed up
the response, and lead to loss of ultrasensitivity and bistability. In essence,
rapid enzyme-substrate rebindings can turn a distributive mechanism into a
processive mechanism. We argue that slow ADP release by the enzymes can protect
the system against these rapid rebindings, thus enabling ultrasensitivity and
bistability.


  We present here a detailed description of the model of ran-driven nuclear
transduction in living cells to be published elswere. The mathematical model
presented is the first to account for the active transport of molecules along
the cytoplasmic microtubules. All parameters entering the models are thoroughly
discussed. The simulations reproduce the behavior observed experimentally.


  As a result of climate change, many populations have to modify their range to
follow the suitable areas - their "climate envelope" - often risking
extinction. During this migration process, they may face absolute boundaries to
dispersal, because of external environmental factors. Consequently, not only
the position, but also the shape of the climate envelope can be modified. We
use a reaction-diffusion model to analyse the effects on population persistence
of simultaneous changes in the climate envelope position and shape. When the
growth term is of logistic type, we show that extinction and persistence are
principally conditioned by the species mobility and the speed of climate
change, but not by the shape of the climate envelope. However, with a growth
term taking an Allee effect into account, we find a high sensitivity to the
variations of the shape of the climate envelope. In this case, the species
which have a high mobility, although they could more easily follow the
migration of the climate envelope, would be at risk of extinction when
encountering a local narrowing of the boundary geometry. This effect can be
attenuated by a progressive opening of the available space at the exit of the
narrowing, even though it transiently leads to a diminished area of the climate
envelope.


  There is a gap between single-species model predictions, and empirical
studies, regarding the effect of habitat fragmentation per se, i.e., a process
involving the breaking apart of habitat without loss of habitat. Empirical
works indicate that fragmentation can have positive as well as negative
effects, whereas, traditionally, single-species models predict a negative
effect of fragmentation. Within the class of reaction-diffusion models, studies
almost unanimously predict such a detrimental effect. In this paper,
considering a single-species reaction-diffusion model with a removal -- or
similarly harvesting -- term, in two dimensions, we find both positive and
negative effects of fragmentation of the reserves, i.e. the protected regions
where no removal occurs. Fragmented reserves lead to higher population sizes
for time-constant removal terms. On the other hand, when the removal term is
proportional to the population density, higher population sizes are obtained on
aggregated reserves, but maximum yields are attained on fragmented
configurations, and for intermediate harvesting intensities.


  We consider the problem of forecasting the regions at higher risk for newly
introduced invasive species. Favourable and unfavourable regions may indeed not
be known a priori, especially for exotic species whose hosts in native range
and newly-colonised areas can be different. Assuming that the species is
modelled by a logistic-like reaction-diffusion equation, we prove that the
spatial arrangement of the favourable and unfavourable regions can
theoretically be determined using only partial measurements of the population
density: 1) a local "spatio-temporal" measurement, during a short time period
and, 2) a "spatial" measurement in the whole region susceptible to
colonisation. We then present a stochastic algorithm which is proved
analytically, and then on several numerical examples, to be effective in
deriving these regions.


  We study the dynamical properties of small regulatory networks treated as non
autonomous dynamical systems called modules when working inside larger networks
or, equivalently when subject to external signal inputs. Particular emphasis is
put on the interplay between the internal properties of the open systems and
the different possible inputs on them to deduce new functionalities of the
modules. We use discrete-time, piecewise-affine and piecewise-contracting
models with interactions of a regulatory nature to perform our study.


  Gene conversion is a mechanism by which a double-strand break in a DNA
molecule is repaired using a homologous DNA molecule as a template. As a
result, one gene is 'copied and pasted' onto the other gene. It was recently
reported that the direction of gene conversion appears to be biased towards G
and C nucleotides. In this paper a stochastic model of the dynamics of the bias
in gene conversion is developed for a finite population of members in a
multigene family. The dual process is the biased voter model, which generates
an ancestral random graph for a given sample. An importance-sampling algorithm
for computing the likelihood of the sample is also given.


  We demonstrate how the stochastic global optimization scheme of Simulated
Annealing can be used to evaluate optimum parameters in the problem of DNA
breathing dynamics. The breathing dynamics is followed in accordance with the
stochastic Gillespie scheme with the denaturation zones in double stranded DNA
studied as a single molecule time series. Simulated Annealing is used to find
the optimum value of the activation energy for which the equilibrium bubble
size distribution matches with a given value. It is demonstrated that the
method overcomes even large noise in the input surrogate data.


  Previous investigation of coding sequence lengths (CDS) in the bacterial
circular chromosome revealed short range correlation in the series of these
data. We have further analyzed the averaged periodograms of these series and we
found that the organization of CDS can be well described by first order
autoregressive processes. This involves interaction between the neighboring
terms. The autoregressive analysis may have great potential in modeling various
physical and biological processes like light emission of galaxies, protein
organization, cell flickering, cognitive processes and perhaps others.


  Networks describing the interaction of the elements that constitute a complex
system grow and develop via a number of different mechanisms, such as the
addition and deletion of nodes, the addition and deletion of edges, as well as
the duplication or fusion of nodes. While each of these mechanisms can have a
different cause depending on whether the network is biological, technological,
or social, their impact on the network's structure, as well as its local and
global properties, is similar. This allows us to study how each of these
mechanisms affects networks either alone or together with the other processes,
and how they shape the characteristics that have been observed. We study how a
network's growth parameters impact the distribution of edges in the network,
how they affect a network's modularity, and point out that some parameters will
give rise to networks that have the opposite tendency, namely to display
anti-modularity. Within the model we are describing, we can search the space of
possible networks for parameter sets that generate networks that are very
similar to well-known and well-studied examples, such as the brain of a worm,
and the network of interactions of the proteins in baker's yeast.


  Motivation: Prediction of ligands for proteins of known 3D structure is
important to understand structure-function relationship, predict molecular
function, or design new drugs. Results: We explore a new approach for ligand
prediction in which binding pockets are represented by atom clouds. Each target
pocket is compared to an ensemble of pockets of known ligands. Pockets are
aligned in 3D space with further use of convolution kernels between clouds of
points. Performance of the new method for ligand prediction is compared to
those of other available measures and to docking programs. We discuss two
criteria to compare the quality of similarity measures: area under ROC curve
(AUC) and classification based scores. We show that the latter is better suited
to evaluate the methods with respect to ligand prediction. Our results on
existing and new benchmarks indicate that the new method outperforms other
approaches, including docking. Availability: The new method is available at
http://cbio.ensmp.fr/paris/ Contact: mikhail.zaslavskiy@mines-paristech.fr


  The "Law of the Minimum" states that growth is controlled by the scarcest
resource (limiting factor). This concept was originally applied to plant or
crop growth (Justus von Liebig, 1840) and quantitatively supported by many
experiments. Some generalizations based on more complicated "dose-response"
curves were proposed. Violations of this law in natural and experimental
ecosystems were also reported. We study models of adaptation in ensembles of
similar organisms under load of environmental factors and prove that violation
of Liebig's law follows from adaptation effects. If the fitness of an organism
in a fixed environment satisfies the Law of the Minimum then adaptation
equalizes the pressure of essential factors and therefore acts against the
Liebig's law. This is the the Law of the Minimum paradox: if for a randomly
chosen pair "organism-environment" the Law of the Minimum typically holds,
then, in a well-adapted system, we have to expect violations of this law.
  For the opposite interaction of factors (a synergistic system of factors
which amplify each other) adaptation leads from factor equivalence to
limitations by a smaller number of factors.
  For analysis of adaptation we develop a system of models based on Selye's
idea of the universal adaptation resource (adaptation energy). These models
predict that under the load of an environmental factor a population separates
into two groups (phases): a less correlated, well adapted group and a highly
correlated group with a larger variance of attributes, which experiences
problems with adaptation. Some empirical data are presented and evidences of
interdisciplinary applications to econometrics are discussed.


  Classification of proteins based on their structure provides a valuable
resource for studying protein structure, function and evolutionary
relationships. With the rapidly increasing number of known protein structures,
manual and semi-automatic classification is becoming ever more difficult and
prohibitively slow. Therefore, there is a growing need for automated, accurate
and efficient classification methods to generate classification databases or
increase the speed and accuracy of semi-automatic techniques. Recognizing this
need, several automated classification methods have been developed. In this
survey, we overview recent developments in this area. We classify different
methods based on their characteristics and compare their methodology, accuracy
and efficiency. We then present a few open problems and explain future
directions.


  Any process in which competing solutions replicate with errors and numbers of
their copies depend on their respective fitnesses is the evolutionary
optimization process. As during carcinogenesis mutated genomes replicate
according to their respective qualities, carcinogenesis obviously qualifies as
the evolutionary optimization process and conforms to common mathematical
basis. The optimization view accents statistical nature of carcinogenesis
proposing that during it the crucial role is actually played by the allocation
of trials. Optimal allocation of trials requires reliable schemas' fitnesses
estimations which necessitate appropriate, fitness landscape dependent,
statistics of population. In the spirit of the applied conceptual framework,
features which are known to decrease efficiency of any evolutionary
optimization procedure (or inhibit it completely) are anticipated as
"therapies" and reviewed. Strict adherence to the evolutionary optimization
framework leads us to some counterintuitive implications which are, however, in
agreement with recent experimental findings, such as sometimes observed more
aggressive and malignant growth of therapy surviving cancer cells.


  The activation/repression of a given gene is typically regulated by multiple
transcription factors (TFs) that bind at the gene regulatory region and recruit
RNA polymerase (RNAP). The interactions between the promoter region and TFs and
between different TFs specify the dynamic responses of the gene under different
physiological conditions. By choosing specific regulatory interactions with up
to three transcription factors, we designed several functional motifs, each of
which is shown to perform a certain function and can be integrated into larger
networks. We analyzed three kinds of networks: (i) Motifs derived from
incoherent feedforward motifs, which behave as `amplitude filters', or
`concentration detectors'. These motifs respond maximally to input
transcription factors with concentrations within a certain range. From these
motifs homeostatic and pulse generating networks are derived. (ii) Tunable
network motifs, which can behave as oscillators or switches for low and high
concentrations of an input transcription factor, respectively. (iii)
Transcription factor controlled adjustable gates, which switch between AND/OR
gate characteristics, depending on the concentration of the input transcription
factor. This study has demonstrated the utility of feedforward networks and the
flexibility of specific transcriptional binding kinetics in generating new
novel behaviors. The flexibility of feedforward networks as dynamic units may
explain the apparent frequency that such motifs are found in real biological
networks.


  Aging is always present, tailoring our interactions with others and
postulating a finite lifespan during which we are able to exercise them. We
consider the prisoner's dilemma game on a square lattice, and examine how
quenched age distributions and different aging protocols influence the
evolution of cooperation when taking the life experience and knowledge
accumulation into account as time passes. In agreement with previous studies,
we find that a quenched assignment of age to players, introducing heterogeneity
to the game, substantially promotes cooperative behavior. Introduction of aging
and subsequent death as a coevolutionary process may act detrimental on
cooperation but enhances it efficiently if the offspring of individuals that
have successfully passed their strategy is considered newborn. We study
resulting age distributions of players, and show that the heterogeneity is
vital yet insufficient for explaining the observed differences in cooperator
abundance on the spatial grid. The unexpected increment of cooperation levels
can be explained by a dynamical effect that has a highly selective impact on
the propagation of cooperator and defector states.


  The distributed genome hypothesis states that the set of genes in a
population of bacteria is distributed over all individuals that belong to the
specific taxon. It implies that certain genes can be gained and lost from
generation to generation. We use the random genealogy given by a Kingman
coalescent in order to superimpose events of gene gain and loss along ancestral
lines. Gene gains occur at a constant rate along ancestral lines. We assume
that gained genes have never been present in the population before. Gene losses
occur at a rate proportional to the number of genes present along the ancestral
line. In this infinitely many genes model we derive moments for several
statistics within a sample: the average number of genes per individual, the
average number of genes differing between individuals, the number of
incongruent pairs of genes, the total number of different genes in the sample
and the gene frequency spectrum. We demonstrate that the model gives a
reasonable fit with gene frequency data from marine cyanobacteria.


  We investigate diffusion-limited reactions between a diffusing particle and a
target site on a semiflexible polymer, a key factor determining the kinetics of
DNA-protein binding and polymerization of cytoskeletal filaments. Our theory
focuses on two competing effects: polymer shape fluctuations, which speed up
association, and the hydrodynamic coupling between the diffusing particle and
the chain, which slows down association. Polymer bending fluctuations are
described using a mean field dynamical theory, while the hydrodynamic coupling
between polymer and particle is incorporated through a simple heuristic
approximation. Both of these we validate through comparison with Brownian
dynamics simulations. Neither of the effects has been fully considered before
in the biophysical context, and we show they are necessary to form accurate
estimates of reaction processes. The association rate depends on the stiffness
of the polymer and the particle size, exhibiting a maximum for intermediate
persistence length and a minimum for intermediate particle radius. In the
parameter range relevant to DNA-protein binding, the rate increase is up to
100% compared to the Smoluchowski result for simple center-of-mass motion. The
quantitative predictions made by the theory can be tested experimentally.


  Characterizing the capabilities, criticalities and response to perturbations
of genome-scale metabolic networks is a basic problem with important
applications. A key question concerns the identification of the potentially
most harmful knockouts. The integration of combinatorial methods with sampling
techniques to explore the space of viable flux states may provide crucial
insights on this issue. We assess the replaceability of every metabolic
conversion in the human red blood cell by enumerating the alternative paths
from substrate to product, obtaining a complete map of the potential damage of
single enzymopathies. Sampling the space of optimal flux states in the healthy
and in the mutated cell reveals both correlations and complementarity between
topologic and dynamical aspects.


  We present a novel approach to investigate the long-time stochastic dynamics
of multi-dimensional classical systems, in contact with a heat-bath. When the
potential energy landscape is rugged, the kinetics displays a decoupling of
short and long time scales and both
  Molecular Dynamics (MD) or Monte Carlo (MC) simulations are generally
inefficient. Using a field theoretic approach, we perform analytically the
average over the short-time stochastic fluctuations. This way, we obtain an
effective theory, which generates the same long-time dynamics of the original
theory, but has a lower time resolution power. Such an approach is used to
develop an improved version of the MC algorithm, which is particularly suitable
to investigate the dynamics of rare conformational transitions. In the specific
case of molecular systems at room temperature, we show that elementary
integration time steps used to simulate the effective theory can be chosen a
factor ~100 larger than those used in the original theory. Our results are
illustrated and tested on a simple system, characterized by a rugged energy
landscape.


  Time series data provided by single-molecule Forster resonance energy
transfer (sm-FRET) experiments offer the opportunity to infer not only model
parameters describing molecular complexes, e.g. rate constants, but also
information about the model itself, e.g. the number of conformational states.
Resolving whether or how many of such states exist requires a careful approach
to the problem of model selection, here meaning discriminating among models
with differing numbers of states. The most straightforward approach to model
selection generalizes the common idea of maximum likelihood-selecting the most
likely parameter values-to maximum evidence: selecting the most likely model.
In either case, such inference presents a tremendous computational challenge,
which we here address by exploiting an approximation technique termed
variational Bayes. We demonstrate how this technique can be applied to temporal
data such as smFRET time series; show superior statistical consistency relative
to the maximum likelihood approach; and illustrate how model selection in such
probabilistic or generative modeling can facilitate analysis of closely related
temporal data currently prevalent in biophysics. Source code used in this
analysis, including a graphical user interface, is available open source via
http://vbFRET.sourceforge.net


  Many biological channels perform highly selective transport without direct
input of metabolic energy and without transitions from a 'closed' to an 'open'
state during transport. Mechanisms of selectivity of such channels serve as an
inspiration for creation of artificial nano-molecular sorting devices and
bio-sensors. To elucidate the transport mechanisms, it is important to
understand the transport on the single molecule level in the experimentally
relevant regime when multiple particles are crowded in the channel. In this
paper we analyze the effects of inter-particle crowding on the non-equilibrium
transport times through a finite-length channel by means of analytical theory
and computer simulations.


  Among the realistic ingredients to be considered in the computational
modeling of infectious diseases, human mobility represents a crucial challenge
both on the theoretical side and in view of the limited availability of
empirical data. In order to study the interplay between small-scale commuting
flows and long-range airline traffic in shaping the spatio-temporal pattern of
a global epidemic we i) analyze mobility data from 29 countries around the
world and find a gravity model able to provide a global description of
commuting patterns up to 300 kms; ii) integrate in a worldwide structured
metapopulation epidemic model a time-scale separation technique for evaluating
the force of infection due to multiscale mobility processes in the disease
dynamics. Commuting flows are found, on average, to be one order of magnitude
larger than airline flows. However, their introduction into the worldwide model
shows that the large scale pattern of the simulated epidemic exhibits only
small variations with respect to the baseline case where only airline traffic
is considered. The presence of short range mobility increases however the
synchronization of subpopulations in close proximity and affects the epidemic
behavior at the periphery of the airline transportation infrastructure. The
present approach outlines the possibility for the definition of layered
computational approaches where different modeling assumptions and granularities
can be used consistently in a unifying multi-scale framework.


  We consider a model of bilayer lipid membrane with bola-lipids. The
bola-lipid is modeled by linking tails of the hydrophobic chains in the
opposite monolayers within bilayer as a first approximation. A number of
thermodynamical characteristics are calculated analytically and compared with
the ones of a regular membrane without chains linkage. Pronounced difference
between lateral pressure profiles at the layers interface for linked and
regular bilayer models is found. In the linked case, the lateral pressure
mid-plane peak disappears, while the free energy per chain increases. We have
also calculated distribution of the orientaional order parameter of linked
chains across the bilayer, and found it is in contrast with the usual lipids
case.


  Signal-processing molecules inside cells are often present at low copy
number, which necessitates probabilistic models to account for intrinsic noise.
Probability distributions have traditionally been found using simulation-based
approaches which then require estimating the distributions from many samples.
Here we present in detail an alternative method for directly calculating a
probability distribution by expanding in the natural eigenfunctions of the
governing equation, which is linear. We apply the resulting spectral method to
three general models of stochastic gene expression: a single gene with multiple
expression states (often used as a model of bursting in the limit of two
states), a gene regulatory cascade, and a combined model of bursting and
regulation. In all cases we find either analytic results or numerical
prescriptions that greatly outperform simulations in efficiency and accuracy.
In the last case, we show that bimodal response in the limit of slow switching
is not only possible but optimal in terms of information transmission.


  The transmission rate of many acute infectious diseases varies significantly
in time, but the underlying mechanisms are usually uncertain. They may include
seasonal changes in the environment, contact rate, immune system response, etc.
The transmission rate has been thought difficult to measure directly. We
present a new algorithm to compute the time-dependent transmission rate
directly from prevalence data, which makes no assumptions about the number of
susceptibles or vital rates. The algorithm follows our complete and explicit
solution of a mathematical inverse problem for SIR-type transmission models. We
prove that almost any infection profile can be perfectly fitted by an SIR model
with variable transmission rate. This clearly shows a serious danger of
over-fitting such transmission models. We illustrate the algorithm with
historic UK measles data and our observations support the common belief that
measles transmission was predominantly driven by school contacts.


  The structure of the genetic code is discussed in formal terms. A rectangular
table of the code ("the code matrix"), whose properties reveal its arithmetical
content tagged with the information symbols in several notations. New
parameters used to analyze of the code matrix, the serial numbers of the
encoded products and coding elements, ordered by molecular mass. The structural
similarity of the amino acid sequences corresponding to two aminoacyl tRNA
synthetases classes is found. The code matrix shows how can be organized the
so-called second genetic code. The symmetrical pattern of the matrix is
supported with the other parameters; it also serves as a basis to construct a
3D model of the genetic code which follows the structure of the simplest Plato
solid, tetrahedron. The reasons for this unusual structure of the genetic code
remains unclear.


  We show that qualitatively different epidemic-like processes from distinct
societal domains (finance, social and commercial blockbusters, epidemiology)
can be quantitatively understood using the same unifying conceptual framework
taking into account the interplay between the timescales of the grouping and
fragmentation of social groups together with typical epidemic transmission
processes. Different domain-specific empirical infection profiles, featuring
multiple resurgences and abnormal decay times, are reproduced simply by varying
the timescales for group formation and individual transmission. Our model
emphasizes the need to account for the dynamic evolution of multi-connected
networks. Our results reveal a new minimally-invasive dynamical method for
controlling such outbreaks, help fill a gap in existing epidemiological theory,
and offer a new understanding of complex system response functions.


  In this work we suggest a statistical mechanics approach to the
classification of high-dimensional data according to a binary label. We propose
an algorithm whose aim is twofold: First it learns a classifier from a
relatively small number of data, second it extracts a sparse signature, {\it
i.e.} a lower-dimensional subspace carrying the information needed for the
classification. In particular the second part of the task is NP-hard, therefore
we propose a statistical-mechanics based message-passing approach. The
resulting algorithm is firstly tested on artificial data to prove its validity,
but also to elucidate possible limitations.
  As an important application, we consider the classification of
gene-expression data measured in various types of cancer tissues. We find that,
despite the currently low quantity and quality of available data (the number of
available samples is much smaller than the number of measured genes, limiting
thus strongly the predictive capacities), the algorithm performs slightly
better than many state-of-the-art approaches in bioinformatics.


  Diffusion is a key element of a large set of phenomena occurring on natural
and social systems modeled in terms of complex weighted networks. Here, we
introduce a general formalism that allows to easily write down mean-field
equations for any diffusive dynamics on weighted networks. We also propose the
concept of annealed weighted networks, in which such equations become exact. We
show the validity of our approach addressing the problem of the random walk
process, pointing out a strong departure of the behavior observed in quenched
real scale-free networks from the mean-field predictions. Additionally, we show
how to employ our formalism for more complex dynamics. Our work sheds light on
mean-field theory on weighted networks and on its range of validity, and warns
about the reliability of mean-field results for complex dynamics.


  We present a mean-field dynamical theory for single semiflexible polymers
which can precisely capture, without fitting parameters, recent fluorescence
correlation spectroscopy results on single monomer kinetics of DNA strands in
solution. Our approach works globally, covering three decades of strand length
and five decades of time: it includes the complex cross-overs occurring between
stiffness-dominated and flexible bending modes, along with larger-scale
rotational and center-of-mass motion. The accuracy of the theory stems in part
from long-range hydrodynamic coupling between the monomers, which makes a
mean-field description more realistic. Its validity extends even to short,
stiff fragments, where we also test the theory through Brownian hydrodynamics
simulations.


  To understand the system-wide organization of metabolism, different lines of
study have devised different categorizations of metabolites. The relationship
and difference between categories can provide new insights for a more detailed
description of the organization of metabolism. In this study, we investigate
the relative organization of three categorizations of metabolites -- pathways,
subcellular localizations and network clusters, by block-model techniques
borrowed from social-network studies and further characterize the categories
from topological point of view. The picture of the metabolism we obtain is that
of peripheral modules, characterized both by being dense network clusters and
localized to organelles, connected by a central, highly connected core.
Pathways typically run through several network clusters and localizations,
connecting them laterally. The strong overlap between organelles and network
clusters suggest that these are natural "modules" -- relatively independent
sub-systems. The different categorizations divide the core metabolism
differently suggesting that this, if possible, should be not be treated as a
module on par with the organelles. Although overlapping more than chance, none
of the pathways correspond very closely to a network cluster or localization.
This, we believe, highlights the benefits of different orthogonal
classifications and future experimental categorizations based on simple
principles.


  Comparing and computing distances between phylogenetic trees are important
biological problems, especially for models where edge lengths play an important
role. The geodesic distance measure between two phylogenetic trees with edge
lengths is the length of the shortest path between them in the continuous tree
space introduced by Billera, Holmes, and Vogtmann. This tree space provides a
powerful tool for studying and comparing phylogenetic trees, both in exhibiting
a natural distance measure and in providing a Euclidean-like structure for
solving optimization problems on trees. An important open problem is to find a
polynomial time algorithm for finding geodesics in tree space. This paper gives
such an algorithm, which starts with a simple initial path and moves through a
series of successively shorter paths until the geodesic is attained.


  Synthetic biology brings together concepts and techniques from engineering
and biology. In this field, computer-aided design (CAD) is necessary in order
to bridge the gap between computational modeling and biological data. An
application named TinkerCell has been created in order to serve as a CAD tool
for synthetic biology. TinkerCell is a visual modeling tool that supports a
hierarchy of biological parts. Each part in this hierarchy consists of a set of
attributes that define the part, such as sequence or rate constants. Models
that are constructed using these parts can be analyzed using various C and
Python programs that are hosted by TinkerCell via an extensive C and Python
API. TinkerCell supports the notion of a module, which are networks with
interfaces. Such modules can be connected to each other, forming larger modular
networks. Because TinkerCell associates parameters and equations in a model
with their respective part, parts can be loaded from databases along with their
parameters and rate equations. The modular network design can be used to
exchange modules as well as test the concept of modularity in biological
systems. The flexible modeling framework along with the C and Python API allows
TinkerCell to serve as a host to numerous third-party algorithms. TinkerCell is
a free and open-source project under the Berkeley Software Distribution
license. Downloads, documentation, and tutorials are available at
www.tinkercell.com.


  In this work, we describe a computational framework for the genome-wide
identification and characterization of mixed
transcriptional/post-transcriptional regulatory circuits in humans. We
concentrated in particular on feed-forward loops (FFL), in which a master
transcription factor regulates a microRNA, and together with it, a set of joint
target protein coding genes. The circuits were assembled with a two step
procedure. We first constructed separately the transcriptional and
post-transcriptional components of the human regulatory network by looking for
conserved over-represented motifs in human and mouse promoters, and 3'-UTRs.
Then, we combined the two subnetworks looking for mixed feed-forward regulatory
interactions, finding a total of 638 putative (merged) FFLs. In order to
investigate their biological relevance, we filtered these circuits using three
selection criteria: (I) GeneOntology enrichment among the joint targets of the
FFL, (II) independent computational evidence for the regulatory interactions of
the FFL, extracted from external databases, and (III) relevance of the FFL in
cancer. Most of the selected FFLs seem to be involved in various aspects of
organism development and differentiation. We finally discuss a few of the most
interesting cases in detail.


  Due to the stochastic nature of biochemical processes, the copy number of any
given type of molecule inside a living cell often exhibits large temporal
fluctuations. Here, we develop analytic methods to investigate how the noise
arising from a bursting input is reshaped by a transport reaction which is
either linear or of the Michaelis-Menten type. A slow transport rate smoothes
out fluctuations at the output end and minimizes the impact of bursting on the
downstream cellular activities. In the context of gene expression in eukaryotic
cells, our results indicate that transcriptional bursting can be substantially
attenuated by the transport of mRNA from nucleus to cytoplasm. Saturation of
the transport mediators or nuclear pores contributes further to the noise
reduction. We suggest that the mRNA transport should be taken into account in
the interpretation of relevant experimental data on transcriptional bursting.


  During embryonic development, differentiating cells respond via gene
expression to positional cues from morphogen gradients. While gene expression
is often highly erratic, embryonic development is precise. We show by theory
and simulations that diffusion of the expressed protein can enhance the
precision of its expression domain. While diffusion lessens the sharpness of
the expression boundary, it also reduces super-Poissonian noise by washing out
bursts of gene expression. Balancing these effects yields an optimal diffusion
constant maximizing the precision of the expression domain.


  Dispersal of species to find a more favorable habitat is important in
population dynamics. Dispersal rates evolve in response to the relative success
of different dispersal strategies. In a simplified deterministic treatment (J.
Dockery, V. Hutson, K. Mischaikow, et al., J. Math. Bio. 37, 61 (1998)) of two
species which differ only in their dispersal rates the slow species always
dominates. We demonstrate that fluctuations can change this conclusion and can
lead to dominance by the fast species or to coexistence, depending on
parameters. We discuss two different effects of fluctuations, and show that our
results are consistent with more complex treatments that find that selected
dispersal rates are not monotonic with the cost of migration.


  We study different configurations of permanent magnets and ferromagnetic
circuit, in order to optimize the magnetic field for the so-called ``magnetic
tweezers'' technique, for studing mechanical properties of DNA molecules. The
magnetic field is used to pull and twist a micron-sized superparamagnetic
bead,tethered to a microscope slide surface by a DNA molecule. The force
applied to the bead must be vertical, pointing upwards, being as strong as
possible, and it must decrease smoothly as the magnets are moved away from the
bead. In order to rotate the bead around the vertical axis, the field must be
horizontal. Moreover, the volume occupied by the magnets is limited by the
optical system. We simulate different configurations by solving the equations
for the static magnetic field; then, we test some of the configurations by
measuring the force acting on a bead tethered by a DNA molecule. One of the
configurations is able to generate a magnetic field ten times stronger than
usually reported.


  The identification of motifs--subgraphs that appear significantly more often
in a particular network than in an ensemble of randomized networks--has become
a ubiquitous method for uncovering potentially important subunits within
networks drawn from a wide variety of fields. We find that the most common
algorithms used to generate the ensemble from the real network change subgraph
counts in a highly correlated manner, so that one subgraph's status as a motif
may not be independent from the statuses of the other subgraphs. We demonstrate
this effect for the problem of 3- and 4-node motif identification in the
transcriptional regulatory networks of E. coli and S. cerevisiae in which
randomized networks are generated via an edge-swapping algorithm (Milo et al.,
Science 298:824, 2002). We show that correlations among 3-node subgraphs are
easily interpreted, and we present an information-theoretic tool that may be
used to identify correlations among subgraphs of any size.


  {\bf Purpose}: Manipulation of the radio-sensitivity of the nucleotide-base
driven by the spin blockade mechanism of diffusive free radicals against
ionizing radiation.
  {\bf Materials and methods}: We theoretically propose a mechanism which uses
the simultaneous application of circularly polarized light and an external
magnetic field to control the polarization of the free radicals and create S=1
electron-hole spin excitations (excitons) on nucleotide-base. We deploy an
ab-initio molecular dynamics model to calculate the characteristic parameters
of the light needed for optical transitions.
  {\bf Results}: As a specific example, we present the numerical results
calculated for a Guanine, in the presence of an OH free radical. To increase
the radio-resistivity of this system, a blue light source for the optical
pumping and induction of excitons on guanine can be used.
  {\bf Conclusions}: The effect of spin-injection on the formation of a free
energy barrier in diffusion controlled chemical reaction pathways leads to the
control of radiation-induced base damage. The proposed method allows us to
manipulate and partially suppress the damage induced by ionizing radiation.


  We apply a simple model system of patchy particles to study monodisperse
self-assembly, using the Platonic solids as target structures. We find marked
differences between the assembly behaviours of the different systems.
Tetrahedra, octahedra and icosahedra assemble easily, while cubes are more
challenging and dodecahedra do not assemble. We relate these differences to the
kinetics and thermodynamics of assembly, with the formation of large disordered
aggregates a particular important competitor to correct assembly. In
particular, the free energy landscapes of those targets that are easy to
assemble are funnel-like, whereas for the dodecahedral system the landscape is
relatively flat with little driving force to facilitate escape from disordered
aggregates.


  We study the self-assembly behaviour of patchy particles with `protein-like'
interactions that can be considered as a minimal model for the assembly of
viral capsids and other shell-like protein complexes. We thoroughly explore the
thermodynamics and dynamics of self assembly as a function of the parameters of
the model and find robust assembly of all target structures considered. Optimal
assembly occurs in the region of parameter space where a free energy barrier
regulates the rate of nucleation, thus preventing the premature exhaustion of
the supply of monomers that can lead to the formation of incomplete shells. The
interactions also need to be specific enough to prevent the assembly of
malformed shells, but whilst maintaining kinetic accessibility. Free-energy
landscapes computed for our model have a funnel-like topography guiding the
system to form the target structure, and show that the torsional component of
the interparticle interactions prevents the formation of disordered aggregates
that would otherwise act as kinetic traps.


  Homeodomain containing proteins are a broad class of DNA binding proteins
that are believed to primarily function as transcription factors.
Electrostatics interactions have been demonstrated to be critical for the
binding of the homeodomain to DNA. An examination of the electrostatic state of
homeodomain residues involved in DNA phosphate binding has demonstrated the
conserved presence of upward shifted pKa values among the basic residue of
lysine and arginine. It is believed that these pKa perturbations work to
facilitate binding to DNA since they ensure that the basic residues always
retain a positive charge.


  In this paper, we consider a simple stochastic epidemic model on large
regular random graphs and the stochastic process that corresponds to this
dynamics in the standard pair approximation. Using the fact that the nodes of a
pair are unlikely to share neighbors, we derive the master equation for this
process and obtain from the system size expansion the power spectrum of the
fluctuations in the quasi-stationary state. We show that whenever the pair
approximation deterministic equations give an accurate description of the
behavior of the system in the thermodynamic limit, the power spectrum of the
fluctuations measured in long simulations is well approximated by the
analytical power spectrum. If this assumption breaks down, then the cluster
approximation must be carried out beyond the level of pairs. We construct an
uncorrelated triplet approximation that captures the behavior of the system in
a region of parameter space where the pair approximation fails to give a good
quantitative or even qualitative agreement. For these parameter values, the
power spectrum of the fluctuations in finite systems can be computed
analytically from the master equation of the corresponding stochastic process.


  We study the effect of migration between coupled populations, or patches, on
the stability properties of multistrain disease dynamics. The epidemic model
used in this work displays a Hopf bifurcation to oscillations in a single well
mixed population. It is shown numerically that migration between two
non-identical patches stabilizes the endemic steady state, delaying the onset
of large amplitude outbreaks and reducing the total number of infections. This
result is motivated by analyzing generic Hopf bifurcations with different
frequencies and with diffusive coupling between them. Stabilization of the
steady state is again seen, indicating that our observation in the full
multistrain model is based on qualitative characteristics of the dynamics
rather than on details of the disease model.


  The spectrum and scale of fluctuations in protein structures affect the range
of cell phenomena, including stability of protein structures or their
fragments, allosteric transitions and energy transfer. The study presents a
statistical-thermodynamic analysis of relationship between the sequence
composition and the distribution of residue fluctuations in protein-protein
complexes. A one-node-per residue elastic network model accounting for the
nonhomogeneous protein mass distribution and the inter-atomic interactions
through the renormalized inter-residue potential is developed. Two factors, a
protein mass distribution and a residue environment, were found to determine
the scale of residue fluctuations. Surface residues undergo larger fluctuations
than core residues, showing agreement with experimental observations. Ranking
residues over the normalized scale of fluctuations yields a distinct
classification of amino acids into three groups. The structural instability in
proteins possibly relates to the high content of the highly fluctuating
residues and a deficiency of the weakly fluctuating residues in irregular
secondary structure elements (loops), chameleon sequences and disordered
proteins. Strong correlation between residue fluctuations and the sequence
composition of protein loops supports this hypothesis. Comparing fluctuations
of binding site residues (interface residues) with other surface residues shows
that, on average, the interface is more rigid than the rest of the protein
surface and Gly, Ala, Ser, Cys, Leu and Trp have a propensity to form more
stable docking patches on the interface. The findings have broad implications
for understanding mechanisms of protein association and stability of protein
structures.


  The population dynamics of predator-prey systems in the presence of
patch-specific predators are explored in a setting where the prey population
has access to both habitats. The emphasis is in situations where patch-prey
abundance drives prey-dispersal between patches, with the fragile prey
populations, that is, populations subject to the Allee effect. The resulting
four-dimensional model's mathematical analysis is carried out via sub-models
that focus in lower dimensional settings. The outcomes depend on, and in fact
they are quite sensitive to, the structure of the system, the range of
parameter values, and initial conditions. We show that the system can support
multi-stability and a diverse set of predator-prey life-history dynamics that
includes rather complex dynamical system outcomes. It is argued that in general
evolution should favor heterogeneous settings including Allee effects,
prey-refuges, and patch-specific predators.


  We study vaccine control for disease spread on an adaptive network modeling
disease avoidance behavior. Control is implemented by adding Poisson
distributed vaccination of susceptibles. We show that vaccine control is much
more effective in adaptive networks than in static networks due to an
interaction between the adaptive network rewiring and the vaccine application.
Disease extinction rates using vaccination are computed, and orders of
magnitude less vaccine application is needed to drive the disease to extinction
in an adaptive network than in a static one.


  Standard deviation measurements of intensity profiles of stationary single
fluorescent molecules are useful for studying axial localization, molecular
orientation, and a fluorescence imaging system's spatial resolution. Here we
report on the analysis of the precision of standard deviation measurements of
intensity profiles of single fluorescent molecules imaged using an EMCCD
camera. We have developed an analytical expression for the standard deviation
measurement error of a single image which is a function of the total number of
detected photons, the background photon noise, and the camera pixel size. The
theoretical results agree well with the experimental, simulation, and numerical
integration results. Using this expression, we show that single-molecule
standard deviation measurements offer nanometer precision for a large range of
experimental parameters.


  In the postgenome era many efforts have been dedicated to systematically
elucidate the complex web of interacting genes and proteins. These efforts
include experimental and computational methods. Microarray technology offers an
opportunity for monitoring gene expression level at the genome scale. By
recourse to information theory, this study proposes a mathematical approach to
reconstruct gene regulatory networks at coarse-grain level from high throughput
gene expression data. The method provides the {\it a posteriori} probability
that a given gene regulates positively, negatively or does not regulate each
one of the network genes. This approach also allows the introduction of prior
knowledge and the quantification of the information gain from experimental data
used in the inference procedure. This information gain can be used to chose
genes to be perturbed in subsequent experiments in order to refine the
knowledge about the architecture of an underlying gene regulatory network. The
performance of the proposed approach has been studied by {\it in numero}
experiments. Our results suggest that the approach is suitable for focusing on
size-limited problems, such as, recovering a small subnetwork of interest by
performing perturbation over selected genes.


  2D display is a fast and economical way of visualizing polymorphism and
comparing genomes, which is based on the separation of DNA fragments in two
steps, according first to their size and then to their sequence composition. In
this paper, we present an exhaustive study of the numerical issues associated
with a model aimed at predicting the final absolute locations of DNA fragments
in 2D display experiments. We show that simple expressions for the mobility of
DNA fragments in both dimensions allow one to reproduce experimental final
absolute locations to better than experimental uncertainties. On the other
hand, our simulations also point out that the results of 2D display experiments
are not sufficient to determine the best set of parameters for the modeling of
fragments separation in the second dimension and that additional detailed
measurements of the mobility of a few sequences are necessary to achieve this
goal. We hope that this work will help in establishing simulations as a
powerful tool to optimize experimental conditions without having to perform a
large number of preliminary experiments and to estimate whether 2D DNA display
is suited to identify a mutation or a genetic difference that is expected to
exist between the genomes of closely related organisms.


  We show that noise-induced oscillations in a gene circuit model display
stochastic coherence, that is, a maximum in the regularity of the oscillations
as a function of noise amplitude. The effect is manifest as a system-size
effect in a purely stochastic molecular reaction description of the circuit
dynamics. We compare the molecular reaction model behavior with that predicted
by a rate equation version of the same system. In addition, we show that
commonly used reduced models that ignore fast operator reactions do not capture
the full stochastic behavior of the gene circuit. Stochastic coherence occurs
under conditions that may be physiologically relevant.


  The usual Euclidean distance may be generalized to extended objects such as
polymers or membranes. Here, this distance is used for the first time as a cost
function to align structures. We examined the alignment of extended strands to
idealized beta-hairpins of various sizes using several cost functions,
including RMSD, MRSD, and the minimal distance. We find that using minimal
distance as a cost function typically results in an aligned structure that is
globally different than that given by an RMSD-based alignment


  We introduce an anisotropic mean-field approach for the dynamics of
semiflexible polymers under intermediate tension, the force range where a chain
is partially extended but not in the asymptotic regime of a nearly straight
contour. The theory is designed to exactly reproduce the lowest order
equilibrium averages of a stretched polymer, and treats the full complexity of
the problem: the resulting dynamics include the coupled effects of long-range
hydrodynamic interactions, backbone stiffness, and large-scale polymer contour
fluctuations. Validated by Brownian hydrodynamics simulations and comparison to
optical tweezer measurements on stretched DNA, the theory is highly accurate in
the intermediate tension regime over a broad dynamical range, without the need
for additional dynamic fitting parameters.


  During the synthesis of integral membrane proteins (IMPs), the hydrophobic
amino acids of the polypeptide sequence are partitioned mostly into the
membrane interior and hydrophilic amino acids mostly into the aqueous exterior.
We analyze the minimum free energy state of polypeptide sequences partitioned
into alpha-helical transmembrane (TM) segments and the role of thermal
fluctuations using a many-body statistical mechanics model. Results suggest
that IMP TM segment partitioning shares important features with general
theories of protein folding. For random polypeptide sequences, the minimum free
energy state at room temperature is characterized by fluctuations in the number
of TM segments with very long relaxation times. Simple assembly scenarios do
not produce a unique number of TM segments and jamming phenomena interfere with
segment placement. For sequences corresponding to IMPs, the minimum free energy
structure with the wildtype number of segments is free of number fluctuations
due to an anomalous gap in the energy spectrum, and simple assembly scenarios
produce this structure. There is a threshold number of random point mutations
beyond which the size of this gap is reduced so that the wildtype groundstate
is destabilized and number fluctuations reappear.


  We analyse a stochastic SIS epidemic amongst a finite population partitioned
into households. Since the population is finite, the epidemic will eventually
go extinct, i.e., have no more infectives in the population. We study the
effects of population size and within household transmission upon the time to
extinction. This is done through two approximations. The first approximation is
suitable for all levels of within household transmission and is based upon an
Ornstein-Uhlenbeck process approximation for the diseases fluctuations about an
endemic level relying on a large population. The second approximation is
suitable for high levels of within household transmission and approximates the
number of infectious households by a simple homogeneously mixing SIS model with
the households replaced by individuals. The analysis, supported by a simulation
study, shows that the mean time to extinction is minimized by moderate levels
of within household transmission.


  Using a variant of parallel tempering, we study the changes in sampling
within a simulation when the all-atom model is coupled to a Go-like potential.
We find that the native structure is not the lowest-energy configuration in the
all-atom force field. Adding a Go-term deforms the energy landscape in a way
that the native configuration becomes the global minimum.


  We study a two-type branching process which provides excellent description of
experimental data on cell dynamics in skin tissue (Clayton et al., 2007). The
model involves only a single type of progenitor cell, and does not require
support from a self-renewed population of stem cells. The progenitor cells
divide and may differentiate into post-mitotic cells. We derive an exact
solution of this model in terms of generating functions for the total number of
cells, and for the number of cells of different types. We also deduce large
time asymptotic behaviors drawing on our exact results, and on an independent
diffusion approximation.


  n-tuple power law widely exists in language, computer program code, DNA and
music. After a vast amount of Zipf analyses of n-tuple power law from empirical
data, we propose a model to explain the n-tuple power law feature existed in
these information translational carriers. Our model is a preferential selection
approach inspired by Simon's model which explained scaling law of single symbol
in a sequence Zipf analysis. The kernel mechanism is neat and simple in our
model. It can be simply described as a randomly copy and paste process, that
is, randomly select a random segment from current sequence and attach it to the
end repeatedly. The simulation of our model shows that n-tuple power law exists
in model generated data. Furthermore, two estimation equations: the Zipf
exponent and the minimal length of n-tuple for power law appears all correspond
to empirical data well. Our model can also reproduce the symmetry breaking
process of ATGC number differences in DNA data.


  We investigate the emergence of target waves in a cyclic predator-prey model
incorporating a periodic current of the three competing species in a small area
situated at the center of the square lattice. The periodic current acts as a
pacemaker, trying to impose its rhythm on the overall spatiotemporal evolution
of the three species. We show that the pacemaker is able to nucleate target
waves that eventually spread across the whole population, whereby three routes
leading to this phenomenon can be distinguished depending on the mobility of
the three species and the oscillation period of the localized current. First,
target waves can emerge due to the synchronization between the periodic current
and oscillations of the density of the three species on the spatial grid. The
second route is similar to the first, the difference being that the
synchronization sets in only intermittently. Finally, the third route towards
target waves is realized when the frequency of the pacemaker is much higher
than that characterizing the oscillations of the overall density of the three
species. By considering mobility and the frequency of the current as variable
parameters, we thus provide insights into the mechanisms of pattern formation
resulting from the interplay between local and global dynamics in systems
governed by cyclically competing species.


  It has been proven that the accessibility of the target sites has a critical
influence for miRNA and siRNA. In this paper, we present a program, rip2.0, not
only the energetically most favorable targets site based on the
hybrid-probability, but also a statistical sampling structure to illustrate the
statistical characterization and representation of the Boltzmann ensemble of
RNA-RNA interaction structures. The outputs are retrieved via backtracing an
improved dynamic programming solution for the partition function based on the
approach of Huang et al. (Bioinformatics). The $O(N^6)$ time and $O(N^4)$ space
algorithm is implemented in C (available from
\url{http://www.combinatorics.cn/cbpc/rip2.html})


  Background: Connectivity mapping is a process to recognize novel
pharmacological and toxicological properties in small molecules by comparing
their gene expression signatures with others in a database. A simple and robust
method for connectivity mapping with increased specificity and sensitivity was
recently developed, and its utility demonstrated using experimentally derived
gene signatures.
  Results: This paper introduces sscMap (statistically significant connections'
map), a Java application designed to undertake connectivity mapping tasks using
the recently published method. The software is bundled with a default
collection of reference gene-expression profiles based on the publicly
available dataset from the Broad Institute Connectivity Map 02, which includes
data from over 7000 Affymetrix microarrays, for over 1000 small-molecule
compounds, and 6100 treatment instances in 5 human cell lines. In addition, the
application allows users to add their custom collections of reference profiles
and is applicable to a wide range of other 'omics technologies.
  Conclusions: The utility of sscMap is two fold. First, it serves to make
statistically significant connections between a user-supplied gene signature
and the 6100 core reference profiles based on the Broad Institute expanded
dataset. Second, it allows users to apply the same improved method to
custom-built reference profiles which can be added to the database for future
referencing. The software can be freely downloaded from
http://purl.oclc.org/NET/sscMap


  We study the stochastic dynamics of growth and shrinkage of single actin
filaments taking into account insertion, removal, and ATP hydrolysis of
subunits either according to the vectorial mechanism or to the random
mechanism. In a previous work, we developed a model for a single actin or
microtubule filament where hydrolysis occurred according to the vectorial
mechanism: the filament could grow only from one end, and was in contact with a
reservoir of monomers. Here we extend this approach in several ways, by
including the dynamics of both ends and by comparing two possible mechanisms of
ATP hydrolysis. Our emphasis is mainly on two possible limiting models for the
mechanism of hydrolysis within a single filament, namely the vectorial or the
random model. We propose a set of experiments to test the nature of the precise
mechanism of hydrolysis within actin filaments.


  The mechanisms for explaining how a stable asymmetric chemical system can be
formed from a symmetric chemical system, in the absence of any asymmetric
influence other than statistical fluctuations, have been developed during the
last decades, focusing on the non-linear kinetic aspects. Besides the absolute
necessity of self-amplification processes, the importance of energetic aspects
is often underestimated. Going down to the most fundamental aspects, the
distinction between a single object -- that can be intrinsically asymmetric --
and a collection of objects -- whose racemic state is the more stable one --
must be emphasized. A system of strongly interacting objects can be described
as one single object retaining its individuality and a single asymmetry; weakly
or non-interacting objects keep their own individuality, and are prone to
racemize towards the equilibrium state. In the presence of energy fluxes,
systems can be maintained in an asymmetric non-equilibrium steady-state. Such
dynamical systems can retain their asymmetry for times longer than their
racemization time.


  The empirical evidence that human color categorization exhibits some
universal patterns beyond superficial discrepancies across different cultures
is a major breakthrough in cognitive science. As observed in the World Color
Survey (WCS), indeed, any two groups of individuals develop quite different
categorization patterns, but some universal properties can be identified by a
statistical analysis over a large number of populations. Here, we reproduce the
WCS in a numerical model in which different populations develop independently
their own categorization systems by playing elementary language games. We find
that a simple perceptual constraint shared by all humans, namely the human Just
Noticeable Difference (JND), is sufficient to trigger the emergence of
universal patterns that unconstrained cultural interaction fails to produce. We
test the results of our experiment against real data by performing the same
statistical analysis proposed to quantify the universal tendencies shown in the
WCS [Kay P and Regier T. (2003) Proc. Natl. Acad. Sci. USA 100: 9085-9089], and
obtain an excellent quantitative agreement. This work confirms that synthetic
modeling has nowadays reached the maturity to contribute significantly to the
ongoing debate in cognitive science.


  Respiration in bacteria involves a sequence of energetically-coupled electron
and proton transfers creating an electrochemical gradient of protons (a
proton-motive force) across the inner bacterial membrane. With a simple kinetic
model we analyze a redox loop mechanism of proton-motive force generation
mediated by a molecular shuttle diffusing inside the membrane. This model,
which includes six electron-binding and two proton-binding sites, reflects the
main features of nitrate respiration in E. coli bacteria. We describe the time
evolution of the proton translocation process. We find that the electron-proton
electrostatic coupling on the shuttle plays a significant role in the process
of energy conversion between electron and proton components. We determine the
conditions where the redox loop mechanism is able to translocate protons
against the transmembrane voltage gradient above 200 mV with a thermodynamic
efficiency of about 37%, in the physiologically important range of temperatures
from 250 to 350 K.


  We study a one-dimensional model of microtubule assembly/disassembly in which
GTP bound to tubulins within the microtubule undergoes stochastic hydrolysis.
In contrast to models that only consider a cap of GTP-bound tubulin, stochastic
hydrolysis allows GTP-bound tubulin remnants to exist within the microtubule.
We find that these buried GTP remnants enable an alternative mechanism of
recovery from shrinkage, and enhances fluctuations of filament lengths. Under
conditions for which this alternative mechanism dominates, an increasing
depolymerization rate leads to a decrease in dissociation rate and thus a net
increase in assembly.


  It has been suggested that microtubules and other cytoskeletal filaments may
act as electrical transmission lines. An electrical circuit model of the
microtubule is constructed incorporating features of its cylindrical structure
with nanopores in its walls. This model is used to study how ionic conductance
along the lumen is affected by flux through the nanopores when an external
potential is applied across its two ends. Based on the results of Brownian
dynamics simulations, the nanopores were found to have asymmetric inner and
outer conductances, manifested as nonlinear IV curves. Our simulations indicate
that a combination of this asymmetry and an internal voltage source arising
from the motion of the C-terminal tails causes a net current to be pumped
across the microtubule wall and propagate down the microtubule through the
lumen. This effect is demonstrated to enhance and add directly to the
longitudinal current through the lumen resulting from an external voltage
source, and could be significant in amplifying low-intensity endogenous
currents within the cellular environment or as a nano-bioelectronic device.


  We study a stochastic process describing the onset of spreading dynamics of
an epidemic in a population composed by individuals of three classes:
susceptible (S), infected (I), and recovered (R). The stochastic process is
defined by local rules and involves the following cyclic process:
S$\to$I$\to$R$\to$S (SIRS). The open process S$\to$I$\to$R (SIR) is studied as
a particular case of the SIRS process. The epidemic process is analyzed at
different levels of description: by a stochastic lattice gas model and by a
birth and death process. By means of Monte Carlo simulations and dynamical
mean-field approximations we show that the SIRS stochastic lattice gas model
exhibit a line of critical points separating two phases: an absorbing phase
where the lattice is completely full of S individuals and an active phase where
S, I and R individuals coexist, which may or may not present population cycles.
The critical line, that corresponds to the onset of epidemic spreading, is
shown to belong in the directed percolation universality class. By considering
the birth and death process we analyse the role of noise in stabilizing the
oscillations.


  The tendencies described in this work were revealed in the course of
examination of adenine and uracil distribution in the mRNA encoding sequence.
The study also discusses the usage of codons occupied by the amino acid
arginine in the table of the universal genetic code. All of the described
tendencies are qualitative, so neither sophisticated methods nor cumbersome
calculations are necessary to reveal and interpret them.


  Phase variation, or stochastic switching between alternative states of gene
expression, is common among microbes, and may be important in coping with
changing environments. We use a theoretical model to assess whether such
switching is a good strategy for growth in environments with occasional
catastrophic events. We find that switching can be advantageous, but only when
the environment is responsive to the microbial population. In our model,
microbes switch randomly between two phenotypic states, with different growth
rates. The environment undergoes sudden "catastrophes", the probability of
which depends on the composition of the population. We derive a simple
analytical result for the population growth rate. For a responsive environment,
two alternative strategies emerge. In the "no switching" strategy, the
population maximises its instantaneous growth rate, regardless of catastrophes.
In the "switching" strategy, the microbial switching rate is tuned to minimise
the environmental response. Which of these strategies is most favourable
depends on the parameters of the model. Previous studies have shown that
microbial switching can be favourable when the environment changes in an
unresponsive fashion between several states. Here, we demonstrate an
alternative role for phase variation in allowing microbes to maximise their
growth in catastrophic responsive environments.


  We mathematically modeled the receptor-activated MAP kinase signaling by
incorporating the regulation through cellular phosphatases. Activation induced
the alignment of a phosphatase cascade in parallel with the MAP kinase pathway.
A novel regulatory motif was thus generated, providing for the combinatorial
control of each MAPK intermediate. This ensured a non-linear mode of signal
transmission with the output being shaped by the balance between the strength
of input signal, and the activity gradient along the phosphatase axis. Shifts
in this balance yielded modulations in topology of the motif, thereby expanding
the repertoire of output responses. Thus we identify an added dimension to
signal processing, wherein the output response to an external stimulus is
additionally filtered through indicators that define the phenotypic status of
the cell.


  We introduce a new way to study molecular evolution within well-established
Hamilton-Jacobi formalism, showing that for a broad class of fitness landscapes
it is possible to derive dynamics analytically within the $1/N$-accuracy, where
$N$ is genome length. For smooth and monotonic fitness function this approach
gives two dynamical phases: smooth dynamics, and discontinuous dynamics. The
latter phase arises naturally with no explicit singular fitness function,
counter-intuitively. The Hamilton-Jacobi method yields straightforward
analytical results for the models that utilize fitness as a function of Hamming
distance from a reference genome sequence. We also show the way in which this
method gives dynamical phase structure for multi-peak fitness.


  The evolution model with parallel mutation-selection scheme is solved for the
case when selection is accompanied by base substitutions, insertions, and
deletions. The fitness is assumed to be either a single-peak function (i.e.,
having one finite discontinuity) or a smooth function of the Hamming distance
from the reference sequence.
  The mean fitness is exactly calculated in large-genome limit. In the case of
insertions and deletions the evolution characteristics depend on the choice of
reference sequence.


  In the coarse grained Brownian Dynamics simulation method the many solvent
molecules are replaced by random thermal kicks and an effective friction acting
on the particles of interest. For Brownian Dynamics the friction has to be so
strong that the particles' velocities are damped much faster than the duration
of an integration timestep. Here we show that this conceptual limit can be
dropped with an analytic integration of the equations of damped motion. In the
resulting Langevin integration scheme our recently proposed approximate form of
the hydrodynamic interactions between the particles can be incorparated
conveniently, leading to a fast multi-particle propagation scheme, which
captures more of the short-time and short-range solvent effects than standard
BD. Comparing the dynamics of a bead-spring model of a short peptide, we
recommend to run simulations of small biological molecules with the Langevin
type finite damping and to include the hydrodynamic interactions.


  Gaussian white noise is frequently used to model fluctuations in physical
systems. In Fokker-Planck theory, this leads to a vanishing probability density
near the absorbing boundary of threshold models. Here we derive the boundary
condition for the stationary density of a first-order stochastic differential
equation for additive finite-grained Poisson noise and show that the response
properties of threshold units are qualitatively altered. Applied to the
integrate-and-fire neuron model, the response turns out to be instantaneous
rather than exhibiting low-pass characteristics, highly non-linear, and
asymmetric for excitation and inhibition. The novel mechanism is exhibited on
the network level and is a generic property of pulse-coupled systems of
threshold units.


  Using methods of statistical physics, we present rigorous theoretical
calculations of Eigen's quasispecies theory with the truncated fitness
landscape which dramatically limits the available sequence space of a
reproducing quasispecies. Depending on the mutation rates, we observe three
phases, a selective one, an intermediate one with some residual order and a
completely randomized phase. Our results are applicable for the general case of
fitness landscape.


  The mechanisms by which adaptive phenotypes spread within an evolving
population after their emergence are understood fairly well. Much less is known
about the factors that influence the evolutionary accessibility of such
phenotypes, a pre-requisite for their emergence in a population. Here, we
investigate the influence of environmental quality on the accessibility of
adaptive phenotypes of Escherichia coli's central metabolic network. We used an
established flux-balance model of metabolism as the basis for a
genotype-phenotype map (GPM). We quantified the effects of seven qualitatively
different environments (corresponding to both carbohydrate and gluconeogenic
metabolic substrates) on the structure of this GPM. We found that the GPM has a
more rugged structure in qualitatively poorer environments, suggesting that
adaptive phenotypes could be intrinsically less accessible in such
environments. Nevertheless, on average ~74% of the genotype can be altered by
neutral drift, in the environment where the GPM is most rugged; this could
allow evolving populations to circumvent such ruggedness. Furthermore, we found
that the normalized mutual information (NMI) of genotype differences relative
to phenotype differences, which measures the GPM's capacity to transmit
information about phenotype differences, is positively correlated with
(simulation-based) estimates of the accessibility of adaptive phenotypes in
different environments. These results are consistent with the predictions of a
simple analytic theory and they suggest an intuitive information-theoretic
principle for evolutionary adaptation; adaptation could be faster in
environments where the GPM has a greater capacity to transmit information about
phenotype differences.