  For ion channel gating, the appearance of two distinct conformational states
and the discrete transitions between them is essential, and therefore of
crucial importance to all living organisms. We show that the physical interplay
between two structural elements that are commonly present in bacterial
mechanosensitive channels, namely a charged vestibule and a hydrophobic
constriction, creates two distinct conformational states, open and closed, as
well as the gating between them. We solve the nonequilibrium
Stokes-Poisson-Nernst-Planck equations, extended to include a molecular
potential of mean force, and show that a first order transition between the
closed and open states arises naturally from the diffusio-osmotic stress caused
by the ions and water inside the channel and the elastic restoring force from
the membrane. Our proposed gating mechanism is likely to be important for a
broad range of ion channels, as well as for biomimetic channels and ion
channel-targeting therapeutics.


  Here we provide a thorough discussion of the model for Min protein dynamics
proposed by Schweizer et al. [11]. The manuscript serves as supplementary
document for our letter to the editor to appear in PNAS. Our analysis is based
on the original COMSOL simulation files that were used for the publication. We
show that all computational data in Schweizer et al. rely on exploitation of
simulation artifacts and various unmentioned modifications of model parameters
that strikingly contradict the experimental setup and experimental data. We
find that the model neither accounts for MinE membrane interactions nor for any
observed MinDE protein patterns. All conclusions drawn from the computational
model are void. There is no evidence at all that persistent MinE membrane
binding has any role in geometry sensing.


  By considering subexponential contributions in large deviation theory, we
determine the fine structure in the probability distribution of the observable
displacement of a bead coupled to a molecular motor. More generally, for any
stochastic motion along a periodic substrate, this approach reveals a discrete
symmetry of this distribution for which hidden degrees of freedom lead to a
periodic modulation of the slope typically associated with the fluctuation
theorem. Contrary to previous interpretations of experimental data, the mean
force exerted by a molecular motor is unrelated to the long-time asymptotics of
this slope and must rather be extracted from its short-time limit.


  Strongly non-Markovian random walks offer a promising modeling framework for
understanding animal and human mobility, yet, few analytical results are
available for these processes. Here we solve exactly a model with long range
memory where a random walker intermittently revisits previously visited sites
according to a reinforced rule. The emergence of frequently visited locations
generates very slow diffusion, logarithmic in time, whereas the walker
probability density tends to a Gaussian. This scaling form does not emerge from
the Central Limit Theorem but from an unusual balance between random and
long-range memory steps. In single trajectories, occupation patterns are
heterogeneous and have a scale-free structure. The model exhibits good
agreement with data of free-ranging capuchin monkeys.


  Dominance and subordinate behaviours are important ingredients in the social
organizations of group living animals. Behavioural observations on the two
eusocial species \textit{Ropalidia marginata} and \textit{Ropalidia
cyathiformis} suggest varying complexities in their social systems. The queen
of R. cyathiformis is an aggressive individual who usually holds the top
position in the dominance hierarchy although she does not necessarily show the
maximum number of acts of dominance, while the R. marginata queen rarely shows
aggression and usually does not hold the top position in the dominance
hierarchy of her colony. These differences are reflected in the distribution of
dominance-subordinate interactions among the hierarchically ranked individuals
in both the species. The percentage of dominance interactions decrease
gradually with hierarchical ranks in R. marginata while in R. cyathiformis it
first increases and then decreases. We use an agent-based model to investigate
the underlying mechanism that could give rise to the observed patterns for both
the species. The model assumes, besides some non-interacting individuals, that
the interaction probabilities of the agents depend on their pre-differentiated
winning abilities. Our simulations show that if the queen takes up a strategy
of being involved in a moderate number of dominance interactions, one could get
the pattern similar to R. cyathiformis, while taking up the strategy of very
low interactions by the queen could lead to the pattern of R. marginata. We
infer that both the species follow a common interaction pattern, while the
differences in their social organization are due to the slight changes in queen
as well as worker strategies. These changes in strategies are expected to
accompany the evolution of more complex societies from simpler ones.


  Models of simple excitable dynamics on graphs are an efficient framework for
studying the interplay between network topology and dynamics. This subject is a
topic of practical relevance to diverse fields, ranging from neuroscience to
engineering. Here we analyze how a single excitation propagates through a
random network as a function of the excitation threshold, that is, the relative
amount of activity in the neighborhood required for an excitation of a node.
Using numerical simulations and analytical considerations, we can understand
the onset of sustained activity as an interplay between topological cycle
statistics and path statistics. Our findings are interpreted in the context of
the theory of network reverberations in neural systems, which is a question of
long-standing interest in computational neuroscience.


  Transitions in the qualitative behavior of chemical reaction dynamics with a
decrease in molecule number have attracted much attention. Here, a method based
on a Markov process with a tridiagonal transition matrix is applied to the
analysis of this transition in reaction dynamics. The transition to bistability
due to the small-number effect and the mean switching time between the bistable
states are analytically calculated in agreement with numerical simulations. In
addition, a novel transition involving the reversal of the chemical reaction
flow is found in the model under an external flow, and also in a
three-component model. The generality of this transition and its correspondence
to biological phenomena are also discussed.


  The interplay between bending of the molecule axis and appearance of
disruptions in circular DNA molecules, with $\sim 100$ base pairs, is
addressed. Three minicircles with different radii and almost equal content of
AT and GC pairs are investigated. The DNA sequences are modeled by a mesoscopic
Hamiltonian which describes the essential interactions in the helix at the
level of the base pair and incorporates twisting and bending degrees of
freedom. Helix unwinding and bubble formation patterns are consistently
computed by a path integral method that sums over a large number of molecule
configurations compatible with the model potential. The path ensembles are
determined, as a function of temperature, by minimizing the free energy of the
system. Fluctuational openings appear along the helix to release the stress due
to the bending of the molecule backbone. In agreement with the experimental
findings, base pair disruptions are found with larger probability in the
smallest minicircle of \textit{66-bps} whose bending angle is $\sim 6^{o} $.
For this minicircle, a sizeable untwisting is obtained with the helical repeat
showing a step-like increase at $\tau =\,315K$. The method can be generalized
to determine the bubble probability profiles of open ends linear sequences.


  Movement is a fundamental behaviour of organisms that brings about beneficial
encounters with resources and mates, but at the same time exposes the organism
to dangerous encounters with predators. The movement patterns adopted by
organisms should reflect a balance between these contrasting processes. This
trade-off can be hypothesized as being evident in the behaviour of plankton,
which inhabit a dilute 3D environment with few refuges or orienting landmarks.
We present an analysis of the swimming path geometries based on a volumetric
Monte Carlo sampling approach, which is particularly adept at revealing such
trade-offs by measuring the self-overlap of the trajectories. Application of
this method to experimentally measured trajectories reveals that swimming
patterns in copepods are shaped to efficiently explore volumes at small scales,
while achieving a large overlap at larger scales. Regularities in the observed
trajectories make the transition between these two regimes always sharper than
in randomized trajectories or as predicted by random walk theory. Thus real
trajectories present a stronger separation between exploration for food and
exposure to predators. The specific scale and features of this transition
depend on species, gender, and local environmental conditions, pointing at
adaptation to state and stage dependent evolutionary trade-offs.


  Risk spreading in bacterial populations is generally regarded as a strategy
to maximize survival. Here, we study its role during range expansion of a
genetically diverse population where growth and motility are two alternative
traits. We find that during the initial expansion phase fast growing cells do
have a selective advantage. By contrast, asymptotically, generalists balancing
motility and reproduction are evolutionarily most successful. These findings
are rationalized by a set of coupled Fisher equations complemented by
stochastic simulations.


  Identifying and quantifying the benefits of sex and recombination is a long
standing problem in evolutionary theory. In particular, contradictory claims
have been made about the existence of a benefit of recombination on high
dimensional fitness landscapes in the presence of sign epistasis. Here we
present a comparative numerical study of sexual and asexual evolutionary
dynamics of haploids on tunably rugged model landscapes under strong selection,
paying special attention to the temporal development of the evolutionary
advantage of recombination and the link between population diversity and the
rate of adaptation. We show that the adaptive advantage of recombination on
static rugged landscapes is strictly transitory. At early times, an advantage
of recombination arises through the possibility to combine individually
occurring beneficial mutations, but this effect is reversed at longer times by
the much more efficient trapping of recombining populations at local fitness
peaks. These findings are explained by means of well established results for a
setup with only two loci. In accordance with the Red Queen hypothesis the
transitory advantage can be prolonged indefinitely in fluctuating environments,
and it is maximal when the environment fluctuates on the same time scale on
which trapping at local optima typically occurs.


  In order to associate complex traits with genetic polymorphisms, genome-wide
association studies process huge datasets involving tens of thousands of
individuals genotyped for millions of polymorphisms. When handling these
datasets, which exceed the main memory of contemporary computers, one faces two
distinct challenges: 1) Millions of polymorphisms and thousands of phenotypes
come at the cost of hundreds of gigabytes of data, which can only be kept in
secondary storage; 2) the relatedness of the test population is represented by
a relationship matrix, which, for large populations, can only fit in the
combined main memory of a distributed architecture. In this paper, by using
distributed resources such as Cloud or clusters, we address both challenges:
The genotype and phenotype data is streamed from secondary storage using a
double buffer- ing technique, while the relationship matrix is kept across the
main memory of a distributed memory system. With the help of these solutions,
we develop separate algorithms for studies involving only one or a multitude of
traits. We show that these algorithms sustain high-performance and allow the
analysis of enormous datasets.


  The epidemic threshold of the susceptible-infected-susceptible (SIS) dynamics
