  Standard models of ion channel voltage gating require substantial movement of
one transmembrane segment, S4, of the voltage sensing domain. Evidence comes
from the accessibility to external methanethiosulfonate (MTS) reagents of the
positively charged arginine residues (R) on this segment. These are first
mutated to cysteines (C), which in turn react with MTS reagents; it is assumed
that the C is passively carried in the S4 movement, becoming accessible on one
side or the other of the membrane. However, the Rs were salt bridged to
negatively charged residues on other transmembrane segments, or hydrogen
bonded, while C reacts as a negative ion. The space available for MTS is fairly
close to the difference in volume between the large R residue and much smaller
C, so the MTS is not severely sterically hindered. A reagent molecule can reach
a C side chain; the C can react if not repelled by a negative charge from the
amino acid to which the R had been salt bridged. Nearby protons may also make
reaction possible unless the C itself is protonated. Therefore interpretation
of the C substitution results requires reconsideration. To test the idea we
have done quantum calculations on part of a mutated S4 and the nearby section
of the channel. The mutation is R300C of the 2A79/3Lut structure, a mutation
that would be done to test MTS reagent access; there is a large cavity where
the R is replaced by C. Two quantum calculations show a substantial difference
in the structure of this cavity with 2 water molecules compared to 4. This
suggests that the structure, and presumably reaction probability, could depend
on water molecules, very likely also protons, in or near the cavity that the
R300C mutation produces.


  The generation of two non-identical membrane compartments via exchange of
vesicles is considered to require two types of vesicles specified by distinct
cytosolic coats that selectively recruit cargo and two membrane-bound SNARE
pairs that specify fusion and differ in their affinities for each type of
vesicles. The mammalian Golgi complex is composed of 6-8 non-identical
cisternae that undergo gradual maturation and replacement yet features only two
SNARE pairs. We present a model that explains how the distinct composition of
Golgi cisternae can be generated with two and even a single SNARE pair and one
vesicle coat. A decay of active SNARE concentration in aging cisternae provides
the seed for a cis > trans SNARE gradient that generates the predominantly
retrograde vesicle flux which further enhances the gradient. This flux in turn
yields the observed inhomogeneous steady-state distribution of Golgi enzymes,
which compete with each other and with the SNAREs for incorporation into
transport vesicles. We show analytically that the steady state SNARE
concentration decays exponentially with the cisterna number. Numerical
solutions of rate equations reproduce the experimentally observed SNARE
gradients, overlapping enzyme peaks in cis, medial and trans and the reported
change in vesicle nature across Golgi: Vesicles originating from younger
cisternae mostly contain Golgi enzymes and SNAREs enriched in these cisternae
and extensively recycle through the Endoplasmic Reticulum (ER), while the other
subpopulation of vesicles contains Golgi proteins prevalent in older cisternae
and hardly reaches the ER.


  Mathematical modeling of biological networks can help to integrate a large
body of information into a consistent framework, which can then be used to
arrive at novel mechanistic insight and predictions. We have previously
developed a detailed, mechanistic model for the regulation of {\sigma}F during
sporulation in Bacillus subtilis. The model was based on a wide range of
quantitative data, and once fitted to the data, the model made predictions that
could be confirmed in experiments. However, the analysis was based on a single
optimal parameter set. We wondered whether the predictions of the model would
be stable for all optimal parameter sets. To that end we conducted a global
parameter screen within the physiological parameter ranges. The screening
approach allowed us to identify sensitive and sloppy parameters, and
highlighted further required datasets during the optimization. Eventually, all
parameter sets that reproduced all available data predicted the physiological
situation correctly.


  Published evidence for the human-mediated extinction of megafauna is examined
and is found to be unsubstantiated. It is shown that the claimed evidence is
not based on data describing the growth of human population but on the
fabricated data. However, even these fabricated data, which were claimed to
support the human-induced extinction of megafauna, contradict this claim. The
belief in the human-induced extinction of megafauna appears to be so strong
that even contradicting evidence based on the fabricated data is interpreted as
the evidence supporting this belief.


  We argue that the stochastic dynamics of interacting agents which replicate,
mutate and die constitutes a non-equilibrium physical process akin to aging in
complex materials. Specifically, our study uses extensive computer simulations
of the Tangled Nature Model (TNM) of biological evolution to show that
punctuated equilibria successively generated by the model's dynamics have
increasing entropy and are separated by increasing entropic barriers. We
further show that these states are organized in a hierarchy and that limiting
the values of possible interactions to a finite interval leads to stationary
fluctuations within a component of the latter. A coarse-grained description
based on the temporal statistics of quakes, the events leading from one
component of the hierarchy to the next, accounts for the logarithmic growth of
the population and the decaying rate of change of macroscopic variables.
Finally, we question the role of fitness in large scale evolution models and
speculate on the possible evolutionary role of rejuvenation and memory effects.


  OPTICS is a density-based clustering algorithm that performs well in a wide
variety of applications. For a set of input objects, the algorithm creates a
so-called reachability plot that can be either used to produce cluster
membership assignments, or interpreted itself as an expressive two-dimensional
representation of the density-based clustering structure of the input set, even
if the input set is embedded in higher dimensions. The main focus of this work
is a visualization method that can be used to assign colours to all entries of
the input database, based on hierarchically represented a-priori knowledge
available for each of these objects. Based on two different,
bioinformatics-related applications we illustrate how the proposed method can
be efficiently used to identify clusters with proven real-life relevance.


  It is demonstrated earlier that the exact Smith-Waterman algorithm yields
more accurate results than the members of the heuristic BLAST family of
algorithms. Unfortunately, the Smith-Waterman algorithm is much slower than the
BLAST and its clones. Here we present a technique and a webserver that uses the
exact Smith-Waterman algorithm, and it is approximately as fast as the BLAST
algorithm. The technique unites earlier methods of extensive preprocessing of
the target sequence database, and CPU-specific coding of the Smith-Waterman
algorithm. The SwissAlign webserver is available at the
http://swissalign.pitgroup.org address.


  Motivated by recent understandings in the stochastic natures of gene
expression, biochemical signaling, and spontaneous reversible epigenetic
switchings, we study a simple deterministic cell population dynamics in which
subpopulations grow with different rates and individual cells can
bi-directionally switch between a small number of different epigenetic
phenotypes. Two theories in the past, the population dynamics and
thermodynamics of master equations, separatedly defined two important concepts
in mathematical terms: the {\em fitness} in the former and the (non-adiabatic)
{\em entropy production} in the latter. Both play important roles in the
evolution of the cell population dynamics. The switching sustains the
variations among the subpopulation growth thus continuous natural selection. As
a form of Price's equation, the fitness increases with ($i$) natural selection
through variations and $(ii)$ a positive covariance between the per capita
growth and switching, which represents a Lamarchian-like behavior. A negative
covariance balances the natural selection in a fitness steady state | "the red
queen" scenario. At the same time the growth keeps the proportions of
subpopulations away from the "intrinsic" switching equilibrium of individual
cells, thus leads to a continous entropy production. A covariance, between the
per capita growth rate and the "chemical potential" of subpopulation,
counter-acts the entropy production. Analytical results are obtained for the
limiting cases of growth dominating switching and vice versa.


  Background: Qualitative frameworks, especially those based on the logical
discrete formalism, are increasingly used to model regulatory and signalling
networks. A major advantage of these frameworks is that they do not require
precise quantitative data, and that they are well-suited for studies of large
networks. While numerous groups have developed specific computational tools
that provide original methods to analyse qualitative models, a standard format
to exchange qualitative models has been missing.
  Results: We present the System Biology Markup Language (SBML) Qualitative
Models Package ("qual"), an extension of the SBML Level 3 standard designed for
computer representation of qualitative models of biological networks. We
demonstrate the interoperability of models via SBML qual through the analysis
of a specific signalling network by three independent software tools.
Furthermore, the cooperative development of the SBML qual format paved the way
for the development of LogicalModel, an open-source model library, which will
facilitate the adoption of the format as well as the collaborative development
of algorithms to analyze qualitative models.
  Conclusion: SBML qual allows the exchange of qualitative models among a
number of complementary software tools. SBML qual has the potential to promote
collaborative work on the development of novel computational approaches, as
well as on the specification and the analysis of comprehensive qualitative
models of regulatory and signalling networks.


  Apoptotic cell death is coordinated through two distinct (type 1 and type 2)
intracellular signaling pathways. How the type 1/type 2 choice is made remains
a fundamental problem in the biology of apoptosis and has implications for
apoptosis related diseases and therapy. We study the problem of type 1/type 2
choice in silico utilizing a kinetic Monte Carlo model of cell death signaling.
Our results show that the type 1/type 2 choice is linked to deterministic
versus stochastic cell death activation, elucidating a unique regulatory
control of the apoptotic pathways. Consistent with previous findings, our
results indicate that caspase 8 activation level is a key regulator of the
choice between deterministic type 1 and stochastic type 2 pathways,
irrespective of cell types. Expression levels of signaling molecules downstream
also regulate the type 1/type 2 choice. A simplified model of DISC clustering
elucidates the mechanism of increased active caspase 8 generation, and type 1
activation, in cancer cells having increased sensitivity to death receptor
activation. We demonstrate that rapid deterministic activation of the type 1
pathway can selectively target those cancer cells, especially if XIAP is also
inhibited; while inherent cell-to-cell variability would allow normal cells
stay protected.


  Assessing the correctness of genome assemblies is an important step in any
genome project. Several methods exist, but most are computationally intensive
and, in some cases, inappropriate. Here I present baa.pl, a fast and
easy-to-use program that uses transcript data to evaluate genomic assemblies.
Through simulations using human chromosome 22, I show that baa.pl excels at
detecting levels of missing sequence and contiguity. The program is freely
available at: https://github.com/josephryan/baa.pl


  Phospholipids are essential building blocks of biological membranes. Despite
