  Understanding factors that shape biodiversity and species coexistence across
scales is of utmost importance in ecology, both theoretically and for
conservation policies. Species-area relationships (SARs), measuring how the
number of observed species increases upon enlarging the sampled area,
constitute a convenient tool for quantifying the spatial structure of
biodiversity. While general features of species-area curves are quite universal
across ecosystems, some quantitative aspects can change significantly. Several
attempts have been made to link these variations to ecological forces. Within
the framework of spatially explicit neutral models, here we scrutinize the
effect of varying the local population size (i.e. the number of individuals per
site) and the level of habitat saturation (allowing for empty sites). We
conclude that species-area curves become shallower when the local population
size increases, while habitat saturation, unless strongly violated, plays a
marginal role. Our findings provide a plausible explanation of why SARs for
microorganisms are flatter than those for larger organisms.


  We consider the application of a popular penalised regression method, Ridge
Regression, to data with very high dimensions and many more covariates than
observations. Our motivation is the problem of out-of-sample prediction and the
setting is high-density genotype data from a genome-wide association or
resequencing study. Ridge regression has previously been shown to offer
improved performance for prediction when compared with other penalised
regression methods. One problem with ridge regression is the choice of an
appropriate parameter for controlling the amount of shrinkage of the
coefficient estimates. Here we propose a method for choosing the ridge
parameter based on controlling the variance of the predicted observations in
the model.
  Using simulated data, we demonstrate that our method outperforms subset
selection based on univariate tests of association and another penalised
regression method, HyperLasso regression, in terms of improved prediction
error. We extend our approach to regression problems when the outcomes are
binary (representing cases and controls, as is typically the setting for
genome-wide association studies) and demonstrate the method on a real data
example consisting of case-control and genotype data on Bipolar Disorder, taken
from the Wellcome Trust Case Control Consortium and the Genetic Association
Information Network.


  Approaches for testing sets of variants, such as a set of rare or common
variants within a gene or pathway, for association with complex traits are
important. In particular, set tests allow for aggregation of weak signal within
a set, can capture interplay among variants, and reduce the burden of multiple
hypothesis testing. Until now, these approaches did not address confounding by
family relatedness and population structure, a problem that is becoming more
important as larger data sets are used to increase power.
  Results: We introduce a new approach for set tests that handles confounders.
Our model is based on the linear mixed model and uses two random effects-one to
capture the set association signal and one to capture confounders. We also
introduce a computational speedup for two-random-effects models that makes this
approach feasible even for extremely large cohorts. Using this model with both
the likelihood ratio test and score test, we find that the former yields more
power while controlling type I error. Application of our approach to richly
structured GAW14 data demonstrates that our method successfully corrects for
population structure and family relatedness, while application of our method to
a 15,000 individual Crohn's disease case-control cohort demonstrates that it
additionally recovers genes not recoverable by univariate analysis.
  Availability: A Python-based library implementing our approach is available
at http://mscompbio.codeplex.com


  Holding on to one's strategy is natural and common if the later warrants
success and satisfaction. This goes against widespread simulation practices of
evolutionary games, where players frequently consider changing their strategy
even though their payoffs may be marginally different than those of the other
players. Inspired by this observation, we introduce an aspiration-based
win-stay-lose-learn strategy updating rule into the spatial prisoner's dilemma
game. The rule is simple and intuitive, foreseeing strategy changes only by
dissatisfied players, who then attempt to adopt the strategy of one of their
nearest neighbors, while the strategies of satisfied players are not subject to
change. We find that the proposed win-stay-lose-learn rule promotes the
evolution of cooperation, and it does so very robustly and independently of the
initial conditions. In fact, we show that even a minute initial fraction of
cooperators may be sufficient to eventually secure a highly cooperative final
state. In addition to extensive simulation results that support our
conclusions, we also present results obtained by means of the pair
approximation of the studied game. Our findings continue the success story of
related win-stay strategy updating rules, and by doing so reveal new ways of
resolving the prisoner's dilemma.


  Is there a functional role for quantum mechanics or coherent quantum effects
in biological processes? While this question is as old as quantum theory, only
recently have measurements on biological systems on ultra-fast time-scales shed
light on a possible answer. In this review we give an overview of the two main
candidates for biological systems which may harness such functional quantum
effects: photosynthesis and magnetoreception. We discuss some of the latest
evidence both for and against room temperature quantum coherence, and consider
whether there is truly a functional role for coherence in these biological
mechanisms. Finally, we give a brief overview of some more speculative examples
of functional quantum biology including the sense of smell, long-range quantum
tunneling in proteins, biological photoreceptors, and the flow of ions across a
cell membrane.


  A class of models is introduced describing the evolution of population
species whose carrying capacities are functionals of these populations. The
functional dependence of the carrying capacities reflects the fact that the
correlations between populations can be realized not merely through direct
interactions, as in the usual predator-prey Lotka-Volterra model, but also
through the influence of species on the carrying capacities of each other. This
includes the self-influence of each kind of species on its own carrying
capacity with delays. Several examples of such evolution equations with
functional carrying capacities are analyzed. The emphasis is given on the
conditions under which the solutions to the equations display extreme events,
such as finite-time death and finite-time singularity. Any destructive action
of populations, whether on their own carrying capacity or on the carrying
capacities of co-existing species, can lead to the instability of the whole
population that is revealed in the form of the appearance of extreme events,
finite-time extinctions or booms followed by crashes.


  In social dilemmas, cooperation among randomly interacting individuals is
often difficult to achieve. The situation changes if interactions take place in
a network where the network structure jointly evolves with the behavioral
strategies of the interacting individuals. In particular, cooperation can be
stabilized if individuals tend to cut interaction links when facing adverse
neighborhoods. Here we consider two different types of reaction to adverse
neighborhoods, and all possible mixtures between these reactions. When faced
with a gloomy outlook, players can either choose to cut and rewire some of
their links to other individuals, or they can migrate to another location and
establish new links in the new local neighborhood. We find that in general
local rewiring is more favorable for the evolution of cooperation than
emigration from adverse neighborhoods. Rewiring helps to maintain the diversity
in the degree distribution of players and favors the spontaneous emergence of
cooperative clusters. Both properties are known to favor the evolution of
cooperation on networks. Interestingly, a mixture of migration and rewiring is
even more favorable for the evolution of cooperation than rewiring on its own.
While most models only consider a single type of reaction to adverse
neighborhoods, the coexistence of several such reactions may actually be an
optimal setting for the evolution of cooperation.


  Recent experiments have illuminated a remarkable growth mechanism of