on random networks having a power law degree distribution with exponent
$\gamma>3$ has been investigated using different mean-field approaches, which
predict different outcomes. We performed extensive simulations in the
quasistationary state for a comparison with these mean-field theories. We
observed concomitant multiple transitions in individual networks presenting
large gaps in the degree distribution and the obtained multiple epidemic
thresholds are well described by different mean-field theories. We observed
that the transitions involving thresholds which vanishes at the thermodynamic
limit involve localized states, in which a vanishing fraction of the network
effectively contribute to epidemic activity, whereas an endemic state, with a
finite density of infected vertices, occurs at a finite threshold. The multiple
transitions are related to the activations of distinct sub-domains of the
network, which are not directly connected.


  Mitogen-activated protein kinase (MAPK) signaling pathways play an essential
role in the transduction of environmental stimuli to the nucleus, thereby
regulating a variety of cellular processes, including cell proliferation,
differentiation and programmed cell death. The components of the MAPK
extracellular activated protein kinase (ERK) cascade represent attractive
targets for cancer therapy as their aberrant activation is a frequent event
among highly prevalent human cancers. MAPK networks are a model for
computational simulation, mostly using Ordinary and Partial Differential
Equations. Key results showed that these networks can have switch-like
behavior, bistability and oscillations. In this work, we consider three
representative ERK networks, one with a negative feedback loop, which present a
binomial steady state ideal under mass-action kinetics. We therefore apply the
theoretical result present in P\'erez Mill\'an et. al (2012) to find a set of
rate constants that allow two significantly different stable steady states in
the same stoichiometric compatibility class for each network. Our approach
makes it possible to study certain aspects of the system, such as
multistationarity, without relying on simulation, since we do not assume a
priori any constant but the topology of the network. As the performed analysis
is general it could be applied to many other important biochemical networks.


  We demonstrate that endonuclease deficient Clustered Regularly Interspaced
Short Palindromic Repeats CRISPR-associated Cas9 protein (dCas9) fused to the
photo-convertible fluorescence protein monomeric mEos3.1 (dCas9-mEos3) can be
used to resolve sub-diffraction limited features of repetitive gene elements,
thus providing a new route to investigate high-order chromatin organization at
these sites.


  At the end of the first larval stage, the C elegans larva chooses between two
developmental pathways, an L2 committed to reproductive development and an L2d,
which has the option of undergoing reproductive development or entering the
dauer diapause. I develop a quantitative model of this choice using
mathematical tools developed for pricing financial options. The model predicts
that the optimal decision must take into account not only the expected
potential for reproductive growth, but also the uncertainty in that expected
potential. Because the L2d has more flexibility than the L2, it is favored in
unpredictable environments. I estimate that the ability to take uncertainty
into account may increase reproductive value by as much as 5%, and discuss
possible experimental tests for this ability.


  Haplotype-based scans to detect natural selection are useful to identify
recent or ongoing positive selection in genomes. As both real and simulated
genomic datasets grow larger, spanning thousands of samples and millions of
markers, there is a need for a fast and efficient implementation of these scans
for general use. Here we present selscan, an efficient multi-threaded
application that implements Extended Haplotype Homozygosity (EHH), Integrated
Haplotype Score (iHS), and Cross-population Extended Haplotype Homozygosity
(XPEHH). selscan accepts phased genotypes in multiple formats, including TPED,
and performs extremely well on both simulated and real data and over an order
of magnitude faster than existing available implementations. It calculates iHS
on chromosome 22 (22,147 loci) across 204 CEU haplotypes in 353s on one thread
(33s on 16 threads) and calculates XPEHH for the same data relative to 210 YRI
haplotypes in 578s on one thread (52s on 16 threads). Source code and binaries
(Windows, OSX and Linux) are available at https://github.com/szpiech/selscan .


  1. Anthropogenic actions cause rapid ecological changes, meaning that animals
have to respond before they have time to adapt. Tools to quantify emergent
spatial patterns from animal-habitat interaction mechanisms are vital for
predicting the population-level effects of such changes.
  2. Environmental perturbations are particularly prevalent in the Amazon
rainforest, and have a profound effect on fragmentation-sensitive insectivorous
bird flocks. Therefore it is important to be able to predict the effects of
such changes on the flocks' space-use patterns.
  3. We use a step selection function (SSF) approach to uncover environmental
drivers behind movement choices. This is used to construct a mechanistic model,
from which we derive predicted utilization distributions (home ranges) of
flocks.
  4. We show that movement decisions are significantly influenced by canopy
height and topography, but not resource depletion and renewal. We quantify the
magnitude of these effects and demonstrate that they are helpful for
understanding various heterogeneous aspects of space use. We compare our
results to recent analytic derivations of space use, demonstrating that they
are only accurate when assuming that there is no persistence in the animals'
movement.
  5. Our model can be translated into other environments or hypothetical
scenarios, such as those given by proposed future anthropogenic actions, to
make predictions of spatial patterns in bird flocks. Furthermore, our approach
is quite general, so could be used to predict the effects of habitat changes on
spatial patterns for a wide variety of animal communities.


  We consider the classical Wright-Fisher model with mutation and selection.
Mutations occur independently in each locus, and selection is performed
according to the sharp peak landscape. In the asymptotic regime studied in [3],
a quasispecies is formed. We find explicitly the distribution of this
quasispecies, which turns out to be the same distribution as for the Moran
model.


  Ecological selection forces for weedy and domesticated traits have influenced
the evolution of seed shape in Setaria resulting in similarity in seed shape
that reflects similarity in ecological function rather than reflecting
phylogenetic relatedness. Seeds from two diploid subspecies of Setaria viridis,
consisting of one weedy subspecies and two races of the domesticated
subspecies, and four other polyploidy weedy species of Setaria. We quantified
seed shape from the silhouettes of the seeds in two separate views. Differences
in shape were compared to ecological role (weed vs. crop) and the evolutionary
trajectory of shape change by phylogenetic grouping from a single reference
species was calculated. Idealized three-dimensional models were created to
examine the differences in shape relative to surface area and volume. All
populations were significantly different in shape, with crops easily
distinguished from weeds, regardless of relatedness between the taxa.
Trajectory of shape change varied by view, but separated crops from weeds and
phylogenetic groupings. Three-dimensional models gave further evidence of
differences in shape reflecting adaptation for environmental exploitation. The
selective forces for weedy and domesticated traits have exceeded phylogenetic
constraints, resulting in seed shape similarity due to ecological role rather
than phylogenetic relatedness. Seed shape and surface-to-volume ratio likely
reflect the importance of the water film that accumulates on the seed surface
when in contact with soil particles. Seed shape may also be a mechanism of
niche separation between taxa.


  The dormancy status of S. faberi seed at abscission was assessed with
reference to tiller and panicle development. Seed from a single genetic line
were grown under field, greenhouse and controlled environment growth chamber
conditions. At abscission, a small fraction (<10%) of S. faberi seed germinated
under favorable conditions. Seed were dissected and germination of caryopses
and embryos also tested. Removal of seed structures exterior to the embryo
increased percentage germination. As the seed rain progressed mean percentage
germination and variation among samples increased. Changes in germination were
correlated with tiller development and relative time of seed maturity within a
panicle. Seed produced on tillers that developed earlier were more likely to be
dormant than seed from later-developing tillers. Seed that matured later on a
panicle were more likely to germinate than seed that matured earlier on the
same panicle. A consistent trend toward later maturing seed having less
dormancy was found for seed grown under different environments which implies an
inherent or parental source for variation in giant foxtail seed dormancy. The
variation in percentage germination at abscission and following stratification
treatments indicates that the S. faberi seed rain consists of individual seeds,
possibly each with a different degree of dormancy.


  This paper uses the reliability polynomial, introduced by Moore and Shannon
in 1956, to analyze the effect of network structure on diffusive dynamics such
as the spread of infectious disease. We exhibit a representation for the
reliability polynomial in terms of what we call {\em structural motifs} that is
well suited for reasoning about the effect of a network's structural properties
on diffusion across the network. We illustrate by deriving several general
results relating graph structure to dynamical phenomena.


  Giant foxtail (Setaria faberi) seeds differ in requirements for germination.
Variable germinability arises during seed development under the influence of
genotype, environment and parent plant. Giant foxtail seed germination has been
shown to be regulated by independent asynchronous or dependent synchronous
action of seed structures. To gain better insight into the process, germination
was divided into axis specific embryo growth categories or states. Three states
were defined for each embryonic axis. The degree of embryo growth (germination
state) after eight to twelve days under germination conditions is believed to
reveal the germinability state (potential for germination) possessed by the
seed before germination. The embryo axes behave independently, which allows any
combination of germination states to occur. In general, the greater the
difference in germination between the axes, the less likely the combination of
states will occur. Photographic evidence of each germination state is shown for
caryopses and seed. Seed with a variety of germinability states is a strategy
for surviving variable environments.


  The limitations of demographic models as well as the opportunities of
evolutionary models are reviewed. Flaws associated with demographic models
include the confounding effects of plant architecture, representation of
heterogeneous individuals in populations, and changes in deme membership
confounding covariance structure. Trait-based evolutionary models include
FoxPatch which represents weedy Setaria spp. seed behavior with explicit life
history process prediction rules and algorithms.


  Variance partitioning methods, which are built upon multivariate statistics,
have been widely applied in different taxa and habitats in community ecology.
Here, I performed a literature review on the development and application of the
methods, and then discussed the limitation of available methods and the
difficulties involved in sampling schemes. The central goal of the work is then
to propose some potential practical methods that might help to overcome
different issues of traditional least-square-based regression modeling. A
variety of regression models has been considered for comparison. In initial
simulations, I identified that generalized additive model (GAM) has the highest
accuracy to predict variation components. Therefore, I argued that other
advanced regression techniques, including the GAM and related models, could be
utilized in variation partitioning for better quantifying the aggregation
scenarios of species distribution.


  Cancer cells can be described as an invasive species that is able to