of vast amount of accurate experimental data the atomistic resolution
structures sampled by the glycerol backbone and choline headgroup in
phoshatidylcholine bilayers are not known. Atomistic resolution molecular
dynamics simulation model would automatically resolve the structures giving an
interpretation of experimental results, if the model would reproduce the
experimental data. In this work we compare the C-H bond vector order parameters
for glycerol backbone and choline headgroup between 14 different atomistic
resolution models and experiments in fully hydrated lipid bilayer. The current
models are not accurately enough to resolve the structure. However, closer
inspection of three best performing models (CHARMM36, GAFFlipid and MacRog)
suggest that improvements in the sampled dihedral angle distributions would
potentilly lead to the model which would resolve the structure. Despite of the
inaccuracy in the fully hydrated structures, the response to the dehydration,
i.e. P-N vector tilting more parallel to membrane normal, is qualitatively
correct in all models. The CHARMM36 and MacRog models describe the interactions
between lipids and cholesterol better than Berger/H\"oltje model. This work has
been, and continues to be, progressed and discussed through the blog:
nmrlipids.blogspot.fi. Everyone is invited to join the discussion and make
contributions through the blog. The manuscript will be eventually submitted to
an appropriate scientific journal. Everyone who has contributed to the work
through the blog will be offered coauthorship. For more details see:
nmrlipids.blogspot.fi.


  Two cases of evolutionary stable strategy with periodic payoffs are studied.
The first is a generalization of Uyttendaele et al. The second is prisoner's
dilemma with periodic payoff. It is shown that reducing the defection payoff by
a periodic term is sufficient to introduce cooperation provided that the
initial fraction of cooperators is greater than 0.6.


  Morphogenesis is a tightly regulated process that has been studied for
decades. We are developing data-based and image-basd mechanistic models for a
range of developmental processes with a view to integrate the available
knowledge and to better understand the underlying regulatory logic. In our
previous papers on simulating organogenesis with COMSOL (German et al COMSOL
Conf Procedings 2011; Menshykau and Iber, COMSOL Conf Proceedings 2012) we
discussed methods to efficiently solve such models on static and growing
domains. A further challenge in modeling morphogenesis is the parameterization
of such models. Here we discuss COMSOL-based methods for parameter
optimization. These routines can be used to determine parameter sets, for which
the simulations reproduce experimental data and constraints. Such data is often
image based, but may also come from classical biochemical or genetic
experiments.


  Ortholog detection (OD) is a critical step for comparative genomic analysis
of protein-coding sequences. In this paper, we begin with a comprehensive
comparison of four popular, methodologically diverse OD methods: MultiParanoid,
Blat, Multiz, and OMA. In head-to-head comparisons, these methods are shown to
significantly outperform one another 12-30% of the time. This high
complementarity motivates the presentation of the first tool for integrating
methodologically diverse OD methods. We term this program MOSAIC, or Multiple
Orthologous Sequence Analysis and Integration by Cluster optimization. Relative
to component and competing methods, we demonstrate that MOSAIC more than
quintuples the number of alignments for which all species are present, while
simultaneously maintaining or improving functional-, phylogenetic-, and
sequence identity-based measures of ortholog quality. Further, we demonstrate
that this improvement in alignment quality yields 40-280% more confidently
aligned sites. Combined, these factors translate to higher estimated levels of
overall conservation, while at the same time allowing for the detection of up
to 180% more positively selected sites. MOSAIC is available as python package.
MOSAIC alignments, source code, and full documentation are available at
http://pythonhosted.org/bio-MOSAIC.


  The stoichiometry of metabolic networks usually gives rise to a family of
conservation laws for the aggregate concentration of specific pools of
metabolites, which not only constrain the dynamics of the network, but also
provide key insight into a cell's production capabilities. When the conserved
quantity identifies with a chemical moiety, extracting all such conservation
laws from the stoichiometry amounts to finding all integer solutions to an
NP-hard programming problem. Here we propose a novel and efficient
computational strategy that combines Monte Carlo, message passing, and
relaxation algorithms to compute the complete set of irreducible integer
conservation laws of a given stoichiometric matrix, also providing a
certificate for correctness and maximality of the solution. The method is
deployed for the analysis of the complete set of irreducible integer pools of
two large-scale reconstructions of the metabolism of the bacterium Escherichia
coli in different growth media. In addition, we uncover a scaling relation that
links the size of the irreducible pool basis to the number of metabolites, for
which we present an analytical explanation.


  Two-dimensional electrophoresis is still a very valuable tool in proteomics,
due to its reproducibility and its ability to analyze complete proteins.
However, due to its sensitivity to dynamic range issues, its most suitable use
in the frame of biomarker discovery is not on very complex fluids such as
plasma, but rather on more proximal, simpler fluids such as CSF, urine, or
secretome samples. Here, we describe the complete workflow for the analysis of
such dilute samples by two-dimensional electrophoresis, starting from sample
concentration, then the two-dimensional electrophoresis step per se, ending
with the protein detection by fluorescence.


  The comparative genomics revolution of the past decade has enabled the
discovery of functional elements in the human genome via sequence comparison.
While that is so, an important class of elements, those specific to humans, is
entirely missed by searching for sequence conservation across species. Here we
present an analysis based on variation data among human genomes that utilizes a
supervised machine learning approach for the identification of human specific
purifying selection in the genome. Using only allele frequency information from
the complete low coverage 1000 Genomes Project dataset in conjunction with a
support vector machine trained from known functional and non-functional
portions of the genome, we are able to accurately identify portions of the
genome constrained by purifying selection. Our method identifies previously
known human-specific gains or losses of function and uncovers many novel
candidates. Candidate targets for gain and loss of function along the human
lineage include numerous putative regulatory regions of genes essential for
normal development of the central nervous system, including a significant
enrichment of gain of function events near neurotransmitter receptor genes.
These results are consistent with regulatory turnover being a key mechanism in
the evolution of human-specific characteristics of brain development. Finally,
we show that the majority of the genome is unconstrained by natural selection
currently, in agreement with what has been estimated from phylogenetic methods
but in sharp contrast to estimates based on transcriptomics or other
high-throughput functional methods.


  This survey focuses on the most important aspects of the mathematical theory
of population genetic models of selection and migration between discrete
niches. Such models are most appropriate if the dispersal distance is short
compared to the scale at which the environment changes, or if the habitat is
fragmented. The general goal of such models is to study the influence of
population subdivision and gene flow among subpopulations on the amount and
pattern of genetic variation maintained. Only deterministic models are treated.
Because space is discrete, they are formulated in terms of systems of nonlinear
difference or differential equations. A central topic is the exploration of the
equilibrium and stability structure under various assumptions on the patterns
of selection and migration. Another important, closely related topic concerns
conditions (necessary or sufficient) for fully polymorphic (internal)
equilibria. First, the theory of one-locus models with two or multiple alleles
is laid out. Then, mostly very recent, developments about multilocus models are
presented. Finally, as an application, analysis and results of an explicit
two-locus model emerging from speciation theory are highlighted.


  Single-cell studies often show significant phenotypic variability due to the
stochastic nature of intra-cellular biochemical reactions. When the numbers of
molecules, e.g., transcription factors and regulatory enzymes, are in low
abundance, fluctuations in biochemical activities become significant and such
"noise" can propagate through regulatory cascades in terms of biochemical
reaction networks. Here we develop an intuitive, yet fully quantitative method
for analyzing how noise affects cellular phenotypes based on identifying a
system's nonlinearities and noise propagations. We observe that such noise can
simultaneously enhance sensitivities in one behavioral region while reducing
sensitivities in another. Employing this novel phenomenon we designed three
biochemical signal processing modules: (a) A gene regulatory network that acts
as a concentration detector with both enhanced amplitude and sensitivity. (b) A
non-cooperative positive feedback system, with a graded dose-response in the
deterministic case, that serves as a bistable switch due to noise-induced
bimodality. (c) A noise-induced linear amplifier for gene regulation that
requires no feedback. The methods developed in the present work allow one to
understand and engineer nonlinear biochemical signal processors based on
fluctuation-induced phenotypes.


  We consider a stochastic Susceptible-Exposed-Infected-Recovered (SEIR)
epidemiological model with a contact rate that fluctuates seasonally. Through
the use of a nonlinear, stochastic projection, we are able to analytically
determine the lower dimensional manifold on which the deterministic and
stochastic dynamics correctly interact. Our method produces a low dimensional
stochastic model that captures the same timing of disease outbreak and the same
amplitude and phase of recurrent behavior seen in the high dimensional model.
Given seasonal epidemic data consisting of the number of infectious
individuals, our method enables a data-based model prediction of the number of
unobserved exposed individuals over very long times.


  Spatial patterns in the early fruit fly embryo emerge from a network of
interactions among transcription factors, the gap genes, driven by maternal
inputs. Such networks can exhibit many qualitatively different behaviors,
separated by critical surfaces. At criticality, we should observe strong
correlations in the fluctuations of different genes around their mean
expression levels, a slowing of the dynamics along some but not all directions
in the space of possible expression levels, correlations of expression
fluctuations over long distances in the embryo, and departures from a Gaussian
distribution of these fluctuations. Analysis of recent experiments on the gap
genes shows that all these signatures are observed, and that the different
signatures are related in ways predicted by theory. While there might be other
explanations for these individual phenomena, the confluence of evidence
suggests that this genetic network is tuned to criticality.


  Dormancy is an essential strategy for microorganisms to cope with