rod-shaped bacteria: proteins associated with cell wall extension move at
constant velocity in circles oriented approximately along the cell
circumference (Garner et al., Science (2011), Dominguez-Escobar et al. Science
(2011), van Teeffelen et al. PNAS (2011). We view these as dislocations in the
partially ordered peptidoglycan structure, activated by glycan strand extension
machinery, and study theoretically the dynamics of these interacting defects on
the surface of a cylinder. Generation and motion of these interacting defects
lead to surprising effects arising from the cylindrical geometry, with
important implications for growth. We also discuss how long range elastic
interactions and turgor pressure affect the dynamics of the fraction of
actively moving dislocations in the bacterial cell wall.


  Most helix-coil transition theories can be characterized by a set of three
parameters: energetic, describing the (free) energy cost of forming a helical
state in one repeating unit; entropic, accounting for the decrease of entropy
due to the helical state formation; and geometric, indicating how many
repeating units are affected by the formation of one helical state. Depending
on their effect on the helix-coil transition, solvents or co-solutes can be
classified with respect to their action on these parameters. Solvent
interactions that alter the entropic cost of helix formation by their osmotic
action can affect both the stability (transition temperature) and the
cooperativity (transition interval) of the helix-coil transition. A consistent
inclusion of osmotic pressure effects in a description of helix-coil transition
for poly(L-glutamic acid) in solution with polyethylene glycol can offer an
explanation of the experimentally observed linear dependence of transition
temperature on osmotic pressure as well as the concurrent changes in the
cooperativity of the transition.


  DNA double-strand breaks (DSBs) represent a serious source of damage for all
living things and thus there have been many quantitative studies of DSBs both
in vivo and in vitro. Despite this fact, the processes that lead to their
production have not yet been clearly understood, and there is no established
theory that can account for the statistics of their production, in particular,
the number of DSBs per base pair per unit Gy, here denoted by P1, which is the
most important parameter for evaluating the degree of risk posed by DSBs. Here,
using the single-molecule observation method with giant DNA molecules (166
kbp), we evaluate the number of DSBs caused by gamma-ray irradiation. We find
that P1 is nearly inversely proportional to the DNA concentration above a
certain threshold DNA concentration. A simple model that accounts for the
marked decrease of P1 shows that it is necessary to consider the
characteristics of giant DNA molecules as semiflexible polymers to interpret
the intrinsic mechanism of DSBs.


  A Systematic Evolution of Ligands by EXponential enrichment (SELEX)
experiment begins in round one with a random pool of oligonucleotides in
equilibrium solution with a target. Over a few rounds, oligonucleotides having
a high affinity for the target are selected. Data from a high throughput SELEX
experiment consists of lists of thousands of oligonucleotides sampled after
each round. Thus far, SELEX experiments have been very good at suggesting the
highest affinity oligonucleotide, but modeling lower affinity recognition site
variants has been difficult. Furthermore, an alignment step has always been
used prior to analyzing SELEX data. We present a novel model, based on a
biochemical parametrization of SELEX, which allows us to use data from all
rounds to estimate the affinities of the oligonucleotides. Most notably, our
model also aligns the oligonucleotides. We use our model to analyze a SELEX
experiment containing double stranded DNA oligonucleotides and the
transcription factor Bicoid as the target. Our SELEX model outperformed other
published methods for predicting putative binding sites for Bicoid as indicated
by the results of an in-vivo ChIP-chip experiment.


  If species abundance distributions are dominated by the simple processes of
individuals in a community giving birth and death independently, the result is
a log series distribution. I calculate this in a number of different ways,
using both master equations and also the combinatoric methods of elementary
statistical mechanics. Considerable light is shed on the nature and role of
'priors' in the language of MaxEnt.


  Genetic Regulatory Networks (GRNs) plays a vital role in the understanding of
complex biological processes. Modeling GRNs is significantly important in order
to reveal fundamental cellular processes, examine gene functions and
understanding their complex relationships. Understanding the interactions
between genes gives rise to develop better method for drug discovery and
diagnosis of the disease since many diseases are characterized by abnormal
behaviour of the genes. In this paper we have reviewed various evolutionary
algorithms-based approach for modeling GRNs and discussed various opportunities
and challenges.


  We consider the problem of learning a high-dimensional multi-task regression
model, under sparsity constraints induced by presence of grouping structures on
the input covariates and on the output predictors. This problem is primarily
motivated by expression quantitative trait locus (eQTL) mapping, of which the
goal is to discover genetic variations in the genome (inputs) that influence
the expression levels of multiple co-expressed genes (outputs), either
epistatically, or pleiotropically, or both. A structured input-output lasso
(SIOL) model based on an intricate l1/l2-norm penalty over the regression
coefficient matrix is employed to enable discovery of complex sparse
input/output relationships; and a highly efficient new optimization algorithm
called hierarchical group thresholding (HiGT) is developed to solve the
resultant non-differentiable, non-separable, and ultra high-dimensional
optimization problem. We show on both simulation and on a yeast eQTL dataset
that our model leads to significantly better recovery of the structured sparse
relationships between the inputs and the outputs, and our algorithm
significantly outperforms other optimization techniques under the same model.
Additionally, we propose a novel approach for efficiently and effectively
detecting input interactions by exploiting the prior knowledge available from
biological experiments.


  In sexual population, recombination reshuffles genetic variation and produces
novel combinations of existing alleles, while selection amplifies the fittest
genotypes in the population. If recombination is more rapid than selection,
populations consist of a diverse mixture of many genotypes, as is observed in
many populations. In the opposite regime, which is realized for example in the
facultatively sexual populations that outcross in only a fraction of
reproductive cycles, selection can amplify individual genotypes into large
clones. Such clones emerge when the fitness advantage of some of the genotypes
is large enough that they grow to a significant fraction of the population
despite being broken down by recombination. The occurrence of this "clonal
condensation" depends, in addition to the outcrossing rate, on the heritability
of fitness. Clonal condensation leads to a strong genetic heterogeneity of the
population which is not adequately described by traditional population genetics
measures, such as Linkage Disequilibrium. Here we point out the similarity
between clonal condensation and the freezing transition in the Random Energy
Model of spin glasses. Guided by this analogy we explicitly calculate the
probability, Y, that two individuals are genetically identical as a function of
the key parameters of the model. While Y is the analog of the spin-glass order
parameter, it is also closely related to rate of coalescence in population
genetics: Two individuals that are part of the same clone have a recent common
ancestor.


  Protein function often involves changes between different conformations.
Central questions are how these conformational changes are coupled to the
binding or catalytic processes during which they occur, and how they affect the
catalytic rates of enzymes. An important model system is the enzyme
dihydrofolate reductase (DHFR) from E. coli, which exhibits characteristic
conformational changes of the active-site loop during the catalytic step and
during unbinding of the product. In this article, we present a general kinetic
framework that can be used (1) to identify the ordering of events in the
coupling of conformational changes, binding and catalysis and (2) to determine
the rates of the substeps of coupled processes from a combined analysis of NMR
R2 relaxation dispersion experiments and traditional enzyme kinetics
measurements. We apply this framework to E. coli DHFR and find that the
conformational change during product unbinding follows a
conformational-selection mechanism, i.e. the conformational change occurs
predominantly prior to unbinding. The conformational change during the
catalytic step, in contrast, is an induced change, i.e. the change occurs after
the chemical reaction. We propose that the reason for these conformational
changes, which are absent in human and other vertebrate DHFRs, is robustness of
the catalytic rate against large pH variations and changes to substrate/product
concentrations in E. coli.


  Water availability is a major environmental driver affecting riparian and
wetland vegetation. The interaction between water table fluctuations and
vegetation in a stochastic environment contributes to the complexity of the
dynamics of these ecosystems. We investigate the possible emergence of spatial
patterns induced by spatio-temporal stochastic resonance in a simple model of
groundwater-dependent ecosystems. These spatio-temporal dynamics are driven by
the combined effect of three components: (i) an additive white Gaussian noise,
accounting for external random disturbances such as fires or fluctuations in
rain water availability, (ii) a weak periodic modulation in time, describing
hydrological drivers such as seasonal fluctuations of water table depth, and
(iii) a spatial coupling term, which takes into account the ability of
vegetation to spread and colonize other parts of the landscape. A suitable
cooperation between these three terms is able to give rise to ordered
structures which show spatial and temporal coherence, and are statistically
steady in time.


  The ability of Gaussian noise to induce ordered states in dynamical systems
is here presented in an overview of the main stochastic mechanisms able to
generate spatial patterns. These mechanisms involve: (i) a deterministic local
dynamics term, accounting for the local rate of variation of the field
variable, (ii) a noise component (additive or multiplicative) accounting for
the unavoidable environmental disturbances, and (iii) a linear spatial coupling
component, which provides spatial coherence and takes into account diffusion
mechanisms. We investigate these dynamics using analytical tools, such as
mean-field theory, linear stability analysis and structure function analysis,
and use numerical simulations to confirm these analytical results.


  We introduce and analyze a minimal model of epigenetic silencing in budding
yeast, built upon known biomolecular interactions in the system. Doing so, we
identify the epigenetic marks essential for the bistability of epigenetic
states. The model explicitly incorporates two key chromatin marks, namely H4K16
acetylation and H3K79 methylation, and explores whether the presence of
multiple marks lead to a qualitatively different systems behavior. We find that
having both modifications is important for the robustness of epigenetic
silencing. Besides the silenced and transcriptionally active fate of chromatin,
our model leads to a novel state with bivalent (i.e., both active and
silencing) marks under certain perturbations (knock-out mutations, inhibition
or enhancement of enzymatic activity). The bivalent state appears under several
perturbations and is shown to result in patchy silencing. We also show that the
titration effect, owing to a limited supply of silencing proteins, can result
in counter-intuitive responses. The design principles of the silencing system
is systematically investigated and disparate experimental observations are
assessed within a single theoretical framework. Specifically, we discuss the
behavior of Sir protein recruitment, spreading and stability of silenced
regions in commonly-studied mutants (e.g., sas2, dot1) illuminating the
controversial role of Dot1 in the systems biology of yeast silencing.


  Spatial reciprocity is a well known tour de force of cooperation promotion. A
thorough understanding of the effects of different population densities is
therefore crucial. Here we study the evolution of cooperation in social
dilemmas on different interaction graphs with a certain fraction of vacant
nodes. We find that sparsity may favor the resolution of social dilemmas,
especially if the population density is close to the percolation threshold of
the underlying graph. Regardless of the type of the governing social dilemma as
well as particularities of the interaction graph, we show that under pairwise
imitation the percolation threshold is a universal indicator of how dense the
occupancy ought to be for cooperation to be optimally promoted. We also
demonstrate that myopic updating, due to the lack of efficient spread of
information via imitation, renders the reported mechanism dysfunctional, which
in turn further strengthens its foundations.


  We review recent work on the statistical mechanics of Von Neumann's growth
model and discuss its application to cellular metabolic networks. In this
context, we present a detailed analysis of the physiological scenario
underlying optimality a la Von Neumann in the metabolism of the bacterium
Escherichia coli, showing that optimal solutions are characterized by a
considerable microscopic flexibility accompanied by a robust emergent picture
for the key physiological functions. This suggests that the ideas behind
optimal economic growth in Von Neumann's model can be helpful in uncovering
functional organization principles of cell energetics.


  Indirect reciprocity in which players cooperate with unacquainted other
players having good reputations is a mechanism for cooperation in relatively
large populations subjected to social dilemma situations. When the population
has group structure, as is often found in social networks, players in
experiments are considered to show behavior that deviates from existing
theoretical models of indirect reciprocity. First, players often show ingroup
favoritism (i.e., cooperation only within the group) rather than full
cooperation (i.e., cooperation within and across groups), even though the
latter is Pareto efficient. Second, in general, humans approximate outgroup
members' personal characteristics, presumably including the reputation used for
indirect reciprocity, by a single value attached to the group. Humans use such
a stereotypic approximation, a phenomenon known as outgroup homogeneity in
social psychology. I propose a model of indirect reciprocity in populations
with group structure to examine the possibility of ingroup favoritism and full
cooperation. In accordance with outgroup homogeneity, I assume that players
approximate outgroup members' personal reputations by a single reputation value
attached to the group. I show that ingroup favoritism and full cooperation are
stable under different social norms (i.e., rules for assigning reputations)
such that they do not coexist in a single model. If players are forced to
consistently use the same social norm for assessing different types of
interactions (i.e., ingroup versus outgroup interactions), only full
cooperation survives. The discovered mechanism is distinct from any form of
group selection. The results also suggest potential methods for reducing
ingroup bias to shift the equilibrium from ingroup favoritism to full
cooperation.


  Biochemical networks are used in computational biology, to model the static
and dynamical details of systems involved in cell signaling, metabolism, and
regulation of gene expression. Parametric and structural uncertainty, as well
as combinatorial explosion are strong obstacles against analyzing the dynamics
of large models of this type. Multi-scaleness is another property of these
networks, that can be used to get past some of these obstacles. Networks with
many well separated time scales, can be reduced to simpler networks, in a way
that depends only on the orders of magnitude and not on the exact values of the
kinetic parameters. The main idea used for such robust simplifications of
networks is the concept of dominance among model elements, allowing
hierarchical organization of these elements according to their effects on the
network dynamics. This concept finds a natural formulation in tropical
geometry. We revisit, in the light of these new ideas, the main approaches to
model reduction of reaction networks, such as quasi-steady state and
quasi-equilibrium approximations, and provide practical recipes for model
reduction of linear and nonlinear networks. We also discuss the application of
model reduction to backward pruning machine learning techniques.


  Following our previous works on the health state of a population and the
related health state function we proceed in developing a method to estimate the
Healthy Life Expectancy in connection to the relative impact of the Mortality
Area in the health state function graph. The resulting tools are applied to the
data sets for 1990, 2000 and 2009 for the Countries of the World Health
Organization (WHO). The application is done in the Excel Chart and it is ready
to be used for other time periods. The results are compared with the estimates
presented in the WHO report for 2000 showing a good relationship between the
estimates of the two methods. However, our proposed method, not based on
collection of data for diseases and other causes affecting a healthy life, is
more flexible. We can estimate the healthy life years from various time periods
when information related to diseases is missing. Keywords: Health state
function, Healthy life expectancy, Deterioration, Loss of healthy years, HALE,
DALE, World Health Organization, WHO


  We have adsorbed plasmid PuC19 DNA on a supported bilayer. The mobility of
the lipids within the bilayer ensured a 2D equilibrium of the DNA molecule. By
varying the fraction of cationic lipids in the membrane, we have tuned the
surface charge. Plasmids conformations were imaged by Atomic Force Microscopy
(AFM).We performed two sets of experiments: deposition from salt free solution
on charged bilayers and deposition from salty solutions on neutral bilayers.
Plasmids can be seen as rings, completely opened structures, or tightly
supercoiled plectonemes, depending on the experimental conditions. The
plectonemic conformation is observed either on charged surfaces (in the absence
of salt) or at 30 mM salt concentration on a neutral bilayer. We demonstrate
the equivalence of surface screening by mobile interfacial charges and bulk
screening from salt ions.


  The study of the interplay between the structure and dynamics of complex
multilevel systems is a pressing challenge nowadays. In this paper, we use a
semi-annealed approximation to study the stability properties of Random Boolean
Networks in multiplex (multi-layered) graphs. Our main finding is that the
multilevel structure provides a mechanism for the stabilization of the dynamics
of the whole system even when individual layers work on the chaotic regime,
therefore identifying new ways of feedback between the structure and the
dynamics of these systems. Our results point out the need for a conceptual
transition from the physics of single layered networks to the physics of
multiplex networks. Finally, the fact that the coupling modifies the phase
diagram and the critical conditions of the isolated layers suggests that
interdependency can be used as a control mechanism.


  Theory purports that animal foraging choices evolve to maximize returns, such
as net energy intake. Empirical research in both human and nonhuman animals
reveals that individuals often attend to the foraging choices of their
competitors while making their own foraging choices. Due to the complications
of gathering field data or constructing experiments, however, broad facts
relating theoretically optimal and empirically realized foraging choices are
only now emerging. Here, we analyze foraging choices of a cohort of
professional day traders who must choose between trading the same stock
multiple times in a row---patch exploitation---or switching to a different
stock---patch exploration---with potentially higher returns. We measure the
difference between a trader's resource intake and the competitors' expected
intake within a short period of time---a difference we call short-term
comparative returns. We find that traders' choices can be explained by foraging
heuristics that maximize their daily short-term comparative returns. However,
we find no one-best relationship between different trading choices and net
income intake. This suggests that traders' choices can be short-term win
oriented and, paradoxically, maybe maladaptive for absolute market returns.


  Disease spreading through human travel networks has been a topic of great
interest in recent years, as witnessed during outbreaks of influenza A (H1N1)
or SARS pandemics. One way to stop spreading over the airline network are
travel restrictions for major airports or network hubs based on the total
number of passengers of an airport. Here, we test alternative strategies using
edge removal, cancelling targeted flight connections rather than restricting
traffic for network hubs, for controlling spreading over the airline network.
We employ a SEIR metapopulation model that takes into account the population of
cities, simulates infection within cities and across the network of the top 500
airports, and tests different flight cancellation methods for limiting the
course of infection. The time required to spread an infection globally, as
simulated by a stochastic global spreading model was used to rank the candidate
control strategies. The model includes both local spreading dynamics at the
level of populations and long-range connectivity obtained from real global
airline travel data. Simulated spreading in this network showed that spreading
infected 37% less individuals after cancelling a quarter of flight connections
between cities, as selected by betweenness centrality. The alternative strategy
of closing down whole airports causing the same number of cancelled connections
only reduced infections by 18%. In conclusion, selecting highly ranked single
connections between cities for cancellation was more effective, resulting in
fewer individuals infected with influenza, compared to shutting down whole
airports. It is also a more efficient strategy, affecting fewer passengers
while producing the same reduction in infections.
  The network of connections between the top 500 airports is available under
the resources link on our website http://www.biological-networks.org.


  Approximate Bayesian computation (ABC) is a likelihood-free approach for
Bayesian inferences based on a rejection algorithm method that applies a
tolerance of dissimilarity between summary statistics from observed and
simulated data. Although several improvements to the algorithm have been
proposed, none of these improvements avoid the following two sources of
approximation: 1) lack of sufficient statistics: sampling is not from the true
posterior density given data but from an approximate posterior density given
summary statistics; and 2) non-zero tolerance: sampling from the posterior
density given summary statistics is achieved only in the limit of zero
tolerance. The first source of approximation can be improved by adding a
summary statistic, but an increase in the number of summary statistics could
introduce additional variance caused by the low acceptance rate. Consequently,
many researchers have attempted to develop techniques to choose informative
summary statistics. The present study evaluated the utility of a kernel-based
ABC method (Fukumizu et al. 2010, arXiv:1009.5736 and 2011, NIPS 24: 1549-1557)
for complex problems that demand many summary statistics. Specifically, kernel
ABC was applied to population genetic inference. We demonstrate that, in
contrast to conventional ABCs, kernel ABC can incorporate a large number of
summary statistics while maintaining high performance of the inference.


  In this work, we inspect the reliability of controlling and quelling an