establish itself in a new environment. The concept of niche construction can be
utilized to describe the process by which cancer cells terraform their
environment, thereby engineering an ecosystem that promotes the genetic fitness
of the species. Ecological dispersion theory can then be utilized to describe
and model the steps and barriers involved in a successful diaspora as the
cancer cells leave the original host organ and migrate to new host organs to
successfully establish a new metastatic community. These ecological concepts
can be further utilized to define new diagnostic and therapeutic areas for
lethal cancers.


  Decisions in a group often result in imitation and aggregation, which are
enhanced in panic, dangerous, stressful or negative situations. Current
explanations of this enhancement are restricted to particular contexts, such as
anti-predatory behavior, deflection of responsibility in humans, or cases in
which the negative situation is associated with an increase in uncertainty. But
this effect is observed across taxa and in very diverse conditions, suggesting
that it may arise from a more general cause, such as a fundamental
characteristic of social decision-making. Current decision-making theories do
not explain it, but we noted that they concentrate on estimating which of the
available options is the best one, implicitly neglecting the cases in which
several options can be good at the same time. We explore a more general model
of decision-making that instead estimates the probability that each option is
good, allowing several options to be good simultaneously. This model predicts
with great generality the enhanced imitation in negative situations. Fish and
human behavioral data showing an increased imitation behavior in negative
circumstances are well described by this type of decisions to choose a good
option.


  Motivation: The availability of thousands of invidual genomes of one species
should boost rapid progress in personalized medicine or understanding of the
interaction between genotype and phenotype, to name a few applications. A key
operation useful in such analyses is aligning sequencing reads against a
collection of genomes, which is costly with the use of existing algorithms due
to their large memory requirements.
  Results: We present MuGI, Multiple Genome Index, which reports all
occurrences of a given pattern, in exact and approximate matching model,
against a collection of thousand(s) genomes. Its unique feature is the small
index size fitting in a standard computer with 16--32\,GB, or even 8\,GB, of
RAM, for the 1000GP collection of 1092 diploid human genomes. The solution is
also fast. For example, the exact matching queries are handled in average time
of 39\,$\mu$s and with up to 3 mismatches in 373\,$\mu$s on the test PC with
the index size of 13.4\,GB. For a smaller index, occupying 7.4\,GB in memory,
the respective times grow to 76\,$\mu$s and 917\,$\mu$s.
  Availability: Software and Suuplementary material:
\url{http://sun.aei.polsl.pl/mugi}.


  To obtain an electron-density map from a macromolecular crystal the
phase-problem needs to be solved, which often involves the use of heavy-atom
derivative crystals and concomitantly the determination of the heavy atom
substructure. This is customarily done by direct methods or Patterson-based
approaches, which however may fail when only poorly diffracting derivative
crystals are available, as often the case for e.g. membrane proteins. Here we
present an approach for heavy atom site identification based on a Molecular
Replacement Parameter Matrix (MRPM) search. It involves an n-dimensional search
to test a wide spectrum of molecular replacement parameters, such as clusters
of different conformations. The result is scored by the ability to identify
heavy-atom positions, from anomalous difference Fourier maps, that allow
meaningful phases to be determined. The strategy was successfully applied in
the determination of a membrane protein structure, the CopA Cu+-ATPase, when
other methods had failed to resolve the heavy atom substructure. MRPM is
particularly suited for proteins undergoing large conformational changes where
multiple search models should be generated, and it enables the identification
of weak but correct molecular replacement solutions with maximum contrast to
prime experimental phasing efforts.


  Motivation: Illumina DNA sequencing is now the predominant source of raw
genomic data, and data volumes are growing rapidly. Bioinformatic analysis
pipelines are having trouble keeping pace. A common bottleneck in such
pipelines is the requirement to read, write, sort and compress large BAM files
multiple times.
  Results: We present SAMBLASTER, a tool that reduces the number of times such
costly operations are performed. SAMBLASTER is designed to mark duplicates in
read-sorted SAM files as a piped post-pass on DNA aligner output before it is
compressed to BAM. In addition, it can simultaneously output into separate
files the discordant read-pairs and/or split-read mappings used for structural
variant calling. As an alignment post-pass, its own runtime overhead is
negligible, while dramatically reducing overall pipeline complexity and
runtime. As a stand-alone duplicate marking tool, it performs significantly
better than PICARD or SAMBAMBA in terms of both speed and memory usage, while
achieving nearly identical results.
  Availability: SAMBLASTER is open source C++ code and freely available from
https://github.com/GregoryFaust/samblaster


  Organisms have increased in complexity through a series of major evolutionary
transitions, in which formerly autonomous entities become parts of a novel
higher-level entity. One intriguing feature of the higher-level entity after
some major transitions is a division of reproductive labor among its
lower-level units. Although it can have clear benefits once established, it is
unknown how such reproductive division of labor originates. We consider a
recent evolution experiment on the yeast Saccharomyces cerevisiae as a unique
platform to address the issue of reproductive differentiation during an
evolutionary transition in individuality. In the experiment, independent yeast
lineages evolved a multicellular "snowflake-like'' cluster form in response to
gravity selection. Shortly after the evolution of clusters, the yeast evolved
higher rates of cell death. While cell death enables clusters to split apart
and form new groups, it also reduces their performance in the face of gravity
selection. To understand the selective value of increased cell death, we create
a mathematical model of the cellular arrangement within snowflake yeast
clusters. The model reveals that the mechanism of cell death and the geometry
of the snowflake interact in complex, evolutionarily important ways. We find
that the organization of snowflake yeast imposes powerful limitations on the
available space for new cell growth. By dying more frequently, cells in
clusters avoid encountering space limitations, and, paradoxically, reach higher
numbers. In addition, selection for particular group sizes can explain the
increased rate of apoptosis both in terms of total cell number and total
numbers of collectives. Thus, by considering the geometry of a primitive
multicellular organism we can gain insight into the initial emergence of
reproductive division of labor during an evolutionary transition in
individuality.


  Understanding protein folding has been one of the great challenges in
biochemistry and molecular biophysics. Over the past 50 years, many
thermodynamic and kinetic studies have been performed addressing the stability
of globular proteins. In comparison, advances in the membrane protein folding
field lag far behind. Although membrane proteins constitute about a third of
the proteins encoded in known genomes, stability studies on membrane proteins
have been impaired due to experimental limitations. Furthermore, no systematic
experimental strategies are available for folding these biomolecules in vitro.
Common denaturing agents such as chaotropes usually do not work on helical
membrane proteins, and ionic detergents have been successful denaturants only
in few cases. Refolding a membrane protein seems to be a craftsman work, which
is relatively straightforward for transmembrane {\beta}-barrel proteins but
challenging for {\alpha}-helical membrane proteins. Additional complexities
emerge in multidomain membrane proteins, data interpretation being one of the
most critical. In this review, we will describe some recent efforts in
understanding the folding mechanism of membrane proteins that have been
reversibly refolded allowing both thermodynamic and kinetic analysis. This
information will be discussed in the context of current paradigms in the
protein folding field.


  Ecologists are interested in understanding and predicting how ecological
communities change through time. While it might seem natural to measure this
through changes in species' abundances, computational limitations mean
transitions between community types are often modelled instead. We present an
approach inspired by DNA substitution models that attempts to estimate historic
interactions between species, and thus estimate turnover rates in ecological
communities. Although our simulations show that the method has some
limitations, our application to butterfly community data shows the method can
detect signal in real data. Open source C++ code implementing the method is
available at http://www.github.com/willpearse/lotto.


  How cytoskeletal filaments collectively undergo growth and shrinkage is an
intriguing question. Collective properties of multiple bio-filaments (actin or
microtubules) undergoing hydrolysis, have not been studied extensively earlier,
within simple theoretical frameworks. In this paper, we show that collective
properties of multiple filaments under force are very distinct from the
properties of a single filament under similar conditions -- these distinctions
manifest as follows: (i) the collapse time during collective catastrophe for a
multifilament system is much larger than that of a single filament with the
same average length, (ii) force-dependence of the cap-size distribution of
multiple filaments are quantitatively different from that of single filament,
(iii) the diffusion constant associated with the system length fluctuations is
distinct for multiple filaments, (iv) switching dynamics of multiple filaments
between capped and uncapped states and the fluctuations therein are also
distinct. We build a unified picture by establishing interconnections among all
these collective phenomena. Additionally, we show that the collapse times
during catastrophes can be sharp indicators of collective stall forces
exceeding the additive contributions of single filaments.


  Genetic variation in a population can sometimes arise so fast as to modify
ecosystem dynamics. Such phenomena have been observed in natural predator-prey
systems, and characterized in the laboratory as showing unusual phase
relationships in population dynamics, including a $\pi$ phase shift between
predator and prey (evolutionary cycles) and even undetectable prey oscillations
compared to those of the predator (cryptic cycles). Here we present a generic
individual-level stochastic model of interacting populations that includes a
subpopulation of low nutritional value to the predator. Using a master equation
formalism, and by mapping to a coherent state path integral solved by a
system-size expansion, we show that evolutionary and cryptic quasicycles can
emerge generically from the combination of intrinsic demographic fluctuations
and clonal mutations alone, without additional biological mechanisms.


  This work presents a modified Kedem-Katchalsky equations for osmosis through
nano-pore. osmotic reflection coefficient of a solute was found to be chiefly
affected by the entrance of the pore while filtration reflection coefficient
can be affected by both the entrance and the internal structure of the pore.
Using an analytical method, we get the quantitative relationship between
osmotic reflection coefficient and the molecule size. The model is verified by
comparing the theoretical results with the reported experimental data of
aquaporin osmosis. Our work is expected to pave the way for a better
understanding of osmosis in bio-system and to give us new ideas in designing
new membranes with better performance.


  We show that disease transmission models in a spatially heterogeneous
environment can have a large number of coexisting endemic equilibria. A general
compartmental model is considered to describe the spread of an infectious
disease in a population distributed over several patches. For disconnected
regions, many boundary equilibria may exist with mixed disease free and endemic
components, but these steady states usually disappear in the presence of
spatial dispersal. However, if backward bifurcations can occur in the regions,
some partially endemic equilibria of the disconnected system move into the
interior of the nonnegative cone and persist with the introduction of mobility
between the patches. We provide a mathematical procedure that precisely
describes in terms of the local reproduction numbers and the connectivity
network of the patches, whether a steady state of the disconnected system is
preserved or ceases to exist for low volumes of travel. Our results are
illustrated on a patchy HIV transmission model with subthreshold endemic
equilibria and backward bifurcation. We demonstrate the rich dynamical behavior
(i.e., creation and destruction of steady states) and the presence of multiple
stable endemic equilibria for various connection networks.


  Molecular docking is a central method in the computer-based screening of
compound libraries as a part of the rational approach to drug design. Although
the method has proved its competence in predicting binding modes correctly, its
inherent complexity puts high demands on computational resources. Moreover the
chemical space to be screened is prohibitively large. Therefore the application
of filtering prior to docking is a promising concept. We implemented a
pre-docking filter based on the tangent distance algorithm originally conceived
for optical character recognition. The challenging transfer of the method from
two-dimensional to three-dimensional data was achieved by representing the
molecular structure by a set of density maps extracted from different views of
the compound. Additionally, our program applies a binary classification using
principal component analysis. Ligand and binding pocket are aligned according
to their centroidal axes, enabling a size-based filtering for the purpose of
enriching the dataset regarding ligands before docking. The evaluation of our
program via redocking produced RMSD values between 8{\AA} and 25{\AA},
indicating that the tangent distance approach is not suited for optimizing the
orientation of a ligand and binding pocket. Investigating probable explanations
lead to the conclusion that a likely cause for these results is the method's
known inability to approximate large transformations. A validation of the
principal component analysis alone performed better: Tests on a dataset of 170
ligands and 6,435 decoys yielded a sensitivity of 0.81, while keeping the
runtime within a reasonable timeframe (1 to 4 seconds). The dataset's
enrichment increased from 2.64% to 2.82%.


  Rod-like bacteria maintain their cylindrical shapes with remarkable precision
during growth. However, they are also capable to adapt their shapes to external
forces and constraints, for example by growing into narrow or curved
confinements. Despite being one of the simplest morphologies, we are still far
from a full understanding of how shape is robustly regulated, and how bacteria
obtain their near-perfect cylindrical shapes with excellent precision. However,
recent experimental and theoretical findings suggest that cell-wall geometry
and mechanical stress play important roles in regulating cell shape in rod-like
bacteria. We review our current understanding of the cell wall architecture and
the growth dynamics, and discuss possible candidates for regulatory cues of
shape regulation in the absence or presence of external constraints. Finally,
we suggest further future experimental and theoretical directions, which may
help to shed light on this fundamental problem.


  Metabolite analysis of peritoneal dialysis (PD) effluent may provide
information regarding onset and progression of complications associated with
prolonged PD therapy. In this context, the NMR detectable small metabolites of
PD effluent samples were characterized using high resolution 1H and 1H-13C NMR
spectroscopy. The various spectra were recorded (at 800 MHz proton frequency)
on PD effluent samples obtained after 4 hour (intraperitoneal) dwell time from
patients with end stage renal failure (ESRF) and continuing normally on PD
therapy. Inspite of devastating spectral feature of PD effluent due to the
presence of intense resonances from glucose and lactate, we were able to
identify about 53 small endogenous metabolites (including many complex coupled
spin systems) and more than 90 % of the total CH cross peaks of 1H-13C HSQC
spectrum were identified specific to various metabolites of PD effluent. We
foresee that the characteristic fingerprints of various metabolites of control
PD effluent samples will be used to identify and distinguish metabolic
differences from PD related complications.


  HIV-1 protease represents an appealing system for directed enzyme re-design,
since it has various different endogenous targets, a relatively simple
structure and it is well studied. Recently Chaudhury and Gray (Structure (2009)
17: 1636 -- 1648) published a computational algorithm to discern the
specificity determining residues of HIV-1 protease. In this paper we present
two computational tools aimed at re-designing HIV-1 protease, derived from the
algorithm of Chaudhuri and Gray. First, we present an energy-only based
methodology to discriminate cleavable and non cleavable peptides for HIV-1
proteases, both wild type and mutant. Secondly, we show an algorithm we
developed to predict mutant HIV-1 proteases capable of cleaving a new target
substrate peptide, different from the natural targets of HIV-1 protease. The
obtained in silico mutant enzymes were analyzed in terms of cleavability and
specificity towards the target peptide using the energy-only methodology. We
found two mutant proteases as best candidates for specificity and cleavability
towards the target sequence.


  How multicellular life forms evolved out from unicellular ones constitutes a
major problem in our understanding of the evolution of our biosphere. A recent
set of experiments involving yeast cell populations has shown that selection
for faster sedimenting cells leads to the appearance of stable aggregates of
cells that are able to split into smaller clusters. It was suggested that the
observed evolutionary patterns could be the result of evolved programs
affecting cell death. Here we show, using a simple model of cell-cell
interactions and evolving adhesion rates, that the observed patterns in cluster
size and localized mortality can be easily interpreted in terms of waste
accumulation and toxicity driven apoptosis. This simple mechanism would have
played a key role in the early evolution of multicellular life forms based on
both aggregative and clonal development. The potential extensions of this work
and its implications for natural and synthetic multicellularity are discussed.


  Energy landscapes provide a valuable means for studying the folding dynamics
of short RNA molecules in detail by modeling all possible structures and their
transitions. Higher abstraction levels based on a macro-state decomposition of
the landscape enable the study of larger systems, however they are still
restricted by huge memory requirements of exact approaches.
  We present a highly parallelizable local enumeration scheme that enables the
computation of exact macro-state transition models with highly reduced memory
requirements. The approach is evaluated on RNA secondary structure landscapes
using a gradient basin definition for macro-states. Furthermore, we demonstrate
the need for exact transition models by comparing two barrier-based appoaches
and perform a detailed investigation of gradient basins in RNA energy
landscapes.
  Source code is part of the C++ Energy Landscape Library available at
http://www.bioinf.uni-freiburg.de/Software/.


  Wells-next-the-Sea and Cromer in Norfolk (England) both rely upon their local
crab populations, since crabbing (gillying) is a major part of their tourist
industry. Compared to a control site with no crabbing, crabs from Wells harbour
and Cromer pier were found to have nearly six times the amount of limb damage.
Crabs caught by the general public had more injuries than crabs caught in
controlled conditions, suggesting the buckets in which the crabs were kept were
to blame. Since there is much evidence that such injuries have negative impacts
on the survival and reproductive success of the shore crab, this is taken as
evidence of non-lethal injury from humans having a population-level effect on
these animals. Questionnaire data demonstrated a public lack of awareness and
want for information, which was then used to obtain funding to produce a
leaflet campaign informing the public of how to crab responsibly. All data
collected is available online at http://dx.doi.org/10.6084/m9.figshare.979288.


  Escherichia coli DNA polymerase V (pol V), a heterotrimeric complex composed
of UmuD'2C, is marginally active. ATP and RecA play essential roles in the
activation of pol V for DNA synthesis including translesion synthesis (TLS). We
have established three features of the roles of ATP and RecA. 1) RecA-activated
DNA polymerase V (pol V Mut), is a DNA-dependent ATPase; 2) bound ATP is
required for DNA synthesis; 3) pol V Mut function is regulated by ATP, with ATP
required to bind primer/template (p/t) DNA and ATP hydrolysis triggering
dissociation from the DNA. Pol V Mut formed with an ATPase-deficient RecA
E38K/K72R mutant hydrolyzes ATP rapidly, establishing the DNA-dependent ATPase
as an intrinsic property of pol V Mut distinct from the ATP hydrolytic activity
of RecA when bound to single-stranded (ss)DNA as a nucleoprotein filament
(RecA*). No similar ATPase activity or autoregulatory mechanism has previously
been found for a DNA polymerase.


  A main challenge of data-driven sciences is how to make maximal use of the
progressively expanding databases of experimental datasets in order to keep
research cumulative. We introduce the idea of a modeling-based dataset
retrieval engine designed for relating a researcher's experimental dataset to
earlier work in the field. The search is (i) data-driven to enable new
findings, going beyond the state of the art of keyword searches in annotations,
(ii) modeling-driven, to both include biological knowledge and insights learned
from data, and (iii) scalable, as it is accomplished without building one
unified grand model of all data. Assuming each dataset has been modeled
beforehand, by the researchers or by database managers, we apply a rapidly
computable and optimizable combination model to decompose a new dataset into
contributions from earlier relevant models. By using the data-driven
decomposition we identify a network of interrelated datasets from a large
annotated human gene expression atlas. While tissue type and disease were major
driving forces for determining relevant datasets, the found relationships were
richer and the model-based search was more accurate than keyword search; it
moreover recovered biologically meaningful relationships that are not
straightforwardly visible from annotations, for instance, between cells in
different developmental stages such as thymocytes and T-cells. Data-driven
links and citations matched to a large extent; the data-driven links even
uncovered corrections to the publication data, as two of the most linked
datasets were not highly cited and turned out to have wrong publication entries
in the database.


  Reciprocity is firmly established as an important mechanism that promotes