environmental stress. However, global ecosystem models typically ignore
microbial dormancy, resulting in major model uncertainties. To facilitate the
consideration of dormancy in these large-scale models, we propose a new
microbial physiology component that works for a wide range of substrate
availabilities. This new model is based on microbial physiological states and
is majorly parameterized with the maximum specific growth and maintenance rates
of active microbes and the ratio of dormant to active maintenance rates. A
major improvement of our model over extant models is that it can explain the
low active microbial fractions commonly observed in undisturbed soils. Our new
model shows that the exponentially-increasing respiration from
substrate-induced respiration experiments can only be used to determine the
maximum specific growth rate and initial active microbial biomass, while the
respiration data representing both exponentially-increasing and
non-exponentially-increasing phases can robustly determine a range of key
parameters including the initial total live biomass, initial active fraction,
the maximum specific growth and maintenance rates, and the half-saturation
constant. Our new model can be incorporated into existing ecosystem models to
account for dormancy in microbially-mediated processes and to provide improved
estimates of microbial activities.


  This paper reviews and clarifies five misunderstandings about cultural
evolution identified by Henrich, Boyd, and Richerson (2008). First, cultural
representations are neither discrete nor continuous; they are distributed
across neurons that respond to microfeatures. This enables associations to be
made, and cultural change to be generated. Second, 'replicator dynamics' do not
ensure natural selection. The replicator notion does not capture the
distinction between actively interpreted self-assembly code and passively
copied self-description, which leads to a fundamental principle of natural
selection: inherited information is transmitted, whereas acquired information
is not. Third, this principle is violated in culture by the ubiquity of
acquired change. Moreover, biased transmission is less important to culture
than the creative processes by which novelty is generated. Fourth, there is no
objective basis for determining cultural fitness. Fifth, the necessity of
randomness is discussed. It is concluded that natural selection is
inappropriate as an explanatory framework for culture.


  A large and diverse array of small hydrophobic molecules induce general
anesthesia. Their efficacy as anesthetics has been shown to correlate both with
their affinity for a hydrophobic environment and with their potency in
inhibiting certain ligand gated ion channels. Here we explore the effects that
n-alcohols and other liquid anesthetics have on the two-dimensional miscibility
critical point observed in cell derived giant plasma membrane vesicles (GPMVs).
We show that anesthetics depress the critical temperature (Tc) of these GPMVs
without strongly altering the ratio of the two liquid phases found below Tc.
The magnitude of this affect is consistent across n-alcohols when their
concentration is rescaled by the median anesthetic concentration (AC50) for
tadpole anesthesia, but not when plotted against the overall concentration in
solution. At AC50 we see a 4{\deg}C downward shift in Tc, much larger than is
typically seen in the main chain transition at these anesthetic concentrations.
GPMV miscibility critical temperatures are also lowered to a similar extent by
propofol, phenylethanol, and isopropanol when added at anesthetic
concentrations, but not by tetradecanol or 2,6 diterbutylphenol, two structural
analogs of general anesthetics that are hydrophobic but have no anesthetic
potency. We propose that liquid general anesthetics provide an experimental
tool for lowering critical temperatures in plasma membranes of intact cells,
which we predict will reduce lipid-mediated heterogeneity in a way that is
complimentary to increasing or decreasing cholesterol. Also, several possible
implications of our results are discussed in the context of current models of
anesthetic action on ligand gated ion channels.


  The problem of targeted network immunization can be defined as the one of
finding a subset of nodes in a network to immunize or vaccinate in order to
minimize a tradeoff between the cost of vaccination and the final (stationary)
expected infection under a given epidemic model. Although computing the
expected infection is a hard computational problem, simple and efficient
mean-field approximations have been put forward in the literature in recent
years. The optimization problem can be recast into a constrained one in which
the constraints enforce local mean-field equations describing the average
stationary state of the epidemic process. For a wide class of epidemic models,
including the susceptible-infected-removed and the
susceptible-infected-susceptible models, we define a message-passing approach
to network immunization that allows us to study the statistical properties of
epidemic outbreaks in the presence of immunized nodes as well as to find
(nearly) optimal immunization sets for a given choice of parameters and costs.
The algorithm scales linearly with the size of the graph and it can be made
efficient even on large networks. We compare its performance with topologically
based heuristics, greedy methods, and simulated annealing.


  Current evolutionary theory suggests that many natural signaling systems
evolved from preexisting cues. In aposematic systems, prey warning signals
benefit both predator and prey. When the signal is highly beneficial, a third
species often evolves to mimic the unpalatable species, exploiting the
signaling system for its own protection. We investigated the evolutionary
development of predator cue utilization and prey signaling in a digital
predator-prey system in which mimicking prey could evolve to alter their
appearance to resemble poison-free or poisonous prey. In predators, we observed
rapid evolution of cue recognition (i.e. active behavioral responses) when
presented with sufficiently poisonous prey. In addition, active signaling (i.e.
mimicry) evolved in prey under all conditions that led to cue utilization. Thus
we show that despite imperfect and dishonest signaling, given a high cost of
consuming poisonous prey, complex systems of interspecific communication can
evolve via predator cue recognition and prey signal manipulation. This provides
evidence supporting hypotheses that cues may serve as stepping-stones in the
evolution of more advanced communication systems and signals incorporating
information about the environment.


  K-mer abundance analysis is widely used for many purposes in nucleotide
sequence analysis, including data preprocessing for de novo assembly, repeat
detection, and sequencing coverage estimation. We present the khmer software
package for fast and memory efficient online counting of k-mers in sequencing
data sets. Unlike previous methods based on data structures such as hash
tables, suffix arrays, and trie structures, khmer relies entirely on a simple
probabilistic data structure, a Count-Min Sketch. The Count-Min Sketch permits
online updating and retrieval of k-mer counts in memory which is necessary to
support online k-mer analysis algorithms. On sparse data sets this data
structure is considerably more memory efficient than any exact data structure.
In exchange, the use of a Count-Min Sketch introduces a systematic overcount
for k-mers; moreover, only the counts, and not the k-mers, are stored. Here we
analyze the speed, the memory usage, and the miscount rate of khmer for
generating k-mer frequency distributions and retrieving k-mer counts for
individual k-mers. We also compare the performance of khmer to several other
k-mer counting packages, including Tallymer, Jellyfish, BFCounter, DSK, KMC,
Turtle and KAnalyze. Finally, we examine the effectiveness of profiling
sequencing error, k-mer abundance trimming, and digital normalization of reads
in the context of high khmer false positive rates. khmer is implemented in C++
wrapped in a Python interface, offers a tested and robust API, and is freely
available under the BSD license at github.com/ged-lab/khmer.


  The application of phylogenetic techniques to the documentation of cultural
history can present a distorted picture due to horizontal transmission and
blending. Moreover, the units of cultural transmission must be communicable
concepts, rather than conveniently measurable attributes, and relatedness
between elements of culture often resides at the conceptual level, something
not captured by phylogenetic methods, which focus on measurable attributes.
(For example, mortars and pestles are as related as two artifacts could be,
despite little similarity at the attribute level.) This paper introduces a new,
cognitively inspired framework for chronicling material cultural history,
building on Lipo's (2005) network-based computational approach. We show that by
incorporating not just superficial attributes of artifact samples (e.g.
fluting) but also conceptual knowledge (e.g. information about function), a
different pattern of cultural ancestry emerges.


  This chapter begins by outlining a promising, new theoretical framework for
the process by which human culture evolves inspired by the views of complexity
theorists on the problem of how life began. Elements of culture, like species,
evolve over time; that is, they exhibit cumulative change that is adaptive in
nature. By studying how biological evolution got started, we gain insight into
not just the specifics of biological evolution, but also general insights into
the initiation of any evolutionary process that may be applicable to culture.
We then explore the implications of this new framework for culture on the
transformative processes of individuals. Specifically, we will address what
this emerging perspective on cultural evolution implies for to go about
attaining a sustainable worldview; that is, a web of habits, understandings,
and ways of approaching situations that is conducive to the development of a
sustainable world.


  Time-dependent distribution of the global extinction of megafauna is compared
with the growth of human population. There is no correlation between the two
processes. Furthermore, the size of human population and its growth rate were
far too small to have any significant impact on the environment and on the life
of megafauna.


  Molecular phylogenetic and phylogeographic reconstructions generally assume
time-homogeneous substitution processes. Motivated by computational
convenience, this assumption sacrifices biological realism and offers little
opportunity to uncover the temporal dynamics in evolutionary histories. Here,
we extend and generalize an evolutionary approach that relaxes the
time-homogeneous process assumption by allowing the specification of different
infinitesimal substitution rate matrices across different time intervals,
called epochs, along the evolutionary history. We focus on an epoch model
implementation in a Bayesian inference framework that offers great modeling
flexibility in drawing inference about any discrete data type characterized as
a continuous-time Markov chain, including phylogeographic traits. To alleviate
the computational burden that the additional temporal heterogeneity imposes, we
adopt a massively parallel approach that achieves both fine- and coarse-grain
parallelization of the computations across branches that accommodate epoch
transitions, making extensive use of graphics processing units. Through
synthetic examples, we assess model performance in recovering evolutionary
parameters from data generated according to different evolutionary scenarios
that comprise different numbers of epochs for both nucleotide and codon
substitution processes. We illustrate the usefulness of our inference framework
in two different applications to empirical data sets: the selection dynamics on
within-host HIV populations throughout infection and the seasonality of global
influenza circulation. In both cases, our epoch model captures key features of
temporal heterogeneity that remained difficult to test using ad hoc procedures.


  We consider an individual based model of phenotypic evolution in
hermaphroditic populations which includes random and assortative mating of
individuals. By increasing the number of individuals to infinity we obtain a
nonlinear transport equation, which describes the evolution of distribution of
phenotypic traits. Existence of an one-dimensional attractor is proved and the
formula for the density of phenotypic traits in the limiting (asymptotic)
population is derived in some particular case.


  Imputation using external reference panels is a widely used approach for