epidemic disease mimicked by a Susceptible-Infected-Susceptible (SIS) model
defined on a complex network by means of current and implementable quarantine
and isolation policies. Specifically, we consider that each individual in the
network is originally linked to two types of individuals: members of the same
household and acquaintances. The topology of this network evolves taking into
account a probability $q$ that aims at representing the quarantine or isolation
process in which the connection with acquaintances is disrupted according to
standard policies of control of epidemics. Within current policies of
self-isolation and standard infection rates, our results show that the
propagation is either only controllable for hypothetical rates of compliance or
uncontrollable at all.


  Adaptation to local environments often occurs through natural selection
acting on a large number of loci, each having a weak phenotypic effect. One way
to detect these loci is to identify genetic polymorphisms that exhibit high
correlation with environmental variables used as proxies for ecological
pressures. Here, we propose new algorithms based on population genetics,
ecological modeling, and statistical learning techniques to screen genomes for
signatures of local adaptation. Implemented in the computer program "latent
factor mixed model" (LFMM), these algorithms employ an approach in which
population structure is introduced using unobserved variables. These fast and
computationally efficient algorithms detect correlations between environmental
and genetic variation while simultaneously inferring background levels of
population structure. Comparing these new algorithms with related methods
provides evidence that LFMM can efficiently estimate random effects due to
population history and isolation-by-distance patterns when computing
gene-environment correlations, and decrease the number of false-positive
associations in genome scans. We then apply these models to plant and human
genetic data, identifying several genes with functions related to development
that exhibit strong correlations with climatic gradients.


  Understanding the causes and effects of spatial aggregation is one of the
most fundamental problems in ecology. Aggregation is an emergent phenomenon
arising from the interactions between the individuals of the population, able
to sense only -at most- local densities of their cohorts. Thus, taking into
account the individual-level interactions and fluctuations is essential to
reach a correct description of the population. Classic deterministic equations
are suitable to describe some aspects of the population, but leave out features
related to the stochasticity inherent to the discreteness of the individuals.
Stochastic equations for the population do account for these
fluctuation-generated effects by means of demographic noise terms but, owing to
their complexity, they can be difficult (or, at times, impossible) to deal
with. Even when they can be written in a simple form, they are still difficult
to numerically integrate due to the presence of the "square-root" intrinsic
noise. In this paper, we discuss a simple way to add the effect of demographic
stochasticity to three classic, deterministic ecological examples where
aggregation plays an important role. We study the resulting equations using a
recently-introduced integration scheme especially devised to integrate
numerically stochastic equations with demographic noise. Aimed at scrutinizing
the ability of these stochastic examples to show aggregation, we find that the
three systems not only show patchy configurations, but also undergo a phase
transition belonging to the directed percolation universality class.


  The nni-distance is a well-known distance measure for phylogenetic trees. We