cooperation. An efficient information exchange is likewise important,
especially on structured populations, where interactions between players are
limited. Motivated by these two facts, we explore the role of facilitators in
social dilemmas on networks. Facilitators are here mirrors to their neighbors
-- they cooperate with cooperators and defect with defectors -- but they do not
participate in the exchange of strategies. As such, in addition to introducing
direct reciprocity, they also obstruct information exchange. In well-mixed
populations, facilitators favor the replacement and invasion of defection by
cooperation as long as their number exceeds a critical value. In structured
populations, on the other hand, there exists a delicate balance between the
benefits of reciprocity and the deterioration of information exchange.
Extensive Monte Carlo simulations of social dilemmas on various interaction
networks reveal that there exists an optimal interplay between reciprocity and
information exchange, which sets in only when a small number of facilitators
occupies the main hubs of the scale-free network. The drawbacks of missing
cooperative hubs are more than compensated by reciprocity and, at the same
time, the compromised information exchange is routed via the auxiliary hubs
with only marginal losses in effectivity. These results indicate that it is not
always optimal for the main hubs to become ''leaders of the masses'', but
rather to exploit their highly connected state to promote tit-for-tat-like
behavior.


  We apply a recently developed model of cytoskeletal force generation to study
a cell intrinsic contractility, as well as its response to external loading.
The model is based on a non-equilibrium thermodynamic treatment of the
mechano-chemistry governing force in the stress fiber-focal adhesion system.
Our computational study suggests that the mechanical coupling between the
stress fibers and focal adhesions leads to a complex, dynamic, mechano-chemical
response. We collect the results in response maps whose regimes are
distinguished by the initial geometry of the stress fiber-focal adhesion
system, and by the external load on the cell. The results from our model
connect qualitatively with recent studies on the force response of smooth
muscle cells on arrays of polymeric microposts (Mann et al., Lab. on a Chip,
12, 731-740, 2012).


  In this communication, we propose a model to study the non-equilibrium
process by which actin stress fibers develop force in contractile cells. The
emphasis here is on the non-equilibrium thermodynamics, which is necessary to
address the mechanics as well as the chemistry of dynamic cell contractility.
In this setting we are able to develop a framework that relates (a) the
dynamics of force generation within the cell and (b) the cell response to
external stimuli to the chemical processes occurring within the cell, as well
as to the mechanics of linkage between the stress fibers, focal adhesions and
extra-cellular matrix.


  The multiple futile cycle is an important building block in networks of
chemical reactions arising in molecular biology. A typical process which it
describes is the addition of n phosphate groups to a protein. It can be
modelled by a system of ordinary differential equations depending on
parameters. The special case n = 2 is called the dual futile cycle. The main
result of this paper is a proof that there are parameter values for which the
system of ODE describing the dual futile cycle has two distinct stable
stationary solutions. The proof is based on bifurcation theory and geometric
singular perturbation theory. An important entity built of three coupled
multiple futile cycles is the MAPK cascade. It is explained how the ideas used
to prove bistability for the dual futile cycle might help to prove the
existence of periodic solutions for the MAPK cascade.


  The current biochemical information processing systems behave in a
predetermined manner because all features are defined during the design phase.
To make such unconventional computing systems reusable and programmable for
biomedical applications, adaptation, learning, and self-modification based on
external stimuli would be highly desirable. However, so far, it has been too
challenging to implement these in wet chemistries. In this paper we extend the
chemical perceptron, a model previously proposed by the authors, to function as
an analog instead of a binary system. The new analog asymmetric signal
perceptron learns through feedback and supports Michaelis-Menten kinetics. The
results show that our perceptron is able to learn linear and nonlinear
(quadratic) functions of two inputs. To the best of our knowledge, it is the
first simulated chemical system capable of doing so. The small number of
species and reactions and their simplicity allows for a mapping to an actual
wet implementation using DNA-strand displacement or deoxyribozymes. Our results
are an important step toward actual biochemical systems that can learn and
adapt.


  ANIMO (Analysis of Networks with Interactive MOdeling) is a software for
modeling biological networks, such as e.g. signaling, metabolic or gene
networks. An ANIMO model is essentially the sum of a network topology and a
number of interaction parameters. The topology describes the interactions
between biological entities in form of a graph, while the parameters determine
the speed of occurrence of such interactions. When a mismatch is observed
between the behavior of an ANIMO model and experimental data, we want to update
the model so that it explains the new data. In general, the topology of a model
can be expanded with new (known or hypothetical) nodes, and enables it to match
experimental data. However, the unrestrained addition of new parts to a model
causes two problems: models can become too complex too fast, to the point of
being intractable, and too many parts marked as "hypothetical" or "not known"
make a model unrealistic. Even if changing the topology is normally the easier
task, these problems push us to try a better parameter fit as a first step, and
resort to modifying the model topology only as a last resource. In this paper
we show the support added in ANIMO to ease the task of expanding the knowledge
on biological networks, concentrating in particular on the parameter settings.


  Proteins belonging to immunoglobulin superfamily(IgSF) show remarkably
conserved nature both in their folded structure and in their folding process,
but they neither originate from very similar sequences nor demonstrate
functional conservation. Treating proteins as fractal objects, without studying
spatial conservation in positioning of particular residues in IgSF, this work
probed the roots structural invariance of immunoglobulins(Ig). Symmetry in
distribution of mass, hydrophobicity, polarizability recorded very similar
extents in Ig and in structurally-closest non-Ig structures. They registered
similar symmetries in dipole-dipole, {\pi}-{\pi}, cation-{\pi} cloud
interactions and also in distribution of active chiral centers, charged
residues and hydrophobic residues. But in contrast to non-Ig proteins, extents
of residual interaction symmetries in Ig.s of largely varying sizes are found
to converge to exactly same magnitude of correlation dimension - these are
named 'structural attractors', who's weightages depend on ensuring exact
convergence of pairwise-interaction symmetries to attractor magnitude. Small
basin of attraction for Ig attractors explained the strict and consistent
quality control in ensuring stability and functionality of IgSF proteins. Low
dependency of attractor weightage on attractor magnitude demonstrated that
residual-interaction symmetry with less pervasive nature can also be crucial in
ensuring Ig stability.


  In this paper we propose a novel theoretical framework for interpreting
long-range dynamical correlations unveiled in proteins through NMR
measurements. The theoretical rationale relies on the hypothesis that
correlated motions in proteins may be reconstructed as large-scale, collective
modes sustained by long-lived nonlinear vibrations known as discrete breathers
(DB) localized at key, hot-spot sites. DBs are spatially localized modes, whose
nonlinear nature hinders resonant coupling with the normal modes, thus
conferring them long lifetimes as compared to normal modes. DBs have been
predicted to exist in proteins, localized at few $hotspot$ residues typically
within the stiffest portions of the structure. We compute DB modes analytically
in the framework of the nonlinear network model, showing that the displacement
patterns of many DBs localized at key sites match to a remarkable extent the
experimentally uncovered correlation blueprint. The computed dispersion
relations prove that it is physically possible for some of these DBs to be
excited out of thermal fluctuations at room temperature.
  Based on our calculations, we speculate that transient energy redistribution
among the vibrational modes in a protein might favor the emergence of DB-like
bursts of long-lived energy at hot-spot sites with lifetimes in the ns range,
able to sustain critical, function-encoding correlated motions. More generally,
our calculations provide a novel quantitative tool to predict fold-spanning
dynamical pathways of correlated residues that may be central to allosteric
cross-talk in proteins.


  We have developed a biosensor based on BSA with the help of metal ions
binding mechanism to detect and remove inorganic As(III), Cu(II), Pb(II) from
water like fishing by hooking system.


  HiRE-RNA is a simplified, coarse-grained RNA model for the prediction of
equilibrium configurations, dynamics and thermodynamics. Using a reduced set of
particles and detailed interactions accounting for base-pairing and stacking we
show that non-canonical and multiple base interactions are necessary to capture
the full physical behavior of complex RNAs. In this paper we give a full
account of the model and we present results on the folding, stability and free
energy surfaces of 16 systems with 12 to 76 nucleotides of increasingly complex
architectures, ranging from monomers to dimers, using a total of 850$\mu$s
simulation time.


  Most approaches to prediction of protein function from primary structure are
based on similarity between the query sequence and sequences of known function.
This approach, however, disregards the occurrence of gene duplication
(paralogy) or convergent evolution of the genes. The analysis of correlated
proteins that share a common domain, taking into consideration the evolutionary
history of genes under study may provide a more reliable annotation of
predicted proteins with unknown function. A computer program that enables
real-time comparison of 'gene trees' with 'species trees' was developed. The
Phylogenetic Genome Annotator (PGA) performs a profile based multiple sequence
alignment of a set of sequences that share a common domain to generate a
phylogenetic gene tree, which is compared to the species phylogeny inferred
from aligned ribossomal RNA data. The correlated protein domains are then
displayed side-by-side with the phylogeny of the corresponding species. The
statistical support of gene clusters (branches) is given by the quartet
puzzling method. This analysis readily discriminates paralogs from orthologs,
enabling the identification of proteins originated by gene duplications and the
prediction of possible functional divergence in groups of similar sequences.
The tool was tested in three distinct subfamilies of the G-protein coupled
receptor superfamily. In the analysed datasets, the paralogy prediction agreed
with the known subfamily grouping, suggesting that subfamily divergence was
facilitated by duplication events in the ancestral nodes.


  Rotavirus (RV) is the major etiologic agent of severe infantile
gastroenteritis; its genome has 11 segments of double stranded RNA, encoding 12
proteins. The non-structural protein 4 (NSP4) encoded by segment 10 is
multifunctional. The aim of this study is to analyze the selective pressure
driving the NSP4 of RV, through the ratio of non-synonymous substitutions per
synonymous substitutions (dN/dS). Our results show that NSP4 is under negative
evolutionary pressure (84.57% of the amino acid sequence) and no site was found
under positive selection. This may support other evolutionary studies of
different RV proteins or viral agents.


  The Phylogenetic Genome Annotator (PGA) is a computer program that enables