increasing power in GWAS and meta-analysis. Existing HMM-based imputation
approaches require individual-level genotypes. Here, we develop a new method
for Gaussian imputation from summary association statistics, a type of data
that is becoming widely available. In simulations using 1000 Genomes (1000G)
data, this method recovers 84% (54%) of the effective sample size for common
(>5%) and low-frequency (1-5%) variants (increasing to 87% (60%) when summary
LD information is available from target samples) versus 89% (67%) for HMM-based
imputation, which cannot be applied to summary statistics. Our approach
accounts for the limited sample size of the reference panel, a crucial step to
eliminate false-positive associations, and is computationally very fast. As an
empirical demonstration, we apply our method to 7 case-control phenotypes from
the WTCCC data and a study of height in the British 1958 birth cohort (1958BC).
Gaussian imputation from summary statistics recovers 95% (105%) of the
effective sample size (as quantified by the ratio of $\chi^2$ association
statistics) compared to HMM-based imputation from individual-level genotypes at
the 227 (176) published SNPs in the WTCCC (1958BC height) data. In addition,
for publicly available summary statistics from large meta-analyses of 4 lipid
traits, we publicly release imputed summary statistics at 1000G SNPs, which
could not have been obtained using previously published methods, and
demonstrate their accuracy by masking subsets of the data. We show that 1000G
imputation using our approach increases the magnitude and statistical evidence
of enrichment at genic vs. non-genic loci for these traits, as compared to an
analysis without 1000G imputation. Thus, imputation of summary statistics will
be a valuable tool in future functional enrichment analyses.


  Molecular traits, such as gene expression levels or protein binding
affinities, are increasingly accessible to quantitative measurement by modern
high-throughput techniques. Such traits measure molecular functions and, from
an evolutionary point of view, are important as targets of natural selection.
We review recent developments in evolutionary theory and experiments that are
expected to become building blocks of a quantitative genetics of molecular
traits. We focus on universal evolutionary characteristics: these are largely
independent of a trait's genetic basis, which is often at least partially
unknown. We show that universal measurements can be used to infer selection on
a quantitative trait, which determines its evolutionary mode of conservation or
adaptation. Furthermore, universality is closely linked to predictability of
trait evolution across lineages. We argue that universal trait statistics
extends over a range of cellular scales and opens new avenues of quantitative
evolutionary systems biology.


  It is well known that behavioral changes in contact patterns may
significantly affect the spread of an epidemic outbreak. Here we focus on
simple endemic models for recurrent epidemics, by modelling the social contact
rate as a function of the available information on the present and past disease
prevalence. We show that social behaviour change alone may trigger sustained
oscillations. This indicates that human behavior might be a critical explaining
factor of oscillations in time-series of endemic diseases. Finally, we briefly
show how the inclusion of seasonal variations in contacts may imply chaos.


  By combining high-density super-resolution imaging with a novel stochastic
analysis, we report here a peculiar nano-structure organization revealed by the
density function of individual AMPA receptors moving on the surface of cultured
hippocampal dendrites. High density regions of hundreds of nanometers for the
trajectories are associated with local molecular assembly generated by direct
molecular interactions due to physical potential wells. We found here that for
some of these regions, the potential wells are organized in ring structures. We
could find up to 3 wells in a single ring. Inside a ring receptors move in a
small band the width of which is of hundreds of nanometers. In addition, rings
are transient structures and can be observed for tens of minutes. Potential
wells located in a ring are also transient and the position of their peaks can
shift with time. We conclude that these rings can trap receptors in a unique
geometrical structure contributing to shape receptor trafficking, a process
that sustains synaptic transmission and plasticity.


  The theory of natural selection cannot describe how early life evolved, in
part because acquired characteristics are passed on through horizontal
exchange. It has been proposed that culture, like life, began with the
emergence of autopoietic form, thus its evolution too cannot be described by
natural selection. The evolution of autopoietic form can be described using a
framework referred to as Context-driven Actualization of Potential (CAP), which
grew out of a generalization of the formalisms of quantum mechanics, and
encompasses nondeterministic as well as deterministic change of state. The
autopoietic structure that evolves through culture is the mind, or more
accurately the conceptual network that yields an individual's internal model of
the world. A branch of CAP research referred to as the state-context-property
(SCOP) formalism provides a mathematical framework for reconciling the
stability of conceptual structure with its susceptibility to context-driven
change. The combination of two or more concepts (an extreme case of contextual
influence), as occurs in insight, is modeled as a state of entanglement.
Theoretical and empirical findings are presented that challenge assumptions
underlying virtually all of cognitive science, such as the notion of spreading
activation and the assumption that cognitive processes can be described with a
Kolmogorovian probability model.


  Conventionally in a PCR reaction, Primers binds to DNA template in an
antiparallel manner and template DNA is amplified as it is. Here we describe an
approach in which First primer binds in a complementary parallel orientation
leading to synthesis in parallel direction. Further reaction happened in usual
way leading to synthesis of final DNA product having opposite polarity then the
template used. We first time proved that synthesis of DNA can happen in a
parallel orientation. We first time have shown that from a single DNA template,
Two different related PCR products can be synthesized.


  The message passing approach of Karrer and Newman [Phys. Rev. E 82, 016101
(2010)] is an exact and practicable representation of
susceptible-infected-recovered dynamics on finite trees. Here we show that,
assuming Poisson contact processes, a pair-based moment closure representation
[Sharkey, J. Math. Biol. 57, 311 (2008)] can be derived from their equations.
We extend the applicability of both representations and discuss their relative
merits. On arbitrary time-independent networks, as was shown for the message
passing formalism, the pair-based moment closure equations also provide a
rigorous lower bound on the expected number of susceptibles at all times.


  Modelling genetic phenomena affecting biological traits is important for the
development of agriculture as it allows breeders to predict the potential of
breeding for certain traits. One such phenomenon is heterosis or hybrid vigor:
crossing individuals from genetically distinct populations often results in
improvements in quantitative traits, such as growth rate, biomass production
and stress resistance. Heterosis has become a very useful tool in global
agriculture, but its genetic basis remains controversial and its effects hard
to predict. We have taken a computational approach to studying heterosis,
developing a simulation of evolution, independent reassortment of alleles and
hybridization of Gene Regulatory Networks (GRNs) in a Boolean framework.
Fitness is measured as the ability of a network to respond to external inputs
in a pre-defined way. Our model reproduced common experimental observations on
heterosis using only biologically justified parameters. Hybrid vigor was
observed and its extent was seen to increase as parental populations diverged,
up until a point of sudden collapse of hybrid fitness. We also reproduce, for
the first time in a model, the fact that hybrid vigor cannot easily be fixed by
within a breeding line, currently an important limitation of the use of hybrid
crops. The simulation allowed us to study the effects of three standard models
for the genetic basis of heterosis and the level of detail in our model allows
us to suggest possible warning signs of the impending collapse of hybrid vigor
in breeding. In addition, the simulation provides a framework that can be
extended to study other aspects of heterosis and alternative evolutionary
scenarios.


  In recent years it is increasingly being recognized that biochemical signals
are not necessarily constant in time and that the temporal dynamics of a signal
can be the information carrier. Moreover, it is now well established that
components are often shared between signaling pathways. Here we show by
mathematical modeling that living cells can multiplex a constant and an
oscillatory signal: they can transmit these two signals through the same
signaling pathway simultaneously, and yet respond to them specifically and
reliably. We find that information transmission is reduced not only by noise
arising from the intrinsic stochasticity of biochemical reactions, but also by
crosstalk between the different channels. Yet, under biologically relevant
conditions more than 2 bits of information can be transmitted per channel, even
when the two signals are transmitted simultaneously. These observations suggest
that oscillatory signals are ideal for multiplexing signals.


  Consider a uniformly mixing population which grows as a super-critical linear
birth and death process. At some time an infectious disease (of SIR or SEIR
type) is introduced by one individual being infected from outside. It is shown
that three different scenarios may occur: 1) an epidemic never takes off, 2) an
epidemic gets going and grows but at a slower rate than the community thus
still being negligible in terms of population fractions, or 3) an epidemic
takes off and grows quicker than the community eventually leading to an endemic
equilibrium. Depending on the parameter values, either scenario 1 is the only
possibility, both scenario 1 and 2 are possible, or scenario 1 and 3 are
possible.


  We consider the problem of identifying an infection source based only on an
observed set of infected nodes in a network, assuming that the infection
process follows a Susceptible-Infected-Susceptible (SIS) model. We derive an
estimator based on estimating the most likely infection source associated with
the most likely infection path. Simulation results on regular trees suggest
that our estimator performs consistently better than the minimum distance
centrality based heuristic.


  The persistence conjecture is a long-standing open problem in chemical
reaction network theory. It concerns the behavior of solutions to coupled ODE
systems that arise from applying mass-action kinetics to a network of chemical
reactions. The idea is that if all reactions are reversible in a weak sense,
then no species can go extinct. A notion that has been found useful in thinking
about persistence is that of "critical siphon." We explore the combinatorics of
critical siphons, with a view towards the persistence conjecture. We introduce
the notions of "drainable" and "self-replicable" (or autocatalytic) siphons. We
show that: every minimal critical siphon is either drainable or
self-replicable; reaction networks without drainable siphons are persistent;
and non-autocatalytic weakly-reversible networks are persistent. Our results
clarify that the difficulties in proving the persistence conjecture are
essentially due to competition between drainable and self-replicable siphons.


  In a seminal paper, Valiant (2006) introduced a computational model for