construct an efficient parallel approximation algorithm for the nni-distance in
the CRCW-PRAM model running in O(log n) time on O(n) processors. Given two
phylogenetic trees T1 and T2 on the same set of taxa and with the same
multi-set of edge-weights, the algorithm constructs a sequence of
nni-operations of weight at most O(log n) \cdot opt, where opt denotes the
minimum weight of a sequence of nni-operations transforming T1 into T2 . This
algorithm is based on the sequential approximation algorithm for the
nni-distance given by DasGupta et al. (2000). Furthermore, we show that the
problem of identifying so called good edge-pairs between two weighted
phylogenies can be computed in O(log n) time on O(n log n) processors.


  Reticulate events play an important role in determining evolutionary
relationships. The problem of computing the minimum number of such events to
explain discordance between two phylogenetic trees is a hard computational
problem. Even for binary trees, exact solvers struggle to solve instances with
reticulation number larger than 40-50. Here we present CycleKiller and
NonbinaryCycleKiller, the first methods to produce solutions verifiably close
to optimality for instances with hundreds or even thousands of reticulations.
Using simulations, we demonstrate that these algorithms run quickly for large
and difficult instances, producing solutions that are very close to optimality.
As a spin-off from our simulations we also present TerminusEst, which is the
fastest exact method currently available that can handle nonbinary trees: this
is used to measure the accuracy of the NonbinaryCycleKiller algorithm. All
three methods are based on extensions of previous theoretical work and are
publicly available. We also apply our methods to real data.


  The goal of this work is to propose a finite population counterpart to
Eigen's model, which incorporates stochastic effects. We consider a Moran model
describing the evolution of a population of size $m$ of chromosomes of length
$\ell$ over an alphabet of cardinality $\kappa$. The mutation probability per
locus is $q$. We deal only with the sharp peak landscape: the replication rate
is $\sigma>1$ for the master sequence and 1 for the other sequences. We study
the equilibrium distribution of the process in the regime where $\ell, m\to
+\infty$, $q\to 0$, $\ell q \to a$, $m/\ell\to\alpha$. We obtain an equation
$\alpha\phi(a)=\ln\kappa$ in the parameter space $(a,\alpha)$ separating the
regime where the equilibrium population is totally random from the regime where
a quasispecies is formed. We observe the existence of a critical population
size necessary for a quasispecies to emerge and we recover the finite
population counterpart of the error threshold. These results are supported by
computer simulations.


  An evolutionary tree (phylogenetic tree) is a binary, rooted, unordered tree
that models the evolutionary history of currently living species in which
leaves are labeled by species. In this paper, we investigate the problem of
finding the maximum consensus evolutionary tree from a set of given rooted
triplets. A rooted triplet is a phylogenetic tree on three leaves and shows the
evolutionary relationship of the corresponding three species. The mentioned
problem is known to be APX-hard. We present two new heuristic algorithms. For a
given set of m triplets on n species, the FastTree algorithm runs in O(mn^2)
which is faster than any other previously known algorithms, although, the
outcome is less satisfactory. The BPMTR algorithm runs in O(mn^3) and in
average performs better than any other previously known approximation
algorithms for this problem.


  Indirect reciprocity is a reputation-based mechanism for cooperation in
social dilemma situations when individuals do not repeatedly meet. The
conditions under which cooperation based on indirect reciprocity occurs have
been examined in great details. Most previous theoretical analysis assumed for
mathematical tractability that an individual possesses a binary reputation
value, i.e., good or bad, which depends on their past actions and other
factors. However, in real situations, reputations of individuals may be
multiple valued. Another puzzling discrepancy between the theory and
experiments is the status of the so-called image scoring, in which cooperation
and defection are judged to be good and bad, respectively, independent of other
factors. Such an assessment rule is found in behavioral experiments, whereas it
is known to be unstable in theory. In the present study, we fill both gaps by
analyzing a trinary reputation model. By an exhaustive search, we identify all
the cooperative and stable equilibria composed of a homogeneous population or a
heterogeneous population containing two types of players. Some results derived
for the trinary reputation model are direct extensions of those for the binary
model. However, we find that the trinary model allows cooperation under image
scoring under some mild conditions.


  Complex networks possess a rich, multi-scale structure reflecting the
dynamical and functional organization of the systems they model. Often there is
a need to analyze multiple networks simultaneously, to model a system by more
than one type of interaction or to go beyond simple pairwise interactions, but
currently there is a lack of theoretical and computational methods to address
these problems. Here we introduce a framework for clustering and community
detection in such systems using hypergraph representations. Our main result is
a generalization of the Perron-Frobenius theorem from which we derive spectral
clustering algorithms for directed and undirected hypergraphs. We illustrate
our approach with applications for local and global alignment of
protein-protein interaction networks between multiple species, for tripartite
community detection in folksonomies, and for detecting clusters of overlapping
regulatory pathways in directed networks.


  We propose a biochemical model providing the kinetic and energetic
descriptions of the processivity dynamics of kinesin and dinein molecular
motors. Our approach is a modified version of a well known model describing
kinesin dynamics and considers the presence of a competitive inhibition
reaction by ADP. We first reconstruct a continuous free-energy landscape of the
cycle catalyst process that allows us to calculate the number of steps given by
a single molecular motor. Then, we calculate an analytical expression
associated to the translational velocity and the stopping time of the molecular
motor in terms of time and ATP concentration. An energetic interpretation of
motor processivity is discussed in quantitative form by using experimental
data. We also predict a time duration of collective processes that agrees with
experimental reports.


  Dynamics in biological networks are in general robust against several
perturbations. We investigate a coupled map network as a model motivated by
gene regulatory networks and design systems which are robust against phenotypic
perturbations (perturbations in dynamics), as well as systems which are robust
against mutation (perturbations in network structure). To achieve such a
design, we apply a multicanonical Monte Carlo method. Analysis based on the
maximum Lyapunov exponent and parameter sensitivity shows that systems with
marginal stability, which are regarded as systems at the edge of chaos, emerge
when robustness against network perturbations is required. This emergence of
the edge of chaos is a self-organization phenomenon and does not need a fine
tuning of parameters.


  In living cells, motor proteins, such as kinesin and dynein can move
processively along microtubule (MT), and also detach from or attach to MT
stochastically. Experiments have found that, the traffic of motor might be
jammed, and various theoretical models have been designed to understand this
traffic jam phenomenon. But previous studies mainly focus on motor
attachment/detachment rate dependent properties. Recent experiment of Leduc
{\it et al.} found that the traffic jam formation of motor protein kinesin
depends also on the length of MT [Proc. Natl. Acad. Sci. U.S.A. {\bf 109},
6100-6105 (2012)]. In this study, the MT length dependent properties of motor
traffic will be analyzed. We found that MT length has one {\it critical value}
$N_c$, traffic jam occurs only when MT length $N>N_c$. The jammed length of MT
increases with total MT length, while the non-jammed MT length might not change
monotonically with the total MT length. The critical value $N_c$ increases with
motor detachment rate from MT, but decreases with motor attachment rate to MT.


  Despite great interest in techniques for stabilizing the dynamics of
biological populations and metapopulations, very few practicable methods have
been developed. We propose an easily implementable method, Adaptive Limiter
Control (ALC), for reducing extinction frequencies and the magnitude of
fluctuation in population sizes and demonstrate its efficacy in stabilizing
laboratory populations and metapopulations of Drosophila melanogaster.
Metapopulation stability was attained through a combination of reduced size
fluctuations and synchrony at the subpopulation level. Simulations indicated
that ALC was effective over a wide range of maximal population growth rates,
migration rates and population dynamics models. Since simulations using widely
applicable, nonspecies-specific models of population dynamics were able to
capture most features of the experimental data, we expect our results to be
applicable to a wide range of species.


  Stochastic thermodynamics as reviewed here systematically provides a
framework for extending the notions of classical thermodynamics like work, heat
and entropy production to the level of individual trajectories of well-defined
non-equilibrium ensembles. It applies whenever a non-equilibrium process is
still coupled to one (or several) heat bath(s) of constant temperature.
Paradigmatic systems are single colloidal particles in time-dependent laser
traps, polymers in external flow, enzymes and molecular motors in single
molecule assays, small biochemical networks and thermoelectric devices
involving single electron transport. For such systems, a first-law like energy
balance can be identified along fluctuating trajectories. Various integral and
detailed fluctuation theorems, which are derived here in a unifying approach
from one master theorem, constrain the probability distributions for work, heat
and entropy production depending on the nature of the system and the choice of
non-equilibrium conditions. For non-equilibrium steady states, particularly
strong results hold like a generalized fluctuation-dissipation theorem
involving entropy production. Ramifications and applications of these concepts
include optimal driving between specified states in finite time, the role of
measurement-based feedback processes and the relation between dissipation and
irreversibility. Efficiency and, in particular, efficiency at maximum power,
can be discussed systematically beyond the linear response regime for two
classes of molecular machines, isothermal ones like molecular motors, and heat
engines like thermoelectric devices, using a common framework based on a cycle
decomposition of entropy production.


  Simulations often involve the use of model parameters which are unknown or
uncertain. For this reason, simulation experiments are often repeated for
multiple combinations of parameter values, often iterating through parameter
values lying on a fixed grid. However, the use of a discrete grid places limits
on the dimension of the parameter space and creates the potential to miss
important parameter combinations which fall in the gaps between grid points.
Here we draw parallels with strategies for numerical integration and describe a
Markov chain Monte-Carlo strategy for exploring parameter values. We illustrate
the approach using examples from phylogenetics, archaeology, and epidemiology.


  We describe a modular modelling approach permitting curation, updating, and