real-time comparison of 'gene trees' versus 'species trees' obtained from
predicted open reading frames of whole genome data. The gene phylogenies are
inferred for each individual genome predicted proteins whereas the species
phylogenies are inferred from rDNA data. The correlated protein domains,
defined by PFAM, are then displayed side-by-side with a phylogeny of the
corresponding species. The statistical support of gene clusters (branches) is
given by the quartet puzzling method. This analysis readily discriminates
paralogs from orthologs, enabling the identification of proteins originated by
gene duplications and the prediction of possible functional divergence in
groups of similar sequences.


  Topological constraints can affect both equilibrium and dynamic properties of
polymer systems, and can play a role in the organization of chromosomes.
Despite many theoretical studies, the effects of topological constraints on the
equilibrium state of a single compact polymer have not been systematically
studied. Here we use simulations to address this longstanding problem. We find
that sufficiently long unknotted polymers differ from knotted ones in the
spatial and topological states of their subchains. The unknotted globule has
subchains that are mostly unknotted and form asymptotically compact $R_G(s)
\sim s^{1/3}$ crumples. However, crumples display high fractal dimension of the
surface $d_b = 2.8$, forming excessive contacts and interpenetrating each
other. We conclude that this topologically constrained equilibrium state
resembles a conjectured crumpled globule [Grosberg et al., Journal de Physique,
1988, 49, 2095], but differs from its idealized hierarchy of self-similar,
isolated and compact crumples.


  Motivation: Whole-genome high-coverage sequencing has been widely used for
personal and cancer genomics as well as in various research areas. However, in
the lack of an unbiased whole-genome truth set, the global error rate of
variant calls and the leading causal artifacts still remain unclear even given
the great efforts in the evaluation of variant calling methods.
  Results: We made ten SNP and INDEL call sets with two read mappers and five
variant callers, both on a haploid human genome and a diploid genome at a
similar coverage. By investigating false heterozygous calls in the haploid
genome, we identified the erroneous realignment in low-complexity regions and
the incomplete reference genome with respect to the sample as the two major
sources of errors, which press for continued improvements in these two areas.
We estimated that the error rate of raw genotype calls is as high as 1 in
10-15kb, but the error rate of post-filtered calls is reduced to 1 in 100-200kb
without significant compromise on the sensitivity.
  Availability: BWA-MEM alignment: http://bit.ly/1g8XqRt; Scripts:
https://github.com/lh3/varcmp; Additional data:
https://figshare.com/articles/Towards_better_understanding_of_artifacts_in_variating_calling_from_high_coverage_samples/981073


  Measures of biodiversity change such as the Living Planet Index describe
proportional change in the abundance of a typical species, which can be thought
of as change in the size of a community. Here, I discuss the orthogonal concept
of change in relative abundances, which I refer to as shape change. To be
logically consistent, a measure of the rate of shape change should be scaling
invariant (have the same value for all data with the same vector of
proportional change over a given time interval), but existing measures do not
have this property. I derive a new, scaling invariant measure. I show that this
new measure and existing measures of biodiversity change such as the Living
Planet Index describe different aspects of dynamics. I show that neither body
size nor environmental variability need affect the rate of shape change. I
extend the measure to deal with colonizations and extinctions, using the
surreal number system. I give examples using data on hoverflies in a garden in
Leicester, UK, and the higher plant community of Surtsey. I hypothesize that
phylogenetically-restricted assemblages will show a higher proportion of size
change than diverse communities.


  Evolutionary game theory is a successful mathematical framework geared
towards understanding the selective pressures that affect the evolution of the
strategies of agents engaged in interactions with potential conflicts. While a
mathematical treatment of the costs and benefits of decisions can predict the
optimal strategy in simple settings, more realistic settings such as finite
populations, non-vanishing mutations rates, stochastic decisions, communication
between agents, and spatial interactions, require agent-based methods where
each agent is modeled as an individual, carries its own genes that determine
its decisions, and where the evolutionary outcome can only be ascertained by
evolving the population of agents forward in time. While highlighting standard
mathematical results, we compare those to agent-based methods that can go
beyond the limitations of equations and simulate the complexity of
heterogeneous populations and an ever-changing set of interactors. We conclude
that agent-based methods can predict evolutionary outcomes where purely
mathematical treatments cannot tread (for example in the weak selection--strong
mutation limit), but that mathematics is crucial to validate the computational
simulations.


  When polygenic traits are under stabilizing selection, many different
combinations of alleles allow close adaptation to the optimum. If alleles have
equal effects, all combinations that result in the same deviation from the
optimum are equivalent. Furthermore, the genetic variance that is maintained by
mutation-selection balance is $2 \mu/S$ per locus, where $\mu$ is the mutation
rate and $S$ the strength of stabilizing selection. In reality, alleles vary in
their effects, making the fitness landscape asymmetric, and complicating
analysis of the equilibria. We show that that the resulting genetic variance
depends on the fraction of alleles near fixation, which contribute by $2
\mu/S$, and on the total mutational effects of alleles that are at intermediate
frequency. The interplay between stabilizing selection and mutation leads to a
sharp transition: alleles with effects smaller than a threshold value of
$2\sqrt{\mu / S}$ remain polymorphic, whereas those with larger effects are
fixed. The genetic load in equilibrium is less than for traits of equal
effects, and the fitness equilibria are more similar. We find that if the
optimum is displaced, alleles with effects close to the threshold value sweep
first, and their rate of increase is bounded by $\sqrt{\mu S}$. Long term
response leads in general to well-adapted traits, unlike the case of equal
effects that often end up at a sub-optimal fitness peak. However, the
particular peaks to which the populations converge are extremely sensitive to
the initial states, and to the speed of the shift of the optimum trait value.


  Biological sensory systems react to changes in their surroundings. They are
characterized by fast response and slow adaptation to varying environmental
cues. Insofar as sensory adaptive systems map environmental changes to changes
of their internal degrees of freedom, they can be regarded as computational
devices manipulating information. Landauer established that information is
ultimately physical, and its manipulation subject to the entropic and energetic
bounds of thermodynamics. Thus the fundamental costs of biological sensory
adaptation can be elucidated by tracking how the information the system has
about its environment is altered. These bounds are particularly relevant for
small organisms, which unlike everyday computers operate at very low energies.
In this paper, we establish a general framework for the thermodynamics of
information processing in sensing. With it, we quantify how during sensory
adaptation information about the past is erased, while information about the
present is gathered. This process produces entropy larger than the amount of
old information erased and has an energetic cost bounded by the amount of new
information written to memory. We apply these principles to the E. coli's
chemotaxis pathway during binary ligand concentration changes. In this regime,
we quantify the amount of information stored by each methyl group and show that
receptors consume energy in the range of the information-theoretic minimum. Our
work provides a basis for further inquiries into more complex phenomena, such
as gradient sensing and frequency response.


  Metagenomics characterizes microbial communities by random shotgun sequencing
of DNA isolated directly from an environment of interest. An essential step in
computational metagenome analysis is taxonomic sequence assignment, which
allows us to identify the sequenced community members and to reconstruct
taxonomic bins with sequence data for the individual taxa. We describe an
algorithm and the accompanying software, taxator-tk, which performs taxonomic
sequence assignments by fast approximate determination of evolutionary
neighbors from sequence similarities. Taxator-tk was precise in its taxonomic
assignment across all ranks and taxa for a range of evolutionary distances and
for short sequences. In addition to the taxonomic binning of metagenomes, it is
well suited for profiling microbial communities from metagenome samples
becauseit identifies bacterial, archaeal and eukaryotic community members
without being affected by varying primer binding strengths, as in marker gene
amplification, or copy number variations of marker genes across different taxa.
Taxator-tk has an efficient, parallelized implementation that allows the
assignment of 6 Gb of sequence data per day on a standard multiprocessor system
with ten CPU cores and microbial RefSeq as the genomic reference data.


  Despite significant progress in structural and functional characterization of
human genome, understanding of mechanisms underlying the genetic basis of human
phenotypic uniqueness remains limited. We report that non-randomly distributed
transposable element-derived sequences, most notably HERV-H/LTR7 and L1HS, are
associated with creation of 99.8% unique to human transcription factor binding
sites in genome of embryonic stem cells (ESC). 4,094 unique to human regulatory
loci display selective and site-specific binding of critical regulators (NANOG,
POU5F1, CTCF, Lamin B1) and are preferentially placed within the matrix of
transcriptionally active DNA segments hyper-methylated in ESC. Unique to human
NANOG-binding sites are enriched near the rapidly evolving in primates
protein-coding genes regulating brain size, pluripotency lncRNAs, hESC
enhancers, and 5-hydroxymethylcytosine-harboring regions immediately adjacent
to binding sites. We propose a proximity placement model explaining how 33-47%
excess of NANOG and POU5F1 proteins immobilized on a DNA scaffold may play a
functional role at distal regulatory elements.


  The genotype-fitness map plays a fundamental role in shaping the dynamics of
evolution. However, it is difficult to directly measure a fitness landscape in
practice, because the number of possible genotypes is astronomical. One
approach is to sample as many genotypes as possible, measure their fitnesses,
and fit a statistical model of the landscape that includes additive and
pairwise interactive effects between loci. Here we elucidate the pitfalls of
using such regressions, by studying artificial but mathematically convenient
fitness landscapes. We identify two sources of bias inherent in these
regression procedures that each tends to under-estimate high fitnesses and
over-estimate low fitnesses. We characterize these biases for random sampling
of genotypes, as well as for samples drawn from a population under selection in
the Wright-Fisher model of evolutionary dynamics. We show that common measures
of epistasis, such as the number of monotonically increasing paths between
ancestral and derived genotypes, the prevalence of sign epistasis, and the
number of local fitness maxima, are distorted in the inferred landscape. As a
result, the inferred landscape will provide systematically biased predictions
for the dynamics of adaptation. We identify the same biases in a computational
RNA-folding landscape, as well as in regulatory sequence binding data, treated
with the same fitting procedure. Finally, we present a method that may
ameliorate these biases in some cases.


  Evolution of cooperation in the prisoner's dilemma and the public goods game