evolution to address the question of complexity that can arise through
Darwinian mechanisms. Valiant views evolution as a restricted form of
computational learning, where the goal is to evolve a hypothesis that is close
to the ideal function. Feldman (2008) showed that (correlational) statistical
query learning algorithms could be framed as evolutionary mechanisms in
Valiant's model. P. Valiant (2012) considered evolvability of real-valued
functions and also showed that weak-optimization algorithms that use
weak-evaluation oracles could be converted to evolutionary mechanisms.
  In this work, we focus on the complexity of representations of evolutionary
mechanisms. In general, the reductions of Feldman and P. Valiant may result in
intermediate representations that are arbitrarily complex (polynomial-sized
circuits). We argue that biological constraints often dictate that the
representations have low complexity, such as constant depth and fan-in
circuits. We give mechanisms for evolving sparse linear functions under a large
class of smooth distributions. These evolutionary algorithms are
attribute-efficient in the sense that the size of the representations and the
number of generations required depend only on the sparsity of the target
function and the accuracy parameter, but have no dependence on the total number
of attributes.


  Several variations of the Watterson estimator of variability for Next
Generation Sequencing (NGS) data have been proposed in the literature. We
present a unified framework for generalized Watterson estimators based on
Maximum Composite Likelihood, which encompasses most of the existing
estimators. We propose this class of unbiased estimators as generalized
Watterson estimators for a large class of NGS data, including pools and trios.
We also discuss the relation with the estimators that have been proposed in the
literature and show that they admit two equivalent but seemingly different
forms, deriving a set of combinatorial identities as a byproduct. Finally, we
give a detailed treatment of Watterson estimators for single or multiple
autopolyploid individuals.


  The complexity of the primary structure of human DNA is explored using
methods from nonequilibrium statistical mechanics, dynamical systems theory and
information theory. The use of chi-square tests shows that DNA cannot be
described as a low order Markov chain of order up to $r=6$. Although detailed
balance seems to hold at the level of purine-pyrimidine notation it fails when
all four basepairs are considered, suggesting spatial asymmetry and
irreversibility. Furthermore, the block entropy does not increase linearly with
the block size, reflecting the long range nature of the correlations in the
human genomic sequences. To probe locally the spatial structure of the chain we
study the exit distances from a specific symbol, the distribution of recurrence
distances and the Hurst exponent, all of which show power law tails and long
range characteristics. These results suggest that human DNA can be viewed as a
non-equilibrium structure maintained in its state through interactions with a
constantly changing environment. Based solely on the exit distance distribution
accounting for the nonequilibrium statistics and using the Monte Carlo
rejection sampling method we construct a model DNA sequence. This method allows
to keep all long range and short range statistical characteristics of the
original sequence. The model sequence presents the same characteristic
exponents as the natural DNA but fails to capture point-to-point details.


  Aquifex aeolicus is a deep-branching hyperthermophilic chemoautotrophic
bacterium restricted to hydrothermal vents and hot springs. These
characteristics make it an excellent model system for studying the early
evolution of metabolism. Here we present the whole-genome metabolic network of
this organism and examine in detail the driving forces that have shaped it. We
make extensive use of phylometabolic analysis, a method we recently introduced
that generates trees of metabolic phenotypes by integrating phylogenetic and
metabolic constraints. We reconstruct the evolution of a range of metabolic
sub-systems, including the reductive citric acid (rTCA) cycle, as well as the
biosynthesis and functional roles of several amino acids and cofactors. We show
that A. aeolicus uses the reconstructed ancestral pathways within many of these
sub-systems, and highlight how the evolutionary interconnections between
sub-systems facilitated several key innovations. Our analyses further highlight
three general classes of driving forces in metabolic evolution. One is the
duplication and divergence of genes for enzymes as these progress from lower to
higher substrate specificity, improving the kinetics of certain sub-systems. A
second is the kinetic optimization of established pathways through fusion of
enzymes, or their organization into larger complexes. The third is the
minimization of the ATP unit cost to synthesize biomass, improving
thermodynamic efficiency. Quantifying the distribution of these classes of
innovations across metabolic sub-systems and across the tree of life will allow
us to assess how a tradeoff between maximizing growth rate and growth
efficiency has shaped the long-term metabolic evolution of the biosphere.


  The genetic code has a high level of error robustness. Using values of
hydrophobicity scales as a proxy for amino acid character, and the Mean Square
measure as a function quantifying error robustness, a value can be obtained for
a genetic code which reflects the error robustness of that code. By comparing
this value with a distribution of values belonging to codes generated by random
permutations of amino acid assignments, the level of error robustness of a
genetic code can be quantified. We present a calculation in which the standard
genetic code is shown to be optimal. We obtain this result by (1) using
recently updated values of polar requirement as input; (2) fixing seven
assignments (Ile, Trp, His, Phe, Tyr, Arg, and Leu) based on aptamer
considerations; and (3) using known biosynthetic relations of the 20 amino
acids. This last point is reflected in an approach of subdivision (restricting
the random reallocation of assignments to amino acid subgroups, the set of 20
being divided in four such subgroups). The three approaches to explain
robustness of the code (specific selection for robustness, amino acid-RNA
interactions leading to assignments, or a slow growth process of assignment
patterns) are reexamined in light of our findings. We offer a comprehensive
hypothesis, stressing the importance of biosynthetic relations, with the code
evolving from an early stage with just glycine and alanine, via intermediate
stages, towards 64 codons carrying todays meaning.


  Paternal age at conception has been found to predict the number of new
genetic mutations. We examined the effect of father's age at birth on offspring
intelligence, head circumference and personality traits. Using the Minnesota
Twin Family Study sample we tested paternal age effects while controlling for
parents' trait levels measured with the same precision as offspring's. From
evolutionary genetic considerations we predicted a negative effect of paternal
age on offspring intelligence, but not on other traits. Controlling for
parental IQ had the effect of turning a positive-zero order association
negative. We found paternal age effects on offspring IQ and MPQ Absorption, but
they were not robustly significant, nor replicable with additional covariates.
No other noteworthy effects were found. Parents' intelligence and personality
correlated with their ages at twin birth, which may have obscured a small
negative effect of advanced paternal age (< 1% of variance explained) on
intelligence. We discuss future avenues for studies of paternal age effects and
suggest that stronger research designs are needed to rule out confounding
factors involving birth order and the Flynn effect.


  Hen egg white lysozyme (HEWL) co-crystallisation conditions of carboplatin
without sodium chloride (NaCl) have been utilised to eliminate partial
conversion of carboplatin to cisplatin observed previously. Tetragonal HEWL
crystals were successfully obtained in 65% MPD with 0.1M citric acid buffer at
pH 4.0 including DMSO. The X-ray diffraction data resolution to be used for the
model refinement was reviewed using several topical criteria together. The
CC1/2 criterion implemented in XDS led to data being significant to 2.0{\AA},
compared to the data only being able to be processed to 3.0{\AA} using the
Bruker software package (SAINT). Then using paired protein model refinements
and DPI values based on the FreeR value, the resolution limit was fine tuned to
be 2.3{\AA}. Interestingly this was compared with results from the EVAL
software package which gave a resolution limit of 2.2{\AA} solely using
<I/sigI> crossing 2, but 2.8{\AA} based on the Rmerge values (60%). The
structural results showed that carboplatin bound to only the N{\delta} binding
site of His-15 one week after crystal growth, whereas five weeks after crystal
growth, two molecules of carboplatin are bound to the His-15 residue. In
summary several new results have emerged: - firstly non-NaCl conditions showed
a carboplatin molecule bound to His-15 of HEWL; secondly binding of one
molecule of carboplatin was seen after one week of crystal growth and two
molecules were bound after five weeks of crystal growth; and thirdly the use of
several criteria to determine the diffraction resolution limit led to the
successful use of data to higher resolution.


  Building a structure using self-assembly of DNA molecules by origami folding
requires finding a route for the scaffolding strand through the desired
structure. When the target structure is a 1-complex (or the geometric
realization of a graph), an optimal route corresponds to an Eulerian circuit
through the graph with minimum turning cost. By showing that it leads to a
solution to the 3-SAT problem, we prove that the general problem of finding an
optimal route for a scaffolding strand for such structures is NP-hard. We then
show that the problem may readily be transformed into a Traveling Salesman
Problem (TSP), so that machinery that has been developed for the TSP may be
applied to find optimal routes for the scaffolding strand in a DNA origami
self-assembly process. We give results for a few special cases, showing for
example that the problem remains intractable for graphs with maximum degree 8,
but is polynomial time for 4-regular plane graphs if the circuit is restricted
to following faces. We conclude with some implications of these results for
related problems, such as biomolecular computing and mill routing problems.


  Picasso's Guernica sketches continue to provide a fruitful testing ground for
examining and assessing the Blind Variation Selective Retention (BVSR) theory
of creativity. Nonmonotonicity--e.g. as indicated by a lack of similarity of
successive sketches--is not evidence of a selectionist process; Darwin's theory
explains adaptive change, not nonmonotonicity. Although the notion of blindness
originally implied randomness, it now encompasses phenomena that bias idea
generation, e.g. the influence of remote associations on sketch ideas. However,
for a selectionist framework is to be applicable, such biases must be
negligible, otherwise evolutionary change is attributed to those biases, not to
selection. The notion of 'variants' should not be applied to creativity;
without a mechanism of inheritance, there is no basis upon which to delineate,
for example, which sketch ideas are or are not variants of a given sketch idea.
The notion of selective retention is also problematic. Selection provides an
explanation when acquired change is not transmitted; it cannot apply to
Picasso's painting (or other creative acts) because his ideas acquired
modifications as he thought them through that were incorporated into paintings
and viewed by others. The generation of one sketch affects the criteria by
which the next is judged, so sequentially generated sketches cannot be treated
as members of a generation, and selected amongst. Although BVSR is
inappropriate as a theoretical framework for creativity, exploring to what
extent selectionism explains the generation of not just biological form but
masterpieces such as Picasso's Guernica is useful for gaining insight into
creativity.


  It has been proposed that, since the origin of life and the ensuing evolution