distributed development of modules through joined community effort overcoming
the problem of keeping a combinatorially exploding number of monolithic models
up to date. For this purpose, the effects of genes and their mutated alleles on
downstream components are modeled by composable, metadata-containing Petri net
models organized in a database with version control, accessible through a web
interface. Gene modules can be coupled to protein modules through mRNA modules
by specific interfaces designed for the automatic, database-assisted
composition. Automatically assembled executable models may then consider cell
type-specific gene expression patterns and the resulting protein
concentrations. Gene modules and allelic interference modules may represent
effects of gene mutation and predict their pleiotropic consequences or uncover
complex genotype/phenotype relationships. Forward and reverse engineered
modules are fully compatible.


  We consider the interplay of vaccination and migration rates on disease
persistence in epidemiological systems. We show that short-term and long-term
migration can inhibit disease persistence. As a result, we show how migration
changes how vaccination rates should be chosen to maintain herd immunity. In a
system of coupled SIR models, we analyze how disease eradication depends
explicitly on vaccine distribution and migration connectivity. The analysis
suggests potentially novel vaccination policies that underscore the importance
of optimal placement of finite resources.


  We propose a new way of looking at global optimization of off-lattice protein
models. We present a dual optimization concept of predicting optimal sequences
as well as optimal folds. We validate the utility of the recently introduced
hidden-force Monte Carlo optimization algorithm by finding significantly lower
energy folds for minimalist protein models than previously reported. Further,
we also find the protein sequence that yields the lowest energy fold amongst
all sequences for a given chain length and residue mixture. In particular, for
protein models with a binary sequence, we show that the sequence-optimized
folds form more compact cores than the lowest energy folds of the historically
fixed, Fibonacci-series sequences of chain lengths of 13, 21, 34, 55, and 89.
We emphasize that while the protein model we used is minimalist, the
methodology is applicable to detailed protein models, and sequence optimization
may yield novel folds and aid de novo protein design.


  We study numerically domain growth and interface fluctuations in one- and
two-dimensional lattice systems composed of four species that interact in a
cyclic way. Particle mobility is implemented through exchanges of particles
located on neighboring lattice sites. For the chain we find that domain growth
strongly depends on the mobility, with a higher mobility yielding a larger
domain growth exponent. In two space dimensions, when also exchanges between
mutually neutral particles are possible, both domain growth and interface
fluctuations display universal regimes that are independent of the predation
and exchange rates.


  The prospect of new biological and industrial applications that require
communication in micro-scale, encourages research on the design of
bio-compatible communication networks using networking primitives already
available in nature. One of the most promising candidates for constructing such
networks is to adapt and engineer specific types of bacteria that are capable
of sensing, actuation, and above all, communication with each other. In this
paper, we describe a new architecture for networks of bacteria to form a data
collecting network, as in traditional sensor networks. The key to this
architecture is the fact that the node in the network itself is a bacterial
colony; as an individual bacterium (biological agent) is a tiny unreliable
element with limited capabilities. We describe such a network under two
different scenarios. We study the data gathering (sensing and multihop
communication) scenario as in sensor networks followed by the consensus problem
in a multi-node network. We will explain as to how the bacteria in the colony
collectively orchestrate their actions as a node to perform sensing and
relaying tasks that would not be possible (at least reliably) by an individual
bacterium. Each single bacterium in the colony forms a belief by sensing
external parameter (e.g., a molecular signal from another node) from the medium
and shares its belief with other bacteria in the colony. Then, after some
interactions, all the bacteria in the colony form a common belief and act as a
single node. We will model the reception process of each individual bacteria
and will study its impact on the overall functionality of a node. We will
present results on the reliability of the multihop communication for data
gathering scenario as well as the speed of convergence in the consensus
scenario.


  MicroRNAs (miRNAs) are a class of non-coding RNAs that regulate gene
expression. Identification of total number of miRNAs even in completely
sequenced organisms is still an open problem. However, researchers have been
using techniques that can predict limited number of miRNA in an organism. In
this paper, we have used homology based approach for comparative analysis of
miRNA of hexapoda group .We have used Apis mellifera, Bombyx mori, Anopholes
gambiae and Drosophila melanogaster miRNA datasets from miRBase repository. We
have done pair wise as well as multiple alignments for the available miRNAs in
the repository to identify and analyse conserved regions among related species.
Unfortunately, to the best of our knowledge, miRNA related literature does not
provide in depth analysis of hexapods. We have made an attempt to derive the
commonality among the miRNAs and to identify the conserved regions which are
still not available in miRNA repositories. The results are good approximation
with a small number of mismatches. However, they are encouraging and may
facilitate miRNA biogenesis for


  In cells, helicase translocation along nucleic acid is essential for many
biological processes. However, so far, the mechanism of this translocation is
not fully understood. Recent studies show that helicase might translocate
through two processes, active process and passive process, with different
translocation rate. In this study, a model including such two processes is
presented. In which, each of these two processes consists of two sub-processes,
chemical sub-process in which needed translocation factors are attached, and
mechanochemical sub-process in which helicase makes a forward translocation
step. Helicase can switch stochastically between these two processes with
external force dependent rates. By this model, ribosome translocation along
message RNA is detailed discussed. We found that, with the increase of external
force, the mean translocation rate of ribosome increases from one lower limit
to one upper limit, and both of these two limits increase with concentrations
of the translocation factors. Under high external force, ribosome translocates
mainly through active process. At saturating concentrations of translocation
factors, the translocation is governed by mechanochemical sub-processes.


  Two models of circulating tumor cell (CTC) dynamics have been proposed to
explain the phenomenon of tumor 'self-seeding', whereby CTCs repopulate the
primary tumor and accelerate growth: Primary Seeding, where cells from a
primary tumor shed into the vasculature and return back to the primary
themselves; and Secondary Seeding, where cells from the primary first
metastasize in a secondary tissue and form microscopic secondary deposits,
which then shed cells into the vasculature returning to the primary. These two
models are difficult to distinguish experimentally, yet the differences between
them is of great importance to both our understanding of the metastatic process
and also for designing methods of intervention. Therefore we developed a
mathematical model to test the relative likelihood of these two phenomena in
the subset of tumours whose shed CTCs first encounter the lung capillary bed,
and show that Secondary Seeding is several orders of magnitude more likely than
Primary seeding. We suggest how this difference could affect tumour evolution,
progression and therapy, and propose several possible methods of experimental
validation.


  I propose a stochastic SIS and SIRS system to include a Poisson measure term
to model anomalies in the dynamics. In particular the positive integrand in the
Poisson term is intended to model quarantine. Conditions are given for the
stability of the disease free equilibrium for both systems.


  We provide an introduction to enumerating and constructing invariants of
group representations via character methods. The problem is contextualised via
two case studies arising from our recent work: entanglement measures, for
characterising the structure of state spaces for composite quantum systems; and
Markov invariants, a robust alternative to parameter-estimation intensive
methods of statistical inference in molecular phylogenetics.


  For more than 50 years we have known that photosynthetic systems harvest
solar energy with almost unit {\it quantum efficiency}. However, recent
experimental evidence of {\it quantum coherence} during the excitonic energy
transport in photosynthetic organisms challenges our understanding of this
fundamental biological function. Currently, and despite numerous efforts, the
causal connection between coherence and efficiency is still a matter of debate.
We show, through the study of extensive simulations of quantum coherent
transport on networks, that three dimensional structures characterized by
centro-symmetric Hamiltonians are statistically more efficient than random
arrangements. Moreover, we demonstrate that the experimental data available for
the electronic Hamiltonians of the Fenna-Mathew-Olson (FMO) complex of sulfur
bacteria and of the crypophyte PC645 complex of marine algae are consistent
with this strong correlation of centro-symmetry with quantum efficiency. These
results show that what appears to be geometrically disordered complexes may
well exhibit a hidden symmetry which enhances the energy transport between
chromophores. We are confident that our results will motivate research to
explore the properties of nearly centro-symmetric Hamiltonians in more
realistic environments, and to unveil the role of symmetries for quantum
effects in biology. The unravelling of such symmetries may open novel
perspectives and suggest new design principles in the development of artificial
devices.


  Complex networks, comprised of individual elements that interact with each
other through reaction channels, are ubiquitous across many scientific and
engineering disciplines. Examples include biochemical, pharmacokinetic,
epidemiological, ecological, social, neural, and multi-agent networks. A common
approach to modeling such networks is by a master equation that governs the
dynamic evolution of the joint probability mass function of the underling
population process and naturally leads to Markovian dynamics for such process.
Due however to the nonlinear nature of most reactions, the computation and
analysis of the resulting stochastic population dynamics is a difficult task.
This review article provides a coherent and comprehensive coverage of recently
developed approaches and methods to tackle this problem. After reviewing a
general framework for modeling Markovian reaction networks and giving specific
examples, the authors present numerical and computational techniques capable of
evaluating or approximating the solution of the master equation, discuss a
recently developed approach for studying the stationary behavior of Markovian
reaction networks using a potential energy landscape perspective, and provide
an introduction to the emerging theory of thermodynamic analysis of such
networks. Three representative problems of opinion formation, transcription
regulation, and neural network dynamics are used as illustrative examples.


  In previous work, we studied host response to a pathogen which uses a cycle
of immunologically distinct stages to establish and maintain infection. We
showed that for generic parameter values, the system has a unique biologically
meaningful stable fixed point. That paper used a simplified model of T-cell
activation, making proliferation depend linearly on antigen-T-cell encounters.
Here we generalize the way in which T-cell proliferation depends on the sizes
of the antigenic populations. In particular, we allow this response to become
saturated at high levels of antigen. As a result, we show that this family of
generalized models shares the same steady-state behavior properties with the
simpler model contemplated in our previous work.


  We obtain conditions for eradication and permanence of infection for a