is studied, where initially players belong to two independent structured
populations. Simultaneously with the strategy evolution, players whose current
utility exceeds a threshold are rewarded by an external link to a player
belonging to the other population. Yet as soon as the utility drops below the
threshold, the external link is terminated. The rewarding of current
evolutionary fitness thus introduces a time-varying interdependence between the
two populations. We show that, regardless of the details of the evolutionary
game and the interaction structure, the self-organization of fitness and reward
gives rise to distinguished players that act as strong catalysts of cooperative
behavior. However, there also exist critical utility thresholds beyond which
distinguished players are no longer able to percolate. The interdependence
between the two populations then vanishes, and cooperators are forced to rely
on traditional network reciprocity alone. We thus demonstrate that a simple
strategy-independent form of rewarding may significantly expand the scope of
cooperation on structured populations. The formation of links outside the
immediate community seems particularly applicable in human societies, where an
individual is typically member in many different social networks.


  The statistical mechanics of polymers grafted on surfaces has been the
subject of intense research activity because of many potential applications. In
this paper, we analytically investigate the conformational changes caused by a
single cross-link on two ideal (Gaussian) chains grafted on a rigid planar
surface. Both the cross-link and the surface reduce the number of allowed
configurations. In the absence of the hard substrate, the sole effect of the
cross-link is a reduction in the effective Kuhn length of a tethered chain. The
cross-link induced shrinkage (collapse) of the grafted chains (mushrooms) turns
out to be a reduction in the variance of the distribution of the height of the
chain rather than a reduction of the height itself.


  Chemistry as an unconventional computing medium presently lacks a systematic
approach to gather, store, and sort data over time. To build more complicated
systems in chemistries, the ability to look at data in the past would be a
valuable tool to perform complex calculations. In this paper we present the
first implementation of a chemical delay line providing information storage in
a chemistry that can reliably capture information over an extended period of
time. The delay line is capable of parallel operations in a single instruction,
multiple data (SIMD) fashion.
  Using Michaelis-Menten kinetics, we describe the chemical delay line
implementation featuring an enzyme acting as a means to reduce copy errors. We
also discuss how information is randomly accessible from any element on the
delay line. Our work shows how the chemical delay line retains and provides a
value from a previous cycle. The system's modularity allows for integration
with existing chemical systems. We exemplify the delay line capabilities by
integration with a threshold asymmetric signal perceptron to demonstrate how it
learns all 14 linearly separable binary functions over a size two sliding
window. The delay line has applications in biomedical diagnosis and treatment,
such as smart drug delivery.


  In this and a companion paper we outline a general framework for the
thermodynamic description of open chemical reaction networks, with special
regard to metabolic networks regulating cellular physiology and biochemical
functions. We first introduce closed networks "in a box", whose thermodynamics
is subjected to strict physical constraints: the mass-action law, elementarity
of processes, and detailed balance. We further digress on the role of solvents
and on the seemingly unacknowledged property of network independence of free
energy landscapes. We then open the system by assuming that the concentrations
of certain substrate species (the chemostats) are fixed, whether because
promptly regulated by the environment via contact with reservoirs, or because
nearly constant in a time window. As a result, the system is driven out of
equilibrium. A rich algebraic and topological structure ensues in the network
of internal species: Emergent irreversible cycles are associated to
nonvanishing affinities, whose symmetries are dictated by the breakage of
conservation laws. These central results are resumed in the relation $a + b =
s^Y$ between the number of fundamental affinities $a$, that of broken
conservation laws $b$ and the number of chemostats $s^Y$. We decompose the
steady state entropy production rate in terms of fundamental fluxes and
affinities in the spirit of Schnakenberg's theory of network thermodynamics,
paving the way for the forthcoming treatment of the linear regime, of
efficiency and tight coupling, of free energy transduction and of thermodynamic
constraints for network reconstruction.


  We consider a three-level food chain in which an epidemics affects the
intermediate population. Two models are presented, respectively either allowing
for unlimited food supply for the bottom prey, or instead assuming for it a
logistic growth. Counterintuitive results related to the paradox of enrichment
are obtained, showing that by providing large amounts of food to the bottom
prey, the top predator and the disease in suitable situations can be
eradicated.


  Evolutionary game theory has become one of the most diverse and far reaching
theories in biology. Applications of this theory range from cell dynamics to
social evolution. However, many applications make it clear that inherent
non-linearities of natural systems need to be taken into account. One way of
introducing such non-linearities into evolutionary games is by the inclusion of
multiple players. An example is of social dilemmas, where group benefits could
e.g.\ increase less than linear with the number of cooperators. Such
multiplayer games can be introduced in all the fields where evolutionary game
theory is already well established. However, the inclusion of non-linearities
can help to advance the analysis of systems which are known to be complex, e.g.
in the case of non-Mendelian inheritance. We review the diachronic theory and
applications of multiplayer evolutionary games and present the current state of
the field. Our aim is a summary of the theoretical results from well-mixed
populations in infinite as well as finite populations. We also discuss examples
from three fields where the theory has been successfully applied, ecology,
social sciences and population genetics. In closing, we probe certain future
directions which can be explored using the complexity of multiplayer games
while preserving the promise of simplicity of evolutionary games.


  Brownian ratchet like stochastic theory for the electrochemical membrane
system of Hodgkin-Huxley (HH) is developed. The system is characterized by a
continuous variable $Q_m(t)$, representing mobile membrane charge density, and
a discrete variable $K_t$ representing ion channel conformational dynamics. A
Nernst-Planck-Nyquist-Johnson type equilibrium is obtained when multiple
conducting ions have a common reversal potential. Detailed balance yields a
previously unknown relation between the channel switching rates and membrane
capacitance, bypassing Eyring-type explicit treatment of gating charge
kinetics. From a molecular structural standpoint, membrane charge $Q_m$ is a
more natural dynamic variable than potential $V_m$; our formalism treats
$Q_m$-dependent conformational transition rates $\lambda_{ij}$ as intrinsic
parameters. Therefore in principle, $\lambda_{ij}$ vs. $V_m$ is experimental
protocol dependent,e.g., different from voltage or charge clamping
measurements. For constant membrane capacitance per unit area $C_m$ and
neglecting membrane potential induced by gating charges, $V_m=Q_m/C_m$, and
HH's formalism is recovered. The presence of two types of ions, with different
channels and reversal potentials, gives rise to a nonequilibrium steady state
with positive entropy production $e_p$. For rapidly fluctuating channels, an
expression for $e_p$ is obtained.


  Tissue cells are in a state of permanent mechanical tension that is
maintained mainly by myosin II minifilaments, which are bipolar assemblies of
tens of myosin II molecular motors contracting actin networks and bundles. Here
we introduce a stochastic model for myosin II minifilaments as two small myosin
II motor ensembles engaging in a stochastic tug-of-war. Each of the two
ensembles is described by the parallel cluster model that allows us to use
exact stochastic simulations and at the same time to keep important molecular
details of the myosin II cross-bridge cycle. Our simulation and analytical
results reveal a strong dependence of myosin II minifilament dynamics on
environmental stiffness that is reminiscent of the cellular response to
substrate stiffness. For small stiffness, minifilaments form transient
crosslinks exerting short spikes of force with negligible mean. For large
stiffness, minifilaments form near permanent crosslinks exerting a mean force
which hardly depends on environmental elasticity. This functional switch arises
because dissociation after the power stroke is suppressed by force (catch
bonding) and because ensembles can no longer perform the power stroke at large
forces. Symmetric myosin II minifilaments perform a random walk with an
effective diffusion constant which decreases with increasing ensemble size, as
demonstrated for rigid substrates with an analytical treatment.


  Molecular dynamics simulations of myoglobin and aspartate aminotransferase,
with explicit solvent, are shown to accurately reproduce the experimentally
measured molar volumes. Single amino-acid substitution at VAL39 of aspartate
aminotransferase is known to produce large volumetric changes in the enzyme,
and this effect is demonstrated in simulation as well. This molecular dynamics
approach, while more computationally expensive that extant computational
methods of determining the apparent volume of biological systems, is quite
feasible with modern computer hardware and is shown to yield accurate
volumetric data with as little as several nanoseconds of dynamics.


  An elastic rod model for semi-flexible polymers is presented. Theory for a
continuum rod is reviewed, and it is shown that a popular discretised model
used in numerical simulations gives the correct continuum limit. Correlation
functions relating to both bending and twisting of the rod are derived for both
continuous and discrete cases, and results are compared with numerical
simulations. Finally, two possible implementations of the discretised model in
the multi-purpose molecular dynamics software package LAMMPS are described.


  We prove that a Moran model converges in probability to Eigen's quasispecies
model in the infinite population limit.


  Starting from common assumptions, we build a rate equation model for
multi-strain disease dynamics in terms of immune repertoire classes. We then
move to a strain-level description where a low-order closure reminiscent of a
pair approximation can be applied. We characterize the endemic equilibrium of
the ensuing model in the absence of mutation and discuss the presence of
degeneracy regarding the prevalence of the different strains. Finally we study
the behavior of the system under the injection of mutant strains.


  The bending of cilia and flagella is driven by forces generated by dynein
motor proteins. These forces slide adjacent microtubule doublets within the
axoneme, the motile cytoskeletal structure. To create regular, oscilla- tory
beating patterns, the activities of the axonemal dyneins must be coordinated
both spatially and temporally. It is thought that coordination is mediated by
stresses or strains, which build up within the moving axoneme, and somehow
regulate dynein activity. While experimenting with axonemes subjected to mild
proteolysis, we observed pairs of doublets associate with each other and form
bends with almost constant curvature. By model- ing the statics of a pair of
filaments, we show that the activity of the motors concentrates at the distal
tips of the doublets. Furthermore, we show that this distribution of motor
activity accords with models in which curvature, or curvature-induced normal
forces, regulates the activity of the motors. These observations, together with
our theoretical analysis, provide evidence that dynein activity can be
regulated by curvature or normal forces, which may, therefore, play a role in
coordinating the beating of cilia and flagella.


  Can one hear the 'sound' of a growing network? We address the problem of