of biological species, a second evolutionary process has appeared on our
planet. It is the evolution of culture-e.g., ideas, beliefs, and artifacts.
Does culture evolve in the same genuine sense as biological life? And if so,
does it evolve through natural selection, or by some other means? Why does no
other species remotely approach the degree of cultural complexity of humans?
These questions lie at the foundation of who we are and what makes our lives
meaningful. Although much research has been done on how selective pressures
operating at the biological level affect cognition and culture, little research
has focused on culture as an evolutionary process in its own right. Like
biological forms, cultural forms-ideas, attitudes, artifacts, mannerisms,
etc.-incrementally adapt to the constraints and affordances of their
environment through descent with modification. In some respects culture appears
to be Darwinian, i.e., a process of differential replication and selection
amongst randomly generated variants. This suggests that knowledge of biological
evolution can be put to use to gain insight into culture. However, attempts to
apply Darwinian theory to culture have not yielded the kind of unifying
framework for the social sciences that it provided for the biological sciences,
largely because of the nonrandom manner in which the mind-the hub of cultural
change-generates and assimilates novelty. This paper investigates how and when
humans became capable of supporting culture, and what previously held it back,
focusing on how we attained the creative powers we now possess. To invent in
the strategic, intuitive manner characteristic of humans requires a cognitive
architecture that supports the capacity to spontaneously adapt concepts to new
circumstances and merge them together to conceptualize new situations.


  In various applications involving complex networks, network measures are
employed to assess the relative importance of network nodes. However, the
robustness of such measures in the presence of link inaccuracies has not been
well characterized. Here we present two simple stochastic models of false and
missing links and study the effect of link errors on three commonly used node
centrality measures: degree centrality, betweenness centrality, and dynamical
importance. We perform numerical simulations to assess robustness of these
three centrality measures. We also develop an analytical theory, which we
compare with our simulations, obtaining very good agreement.


  The "developmental hourglass" describes a pattern of increasing morphological