nonautonomous SIQR model with time-dependent parameters, that are not assumed
to be periodic. The incidence is given by functions of all compartments and the
threshold conditions are given by some numbers that play the role of the basic
reproduction number. We obtain simple threshold conditions in the autonomous,
asymptotically autonomous and periodic settings and show that our thresholds
coincide with the ones already established. Additionally, we obtain threshold
conditions for the general nonautonomous model with mass-action, standard and
quarantine-adjusted incidence.


  Wolbachia are maternally-inherited symbiotic bacteria commonly found in
arthropods, which are able to manipulate the reproduction of their host in
order to maximise their transmission. Here we use whole genome resequencing
data from 290 lines of Drosophila melanogaster from North America, Europe and
Africa to predict Wolbachia infection status, estimate cytoplasmic genome copy
number, and reconstruct Wolbachia and mtDNA genome sequences. Complete
Wolbachia and mitochondrial genomes show congruent phylogenies, consistent with
strict vertical transmission through the maternal cytoplasm and imperfect
transmission of Wolbachia. Bayesian phylogenetic analysis reveals that the most
recent common ancestor of all Wolbachia and mitochondrial genomes in D.
melanogaster dates to around 8,000 years ago. We find evidence for a recent
incomplete global replacement of ancestral Wolbachia and mtDNA lineages, which
is likely to be one of several similar incomplete replacement events that have
occurred since the out-of-Africa migration that allowed D. melanogaster to
colonize worldwide habitats.


  Bidirectional cargo transport by molecular motors in cells is a complex