recognizing the topology of evolving biological or social networks. Starting
from percolation theory, we analytically prove a linear inverse relationship
between two simple graph parameters--the logarithm of the average cluster size
and logarithm of the ratio of the edges of the graph to the theoretically
maximum number of edges for that graph--that holds for all growing power law
graphs. The result establishes a novel property of evolving power-law networks
in the asymptotic limit of network size. Numerical simulations as well as
fitting to real-world citation co-authorship networks demonstrate that the
result holds for networks of finite sizes, and provides a convenient measure of
the extent to which an evolving family of networks belongs to the same
power-law class.


  One of the most popular models for quantitatively understanding the emergence
of drug resistance both in bacterial colonies and in malignant tumors was
introduced long ago by Luria and Delbr\"uck. Here, individual resistant mutants
emerge randomly during the birth events of an exponentially growing sensitive
population. A most interesting limit of this process occurs when the population
size $N$ is large and mutation rates are low, but not necessarily small
compared to $1/N$. Here we provide a scaling solution valid in this limit,
making contact with the theory of Levy $\alpha$-stable distributions, in
particular one discussed long ago by Landau. One consequence of this
association is that moments of the distribution are highly misleading as far as
characterizing typical behavior. A key insight that enables our solution is
that working in the fixed population size ensemble is not the same as working
in a fixed time ensemble. Some of our results have been presented previously in
shortened form.


  The average size of internal translated exons, ranging from 120 to 165 nt
across metazoans, is approximately the size of the typical mononucleosome (147
nt). Genome-wide study has also shown that nucleosome occupancy is
significantly higher in exons than in introns, which might indicate that the
evolution of exon size is related to its nucleosome occupancy. By grouping
exons by the GC contents of their flanking introns, we show that the average
exon size is positively correlated with its GC content. Using the sequencing
data from direct mapping of Homo sapiens nucleosomes with limited nuclease
digestion, we show that the level of nucleosome occupancy is also positively
correlated with the exon GC content in a similar fashion. We then demonstrated
that exon size is positively correlated with their nucleosome occupancy. The
strong correlation between exon size and the nucleosome occupancy suggests that
chromatin organization may be related to the evolution of exon sizes.


  All genes on the human Y-chromosome were studied using fractal dimension and
Shannon entropy. Clear outlier clusters were identified. Among these were 6
sequences that have since been withdrawn as CDSs and 1 additional sequence that
is not in the current assembly. A methodology for ranking the sequences based
on deviation from average values of FD and SE was developed. The group of
sequences scored among the 10% largest deviations had abnormally high
likelihood to be from centromeric or pseudoautosomal regions and low likelihood
to be from X-chromosome transposed regions. lncRNA sequences were also enriched
among the outliers. In addition, the number of expressed genes previously
identified for evolutionary study tended to not have large deviations from the
average.
  Keywords: Y-chromosome; Shannon di-nucleotide entropy; fractal dimension;
centromeric genes; gene degredation; lncRNA


  The fractal globule state is a popular model for describing chromatin packing
in eukaryotic nuclei. Here we provide a scaling theory and dissipative particle
dynamics (DPD) computer simulation for the thermal motion of monomers in the
fractal globule state. Simulations starting from different entanglement-free
initial states show good convergence which provides evidence supporting the
existence of unique metastable fractal globule state. We show monomer motion in
this state to be sub-diffusive described by $\langle X^2 (t)\rangle \sim
t^{\alpha_F}$ with $\alpha_F$ close to 0.4. This result is in good agreement
with existing experimental data on the chromatin dynamics which makes an
additional argument in support of the fractal globule model of chromatin
packing.


  This letter introduced a new R package 'coexist' which can perform species
coexistence simulation and analysis. The package was initially developed for
understanding the role of different combinations of varying species growth
rates, dispersal ability, habitat heterogeneity, and competition interaction on
influencing species coexistence under fluctuating environment. Further versions
could introduce more coexistence models under different ecological conditions.
The source package is available at the Google Code
(https://code.google.com/p/coexist/) and R CRAN repository
(http://cran.r-project.org/web/packages/coexist/).


  We present a model for the uncoating of the HIV virus driven by forces
exerted on the protein shell of HIV generated by DNA confinement.


  Inferring the ancestral state at the root of a phylogenetic tree from states
observed at the leaves is a problem arising in evolutionary biology. The
simplest technique -- majority rule -- estimates the root state by the most
frequently occurring state at the leaves. Alternative methods -- such as
maximum parsimony - explicitly take the tree structure into account. Since
either method can outperform the other on particular trees, it is useful to
consider the accuracy of the methods on trees generated under some evolutionary
null model, such as a Yule pure-birth model. In this short note, we answer a
recently posed question concerning the performance of majority rule on Yule
trees under a symmetric 2-state Markovian substitution model of character state
change. We show that majority rule is accurate precisely when the ratio of the
birth (speciation) rate of the Yule process to the substitution rate exceeds
the value $4$. By contrast, maximum parsimony has been shown to be accurate
only when this ratio is at least 6. Our proof relies on a second moment
calculation, coupling, and a novel application of a reflection principle.


  Predator-prey relationships are one of the most studied interactions in
population ecology. However, little attention has been paid to the possibility
of role exchange between species once determined as predators and preys,
despite firm field evidence of such phenomena in the nature. In this paper, we
build a model capable of reproducing the main phenomenological features of one
reported predator-prey role-reversal system, and present results for both the
homogeneous and the space explicit cases. We find that, depending on the choice
of parameters, our role-reversal dynamical system exhibits excitable-like
behaviour, generating waves of species' concentrations that propagate through
space.


  Biological information is not only stored in the digital chemical sequence of
double helical DNA, but is also encoded in the mechanical properties of the DNA
strands, which can influence biochemical processes involving its readout. For
example, loop formation in the Lac operon can regulate the expression of key
genes, and DNA supercoiling is closely correlated to rhythmic circardian gene
expression in cyanobacteria. Supercoiling is also important for large scale
organisation of the genome in both eukaryotic and prokaryotic cells. DNA can
respond to torsional stress by writhing to form looped structures called
plectonemes, thus transferring energy stored as twist into energy stored in
bending. Denaturation bubbles can also relax torsional stress, with the
enthalpic cost of breaking bonds being compensated by their ability to absorb
undertwist. Here we predict a novel regime where bubbles form at the tips of
plectonemes, and study its properties using coarse-grained simulations. These
tip bubbles can occur for both positive and negative supercoiling and greatly
reduce plectoneme diffusion by a pinning mechanism. They can cause plectonemes
to preferentially localise to AT rich regions, because bubbles more easily form
there. The tip-bubble regime occurs for supercoiling densities and forces that
are typically encountered for DNA in vivo, and may be exploited for biological
control of genomic processes.


  Direct reciprocity and conditional cooperation are important mechanisms to
prevent free riding in social dilemmas. But in large groups these mechanisms
may become ineffective, because they require single individuals to have a
substantial influence on their peers. However, the recent discovery of the
powerful class of zero-determinant strategies in the iterated prisoner's
dilemma suggests that we may have underestimated the degree of control that a
single player can exert. Here, we develop a theory for zero-determinant
strategies for multiplayer social dilemmas, with any number of involved
players. We distinguish several particularly interesting subclasses of
strategies: fair strategies ensure that the own payoff matches the average
payoff of the group; extortionate strategies allow a player to perform above
average; and generous strategies let a player perform below average. We use
this theory to explore how individuals can enhance their strategic options by
forming alliances. The effects of an alliance depend on the size of the
alliance, the type of the social dilemma, and on the strategy of the allies:
fair alliances reduce the inequality within their group; extortionate alliances
outperform the remaining group members; but generous alliances increase
welfare. Our results highlight the critical interplay of individual control and
alliance formation to succeed in large groups.


  The universality of many pathways of core metabolism suggests a strong role
for evolutionary selection, but it remains unclear whether existing pathways
have been selected from a large or small set of biochemical possibilities. To
address this question, we construct "in silico" all possible biochemically
feasible alternatives to the trunk pathway of glycolysis and gluconeogenesis,
one of the most highly conserved pathways in metabolism. We show that, even
though a large number of alternative pathways exist, the alternatives carry
lower flux than the real pathway under typical physiological conditions.
Alternative pathways that could potentially carry higher flux often lead to
infeasible intermediate metabolite concentrations. We also find that if
physiological conditions were different, different pathways could outperform
those found in nature. Our results demonstrate how the rules of biochemistry
restrict the alternatives that are open to evolution, and suggest that the
existing trunk pathway of glycolysis and gluconeogenesis represents a maximal
flux solution.


  Coevolution between strategy and network structure is established as a means
to arrive at optimal conditions for resolving social dilemmas. Yet recent
research highlights that the interdependence between networks may be just as
important as the structure of an individual network. We therefore introduce
coevolution of strategy and network interdependence to study whether it can
give rise to elevated levels of cooperation in the prisoner's dilemma game. We
show that the interdependence between networks self-organizes so as to yield
optimal conditions for the evolution of cooperation. Even under extremely
adverse conditions cooperators can prevail where on isolated networks they
would perish. This is due to the spontaneous emergence of a two-class society,
with only the upper class being allowed to control and take advantage of the
interdependence. Spatial patterns reveal that cooperators, once arriving to the
upper class, are much more competent than defectors in sustaining compact
clusters of followers. Indeed, the asymmetric exploitation of interdependence
confers to them a strong evolutionary advantage that may resolve even the
toughest of social dilemmas.