divergence towards earlier and later embryonic development, separated by a
period of significant conservation across distant species (the "phylotypic
stage"). Recent studies have found evidence in support of the hourglass effect
at the genomic level. For instance, the phylotypic stage expresses the oldest
and most conserved transcriptomes. However, the regulatory mechanism that
causes the hourglass pattern remains an open question. Here, we use an
evolutionary model of regulatory gene interactions during development to
identify the conditions under which the hourglass effect can emerge in a
general setting. The model focuses on the hierarchical gene regulatory network
that controls the developmental process, and on the evolution of a population
under random perturbations in the structure of that network. The model
predicts, under fairly general assumptions, the emergence of an hourglass
pattern in the structure of a temporal representation of the underlying gene
regulatory network. The evolutionary age of the corresponding genes also
follows an hourglass pattern, with the oldest genes concentrated at the
hourglass waist. The key behind the hourglass effect is that developmental
regulators should have an increasingly specific function as development
progresses. Analysis of developmental gene expression profiles from Drosophila
melanogaster and Arabidopsis thaliana provide consistent results with our
theoretical predictions.


  We perform an exhaustive analysis of genome statistics for organisms,
particularly extremophiles, growing in a wide range of physicochemical
conditions. Specifically, we demonstrate how the correlation between the
frequency of amino acids and their molecular weight, preserved on average,
typically decreases as optimal growth temperature increases. We show how the
relation between codon degeneracy and amino acid mass is enforced across these
organisms. We assess the occurrence of contiguous amino acids, finding several
significant short words, often containing cysteine, histidine or proline.
Typically, the significance of these words is independent of growth
temperature. In a novel approach, first-passage distributions are used to
capture correlations between discontiguous residues. We find a nearly universal
exponential background that we relate to properties of the aforementioned
individual amino acid frequencies. We find this approach reliably extracts
correlations that depend on growth temperature, some of which have not been
previously characterized.


  The regulation of metabolic activity by tuning enzyme expression levels is
crucial to sustain cellular growth in changing environments. Metabolic networks
are often studied at steady state using constraint-based models and
optimization techniques. However, metabolic adaptations driven by changes in
gene expression cannot be analyzed by steady state models, as these do not
account for temporal changes in biomass composition. Here we present a dynamic
optimization framework that integrates the metabolic network with the dynamics
of biomass production and composition, explicitly taking into account enzyme
production costs and enzymatic capacity. In contrast to the established dynamic
flux balance analysis, our approach allows predicting dynamic changes in both
the metabolic fluxes and the biomass composition during metabolic adaptations.
We applied our algorithm in two case studies: a minimal nutrient uptake
network, and an abstraction of core metabolic processes in bacteria. In the
minimal model, we show that the optimized uptake rates reproduce the empirical
Monod growth for bacterial cultures. For the network of core metabolic
processes, the dynamic optimization algorithm predicted commonly observed
metabolic adaptations, such as a diauxic switch with a preference ranking for
different nutrients, re-utilization of waste products after depletion of the
original substrate, and metabolic adaptation to an impending nutrient
depletion. These examples illustrate how dynamic adaptations of enzyme
expression can be predicted solely from an optimization principle.


  The Weismann barrier, or the impossibility of inheritance of acquired traits,
comprises a foundation of modern biology, and it has been a major obstacle in
establishing the connection between evolution and ontogenesis. We propose the
cooperative model based on the assumption that evolution is achieved by a
cooperation between genetic mutations and acquired changes (phenotypic
plasticity). It is also assumed in this model that natural selection operates
on phenotypes, rather than genotypes, of individuals, and that the relationship
between phenotypes and genotypes is one-to-many. In the simulations based on
these assumptions, individuals exhibited phenotypic changes in response to an
environmental change, corresponding multiple genetic mutations were
increasingly accumulated in individuals in the population, and phenotypic
plasticity was gradually replaced with genetic mutations. This result suggests
that Lamarck's law of use and disuse can effectively hold without conflicting
the Weismann barrier, and thus evolution can be logically connected with
ontogenesis.


  Bacterial spores in a metabolically dormant state can survive long periods
without nutrients under extreme environmental conditions. The molecular basis
of spore dormancy is not well understood, but the distribution and physical
state of water within the spore is thought to play an important role. Two
scenarios have been proposed for the spore's core region, containing the DNA
and most enzymes. In the gel scenario, the core is a structured macromolecular
framework permeated by mobile water. In the glass scenario, the entire core,
including the water, is an amorphous solid and the quenched molecular diffusion
accounts for the spore's dormancy and thermal stability. Here, we use $^2$H
magnetic relaxation dispersion to selectively monitor water mobility in the
core of Bacillus subtilis spores in the presence and absence of core Mn$^{2+}$
ions. We also report and analyze the solid-state $^2$H NMR spectrum from these
spores. Our NMR data clearly support the gel scenario with highly mobile core
water (~ 25 ps average rotational correlation time). Furthermore, we find that
the large depot of manganese in the core is nearly anhydrous, with merely 1.7 %
on average of the maximum sixfold water coordination.


  The combination of multiple classifiers using ensemble methods is
increasingly important for making progress in a variety of difficult prediction
problems. We present a comparative analysis of several ensemble methods through
two case studies in genomics, namely the prediction of genetic interactions and
protein functions, to demonstrate their efficacy on real-world datasets and
draw useful conclusions about their behavior. These methods include simple
aggregation, meta-learning, cluster-based meta-learning, and ensemble selection
using heterogeneous classifiers trained on resampled data to improve the
diversity of their predictions. We present a detailed analysis of these methods
across 4 genomics datasets and find the best of these methods offer
statistically significant improvements over the state of the art in their
respective domains. In addition, we establish a novel connection between
ensemble selection and meta-learning, demonstrating how both of these disparate
methods establish a balance between ensemble diversity and performance.


  The sample frequency spectrum (SFS) is a widely-used summary statistic of
genomic variation in a sample of homologous DNA sequences. It provides a highly
efficient dimensional reduction of large-scale population genomic data and its
mathematical dependence on the underlying population demography is well
understood, thus enabling the development of efficient inference algorithms.
However, it has been recently shown that very different population demographies
can actually generate the same SFS for arbitrarily large sample sizes. Although
in principle this nonidentifiability issue poses a thorny challenge to
statistical inference, the population size functions involved in the
counterexamples are arguably not so biologically realistic. Here, we revisit
this problem and examine the identifiability of demographic models under the
restriction that the population sizes are piecewise-defined where each piece
belongs to some family of biologically-motivated functions. Under this
assumption, we prove that the expected SFS of a sample uniquely determines the
underlying demographic model, provided that the sample is sufficiently large.
We obtain a general bound on the sample size sufficient for identifiability;
the bound depends on the number of pieces in the demographic model and also on
the type of population size function in each piece. In the cases of
piecewise-constant, piecewise-exponential and piecewise-generalized-exponential
models, which are often assumed in population genomic inferences, we provide
explicit formulas for the bounds as simple functions of the number of pieces.
Lastly, we obtain analogous results for the "folded" SFS, which is often used
when there is ambiguity as to which allelic type is ancestral. Our results are
proved using a generalization of Descartes' rule of signs for polynomials to
the Laplace transform of piecewise continuous functions.


  The domestication and development of cattle has considerably impacted human
societies, but the histories of cattle breeds have been poorly understood
especially for African, Asian, and American breeds. Using genotypes from 43,043
autosomal single nucleotide polymorphism markers scored in 1,543 animals, we
evaluate the population structure of 134 domesticated bovid breeds. Regardless
of the analytical method or sample subset, the three major groups of Asian
indicine, Eurasian taurine, and African taurine were consistently observed.
Patterns of geographic dispersal resulting from co-migration with humans and
exportation are recognizable in phylogenetic networks. All analytical methods
reveal patterns of hybridization which occurred after divergence. Using 19
breeds, we map the cline of indicine introgression into Africa. We infer that
African taurine possess a large portion of wild African auroch ancestry,
causing their divergence from Eurasian taurine. We detect exportation patterns
in Asia and identify a cline of Eurasian taurine/indicine hybridization in
Asia. We also identify the influence of species other than Bos taurus in the
formation of Asian breeds. We detect the pronounced influence of Shorthorn
cattle in the formation of European breeds. Iberian and Italian cattle possess
introgression from African taurine. American Criollo cattle are shown to be of
Iberian, and not African, decent. Indicine introgression into American cattle
occurred in the Americas, and not Europe. We argue that cattle migration,
movement and trading followed by admixture have been important forces in
shaping modern bovine genomic variation.


  In many sensory systems, transmembrane receptors are spatially organized in
large clusters. Such arrangement may facilitate signal amplification and the
integration of multiple stimuli. However, this organization likely also affects
the kinetics of signaling since the cytoplasmic enzymes that modulate the
activity of the receptors must localize to the cluster prior to receptor
modification. Here we examine how these spatial considerations shape signaling
dynamics at rest and in response to stimuli. As a model, we use the chemotaxis
pathway of Escherichia coli, a canonical system for the study of how organisms
sense, respond, and adapt to environmental stimuli. In bacterial chemotaxis,
adaptation is mediated by two enzymes that localize to the clustered receptors
and modulate their activity through methylation-demethylation. Using a novel
stochastic simulation, we show that distributive receptor methylation is
necessary for successful adaptation to stimulus and also leads to large
fluctuations in receptor activity in the steady state. These fluctuations arise
from noise in the number of localized enzymes combined with saturated
modification kinetics between localized enzymes and receptor substrate. An
analytical model explains how saturated enzyme kinetics and large fluctuations
can coexist with an adapted state robust to variation in the expression level
of the pathway constituents, a key requirement to ensure the functionality of
individual cells within a population. This contrasts with the well-mixed
covalent modification system studied by Goldbeter and Koshland in which mean
activity becomes ultrasensitive to protein abundances when the enzymes operate
at saturation. Large fluctuations in receptor activity have been quantified
experimentally. Here we clarify their mechanistic relationship with
well-studied aspects of the chemotaxis system, precise adaptation and
functional robustness.


  How ATP binding initiates the docking process of kinesin's neck linker is a
key question in understanding kinesin mechanism. It is believed that the
formation of an extra turn structure by the first three amino acids of neck
linker (LYS325, THR326, ILE327 in 2KIN) is crucial for initiating the docking
process. But the initial conformation of neck linker (specially the three amino
acids of the extra turn) and the neck linker docking initiation mechanism
remain unclear. By using molecular dynamics method, we investigate the initial
conformation of kinesin's neck linker in the docking process. We find that, in
the initial state of NL docking process, NL still has interactions with
{\beta}0 and forms a conformation similar to the "cover-neck bundle" structure
proposed by Hwang et al. [Structure 2008, 16(1): 62-71]. From this initial
structure, the docking of the "cover-neck bundle" structure can be achieved.
The motor head provides a forward force on the initial cover-neck bundle
structure through ATP-induced rotation. This force, together with the
hydrophobic interaction of ILE327 with the hydrophobic pocket on the motor
head, drives the formation of the extra turn and initiates the neck linker
docking process.


  Nanodiamonds (NDs) containing silicon vacancy (SiV) defects were evaluated as
a potential biomarker for the labeling and fluorescent imaging of neural
precursor cells (NPCs). SiV-containing NDs were synthesized using chemical
vapor deposition and silicon ion implantation. Spectrally, SiV-containing NDs
exhibited extremely stable fluorescence and narrow bandwidth emission with an
excellent signal to noise ratio exceeding that of NDs containing
nitrogen-vacancy (NV) centers. NPCs labeled with NDs exhibited normal cell
viability and proliferative properties consistent with biocompatibility. We
conclude that SiVcontaining NDs are a promising biomedical research tool for
cellular labeling and optical imaging in stem cell research.


  We unravel how functional plasticity and redundancy are essential mechanisms
underlying the ability to survive of metabolic networks. We perform an
exhaustive computational screening of synthetic lethal reaction pairs in
Escherichia coli in a minimal medium and we find that synthetic lethal pairs
divide in two different groups depending on whether the synthetic lethal
interaction works as a backup or as a parallel use mechanism, the first
corresponding to essential plasticity and the second to essential redundancy.
In E. coli, the analysis of pathways entanglement through essential redundancy
supports the view that synthetic lethality affects preferentially a single
function or pathway. In contrast, essential plasticity, the dominant class,
tends to be inter-pathway but strongly localized and unveils Cell Envelope
Biosynthesis as an essential backup for Membrane Lipid Metabolism. When
comparing E. coli and Mycoplasma pneumoniae, we find that the metabolic
networks of the two organisms exhibit a large difference in the relative
importance of plasticity and redundancy which is consistent with the conjecture
that plasticity is a sophisticated mechanism that requires a complex
organization. Finally, coessential reaction pairs are explored in different
environmental conditions to uncover the interplay between the two mechanisms.
We find that synthetic lethal interactions and their classification in
plasticity and redundancy are basically insensitive to medium composition, and
are highly conserved even when the environment is enriched with nonessential
compounds or overconstrained to decrease maximum biomass formation.


  Recent advances in high-throughput cDNA sequencing (RNA-Seq) technology have
revolutionized transcriptome studies. A major motivation for RNA-Seq is to map
the structure of expressed transcripts at nucleotide resolution. With accurate
computational tools for transcript reconstruction, this technology may also
become useful for genome (re-)annotation, which has mostly relied on de novo
gene finding where gene structures are primarily inferred from the genome
sequence. We developed a machine-learning method, called mTim (margin-based
transcript inference method) for transcript reconstruction from RNA-Seq read
alignments that is based on discriminatively trained hidden Markov support
vector machines. In addition to features derived from read alignments, it
utilizes characteristic genomic sequences, e.g. around splice sites, to improve
transcript predictions. mTim inferred transcripts that were highly accurate and
relatively robust to alignment errors in comparison to those from Cufflinks, a
widely used transcript assembly method.


  We study evolutionary games in a spatial diluted grid environment in which
agents strategically interact locally but can also opportunistically move to
other positions within a given migration radius. Using the imitation of the
best rule for strategy revision, it is shown that cooperation may evolve and be
stable in the Prisoner's Dilemma game space for several migration distances but
only for small game interaction radius while the Stag Hunt class of games
become fully cooperative. We also show that only a few trials are needed for
cooperation to evolve, i.e. searching costs are not an issue. When the
stochastic Fermi strategy update protocol is used cooperation cannot evolve in
the Prisoner's Dilemma if the selection intensity is high in spite of
opportunistic migration. However, when imitation becomes more random, fully or
partially cooperative states are reached in all games for all migration
distances tested and for short to intermediate interaction radii.


  We adapted the highly successful scattering density profile (SDP) model,
which is used to jointly analyze small angle x-ray and neutron scattering data
from unilamellar vesicles, for use with data from fully hydrated, liquid
crystalline multilamellar vesicles (MLVs). By using a genetic algorithm, this
new method is capable of providing high resolution structural information, as
well as determining bilayer elastic bending fluctuations from standalone x-ray
data. Important structural parameters, such as bilayer thickness and area per
lipid were determined for a series of saturated and unsaturated lipids, as well
as binary mixtures with cholesterol. Results are in good agreement with
previously reported SDP data, which used both neutron and x-ray data. The
addition of deuterated and non-deuterated MLV neutron data to the analysis
improved lipid backbone information, but did not improve, within experimental
uncertainties, the structural information regarding bilayer thickness and area
per lipid.


  Metadynamics is a commonly used and successful enhanced sampling method. By
the introduction of a history dependent bias which depends on a restricted
number of collective variables(CVs) it can explore complex free energy surfaces
characterized by several metastable states separated by large free energy
barriers. Here we extend its scope by introducing a simple yet powerful method
for calculating the rates of transition between different metastable states.
The method does not rely on a previous knowledge of the transition states or
reaction co-ordinates, as long as CVs are known that can distinguish between
the various stable minima in free energy space. We demonstrate that our method
recovers the correct escape rates out of these stable states and also preserves
the correct sequence of state-to-state transitions, with minimal extra
computational effort needed over ordinary metadynamics. We apply the formalism
to three different problems and in each case find excellent agreement with the
results of long unbiased molecular dynamics runs.


  Histone modification is a vital epigenetic mechanism for transcriptional
control in eukaryotes. High-throughput techniques have enabled whole-genome
analysis of histone modifications in recent years. However, most studies assume
one combination of histone modification invariantly translates to one
transcriptional output regardless of local chromatin environment. In this study
we hypothesize that, the genome is organized into local domains that manifest
similar enrichment pattern of histone modification, which leads to orchestrated
regulation of expression of genes with relevant bio- logical functions. We
propose a multivariate Bayesian Change Point (BCP) model to segment the
Drosophila melanogaster genome into consecutive blocks on the basis of
combinatorial patterns of histone marks. By modeling the sparse distribution of
histone marks across the chromosome with a zero-inflated Gaussian mixture, our
partitions capture local BLOCKs that manifest relatively homogeneous enrichment
pattern of histone modifications. We further characterized BLOCKs by their
transcription levels, distribution of genes, degree of co-regulation and GO
enrichment. Our results demonstrate that these BLOCKs, although inferred merely
from histone modifications, reveal strong relevance with physical domains,
which suggest their important roles in chromatin organization and coordinated
gene regulation.


  We introduce a multi-agent model for exploring how selection of neighbours
determines some aspects of order and cohesion in swarms. The model algorithm
states that every agents' motion seeks for an optimal distance from the nearest
topological neighbour encompassed in a limited attention field. Despite the
great simplicity of the implementation, varying the amplitude of the attention
landscape, swarms pass from cohesive and regular structures towards fragmented
and irregular configurations. Interestingly, this movement rule is an ideal
candidate for implementing the selfish herd hypothesis which explains
aggregation of alarmed group of social animals.


  Recently, Bogu\~{n}\'{a} {\it et. al.} [Phys. Rev. Lett. {\bf 111}, 068701
(2013), arXiv:1305.4819] claimed that the epidemic threshold of the
susceptible-infected-susceptible (SIS) model is zero on random scale-free (SF)
networks with the small-world (SW) property. They drew the conclusion by taking
into account a long-range reinfection mechanism. In this Comment, we will show
that such an effect is too weak to guarantee an endemic phase with a {\em
finite} fraction of infected nodes. Consequently, the epidemic threshold can be
finite in random SF networks.


  The original ideas of Cooper and Dryden, that allosteric signalling can be