phenomenon, in which the cargo (usually a vesicle) alternately moves in
retrograde and anterograde directions. In this case, teams of oppositely
pulling motors (eg., kinesin and dynein) bind to the cargo simultaneously, and
`coordinate' their activity such that the motion consists of spells of
positively and negatively directed segments, separated by pauses of varying
duration. A set of recent experiments have analyzed the bidirectional motion of
endosomes in the amoeba D. discoideum in detail. It was found that in between
directional switches, a team of 5-6 dyneins stall a cargo against a stronger
kinesin in tug of war, which lasts for almost a second. As the mean detachment
time of a kinesin under its stall load was also observed to be ~ 1s, we infer
that the collective detachment time of the dynein assembly must also be
similar. Here, we analyze this inference from a modeling perspective, using
experimentally measured single-molecule parameters as inputs. We find that the
commonly assumed exponential load-dependent detachment rate is inconsistent
with observations, as it predicts that a 5-dynein assembly will detach under
its combined stall load in less than a hundredth of a second. A modified model
where the load-dependent unbinding rate is assumed to saturate at stall-force
level for super-stall loads gives results which are in agreement with
experimental data. Our analysis suggests that the load-dependent detachment of
a dynein in a team is qualitatively different at sub-stall and super-stall
loads, a conclusion which is likely to have implications in other situations
involving collective effects of many motors.


  Nucleosome organization in eukaryotic genomes has a deep impact on gene
function. Although progress has been recently made in the identification of
various concurring factors influencing nucleosome positioning, it is still
unclear whether nucleosome positions are sequence dictated or determined by a
random process. It has been postulated for a long time that,in the proximity of
TSS, a barrier determines the position of the +1 nucleosome and then geometric
constraints alter the random positioning process determining nucleosomal
phasing. Such a pattern fades out as one moves away from the barrier to become
again a random positioning process. Although this statistical model is widely
accepted,the molecular nature of the barrier is still unknown. Moreover,we are
far from the identification of a set of sequence rules able:to account for the
genome-wide nucleosome organization;to explain the nature of the barriers on
which the statistical mechanism hinges;to allow for a smooth transition from
sequence-dictated to statistical positioning and back. We show that sequence
complexity,quantified via various methods, can be the rule able to at least
partially account for all the above.In particular, we have conducted our
analyses on 4 high resolution nucleosomal maps of the model eukaryotes and
found that nucleosome depleted regions can be well distinguished from
nucleosome enriched regions by sequence complexity measures.In particular, (a)
the depleted regions are less complex than the enriched ones, (b) around TSS
complexity measures alone are in striking agreement with in vivo nucleosome
occupancy,in particular precisely indicating the positions of the +1 and -1
nucleosomes. Those findings indicate that the intrinsic richness of
subsequences within sequences plays a role in nucleosomal formation in genomes,
and that sequence complexity constitutes the molecular nature of nucleosome
barrier.


  We attempt to set a mathematical foundation of immunology and amino acid
chains. To measure the similarities of these chains, a kernel on strings is
defined using only the sequence of the chains and a good amino acid
substitution matrix (e.g. BLOSUM62). The kernel is used in learning machines to
predict binding affinities of peptides to human leukocyte antigens DR (HLA-DR)
molecules. On both fixed allele (Nielsen and Lund 2009) and pan-allele (Nielsen
et.al. 2010) benchmark databases, our algorithm achieves the state-of-the-art
performance. The kernel is also used to define a distance on an HLA-DR allele
set based on which a clustering analysis precisely recovers the serotype
classifications assigned by WHO (Nielsen and Lund 2009, and Marsh et.al. 2010).
These results suggest that our kernel relates well the chain structure of both
peptides and HLA-DR molecules to their biological functions, and that it offers
a simple, powerful and promising methodology to immunology and amino acid chain
studies.


  We study the evolutionary dynamics of an asexual population of nonmutators
and mutators on a class of epistatic fitness landscapes. We consider the
situation in which all mutations are deleterious and mutators are produced from
nonmutators continually at a constant rate. We find that in an infinitely large
population, a minimum nonmutator-to-mutator conversion rate is required to fix
the mutators but an arbitrarily small conversion rate results in the fixation
of mutators in a finite population. We calculate analytical expressions for the
mutator fraction at mutation-selection balance and fixation time for mutators
in a finite population when mutational effects are weaker (regime I) and
stronger (regime II) than the selective effects. Our main result is that in
regime I, the mutator fraction and the fixation time are independent of
epistasis but in regime II, mutators are rarer and take longer to fix when the
decrease in fitness with the number of deleterious mutations occurs at an
accelerating rate (synergistic epistasis) than at a diminishing rate
(antagonistic epistasis). Our analytical results are compared with numerics and
their implications are discussed.


  A present challenge in complex systems is to identify mesoscale structures
that have distinct dynamical implications. In this paper we present a detailed
investigation of a previously observed dynamical equivalence of certian
ecological food webs. We show that this equivalence is rooted in mesoscale
symmetries that exist in these webs. Certain eigenvectors of the Jacobian
describing dynamical modes of the system, such as specific instabilities or
responses to perturbations, localize on these symmetric motifs. On the one hand
this means that by removing a symmetry from the network one obtains a system
which has identical dynamics except for the removal of the localized mode. This
explains the previously observed equivalence. On the other hand it means that
we can identify dynamical modes that only depend on the symmetric motif.
Symmetric structures thus provide an example for mesoscale network motifs
having distinct and exact implications for the dynamics.


  We investigate the population dynamics in generalized Rock-Paper-Scissors
models with an arbitrary number of species $N$. We show, for the first time,
that spiral patterns with $N$-arms may develop both for odd and even $N$, in
particular in models where a bidirectional predation interaction of equal
strength between all species is modified to include one N-cyclic predator-prey
rule. While the former case gives rise to an interface network with Y-type
junctions obeying the scaling law $L \propto t^{1/2}$, where $L$ is the
characteristic length of the network and $t$ is the time, the later can lead to
a population network with $N$-armed spiral patterns, having a roughly constant
characteristic length scale. We explicitly demonstrate the connection between
interface junctions and spiral patterns in these models and compute the
corresponding scaling laws. This work significantly extends the results of
previous studies of population dynamics and could have profound implications
for the understanding of biological complexity in systems with a large number
of species.


  We present a phenomenological dynamical model able to describe the stretching
features of the curve of DNA length vs applied force. As concerns the chain,
the model is based on the discrete wormlike chain model with elastic
modifications, which properly describes the elongation features at low and
intermediate forces. The dynamics is developed under a double-well potential
with a linear term, which, at high forces, accounts for the narrow transition
present in the DNA elongation (overstretching). A quite good agreement between
simulation and experiment is obtained.


  This paper addresses the question of optimal phenotypic plasticity as a
response to environmental fluctuations while optimizing the cost/benefit ratio,
where the cost is energetic expense of plasticity, and benefit is fitness. The
dispersion matrix \Sigma of the genes' response (H = ln|\Sigma|) is used: (i)
in a numerical model as a metric of the phenotypic variance reduction in the
course of fitness optimization, then (ii) in an analytical model, in order to
optimize parameters under the constraint of limited energy availability.
Results lead to speculate that such optimized organisms should maximize their
exergy and thus the direct/indirect work they exert on the habitat. It is shown
that the optimal cost/benefit ratio belongs to an interval in which differences
between individuals should not substantially modify their fitness.
Consequently, even in the case of an ideal population, close to the optimal
plasticity, a certain level of genetic diversity should be long conserved, and
a part, still to be determined, of intra-populations genetic diversity probably
stem from environment fluctuations. Species confronted to monotonous factors
should be less plastic than vicariant species experiencing heterogeneous
environments. Analogies with the MaxEnt algorithm of E.T. Jaynes (1957) are
discussed, leading to the conjecture that this method may be applied even in
case of multivariate but non multinormal distributions of the responses.


  It is well known since the pioneering work of Goldbeter and Koshland [Proc.
Natl. Acad. Sci. USA, vol. 78, pp. 6840-6844 (1981)] that cellular
phosphorylation- dephosphorylation cycle (PdPC), catalyzed by kinase and
phosphatase under saturated condition with zeroth order enzyme kinetics,
exhibits ultrasensitivity, sharp transition. We analyse the dynamics aspects of
the zeroth order PdPC kinetics and show a critical slowdown akin to the phase
transition in condensed matter physics. We demonstrate that an extremely
simple, though somewhat mathematically "singular" model is a faithful
representation of the ultrasentivity phenomenon. The simplified mathematical
model will be valuable, as a component, in developing complex cellular
signaling network theory as well as having a pedagogic value.


  A multi-labeled tree, or MUL-tree, is a phylogenetic tree where two or more
leaves share a label, e.g., a species name. A MUL-tree can imply multiple
conflicting phylogenetic relationships for the same set of taxa, but can also
contain conflict-free information that is of interest and yet is not obvious.
We define the information content of a MUL-tree T as the set of all
conflict-free quartet topologies implied by T, and define the maximal reduced
form of T as the smallest tree that can be obtained from T by pruning leaves
and contracting edges while retaining the same information content. We show
that any two MUL-trees with the same information content exhibit the same
reduced form. This introduces an equivalence relation in MUL-trees with
potential applications to comparing MUL-trees. We present an efficient
algorithm to reduce a MUL-tree to its maximally reduced form and evaluate its
performance on empirical datasets in terms of both quality of the reduced tree
and the degree of data reduction achieved.


  We study the stochastic evolution of four species in cyclic competition in a
well mixed environment. In systems composed of a finite number $N$ of particles
these simple interaction rules result in a rich variety of extinction
scenarios, from single species domination to coexistence between
non-interacting species. Using exact results and numerical simulations we
discuss the temporal evolution of the system for different values of $N$, for
different values of the reaction rates, as well as for different initial
conditions. As expected, the stochastic evolution is found to closely follow
the mean-field result for large $N$, with notable deviations appearing in
proximity of extinction events. Different ways of characterizing and predicting
extinction events are discussed.


  Intransitivity is a property of connected, oriented graphs representing
species interactions that may drive their coexistence even in the presence of
competition, the standard example being the three species Rock-Paper-Scissors
game. We consider here a generalization with four species, the minimum number
of species allowing other interactions beyond the single loop (one predator,
one prey). We show that, contrary to the mean field prediction, on a square
lattice the model presents a transition, as the parameter setting the rate at
which one species invades another changes, from a coexistence to a state in
which one species gets extinct. Such a dependence on the invasion rates shows
that the interaction graph structure alone is not enough to predict the outcome
of such models. In addition, different invasion rates permit to tune the level
of transitiveness, indicating that for the coexistence of all species to
persist, there must be a minimum amount of intransitivity.


  The composition of narrow distribution poly ethylene oxide-propylene oxide
copolymer (Mw ~ 8700 Da) was studied using matrix assisted laser desorption
ionization (MALDI) mass spectrometry. The ethylene oxide-propylene oxide
copolymer produced oligomers separated by 14 Da. The average resolving power
over the entire spectrum was 28,000. Approximately 448 isotopically resolved
peaks representing about 56 oligomers are identified. Although agreement
between experimental and calculated isotopic distributions was strong, the
compositional assignment was difficult. This is due to the large number of
possible isobaric components. The purpose of this research is to resolve and
study the composition of high mass copolymer such as ethylene oxide-propylene
oxide.


  Data on the number of Open Reading Frames (ORFs) coded by genomes from the 3
domains of Life show some notable general features including essential
differences between the Prokaryotes and Eukaryotes, with the number of ORFs
growing linearly with total genome size for the former, but only
logarithmically for the latter. Assuming that the (protein) coding and
non-coding fractions of the genome must have different dynamics and that the
non-coding fraction must be controlled by a variety of (unspecified)
probability distribution functions, we are able to predict that the number of
ORFs for Eukaryotes follows a Benford distribution and has a specific
logarithmic form. Using the data for 1000+ genomes available to us in early
2010, we find excellent fits to the data over several orders of magnitude, in
the linear regime for the Prokaryote data, and the full non-linear form for the
Eukaryote data. In their region of overlap the salient features are
statistically congruent, which allows us to: interpret the difference between
Prokaryotes and Eukaryotes as the manifestation of the increased demand in the
biological functions required for the larger Eukaryotes, estimate some minimal
genome sizes, and predict a maximal Prokaryote genome size on the order of 8-12
megabasepairs. These results naturally allow a mathematical interpretation in
terms of maximal entropy and, therefore, most efficient information
transmission.


  We have studied the effects of an external sinusoidal force in protein
folding kinetics. The externally applied force field acts on the each amino
acid residues of polypeptide chains. Our simulation results show that mean
protein folding time first increases with driving frequency and then decreases
passing through a maximum. With further increase of the driving frequency the
mean folding time starts increasing as the noise-induced hoping event (from the
denatured state to the native state) begins to experience many oscillations
over the mean barrier crossing time period. Thus unlike one-dimensional barrier
crossing problems, the external oscillating force field induces both
\emph{stabilization or destabilization of the denatured state} of a protein. We
have also studied the parametric dependence of the folding dynamics on
temperature, viscosity, non-Markovian character of bath in presence of the
external field.


  Force extension curves (FECs), which quantify the response of a variety of
biomolecules subject to mechanical force ($f$), are often quantitatively fit
using worm-like chain (WLC) or freely-jointed chain (FJC) models. These models
predict that the chain extension, $x$, normalized by the contour length
increases linearly at small $f$ and at high forces scale as $x \sim (1 -
f^{-\alpha})$ where $\alpha$= 0.5 for WLC and unity for FJC. In contrast,
experiments on ssDNA show that over a range of $f$ and ionic concentration, $x$
scales as $x\sim\ln f$, which cannot be explained using WLC or FJC models.
Using theory and simulations we show that this unusual behavior in FEC in ssDNA
is due to sequence-independent polyelectrolyte effects. We show that the $x\sim
\ln f$ arises because in the absence of force the tangent correlation function,
quantifying chain persistence, decays algebraically on length scales on the
order of the Debye length. Our theory, which is most appropriate for monovalent
salts, quantitatively fits the experimental data and further predicts that such
a regime is not discernible in double stranded DNA.


  A scheme for detection of abnormality in molecular nano-networks is proposed.
This is motivated by the fact that early diagnosis, classification and
detection of diseases such as cancer play a crucial role in their successful
treatment. The proposed nano-abnormality detection scheme (NADS) comprises of a
two-tier network of sensor nano-machines (SNMs) in the first tier and a data
gathering node (DGN) at the sink. The SNMs detect the presence of competitor
cells as abnormality that is captured by variations in parameters of a
nano-communications channel. In the second step, the SNMs transmit micro-scale
messages over a noisy micro communications channel (MCC) to the DGN, where a
decision is made upon fusing the received signals. The detection performance of
each SNM is analyzed by setting up a Neyman-Pearson test. Next, taking into
account the effect of the MCC, the overall performance of the proposed NADS is
quantified in terms of probabilities of misdetection and false alarm. A design
problem is formulated, when the optimized concentration of SNMs in a sample is
obtained for a high probability of detection and a limited probability of false
alarm.


  When performing an analysis on a collection of molecular sequences, it can be
convenient to reduce the number of sequences under consideration while
maintaining some characteristic of a larger collection of sequences. For
example, one may wish to select a subset of high-quality sequences that
represent the diversity of a larger collection of sequences. One may also wish
to specialize a large database of characterized "reference sequences" to a
smaller subset that is as close as possible on average to a collection of
"query sequences" of interest. Such a representative subset can be useful
whenever one wishes to find a set of reference sequences that is appropriate to
use for comparative analysis of environmentally-derived sequences, such as for
selecting "reference tree" sequences for phylogenetic placement of metagenomic
reads. In this paper we formalize these problems in terms of the minimization
of the Average Distance to the Closest Leaf (ADCL) and investigate algorithms
to perform the relevant minimization. We show that the greedy algorithm is not
effective, show that a variant of the Partitioning Among Medoids (PAM)
heuristic gets stuck in local minima, and develop an exact dynamic programming
approach. Using this exact program we note that the performance of PAM appears
to be good for simulated trees, and is faster than the exact algorithm for
small trees. On the other hand, the exact program gives solutions for all
numbers of leaves less than or equal to the given desired number of leaves,
while PAM only gives a solution for the pre-specified number of leaves. Via
application to real data, we show that the ADCL criterion chooses chimeric
sequences less often than random subsets, while the maximization of
phylogenetic diversity chooses them more often than random. These algorithms
have been implemented in publicly available software.


  Exploring the genetic basis of heritable traits remains one of the central
challenges in biomedical research. In simple cases, single polymorphic loci
explain a significant fraction of the phenotype variability. However, many
traits of interest appear to be subject to multifactorial control by groups of
genetic loci instead. Accurate detection of such multivariate associations is
nontrivial and often hindered by limited power. At the same time, confounding
influences such as population structure cause spurious association signals that
result in false positive findings if they are not accounted for in the model.
Here, we propose LMM-Lasso, a mixed model that allows for both, multi-locus
mapping and correction for confounding effects. Our approach is simple and free
of tuning parameters, effectively controls for population structure and scales
to genome-wide datasets. We show practical use in genome-wide association
studies and linkage mapping through retrospective analyses. In data from
Arabidopsis thaliana and mouse, our method is able to find a genetic cause for
significantly greater fractions of phenotype variation in 91% of the phenotypes
considered. At the same time, our model dissects this variability into
components that result from individual SNP effects and population structure. In
addition to this increase of genetic heritability, enrichment of known
candidate genes suggests that the associations retrieved by LMM-Lasso are more
likely to be genuine.


  We study a stochastic epidemic model consisting of elements (organisms in a
community or cells in tissue) with fixed positions, in which damage or disease
is transmitted by diffusing agents ("signals") emitted by infected individuals.
The signals decay as well as diffuse; since they are assumed to be produced in
large numbers, the signal concentration is treated deterministically. The
model, which includes four cellular states (susceptible, transformed, depleted,
and removed), admits various interpretations: spread of an infection or
infectious disease, or of damage in a tissue in which injured cells may
themselves provoke further damage, and as a description of the so-called
radiation-induced bystander effect, in which the signals are molecules capable
of inducing cell damage and/or death in unirradiated cells. The model exhibits
a continuous phase transition between spreading and nonspreading phases. We
formulate two mean-field theory (MFT) descriptions of the model, one of which
ignores correlations between the cellular state and the signal concentration,
and another that treats such correlations in an approximate manner. Monte Carlo
simulations of the spread of infection on the square lattice yield values for
the critical exponents and the fractal dimension consistent with the dynamic
percolation universality class.


  During the last decade, network approaches became a powerful tool to describe
protein structure and dynamics. Here, we describe first the protein structure
networks of molecular chaperones, then characterize chaperone containing
sub-networks of interactomes called as chaperone-networks or chaperomes. We
review the role of molecular chaperones in short-term adaptation of cellular
networks in response to stress, and in long-term adaptation discussing their
putative functions in the regulation of evolvability. We provide a general
overview of possible network mechanisms of adaptation, learning and memory
formation. We propose that changes of network rigidity play a key role in
learning and memory formation processes. Flexible network topology provides
"learning competent" state. Here, networks may have much less modular
boundaries than locally rigid, highly modular networks, where the learnt
information has already been consolidated in a memory formation process. Since
modular boundaries are efficient filters of information, in the "learning
competent" state information filtering may be much smaller, than after memory
formation. This mechanism restricts high information transfer to the "learning
competent" state. After memory formation, modular boundary-induced segregation
and information filtering protect the stored information. The flexible networks
of young organisms are generally in a "learning competent" state. On the
contrary, locally rigid networks of old organisms have lost their "learning
competent" state, but store and protect their learnt information efficiently.
We anticipate that the above mechanism may operate at the level of both
protein-protein interaction and neuronal networks.


  The ABC (ATP Binding Cassette) transporter protein superfamily comprises a
large number of ubiquitous and functionally versatile proteins conserved from
archaea to humans. ABC transporters have a key role in many human diseases and
also in the development of multidrug resistance in cancer and in parasites.
Although a dramatic progress has been achieved in ABC protein studies in the
last decades, we are still far from a detailed understanding of their molecular
functions. Several aspects of pharmacological ABC transporter targeting also
remain unclear. Here we summarize the conformational and protonation changes of
ABC transporters and the potential use of this information in pharmacological
design. Network related methods, which recently became useful tools to describe
protein structure and dynamics, have not been applied to study allosteric
coupling in ABC proteins as yet. A detailed description of the strengths and
limitations of these methods is given, and their potential use in describing
ABC transporter dynamics is outlined. Finally, we highlight possible future
aspects of pharmacological utilization of network methods and outline the
future trends of this exciting field.


  In equilibrium systems with short-ranged interactions, the relative stability
of different thermodynamic states generally does not depend on system size (as
long as this size is larger than the interaction range). Here, we use a large
deviations approach to show that, in contrast, different states can exchange
stability as system size is varied in a driven, bistable reaction-diffusion
system. This striking effect is related to a shift from a spatially uniform to
a nonuniform transition state and should generically be possible in a wide
range of nonequilibrium physical and biological systems.


  We use the Allen Gene Expression Atlas (AGEA) and the OMA ortholog dataset to
investigate the evolution of mouse-brain neuroanatomy from the standpoint of
the molecular evolution of brain-specific genes. For each such gene, using the
phylogenetic tree for all fully sequenced species and the presence of orthologs
of the gene in these species, we construct and assign a discrete measure of
evolutionary age. The gene expression profile of all gene of similar age,
relative to the average gene expression profile, distinguish regions of the
brain that are over-represented in the corresponding evolutionary timescale. We
argue that the conclusions one can draw on evolution of twelve major brain
regions from such a molecular level analysis supplements existing knowledge of
mouse brain evolution and introduces new quantitative tools, especially for
comparative studies, when AGEA-like data sets for other species become
available. Using the functional role of the genes representational of a certain
evolutionary timescale and brain region we compare and contrast, wherever
possible, our observations with existing knowledge in evolutionary
neuroanatomy.


  A further generalization of the stochastic replicator dynamic derived by
Fudenberg and Harris \cite{FH92} is considered. In particular, a Poissonian
integral is introduced to the fitness to simulate the affects of anomalous
events. For the two strategy population, an estimation of the long run behavior
of the dynamic is derived. For the population with many strategies, conditions
for stability to pure strict Nash equilibria, extinction of dominated pure
strategies, and recurrence in a neighborhood of an internal evolutionary stable
strategy are derived. This extends the results given by Imhof \cite{I05}.


  Malaria parasites are known to mediate the induction of inflammatory immune
responses at the maternal-foetal interface during placental malaria (PM)
leading to adverse consequences like pre-term deliveries and abortions.
Immunological events that take place within the malaria-infected placental
micro-environment leading to retarded foetal growth and disruption of
pregnancies are among the critical parameters that are still in need of further
elucidation. The establishment of more animal models for studying placental
malaria can provide novel ways of circumventing problems experienced during
placental malaria research in humans such as inaccurate estimation of
gestational ages. Using the newly established olive baboon (Papio
anubis)-Plasmodium knowlesi (P. knowlesi) H strain model of placental malaria,
experiments were carried out to determine placental cytokine profiles
underlying the immunopathogenesis of placental malaria. Four pregnant olive
baboons were infected with blood stage P. knowlesi H strain parasites on the
one fiftieth day of gestation while four other uninfected pregnant olive
baboons were maintained as uninfected controls. After nine days of infection,
placentas were extracted from all the eight baboons through cesarean surgery
and used for the processing of placental plasma and sera samples for cytokine
sandwich enzyme linked immunosorbent assays (ELISA). Results indicated that the
occurrence of placental malaria was associated with elevated concentrations of
tumour necrosis factor alpha (TNF-{\alpha}) and interleukin 12 (IL-12).
Increased levels of IL-4, IL-6 and IL-10 and interferon gamma (IFN-{\gamma})
levels were detected in uninfected placentas. These findings match previous
reports regarding immunity during PM thereby demonstrating the reliability of
the olive baboon-P. knowlesi model for use in further studies.


  One of the challenges in analyzing high-dimensional expression data is the
detection of important biological signals. A common approach is to apply a
dimension reduction method, such as principal component analysis. Typically,
after application of such a method the data is projected and visualized in the
new coordinate system, using scatter plots or profile plots. These methods
provide good results if the data have certain properties which become visible
in the new coordinate system and which were hard to detect in the original
coordinate system. Often however, the application of only one method does not
suffice to capture all important signals. Therefore several methods addressing
different aspects of the data need to be applied. We have developed a framework
for linear and non-linear dimension reduction methods within our visual
analytics pipeline SpRay. This includes measures that assist the interpretation
of the factorization result. Different visualizations of these measures can be
combined with functional annotations that support the interpretation of the
results. We show an application to high-resolution time series microarray data
in the antibiotic-producing organism Streptomyces coelicolor as well as to
microarray data measuring expression of cells with normal karyotype and cells
with trisomies of human chromosomes 13 and 21.


  Background: Nowadays, the reconstruction of genome scale metabolic models is
a non-automatized and interactive process based on decision taking. This
lengthy process usually requires a full year of one person's work in order to
satisfactory collect, analyze and validate the list of all metabolic reactions
present in a specific organism. In order to write this list, one manually has
to go through a huge amount of genomic, metabolomic and physiological
information. Currently, there is no optimal algorithm that allows one to
automatically go through all this information and generate the models taking
into account probabilistic criteria of unicity and completeness that a
biologist would consider. Results: This work presents the automation of a
methodology for the reconstruction of genome scale metabolic models for any
organism. The methodology that follows is the automatized version of the steps
implemented manually for the reconstruction of the genome scale metabolic model
of a photosynthetic organism, {\it Synechocystis sp. PCC6803}. The steps for
the reconstruction are implemented in a computational platform (COPABI) that
generates the models from the probabilistic algorithms that have been
developed. Conclusions: For validation of the developed algorithm robustness,
the metabolic models of several organisms generated by the platform have been
studied together with published models that have been manually curated. Network
properties of the models like connectivity and average shortest mean path of
the different models have been compared and analyzed.


  The metabolic network of a living cell involves several hundreds or thousands
of interconnected biochemical reactions. Previous research has shown that under
realistic conditions only a fraction of these reactions is concurrently active
in any given cell. This is partially determined by nutrient availability, but
is also strongly dependent on the metabolic function and network structure.
Here, we establish rigorous bounds showing that the fraction of active
reactions is smaller (rather than larger) in metabolic networks evolved or
engineered to optimize a specific metabolic task, and we show that this is
largely determined by the presence of thermodynamically irreversible reactions
in the network. We also show that the inactivation of a certain number of
reactions determined by irreversibility can generate a cascade of secondary
reaction inactivations that propagates through the network. The mathematical
results are complemented with numerical simulations of the metabolic networks
of the bacterium Escherichia coli and of human cells, which show,
counterintuitively, that even the maximization of the total reaction flux in
the network leads to a reduced number of active reactions.