induced between distant binding sites on proteins without any change in mean
structural conformation, has proved to be a remarkably prescient insight into
the rich structure of protein dynamics. It represents an alternative to the
celebrated Monod-Wyman-Changeux mechanism and proposes that modulation of the
amplitude of thermal fluctuations around a mean structure, rather than shifts
in the structure itself, give rise to allostery in ligand binding. In a
complementary approach to experiments on real proteins, here we take a
theoretical route to identify the necessary structural components of this
mechanism. By reviewing and extending an approach that moves from very
coarse-grained to more detailed models, we show that, a fundamental requirement
for a body supporting fluctuation-induced allostery is a strongly inhomogeneous
elastic modulus. This requirement is reflected in many real proteins, where a
good approximation of the elastic structure maps strongly coherent domains onto
rigid blocks connected by more flexible interface regions.


  We modify a mathematical model of photosynthesis to quantify the
perturbations that high energy muons could make on aquatic primary
productivity. Then we apply this in the context of the extragalactic shock
model, according to which Earth receives an enhanced dose of high-energy cosmic
rays when it is at the galactic north. We obtain considerable reduction in the
photosynthesis rates, consistent with potential drops in biodiversity.


  A variant of microcephalin, MCPH1 gene, was introgressed about 37,000 years
ago into Homo sapiens genetic pool from an archaic (Homo erectus) lineage and
rose to exceptionally high frequency of around 70 percent worldwide today. It
is involved in regulating neuroblast proliferation and its changes alter the
rate of division and/or differentiation of neuroblasts during the neurogenic
phase of embriogenesis, which could alter the size and structure of the
resulting brain.
  At the time of introgression, images had already been painted on the walls of
caves and speech has been in use for over 100,000 years, as had been abstract
thinking. Like today, reasoning and thinking were the primary faculties of
individuals. Homo erectus either did not possess those faculties or was
markedly inferior to Homo sapiens in them. Its brain was smaller and the cortex
was apparently less convoluted. Thus, introgressed microcephalin allele
directed neurogenesis evolutionary back to less complicated brain structure
typical for our evolutionary forefathers, slightly decreasing the level of
complexity already achieved by Homo sapiens 37,000 years ago. Despite that, it
proliferated at a rapid pace.
  It yields a supposition: 37,000 years ago the brains of Homo sapiens were too
big and too complicated for the kind of thinking needed for the highest fitness
of individuals. Since adaptation cannot by definition surpass selection
requirements, the volume and complication of the human brain did not originate
under selective pressure to improve effective thinking and they cannot be
explained in terms of such selection.
  A proposal to solve this quandary is presented, claiming that Homo sapiens
originated just by chance. Endurance running led to the emergence of Homo
sapiens. The human mind and larynx used for speech are side-effects of more
than a million years of endurance running by pre-human hunters.


  It is important to be clear as to whether a theory such as evolutionary
archaeology pertains to biological evolution, in which acquired change is
obliterated at the end of each generation, or cultural change, in which
acquired change is retained. In evolutionary archaeology, (1) the population is
said to consist of artifacts, yet (2) artifacts are said to be phenotypic.
Neither (1) nor (2) is necessarily problematic in and of itself, but the two
are inconsistent, as the first pertains to cultural change whereas the second
to the biological evolution of humans. A first step to avoiding this problem is
to recognize that there is a need for a theory of change specific to human
culture. Referring to ongoing work using a related approach to cultural change,
it is suggested that the inconsistencies in evolutionary archaeology, though
problematic, are not insurmountable.


  Computer simulations carried out by Alroy (2001) are examined. Contrary to
the claim of their Author, there is no convincing evidence that the extinction
of megafauna was caused by humans. Intentionally or unintentionally, attempt
was made to force the human-induced extinction of megafauna by assuming an
absurdly fast growth of the hypothetical human population in North America. The
assumed growth rate was around two-orders of magnitude larger than normally
expected. It is well known that the past growth of human population was slow,
but in these simulations, the growth of the hypothetical human population was
unreasonably fast. However, even under this unreasonable assumption about the
growth of human population, computer simulations do not support the postulate
of human-assisted extinction of megafauna because there is no clear correlation
between the growth of the hypothetical human population and the distributions
describing the decline in the population of megafauna.


  Riboswitches are structured mRNA elements that modulate gene expression. They
undergo conformational changes triggered by highly specific interactions with
sensed metabolites. Among the structural rearrangements engaged by
riboswitches, the forming and melting of the aptamer terminal helix, the
so-called P1 stem, is essential for genetic control. The structural mechanisms
by which this conformational change is modulated upon ligand binding mostly
remain to be elucidated. Here, we used pulling molecular dynamics simulations
to study the thermodynamics of the P1 stem in the add adenine riboswitch. The
P1 ligand-dependent stabilization was quantified in terms of free energy and
compared with thermodynamic data. This comparison suggests a model for the
aptamer folding in which direct P1-ligand interactions play a minor role on the
conformational switch when compared with those related to the ligand-induced
aptamer preorganization.


  Empirical growth rates of the ancient human population in Australia have been
calculated using the data for the number of rock-shelter sites. They confirm
the earlier conclusion that there was no intensification of growth between
10,000 and 1000 years BP.


  Collective dynamics and force generation by cytoskeletal filaments are
crucial in many cellular processes. Investigating growth dynamics of a bundle
of N independent cytoskeletal filaments pushing against a wall, we show that
chemical switching (ATP/GTP hydrolysis) leads to a collective phenomenon that
is currently unknown. Obtaining force-velocity relations for different models
that capture chemical switching, we show, analytically and numerically, that
the collective stall force of N filaments is greater than N times the stall
force of a single filament. Employing an exactly solvable toy model, we
analytically prove the above result for N=2. We, further, numerically show the
existence of this collective phenomenon, for N>=2, in realistic models (with
random and sequential hydrolysis) that simulate actin and microtubule bundle
growth. We make quantitative predictions for the excess forces, and argue that
this collective effect is related to the non-equilibrium nature of chemical
switching.


  A major challenge in biology is to understand how molecular processes
determine phenotypic features. We address this fundamental problem in a class
of model systems by developing a general mathematical framework that allows the
calculation of mesoscopic properties from the knowledge of microscopic
Markovian transition probabilities. We show how exact analytic formulae for the
first and second moments of resident time distributions in mesostates can be
derived from microscopic resident times and transition probabilities even for
systems with a large number of microstates. We apply our formalism to models of
the inositol trisphosphate receptor, which plays a key role in generating
calcium signals triggering a wide variety of cellular responses. We demonstrate
how experimentally accessible quantities, such as opening and closing times and
the coefficient of variation of inter-spike intervals, and other, more
elaborated, quantities can be analytically calculated from the underlying
microscopic Markovian dynamics. A virtue of our approach is that we do not need
to follow the detailed time evolution of the whole system, as we derive the
relevant properties of its steady state without having to take into account the
often extremely complicated transient features. We emphasize that our formulae
fully agree with results obtained by stochastic simulations and approaches
based on a full determination of the microscopic system's time evolution. We
also illustrate how experiments can be devised to discriminate between
alternative molecular models of the inositol trisphosphate receptor. The
developed approach is applicable to any system described by a Markov process
and, owing to the analytic nature of the resulting formulae, provides an easy
way to characterize also rare events that are of particular importance to
understand the intermittency properties of complex dynamic systems.


  Quite generally, constraint-based metabolic flux analysis describes the space
of viable flux configurations for a metabolic network as a high-dimensional
polytope defined by the linear constraints that enforce the balancing of
production and consumption fluxes for each chemical species in the system. In
some cases, the complexity of the solution space can be reduced by performing
an additional optimization, while in other cases, knowing the range of
variability of fluxes over the polytope provides a sufficient characterization
of the allowed configurations. There are cases, however, in which the thorough
information encoded in the individual distributions of viable fluxes over the
polytope is required. Obtaining such distributions is known to be a highly
challenging computational task when the dimensionality of the polytope is
sufficiently large, and the problem of developing cost-effective {\it ad hoc}
algorithms has recently seen a major surge of interest. Here, we propose a
method that allows us to perform the required computation heuristically in a
time scaling {\it linearly} with the number of reactions in the network,
overcoming some limitations of similar techniques employed in recent years. As
a case study, we apply it to the analysis of the human red blood cell metabolic
network, whose solution space can be sampled by different exact techniques,
like Hit-and-Run Monte Carlo (scaling roughly like the third power of the
system size). Remarkably accurate estimates for the true distributions of
viable reaction fluxes are obtained, suggesting that, although further
improvements are desirable, our method enhances our ability to analyze the
space of allowed configurations for large biochemical reaction networks.