  According to the `ceRNA hypothesis', microRNAs (miRNAs) may act as mediators
of an effective positive interaction between long coding or non-coding RNA
molecules, carrying significant potential implications for a variety of
biological processes. Here, inspired by recent work providing a quantitative
description of small regulatory elements as information-conveying channels, we
characterize the effectiveness of miRNA-mediated regulation in terms of the
optimal information flow achievable between modulator (transcription factors)
and target nodes (long RNAs). Our findings show that, while a sufficiently
large degree of target derepression is needed to activate miRNA-mediated
transmission, (a) in case of differential mechanisms of complex processing
and/or transcriptional capabilities, regulation by a post-transcriptional
miRNA-channel can outperform that achieved through direct transcriptional
control; moreover, (b) in the presence of large populations of weakly
interacting miRNA molecules the extra noise coming from titration disappears,
allowing the miRNA-channel to process information as effectively as the direct
channel. These observations establish the limits of miRNA-mediated
post-transcriptional cross-talk and suggest that, besides providing a degree of
noise buffering, this type of control may be effectively employed in cells both
as a failsafe mechanism and as a preferential fine tuner of gene expression,
pointing to the specific situations in which each of these functionalities is
maximized.


  In this paper, a network of hypertension-related genes is constructed by
analyzing the correlations of gene expression data among the Dahl
salt-sensitive rat and two consomic rat strains. The numerical calculations
show that this sparse and assortative network has small-world and scale-free
properties. Further, 16 key hub genes (Col4a1, Lcn2, Cdk4, etc.) are determined
by introducing an integrated centrality and have been confirmed by
biological/medical research to play important roles in hypertension.


  Although dengue fever is an ancient disease that should have been eradicated
since time ago, in Mexico as well as in other regions worldwide the cases of
classic and hemorrhagic dengue continue to be reported presently. Since 2002
(and before), a worrying outbreak of dengue was observed in several states of
Mexico, notably the case of Colima for which a "descacharrizacion" campaign was
implemented in 2009 to encourage the population to get rid of old objects
accumulated in their backyards and outdoors so as to prevent proliferation of
mosquitoes Aedes Aegypti. To understand the effect of such campaign, we studied
the incidence of dengue in four municipalities of the mentioned state, namely
Manzanillo, Armeria, Tecoman and Ixtlahuacan. In particular, we observed that
the incidence of dengue in these four municipalities had a strong relationship
with climatological normals, specifically precipitation and minimum
temperature. Furthermore, the results showed that the mentioned campaign helped
prevent mosquito eggs incubation and their development, thus leading to
reduction of the incidence of dengue in the human population.


  Data describing the historical growth of human population global and regional
(Western Europe, Eastern Europe, Asia, former USSR, Africa and Latin America)
are analysed. Results are in harmony with the earlier analysis of the
historical growth of the world population in the past 12,000 years and with a
similar but limited study carried out over 50 years ago. This analysis is also
in harmony with the study of the historical economic growth. Within the range
of analysable data, there was no Malthusian stagnation. Takeoffs from
stagnation to growth, postulated by the Unified Growth Theory never happened.
There were no escapes from the Malthusian trap because there was no trap. This
analysis and the earlier studies of the Gross Domestic Product lead to the
conclusion that there were also no takeoffs in the income per capita
distributions, claimed by the Unified Growth Theory. Consequently, the claimed
in this theory differential timing in takeoffs never happened. Unified Growth
Theory is contradicted yet again by the mathematical analysis of the same data,
which were used, but never analysed, during the formulation of this theory.
However, this study, as well as the earlier publications on the related topics,
shows also that some fundamental postulates used in the economic and
demographic research are repeatedly contradicted by the mathematical analysis
of data.


  Recent years have seen a large amount of interest in epidemics on networks as
a way of representing the complex structure of contacts capable of spreading
infections through the modern human population. The configuration model is a
popular choice in theoretical studies since it combines the ability to specify
the distribution of the number of contacts (degree) with analytical
tractability. Here we consider the early real-time behaviour of the Markovian
SIR epidemic model on a configuration model network using a multi-type
branching process. We find closed-form analytic expressions for the mean and
variance of the number of infectious individuals as a function of time and the
degree of the initially infected individual(s), and write down a system of
differential equations for the probability of extinction that are numerically
fast compared to Monte Carlo simulation. We show that these quantities are all
sensitive to the degree distribution - in particular we confirm that the mean
prevalence of infection depends on the first two moments of the degree
distribution and the variance in prevalence depends on the first three moments
of the degree distribution. In contrast to most existing analytic approaches,
the accuracy of these results does not depend on having a large number of
infectious individuals, meaning that in the large population limit they would
be asymptotically exact even for one initial infectious individual.


  Recent advances have allowed for both morphological fossil evidence and
molecular sequences to be integrated into a single combined inference of
divergence dates under the rule of Bayesian probability. In particular the
fossilized birth-death tree prior and the Lewis-Mk model of discrete
morphological evolution allow for the estimation of both divergence times and
phylogenetic relationships between fossil and extant taxa. We exploit this
statistical framework to investigate the internal consistency of these models
by producing phylogenetic estimates of the age of each fossil in turn, within
two rich and well-characterized data sets of fossil and extant species
(penguins and canids). We find that the estimation accuracy of fossil ages is
generally high with credible intervals seldom excluding the true age and median
relative error in the two data sets of 5.7% and 13.2% respectively. The median
relative standard error (RSD) was 9.2% and 7.2% respectively, suggesting good
precision, although with some outliers. In fact in the two data sets we analyze
the phylogenetic estimates of fossil age is on average < 2 My from the midpoint
age of the geological strata from which it was excavated. The high level of
internal consistency found in our analyses suggests that the Bayesian
statistical model employed is an adequate fit for both the geological and
morphological data, and provides evidence from real data that the framework
used can accurately model the evolution of discrete morphological traits coded
from fossil and extant taxa. We anticipate that this approach will have diverse
applications beyond divergence time dating, including dating fossils that are
temporally unconstrained, testing of the "morphological clock", and for
uncovering potential model misspecification and/or data errors when
controversial phylogenetic hypotheses are obtained based on combined divergence
dating analyses.


  We use a multivariate formulation of sequential Monte Carlo filter that
utilizes mechanistic models for Ebola virus propagation and available incidence
data to simultaneously estimate the disease progression states and the model
parameters. This method has the advantage of performing the inference online as
the new data becomes available and estimates the evolution of basic
reproductive ratio $R_0(t)$ of the Ebola outbreak through time. Our analysis
identifies a peak in the basic reproductive ratio close to the time when Ebola
cases were reported in Europe and the USA.


  Complex network approach has become an effective way to describe
interrelationships among large amounts of biological data, which is especially
useful in finding core functions and global behavior of biological systems.
Hypertension is a complex disease caused by many reasons including genetic,
physiological, psychological and even social factors. In this paper, based on
the information of biological pathways, we construct a network model of
hypertension-related genes of the salt-sensitive rat to explore the
interrelationship between genes. Statistical and topological characteristics
show that the network has the small-world but not scale-free property, and
exhibits a modular structure, revealing compact and complex connections among
these genes. By the threshold of integrated centrality larger than 0.71, seven
key hub genes are found: Jun, Rps6kb1, Cycs, Creb3l2, Cdk4, Actg1 and RT1-Da.
These genes should play an important role in hypertension, suggesting that the
treatment of hypertension should focus on the combination of drugs on multiple
genes.


  From footpaths to flight routes, human mobility networks facilitate the
spread of communicable diseases. Control and elimination efforts depend on
characterizing these networks in terms of connections and flux rates of
individuals between contact nodes. In some cases, transport can be
parameterized with gravity-type models or approximated by a diffusive random
walk. As a alternative, we have isolated intranational commercial air traffic
as a case study for the utility of non-diffusive, heavy-tailed transport
models. We implemented new stochastic simulations of a prototypical
influenza-like infection, focusing on the dense, highly-connected United States
air travel network. We show that mobility on this network can be described
mainly by a power law, in agreement with previous studies. Remarkably, we find
that the global evolution of an outbreak on this network is accurately
reproduced by a two-parameter space-fractional diffusion equation, such that
those parameters are determined by the air travel network.


  We propose a faster algorithm for individual based simulations for adaptive
dynamics based on a simple modification to the standard Gillespie Algorithm for
simulating stochastic birth-death processes. We provide an analytical
explanation that shows that simulations based on the modified algorithm, in the
deterministic limit, lead to the same equations of adaptive dynamics as well as
same conditions for evolutionary branching as those obtained from the standard
Gillespie algorithm. Based on this algorithm, we provide an intuitive and
simple interpretation of the canonical equation of adaptive dynamics. With the
help of examples we compare the performance of this algorithm to the standard
Gillespie algorithm and demonstrate its efficiency. We also study an example
using this algorithm to study evolutionary dynamics in a multi-dimensional
phenotypic space and study the question of predictability of evolution.


  Lateral transfer, a process whereby species exchange evolutionary traits
through non-ancestral relationships, is a frequent source of model
misspecification in phylogenetic inference. Lateral transfer obscures the
phylogenetic signal in the data as the histories of affected traits are mosaics
of the overall phylogeny. We control for the effect of lateral transfer in a
Stochastic Dollo model and a Bayesian setting. Our likelihood is highly
intractable as the parameters are the solution of a sequence of large systems
of differential equations representing the expected evolution of traits along a
tree. We illustrate our method on a data set of lexical traits in Eastern
Polynesian languages and obtain an improved fit over the corresponding model
without lateral transfer.


  Direct reciprocity is a mechanism for the evolution of cooperation based on
repeated interactions. When individuals meet repeatedly, they can use
conditional strategies to enforce cooperative outcomes that would not be
feasible in one-shot social dilemmas. Direct reciprocity requires that
individuals keep track of their past interactions and find the right response.
However, there are natural bounds on strategic complexity: Humans find it
difficult to remember past interactions accurately, especially over long
timespans. Given these limitations, it is natural to ask how complex strategies
need to be for cooperation to evolve. Here, we study stochastic evolutionary
game dynamics in finite populations to systematically compare the evolutionary
performance of reactive strategies, which only respond to the co-player's
previous move, and memory-one strategies, which take into account the own and
the co-player's previous move. In both cases, we compare deterministic strategy
and stochastic strategy spaces. For reactive strategies and small costs, we
find that stochasticity benefits cooperation, because it allows for
generous-tit-for-tat. For memory one strategies and small costs, we find that
stochasticity does not increase the propensity for cooperation, because the
deterministic rule of win-stay, lose-shift works best. For memory one
strategies and large costs, however, stochasticity can augment cooperation.


  F1F0 ATP synthase (ATPase) either facilitates the synthesis of ATP in the
mitochondrial membranes and bacterial inner membranes in a process driven by
the proton moving force (pmf), or uses the energy from ATP hydrolysis to pump
protons against the concentration gradient across the membrane. ATPase is
composed of two rotary motors, F0 and F1, which generate the opposing rotation
and compete for control of their shared central gamma-shaft. Here we present a
self-consistent physical model of the F1 motor as a simplified two-state
Brownian ratchet based on the asymmetry of torsional elastic energy of the
coiled-coil gamma-shaft. This stochastic model unifies the physical description
of linear and rotary motors and explains the stepped unidirectional rotation of
the $\gamma$-shaft, in agreement with the `binding-change' ideas of Boyer.
Substituting the model parameters, all independently known from recent
experiments, our model quantitatively reproduces the ATPase operation, e.g. the
`no-load' angular velocity is ca. 400~rad/s anticlockwise at 4 mM ATP, in close
agreement with experiment. Increasing the pmf torque exerted by F0 can slow,
stop and overcome the torque generated by F1, switching from ATP hydrolysis to
synthesis at a very low value of `stall torque'. We discuss the matters of the
motor efficiency, which is very low if calculated from the useful mechanical
work it produces - but is quite high when the `useful outcome' is measured in
the number of H+ pushed against the chemical gradient in the F1 ATP-driven
operation.


  Recent studies show differences in individual motion and shoaling tendency
between strains of the same species. Here, we analyse the collective motion and
the response to visual stimuli in two morphologically different strains (TL and
AB) of zebrafish. For both strains, we observe 10 groups of 5 and 10 zebrafish
swimming freely in a large experimental tank with two identical attractive
landmarks (cylinders or disks) for one hour. We track the positions of the fish
by an automated tracking method and compute several metrics at the group level.
First, the probability of presence shows that both strains avoid free space and
are more likely to swim in the vicinity of the walls of the tank and the
attractive landmarks. Second, the analysis of landmarks occupancy shows that AB
zebrafish are more present in their vicinity than TL zebrafish and that both
strains regularly transit from one landmark to the other with no preference on
the long duration. Finally, TL zebrafish show a higher cohesion than AB
zebrafish. Thus, landmarks and duration of the repicates allow to reveal
collective behavioural variabilities among different strains of zebrafish.
These results provide a new insight into the need to take into account
individual variability of zebrafish strains for studying collective behaviour.


  Synthetic lethality refers to a combination of two or more genetic events
(typically affecting different genes) in which the co-occurrence of the events
results in cell or organismal lethality, but the cell or organism remains
viable when only one of the events occurs. Synthetic lethality has gained
attention in the last few years for its value in selective killing of cancer
cells: by targeting the synthetic lethal partner of an altered gene in cancer,
only the cancer cells can be killed while sparing normal cells. In a recent
study, we showed that mutual exclusive combinations of genetic events in cancer
hint at naturally occurring synthetic lethal combinations, and therefore by
systematically mining for these combinations we can identify novel therapeutic
targets for cancer. Based on this, we had identified a list of 718 genes that
are mutually exclusive to six DNA-damage response genes in cancer. Here, we
extend these results to identify a subset of 43 genes whose over-expression
correlates with significantly poor survival in estrogen receptor-negative
breast cancers, and thus provide a promising list of potential therapeutic
targets and/or biomarkers.


  The mutation rate of a well adapted population is prone to reduction so as to
have a lower mutational load. We aim to understand the role of epistatic
interactions between the fitness affecting mutations in this process. Using a
multitype branching process, the fixation probability of a single non-mutator
emerging in a large asexual mutator population is analytically calculated here.
The mutator population undergoes deleterious mutations at constant, but at a
much higher rate than that of the non-mutator. We find that antagonistic
epistasis lowers the chances of mutation rate reduction, while synergistic
epistasis enhances it. Below a critical value of epistasis, the fixation
probability behaves non-monotonically with variation in mutation rate of the
background population. Moreover, the variation of this critical value of the
epistasis parameter with the strength of the mutator is discussed in the
Appendix. For synergistic epistasis, when selection is varied, the fixation
probability reduces overall, with damped oscillations.


  Reducing the fuel load in fire-prone landscapes is aimed at mitigating the
risk of catastrophic wildfires but there are ecological consequences.
Maintaining habitat for fauna of both sufficient extent and connectivity while
fragmenting areas of high fuel loads presents land managers with seemingly
contrasting objectives. Faced with this dichotomy, we propose a Mixed Integer
Programming (MIP) model that can optimally schedule fuel treatments to reduce
fuel hazards by fragmenting high fuel load regions while considering critical
ecological requirements over time and space. The model takes into account both
the frequency of fire that vegetation can tolerate and the frequency of fire
necessary for fire-dependent species. Our approach also ensures that suitable
alternate habitat is available and accessible to fauna affected by a treated
area. More importantly, to conserve fauna the model sets a minimum acceptable
target for the connectivity of habitat at any time . These factors are all
included in the formulation of a model that yields a multi-period
spatially-explicit schedule for treatment planning. Our approach is then
demonstrated in a series of computational experiments with hypothetical
landscapes, a single vegetation type and a group of faunal species with the
same habitat requirements. Our experiments show that it is possible to reduce
the risk of wildfires while ensuring sufficient connectivity of habitat over
both space and time. Furthermore, it is demonstrated that the habitat
connectivity constraint is more effective than neighbourhood habitat
constraints. This is critical for the conservation of fauna and of special
concern for vulnerable or endangered species.


  Based on the glocal HIM metric and its induced graph kernel, we propose a
novel solution in differential network analysis that integrates network
comparison and classification tasks. The HIM distance is defined as the
one-parameter family of product metrics linearly combining the normalised
Hamming distance H and the normalised Ipsen-Mikhailov spectral distance IM. The
combination of the two components within a single metric allows overcoming
their drawbacks and obtaining a measure that is simultaneously global and
local. Furthermore, plugging the HIM kernel into a Support Vector Machine gives
us a classification algorithm based on the HIM distance. First, we outline the
theory underlying the metric construction. We introduce two diverse
applications of the HIM distance and the HIM kernel to biological datasets.
This versatility supports the adoption of the HIM family as a general tool for
information extraction, quantifying difference among diverse in- stances of a
complex system. An Open Source implementation of the HIM metrics is provided by
the R package nettols and in its web interface ReNette.


  This work is a systematic study of discrete Markov chains that are used to
describe the evolution of a two-types population. Motivated by results valid
for the well-known Moran (M) and Wright-Fisher (WF) processes, we define a
general class of Markov chains models which we term the Kimura class. It
comprises the majority of the models used in population genetics, and we show
that many well-known results valid for M and WF processes are still valid in
this class. In all Kimura processes, a mutant gene will either fixate or become
extinct, and we present a necessary and sufficient condition for such processes
to have the probability of fixation strictly increasing in the initial
frequency of mutants. This condition implies that there are WF processes with
decreasing fixation probability --- in contradistinction to M processes which
always have strictly increasing fixation probability. As a by-product, we show
that an increasing fixation probability defines uniquely an M or WF process
which realises it, and that any fixation probability with no state having
trivial fixation can be realised by at least some WF process. These results are
extended to a subclass of processes that are suitable for describing
time-inhomogeneous dynamics. We also discuss the traditional identification of
frequency dependent fitnesses and pay-offs, extensively used in evolutionary
game theory, the role of weak selection when the population is finite, and the
relations between jumps in evolutionary processes and frequency dependent
fitnesses.


  We consider random search processes alternating stochastically between
diffusion and ballistic motion, in which the distribution function of ballistic
motion directions varies from point to point in space. The specific space
dependence of the directional distribution together with the switching rates
between the two modes of motion establishes a spatially inhomogeneous search
strategy. We show that the mean first passage times for several standard search
problems - narrow escape, reaction partner finding, reaction-escape - can be
minimized with a directional distribution that is reminiscent of the spatial
organization of the cytoskeleton filaments of cells with a centrosome: radial
ballistic transport from center to periphery and back, and ballistic transport
in random directions within a concentric shell of thickness $\Delta_{\rm opt}$
along the domain boundary. The results suggest that living cells realize
efficient search strategies for various intracellular transport problems
economically through a spatial cytoskeleton organization that involves radial
microtubules in the central region and only a narrow actin cortex rather than a
cell body filled with randomly oriented actin filaments.


  This work presents a comparative analysis of Influenzanet data for influenza
itself and common cold in the Netherlands during the last 5 years, from the
point of view of modelling by linearised SIRS equations parametrically driven
by the ambient temperature. It is argued that this approach allows for the
forecast of common cold, but not of influenza in a strict sense. The difference
in their kinetic models is discussed with reference to the clinical background.


  Neonatal sepsis (NS) is responsible for over a 1 million yearly deaths
worldwide. In the developing world NS is often treated without an identified
microbial pathogen. Amplicon sequencing of the bacterial 16S rRNA gene can be
used to identify organisms that are difficult to detect by routine
microbiological methods. However, contaminating bacteria are ubiquitous in both
hospital settings and research reagents, and must be accounted for to make
effective use of these data. In the present study, we sequenced the bacterial
16S rRNA gene obtained from blood and cerebrospinal fluid (CSF) of 80 neonates
presenting with NS to the Mbarara Regional Hospital in Uganda. Assuming that
patterns of background contamination would be independent of pathogenic
microorganism DNA, we applied a novel quantitative approach using principal
orthogonal decomposition to separate background contamination from potential
pathogens in sequencing data. We designed our quantitative approach contrasting
blood, CSF, and control specimens, and employed a variety of statistical random
matrix bootstrap hypotheses to estimate statistical significance. These
analyses demonstrate that Leptospira appears present in some infants presenting
within 48 hr of birth, indicative of infection in utero, and up to 28 days of
age, suggesting environmental exposure. This organism cannot be cultured in
routine bacteriological settings, and is enzootic in the cattle that the rural
peoples of western Uganda often live in close proximity. Our findings
demonstrate that statistical approaches to remove background organisms common
in 16S sequence data can reveal putative pathogens in small volume biological
samples from newborns. This computational analysis thus reveals an important
medical finding that has the potential to alter therapy and prevention efforts
in a critically ill population.


  Residential mobility is deeply entangled with all aspects of hunter-gatherer
life ways, and is therefore an issue of central importance in hunter-gatherer
studies. Hunter-gatherers vary widely in annual rates of residential mobility,
and understanding the sources of this variation has long been of interest to
anthropologists and archaeologists. Since mobility is, to a large extent,
driven by the need for a continuous supply of food, a natural framework for
addressing this question is provided by the metabolic theory of ecology. This
provides a powerful framework for formulating formal testable hypotheses
concerning evolutionary and ecological constraints on the scale and variation
of hunter-gatherer residential mobility. We evaluate these predictions using
extant data and show strong support for the hypotheses. We show that the
overall scale of hunter-gatherer residential mobility is predicted by average
human body size, and the limited capacity of mobile hunter-gatherers to store
energy internally. We then show that the majority of variation in residential
mobility observed across cultures is predicted by energy availability in local
ecosystems. Our results demonstrate that large-scale evolutionary and
ecological processes, common to all plants and animals, constrain
hunter-gatherers in predictable ways as they move through territories to
effectively exploit resources over the course of a year. Moreover, our results
extend the scope of the metabolic theory of ecology by showing how it
successfully predicts variation in the behavioral ecology of populations within
a species.


  The amplification of the gene MYCN (V-myc myelocytomatosis viral-valeted
oncogene, neuroblastoma derived) has been a well-documented indicator for poor
prognosis in neuroblastoma, a childhood cancer. Unfortunately, there has been
limited success in understanding MYCN functionality in the landscape of
neuroblastoma and more importantly, given that MYCN has been deemed
undruggable, the need to potentially illuminate key opportunities that
indirectly target MYCN is of great interest. To this end, this work employs an
emerging quantitative technique from network science, namely network curvature,
to quantify the biological robustness of MYCN and its surrounding neighborhood.
In particular, when amplified in Stage IV cancer, MYCN exhibits higher
curvature (more robust) than those samples with under expressed MYCN levels.
When examining the surrounding neighborhood, the above argument still holds for
network curvature, but is lost when only analyzing differential expression - a
common technique amongst oncologists and computational/molecular biologists.
This finding points to the problem (and possible solution) of drug targeting in
the context of complexity and indirect cell signaling affects that have often
been obfuscated through traditional techniques.


  The alignment of protein-protein interaction (PPI) networks has many
applications, such as the detection of conserved biological network motifs, the
prediction of protein interactions, and the reconstruction of phylogenetic
trees [1, 2, 3]. IsoRank is one of the first global network alignment
algorithms [4, 5, 6], where the goal is to match all (or most) of the nodes of
two PPI networks. The IsoRank algorithm first computes a pairwise node
similarity metric, and then generates a matching between the two node sets
based on this metric. The metric is a convex combination of a structural
similarity score (with weight $ \alpha $) and an extraneous amino-acid sequence
similarity score for two proteins (with weight $ 1 - \alpha $). In this short
paper, we make two contributions. First, we show that when IsoRank similarity
depends only on network structure ($\alpha = 1$), the similarity of two nodes
is only a function of their degrees. In other words, IsoRank similarity is
invariant to any network rewiring that does not affect the node degrees. This
result suggests a reason for the poor performance of IsoRank in structure-only
($ \alpha = 1 $) alignment. Second, using ideas from [7, 8], we develop an
approximation algorithm that outperforms IsoRank (including recent versions
with better scaling, e.g., [9]) by several orders of magnitude in time and
memory complexity, despite only a negligible loss in precision.


  Biological aging is characterized by an age-dependent increase in the
probability of death and by a decrease in the reproductive capacity. Individual
age-dependent rates of survival and reproduction have a strong impact on
population dynamics, and the genetic elements determining survival and
reproduction are under different selective forces throughout an organism
lifespan. Here we develop a highly versatile numerical model of genome
evolution --- both asexual and sexual --- for a population of virtual
individuals with overlapping generations, where the genetic elements affecting
survival and reproduction rate at different life stages are free to evolve due
to mutation and selection. Our model recapitulates several emerging properties
of natural populations, developing longer reproductive lifespan under stable
conditions and shorter survival and reproduction in unstable environments.
Faster aging results as the consequence of the reduced strength of purifying
selection in more unstable populations, which have large portions of the genome
that accumulate detrimental mutations. Unlike sexually reproducing populations
under constant resources, asexually reproducing populations fail to develop an
age-dependent increase in death rates and decrease in reproduction rates,
therefore escaping senescence. Our model provides a powerful in silico
framework to simulate how populations and genomes change in the context of
biological aging and opens a novel analytical opportunity to characterize how
real populations evolve their specific aging dynamics.


  A wide range of organisms features molecular machines, circadian clocks,
which generate endogenous oscillations with ~24 h periodicity and thereby
synchronize biological processes to diurnal environmental fluctuations.
Recently, it has become clear that plants harbor more complex gene regulatory
circuits within the core circadian clocks than other organisms, inspiring a
fundamental question: are all these regulatory interactions between clock genes
equally crucial for the establishment and maintenance of circadian rhythms? Our
mechanistic simulation for Arabidopsis thaliana demonstrates that at least half
of the total regulatory interactions must be present to express the circadian
molecular profiles observed in wild-type plants. A set of those essential
interactions is called herein a kernel of the circadian system. The kernel
structure unbiasedly reveals four interlocked negative feedback loops
contributing to circadian rhythms, and three feedback loops among them drive
the autonomous oscillation itself. Strikingly, the kernel structure, as well as
the whole clock circuitry, is overwhelmingly composed of inhibitory, rather
than activating, interactions between genes. We found that this tendency
underlies plant circadian molecular profiles which often exhibit
sharply-shaped, cuspidate waveforms. Through the generation of these cuspidate
profiles, inhibitory interactions may facilitate the global coordination of
temporally-distant clock events that are markedly peaked at very specific times
of day. Our systematic approach resulting in experimentally-testable
predictions provides insights into a design principle of biological clockwork,
with implications for synthetic biology.


  Understanding protein self-assembly is important for many biological and
industrial processes. Proteins can self-assemble into crystals, filaments,
gels, and other amorphous aggregates. The final forms include virus capsids and
condensed phases associated with diseases, such as amyloid fibrils. Although
seemingly different, these assemblies all originate from fundamental protein
interactions and are driven by similar thermodynamic and kinetic factors. Here
we review recent advances in understanding protein self-assembly through a soft
condensed matter perspective with an emphasis on three specific systems:
globular proteins, viruses and amyloid fibers. We conclude with a discussion of
unanswered questions in the field.


  An example of a nonspecific molecular bond is the affinity of any positive
charge for any negative charge (like-unlike), or of nonpolar material for
itself when in aqueous solution (like-like). This contrasts specific bonds such
as the affinity of the DNA base A for T, but not for C, G, or another A. Recent
experimental breakthroughs in DNA nanotechnology demonstrate that a particular
nonspecific like-like bond ("blunt-end DNA stacking" that occurs between the
ends of any pair of DNA double-helices) can be used to create specific
"macrobonds" by careful geometric arrangement of many nonspecific blunt ends,
motivating the need for sets of macrobonds that are orthogonal: two macrobonds
not intended to bind should have relatively low binding strength, even when
misaligned.
  To address this need, we introduce geometric orthogonal codes that abstractly
model the engineered DNA macrobonds as two-dimensional binary codewords. While
motivated by completely different applications, geometric orthogonal codes
share similar features to the optical orthogonal codes studied by Chung,
Salehi, and Wei. The main technical difference is the importance of 2D geometry
in defining codeword orthogonality.


  Since the sequencing of large genomes, many statistical features of their
sequences have been found. One intriguing feature is that certain subsequences
are much more abundant than others. In fact, abundances of subsequences of a
given length are distributed with a scale-free power-law tail, resembling
properties of human texts, such as the Zipf's law. Despite recent efforts, the
understanding of this phenomenon is still lacking. Here we find that selfish
DNA elements, such as those belonging to the Alu family of repeats, dominate
the power-law tail. Interestingly, for the Alu elements the power-law exponent
increases with the length of the considered subsequences. Motivated by these
observations, we develop a model of selfish DNA expansion. The predictions of
this model qualitatively and quantitatively agree with the empirical
observations. This allows us to estimate parameters for the process of selfish
DNA spreading in a genome during its evolution. The obtained results shed light
on how evolution of selfish DNA elements shapes non-trivial statistical
properties of genomes.


  Statistical inference about the speciation process has often been based on
the isolation-with-migration (IM) model, especially when the research aim is to
learn about the presence or absence of gene flow during divergence. The
generalised IM model introduced in this paper extends both the standard
two-population IM model and the isolation-with-initial-migration (IIM) model,
and encompasses both these models as special cases. It can be described as a
two-population IM model in which migration rates and population sizes are
allowed to change at some point in the past. By developing a maximum-likelihood
implementation of this GIM model, we enable inference on both historical and
contemporary rates of gene flow between two closely related species. Our method
relies on the spectral decomposition of the coalescent generator matrix and is
applicable to data sets consisting of the numbers of nucleotide differences
between one pair of DNA sequences at each of a large number of independent
loci.


  We prove that, for Poisson transmission and recovery processes, the classic
Susceptible $\to$ Infected $\to$ Recovered (SIR) epidemic model of Kermack and
McKendrick provides, for any given time $t>0$, a strict lower bound on the
expected number of suscpetibles and a strict upper bound on the expected number
of recoveries in the general stochastic SIR epidemic. The proof is based on the
recent message passing representation of SIR epidemics applied to a complete
graph.


  The Gene Ontology (GO) is a formidable resource but there are several
considerations about it that are essential to understand the data and interpret
it correctly. The GO is sufficiently simple that it can be used without deep
understanding of its structure or how it is developed, which is both a strength
and a weakness. In this chapter, we discuss some common misinterpretations of
the ontology and the annotations. A better understanding of the pitfalls and
the biases in the GO should help users make the most of this very rich
resource. We also review some of the misconceptions and misleading assumptions
commonly made about GO, including the effect of data incompleteness, the
importance of annotation qualifiers, and the transitivity or lack thereof
associated with different ontology relations. We also discuss several biases
that can confound aggregate analyses such as gene enrichment analyses. For each
of these pitfalls and biases, we suggest remedies and best practices.


  The Gene Ontology (GO) project is the largest resource for cataloguing gene
function. The combination of solid conceptual underpinnings and a practical set
of features have made the GO a widely adopted resource in the research
community and an essential resource for data analysis. In this chapter, we
provide a concise primer for all users of the GO. We briefly introduce the
structure of the ontology and explain how to interpret annotations associated
with the GO.


  Coronary stents have changed the way in which coronary artery diseases are
managed. Although bare-metal stents can be traced back to 1994, long term
efficacy of these stents has been shattered by ISR (In-Stent restenosis) and
late stent thrombosis. Research on finding a solution to these issues has led
to the development of DESs (Drug eluting stents). However long term
effectiveness of DESs and various drug delivery systems have raised concerns.
Also current DESs does not have the capability to adjust the drug dose or
release kinetics that corresponds to the diseased status of the affected
vessel. Through the use of a drug delivery system which employs controlled drug
release, it may be possible to control the release rate of the pharmacological
drugs or substances of interest, and therefore one might be able to circumvent
the range of events that leads to ISR, thus preventing the need for further
invasive interventions. To overcome the limitations of drug delivery systems,
the administration of drugs using nanoparticle based microcapsules as carriers
is being researched for their ability to facilitate a sustained, prolonged and
controlled release of drugs. But still much research is needed to evaluate and
fix issues relating toxicity, chemical and mechanical properties of the
nanoparticles. This project combines a range of nanobiotechnological techniques
to synthesise nano-components which are finally assembled together to get the
microcapsules which could be used in drug delivery systems. The microcapsules
were engineered keeping in mind their potential use in an ultrasound based drug
delivery system.


  The ribosomal density along the coding region of the mRNA molecule affect
various fundamental intracellular phenomena including: protein production
rates, organismal fitness, ribosomal drop off, and co-translational protein
folding. Thus, regulating translation in order to obtain a desired ribosomal
profile along the mRNA molecule is an important biological problem. We study
this problem using a model for mRNA translation, called the ribosome flow model
(RFM). In the RFM, the mRNA molecule is modeled as chain of n sites. The n
state-variables describe the ribosomal density profile along the mRNA molecule,
whereas the transition rates from each site to the next are controlled by n+1
positive constants. To study the problem of controlling the density profile, we
consider some or all of the transition rates as time-varying controls. We
consider the following problem: given an initial and a desired ribosomal
density profile, determine the time-varying values of the transition rates that
steer the RFM to this density profile, if they exist. Specifically, we consider
two control problems. In the first, all transition rates can be regulated and
the goal is to steer the ribosomal density profile and the protein production
rate from a given initial value to a desired value. In the second, a single
transition rate is controlled and the goal is to steer the production rate to a
desired value. In the first case, we show that the system is controllable, i.e.
the control is powerful enough to steer the RFM to any desired value, and we
provide closed-form expressions for constant control functions (or transition
rates) asymptotically steering the RFM to the desired value. For the second
problem, we show that the production rate can be steered to any desired value
in a feasible region determined by the other, constant transition rates. We
discuss some of the biological implications of these results.


  The Landau-Ginzburg-Wilson paradigm is applied to model the low-temperature
crystallographic C$\alpha$ backbone structure of sperm whale myoglobin. The
Glauber protocol is employed to simulate its response to an increase in ambient
temperature. The myoglobin is found to unfold from its native state by a
succession of $\alpha$-helical intermediates, fully in line with the observed
folding and unfolding patterns in denaturation experiments. In particular, a
molten globule intermediate is identified with experimentally correct
attributes. A detailed, experimentally testable contact map is constructed to
characterise the specifics of the unfolding pathway, including the formation of
long range interactions. The results reveal how the unfolding process of a
protein is driven by the interplay between, and a successive melting of, its
modular secondary structure components.


  Recent experiments have demonstrated that dynein motor exhibits catch bonding
behaviour, in which the unbinding rate of a single dynein decreases with
increasing force, for a certain range of force. Motivated by these experiments,
we propose a model for catch bonding in dynein using a threshold force bond
deformation (TFBD) model wherein catch bonding sets in beyond a critical
applied load force. We study the effect of catch bonding on unidirectional
transport properties of cellular cargo carried by multiple dynein motors within
the framework of this model. We find catch bonding can result in dramatic
changes in the transport properties, which are in sharp contrast to kinesin
driven unidirectional transport, where catch bonding is absent. We predict
that, under certain conditions, the average velocity of the cellular cargo can
actually increase as applied load is increased. We characterize the transport
properties in terms of a velocity profile phase plot in the parameter space of
the catch bond strength and the stall force of the motor. This phase plot
yields predictions that may be experimentally accessed by suitable
modifications of motor transport and binding properties. Our work necessitates
a reexamination of existing theories of collective bidirectional transport of
cellular cargo where the catch bond effect of dynein described in this paper is
expected to play a crucial role.


  Traditional Indian (Ayurvedic) and Chinese herbal medicines have been used in
the treatment of a variety of diseases for thousands of years because of their
natural origin and lesser side effects. However, the safety and efficacy data
(including dose and quality parameters) on most of these traditional medicines
are far from sufficient to meet the criteria needed to support their world-wide
therapeutic use. Also, the mechanistic understanding of most of these herbal
medicines is still lacking due to their complex components which further limits
their wider application and acceptance. Metabolomics -a novel approach to
reveal altered metabolism (biochemical effects) produced in response to a
disease or its therapeutic intervention- has huge potential to assess the
pharmacology and toxicology of traditional herbal medicines (THMs). Therefore,
it is gradually becoming a mutually complementary technique to genomics,
transcriptomics and proteomics for therapeutic evaluation of pharmaceutical
products (including THMs); the approach is so called pharmaco-metabolomics. The
whole paradigm is based on its ability to provide metabolic signatures to
confirm the diseased condition and then to use the concentration profiles of
these biomarkers to assess the therapeutic response. Nuclear magnetic resonance
(NMR) spectroscopy coupled with multivariate data analysis is currently the
method of choice for pharmaco-metabolomics studies owing to its unbiased,
non-destructive nature and minimal sample preparation requirement. In recent
past, dozens of NMR based pharmaco-metabolomic studies have been devoted to
prove the therapeutic efficacy/safety and to explore the underlying mechanisms
of THMs, with promising results. The current perspective article summarizes
various such studies in addition to describing the technical and conceptual
aspects involved in NMR based pharmaco-metabolomics.


  We use the optimal foraging theory to study coexistence between two plant
species and a generalist pollinator. We compare conditions for plant
coexistence for non-adaptive vs. adaptive pollinators that adjust their
foraging strategy to maximize fitness. When pollinators have fixed preferences,
we show that plant coexistence typically requires both weak competition between
plants for resources (e.g., space or nutrients) and pollinator preferences that
are not too biased in favour of either plant. We also show how plant
coexistence is promoted by indirect facilitation via the pollinator. When
pollinators are adaptive foragers, pollinator's diet maximizes pollinator's
fitness measured as the per capita population growth rate. Simulations show
that this has two conflicting consequences for plant coexistence. On the one
hand, when competition between pollinators is weak, adaptation favours
pollinator specialization on the more profitable plant which increases
asymmetries in plant competition and makes their coexistence less likely. On
the other hand, when competition between pollinators is strong, adaptation
promotes generalism, which facilitates plant coexistence. In addition, adaptive
foraging allows pollinators to survive sudden loss of the preferred plant host,
thus preventing further collapse of the entire community.
  Keywords: mutualism, competition, optimal foraging, evolutionarily stable
strategy, coexistence, adaptation rate


  We address an open question of Francis and Steel about phylogenetic networks
and trees. They give a polynomial time algorithm to decide if a phylogenetic
network, N, is tree-based and pose the problem: given a fixed tree T and
network N, is N based on T? We show that it is NP-hard to decide, by reduction
from 3-Dimensional Matching (3DM), and further, that the problem is fixed
parameter tractable.


  Repeated games have a long tradition in the behavioral sciences and
evolutionary biology. Recently, strategies were discovered that permit an
unprecedented level of control over repeated interactions by enabling a player
to unilaterally enforce linear constraints on payoffs. Here, we extend this
theory of "zero-determinant" (or, more generally, "autocratic") strategies to
alternating games, which are often biologically more relevant than traditional
synchronous games. Alternating games naturally result in asymmetries between
players because the first move matters or because players might not move with
equal probabilities. In a strictly-alternating game with two players, $X$ and
$Y$, we give conditions for the existence of autocratic strategies for player
$X$ when (i) $X$ moves first and (ii) $Y$ moves first. Furthermore, we show
that autocratic strategies exist even for (iii) games with randomly-alternating
moves. Particularly important categories of autocratic strategies are
extortionate and generous strategies, which enforce unfavorable and favorable
outcomes for the opponent, respectively. We illustrate these strategies using
the continuous Donation Game, in which a player pays a cost to provide a
benefit to the opponent according to a continuous cooperative investment level.
Asymmetries due to alternating moves could easily arise from dominance
hierarchies, and we show that they can endow subordinate players with more
autocratic strategies than dominant players.


  Construction of phylogenetic trees and networks for extant species from their
characters represents one of the key problems in phylogenomics. While solution
to this problem is not always uniquely defined and there exist multiple methods
for tree/network construction, it becomes important to measure how well the
constructed networks capture the given character relationship across the
species.
  In the current study, we propose a novel method for measuring the specificity
of a given phylogenetic network in terms of the total number of distributions
of character states at the leaves that the network may impose. While for binary
phylogenetic trees, this number has an exact formula and depends only on the
number of leaves and character states but not on the tree topology, the
situation is much more complicated for non-binary trees or networks.
Nevertheless, we develop an algorithm for combinatorial enumeration of such
distributions, which is applicable for arbitrary trees and networks under some
reasonable assumptions.


  We propose a model for membrane-cortex adhesion which couples membrane
deformations, hydrodynamics and kinetics of membrane-cortex ligands. In its
simplest form, the model gives explicit predictions for the critical pressure
for membrane detachment and for the value of adhesion energy. We show that
these quantities exhibit a significant dependence on the active acto-myosin
stresses. The model provides a simple framework to access quantitative
information on cortical activity by means of micropipette experiments. We also
extend the model to incorporate fluctuations and show that detailed information
on the stability of membrane-cortex coupling can be obtained by a combination
of micropipette aspiration and fluctuation spectroscopy measurements.


  Rate of growth of the amount of data generated using the high throughput
sequencing (HTS) platforms now exceeds the growth stipulated by Moore's Law.
The HTS data is expected to surpass those of other "big data" domains such as
astronomy, before the year 2025. In addition to sequencing genomes for research
purposes, genome and exome sequencing in clinical settings will be a routine
part of health care. The analysis of such large amounts of data, however, is
not without computational challenges. This burden is even more increased due to
the periodic updates to reference genomes, which typically require re-analysis
of existing data. Here we propose Coin-Application Mediator Interface (Coinami)
to distribute the workload for mapping reads to reference genomes using a
volunteer grid computer approach similar to Berkeley Open Infrastructure for
Network Computing (BOINC). However, since HTS read mapping requires substantial
computational resources and fast analysis turnout is desired, Coinami uses the
HTS read mapping as proof-of-work to generate valid blocks to main its own
cryptocurrency system, which may help motivate volunteers to dedicate more
resources. The Coinami protocol includes mechanisms to ensure that jobs
performed by volunteers are correct, and provides genomic data privacy. The
prototype implementation of Coinami is available at http://coinami.github.io/.


  In this article, we propose an approach to breeding which focuses on mating
instead of truncation selection, our method uses genome-wide marker information
in a similar fashion to genomic selection so we refer it to as genomic mating.
Using concepts of estimated breeding values, risk (usefulness) and inbreeding,
an efficient mating approach is formulated for improvement of breeding values
in the long run. We have used a genetic algorithm to find solutions to this
optimization problem. Results from our simulations point to the efficiency of
genomic mating for breeding complex traits compared to truncation selection.


  The diversity of T-cell receptors recognizing foreign pathogens is generated
through a highly stochastic recombination process, making the independent
production of the same sequence rare. Yet unrelated individuals do share
receptors, which together constitute a "public" repertoire of abundant
clonotypes. The TCR repertoire is initially formed prenatally, when the enzyme
inserting random nucleotides is downregulated, producing a limited diversity
subset. By statistically analyzing deep sequencing T-cell repertoire data from
twins, unrelated individuals of various ages, and cord blood, we show that
T-cell clones generated before birth persist and maintain high abundances in
adult organisms for decades, slowly decaying with age. Our results suggest that
large, low-diversity public clones are created during pregnancy, and survive
over long periods, providing the basis of the public repertoire.


  Motivation: Next-generation sequencing tools have enabled producing of huge
amount of genomic information at low cost. Unfortunately, presence of
sequencing errors in such data affects quality of downstream analyzes. Accuracy
of them can be improved by performing error correction. Because of huge amount
of such data correction algorithms have to: be fast, memory-frugal, and provide
high accuracy of error detection and elimination for variously-sized organisms.
  Results: We introduce a new algorithm for genomic data correction, capable of
processing eucaryotic 300 Mbp-genome-size, high error-rated data using less
than 4 GB of RAM in less than 40 minutes on 16-core CPU. The algorithm allows
to correct sequencing data at better or comparable level than competitors. This
was achieved by using very robust KMC~2 $k$-mer counter, new method of
erroneous regions correction based on both $k$-mer counts and FASTQ quality
indicators as well as careful optimization. Availability: Program is freely
available at http://sun.aei.posl.pl/REFRESH/reckoner. Contact:
sebastian.deorowicz@polsl.pl


  Environmental stochasticity is known to be a destabilizing factor, increasing
abundance fluctuations and extinction rates of populations. However, the
stability of a community may benefit from the differential response of species
to environmental variations due to the storage effect. This paper provides a
systematic and comprehensive discussion of these two contradicting tendencies,
using the metacommunity version of the recently proposed time-average neutral
model of biodiversity which incorporates environmental stochasticity and
demographic noise and allows for extinction and speciation. We show that the
incorporation of demographic noise into the model is essential to its
applicability, yielding realistic behavior of the system when fitness
variations are relatively weak. The dependence of species richness on the
strength of environmental stochasticity changes sign when the correlation time
of the environmental variations increases. This transition marks the point at
which the storage effect no longer succeeds in stabilizing the community.


  Given a digraph $G$, a lot of attention has been deserved on the maximum
number $\phi(G)$ of fixed points in a Boolean network $f:\{0,1\}^n\to\{0,1\}^n$
with $G$ as interaction graph. In particular, a central problem in network
coding consists in studying the optimality of the classical upper bound
$\phi(G)\leq 2^{\tau}$, where $\tau$ is the minimum size of a feedback vertex
set of $G$. In this paper, we study the maximum number $\phi_m(G)$ of fixed
points in a {\em monotone} Boolean network with interaction graph $G$. We
establish new upper and lower bounds on $\phi_m(G)$ that depends on the cycle
structure of $G$. In addition to $\tau$, the involved parameters are the
maximum number $\nu$ of vertex-disjoint cycles, and the maximum number
$\nu^{*}$ of vertex-disjoint cycles verifying some additional technical
conditions. We improve the classical upper bound $2^\tau$ by proving that
$\phi_m(G)$ is at most the largest sub-lattice of $\{0,1\}^\tau$ without chain
of size $\nu+1$, and without another forbidden-pattern of size $2\nu^{*}$.
Then, we prove two optimal lower bounds: $\phi_m(G)\geq \nu+1$ and
$\phi_m(G)\geq 2^{\nu^{*}}$. As a consequence, we get the following
characterization: $\phi_m(G)=2^\tau$ if and only if $\nu^{*}=\tau$. As another
consequence, we get that if $c$ is the maximum length of a chordless cycle of
$G$ then $2^{\nu/3^c}\leq\phi_m(G)\leq 2^{c\nu}$. Finally, with the technics
introduced, we establish an upper bound on the number of fixed points of any
Boolean network according to its signed interaction graph.


  Within bioinformatics, the textual alignment of amino acid sequences has long
dominated the determination of similarity between proteins, with all that
implies for shared structure, function and evolutionary descent. Despite the
relative success of modern-day sequence alignment algorithms, so-called
alignment-free approaches offer a complementary means of determining and
expressing similarity, with potential benefits in certain key applications,
such as regression analysis of protein structure-function studies, where
alignment-base similarity has performed poorly. Here, we offer a fresh,
statistical physics-based perspective focusing on the question of
alignment-free comparison, in the process adapting results from 'first passage
probability distribution' to summarize statistics of ensemble averaged amino
acid propensity values. In this article, we introduce and elaborate this
approach.


  Real-time vaccination following an outbreak can effectively mitigate the
damage caused by an infectious disease. However, in many cases, available
resources are insufficient to vaccinate the entire at-risk population,
logistics result in delayed vaccine deployment, and the interaction between
members of different cities facilitates a wide spatial spread of infection.
Limited vaccine, time delays, and interaction (or coupling) of cities lead to
tradeoffs that impact the overall magnitude of the epidemic. These tradeoffs
mandate investigation of optimal strategies that minimize the severity of the
epidemic by prioritizing allocation of vaccine to specific subpopulations. We
use an SIR model to describe the disease dynamics of an epidemic which breaks
out in one city and spreads to another. We solve a master equation to determine
the resulting probability distribution of the final epidemic size. We then
identify tradeoffs between vaccine, time delay, and coupling, and we determine
the optimal vaccination protocols resulting from these tradeoffs.


  We study the nucleation of blebs, i.e., protrusions arising from a local
detachment of the membrane from the cortex of a cell. Based on a simple model
of elastic linkers with force-dependent kinetics, we show that bleb nucleation
is governed by membrane peeling. By this mechanism, the growth or shrinkage of
a detached membrane patch is completely determined by the linker kinetics,
regardless of the energetic cost of the detachment. We predict the critical
nucleation radius for membrane peeling and the corresponding effective energy
barrier. These may be typically smaller than those predicted by classical
nucleation theory, implying a much faster nucleation. We also perform
simulations of a continuum stochastic model of membrane-cortex adhesion to
obtain the statistics of bleb nucleation times as a function of the stress on
the membrane. The determinant role of membrane peeling changes our
understanding of bleb nucleation and opens new directions in the study of
blebs.


  The sustainability of structured biological, social, economic and ecological
communities are often determined by the outcome of social conflicts between
cooperative and selfish individuals (cheaters). Cheaters avoid the cost of
contributing to the community and can occasionally spread in the population
leading to the complete collapse of cooperation. Although such a collapse often
unfolds unexpectedly bearing the traits of a critical transition, it is unclear
whether one can detect the rising risk of cheater's invasions and loss of
cooperation in an evolving community. Here, we combine dynamical networks and
evolutionary game theory to study the abrupt loss of cooperation as a critical
transition. We estimate the risk of collapse of cooperation after the
introduction of a single cheater under gradually changing conditions. We
observe a systematic increase in the average time it takes for cheaters to be
eliminated from the community as the risk of collapse increases. We detect this
risk based on changes in community structure and composition. Nonetheless,
reliable detection depends on the mechanism that governs how cheaters evolve in
the community. Our results suggest possible avenues for detecting the loss of
cooperation in evolving communities


  Molecular simulations can provide microscopic insight into the physical and
chemical driving forces of complex molecular processes. Despite continued
advancement of simulation methodology, model errors may lead to inconsistencies
between simulated and reference (e.g., from experiments or higher-level
simulations) observables. To bound the microscopic information generated by
computer simulations within reference measurements, we propose a method that
reweights the microscopic transitions of the system to improve consistency with
a set of coarse kinetic observables. The method employs the well-developed
Markov state modeling framework to efficiently link microscopic dynamics with
long-time scale constraints, thereby consistently addressing a wide range of
time scales. To emphasize the robustness of the method, we consider two
distinct coarse-grained models with significant kinetic inconsistencies. When
applied to the simulated conformational dynamics of small peptides, the
reweighting procedure systematically improves the time scale separation of the
slowest processes. Additionally, constraining the forward and backward rates
between metastable states leads to slight improvement of their relative
stabilities and, thus, refined equilibrium properties of the resulting model.
Finally, we find that difficulties in simultaneously describing both the
simulated data and the provided constraints can help identify specific
limitations of the underlying simulation approach.


  Motivated by applications in domains such as social networks and
computational biology, we study the problem of community recovery in graphs
with locality. In this problem, pairwise noisy measurements of whether two
nodes are in the same community or different communities come mainly or
exclusively from nearby nodes rather than uniformly sampled between all nodes
pairs, as in most existing models. We present an algorithm that runs nearly
linearly in the number of measurements and which achieves the information
theoretic limit for exact recovery.


  The galactose network is a complex system responsible for galactose
metabolism. It has been extensively studied experimentally and mathematically
at the unicellular level to broaden our understanding of its regulatory
mechanisms at higher order species. Although the key molecular players involved
in the metabolic and regulatory processes underlying this system have been
known for decades, their interactions and chemical kinetics remain incompletely
understood. Mathematical models can provide an alternative method to study the
dynamics of this network from a quantitative and a qualitative perspective.
Here, we employ such approaches to unravel the main properties of the galactose
network, including equilibrium binary and temporal responses, as a way to
decipher its adaptation to actively-changing inputs. We combine the two main
components of the network; namely, the genetic branch which allows for bistable
responses, and a metabolic branch, encompassing the relevant metabolic
processes and glucose repressive reactions. We use both computational tools to
estimate model parameters based on published experimental data, as well as
bifurcation analysis to decipher the properties of the system in various
parameter regimes. Our model analysis reveals that the interplay between the
inducer (galactose) and the repressor (glucose) creates the bistability regime
which dictates the temporal responses of the system. Based on the same
bifurcation techniques, we can also explain why the system is robust to genetic
mutations and molecular instabilities. These findings may provide
experimentalists with a theoretical framework upon which they can determine how
the galactose network functions under various conditions.


  Accurate estimation of evolutionary distances between taxa is important for
many phylogenetic reconstruction methods. In the case of bacteria, distances
can be estimated using a range of different evolutionary models, from single
nucleotide polymorphisms to large-scale genome rearrangements. In the case of
sequence evolution models (such as the Jukes-Cantor model and associated
metric) have been used to correct pairwise distances. Similar correction
methods for genome rearrangement processes are required to improve inference.
Current attempts at correction fall into 3 categories: Empirical computational
studies, Bayesian/MCMC approaches, and combinatorial approaches. Here we
introduce a maximum likelihood estimator for the inversion distance between a
pair of genomes, using the group-theoretic approach to modelling inversions
introduced recently. This MLE functions as a corrected distance: in particular,
we show that because of the way sequences of inversions interact with each
other, it is quite possible for minimal distance and MLE distance to
differently order the distances of two genomes from a third. This has obvious
implications for the use of minimal distance in phylogeny reconstruction. The
work also tackles the above problem allowing free rotation of the genome.
Generally a frame of reference is locked, and all computation made accordingly.
This work incorporates the action of the dihedral group so that distance
estimates are free from any a priori frame of reference.


  Forest regeneration is a spatial multifaceted process with numerous unknown
variables. The individual area or the area potentially available (APA) to an
individual plant embodies an idea widely used in population ecology but it has
fewer applications in forest research. It was used a Voronoi/ Thiessen
tessellation in order to determine the individual area of each sapling. The
study was conducted in a naturally regenerated area; using data collected from
over seven thousands saplings positioned in a network of permanent rectangular
sampling plots. The Voronoi/ Thiessen polygons were used to characterize
spatial pattern of sapling distribution as well as the competition relations
between the individuals. It is obvious that, at least from the mathematically
point of view, the Voronoi tessellation represents one of the best solutions to
determine neighbouring competitors of a tree. There were studied the
correlations between APA values and the main biometrical attributes, height
growth and competition indices. Furthermore, it is shown that APA coefficient
of variation is a straight-forward indicator with positive results as an
indicator of spatial pattern. The statistical significance of this indicator
was evaluated by comparing the results with the values of a 95% confidence
envelope generated by Monte-Carlo simulations. Two practical software tools
were produced using Visual Basic (VORONOI and ARIA VORONOI) in order to
simplify the analyses.


  Ranaspumin-2 (Rsn-2) is a surfactant protein found in the foam nests of the
t\'{u}ngara frog. Previous experimental work has led to a proposed model of
adsorption which involves an unusual clam shell-like `unhinging' of the protein
at an interface. Interestingly, there is no concomitant denaturation of the
secondary structural elements of Rsn-2 with the large scale transformation of
its tertiary structure. In this work we use both experiment and simulation to
better understand the driving forces underpinning this unusual process. We
develop a modified G\={o}-model approach where we have included explicit
representation of the side-chains in order to realistically model the
interaction between the secondary structure elements of the protein and the
interface. Doing so allows for the study of the underlying energy landscape
which governs the mechanism of Rsn-2 interfacial adsorption. Experimentally, we
study targeted mutants of Rsn-2, using the Langmuir trough, pendant drop
tensiometry and circular dichroism, to demonstrate that the clam-shell model is
correct. We find that Rsn-2 adsorption is in fact a two-step process: the
hydrophobic N-terminal tail recruits the protein to the interface after which
Rsn-2 undergoes an unfolding transition which maintains its secondary
structure. Intriguingly, our simulations show that the conformation Rsn-2
adopts at an interface is an arrested state along the denaturation pathway.
More generally, our computational model should prove a useful, and
computationally efficient, tool in studying the dynamics and energetics of
protein-interface interactions.


  It is possible to model vector-borne infection using the classical
Ross-Macdonald model. This attempt, however fails in several respects. First,
using measured (or estimated) parameters, the model predicts a much greater
number of cases than what is usually observed. Second, the model predicts a
single huge outbreaks that is followed after decades of much smaller outbreaks.
This is not what is observed. Usually towns or cities report a number of cases
that recur for many years, even when environmental changes cannot explain the
disappearance of the infection in-between the peaks. In this paper we continue
to examine the pitfalls in modeling this class of infections, and explain that,
in fact, if properly used, the Ross-Macdonald model works, can be used to
understand the patterns of epidemics and even, to some extents, to make some
predictions. We model several outbreaks of dengue fever and show that the
variable pattern of year recurrence (or absence of it) can be understood and
explained by a simple Ross-Macdonald model modified to take into account human
movement across a range of neighborhoods inside a city. In addition, we analyze
the effect of seasonal variations in the parameters determining the number,
longevity and biting behavior of mosquitoes. Based on the size of the first
outbreak, we show that it is possible to estimate the proportion of the
remaining susceptibles and predict the likelihood and magnitude of eventual
subsequent outbreaks. The approach is exemplified by actual dengue outbreaks
with different recurrence patterns from some Brazilian regions.


  Adenosine triphosphate (ATP) is the universal chemical energy currency for
cellular activities provided mainly by the membrane enzyme FoF1-ATP synthase in
bacteria, chloroplasts and mitochondria. Synthesis of ATP is accompanied by
subunit rotation within the enzyme. Over the past 15 years we have developed a
variety of single-molecule FRET (smFRET) experiments to monitor catalytic
action of individual bacterial enzymes in vitro. By specifically labeling
rotating and static subunits within a single enzyme we were able to observe
three-stepped rotation in the F1 motor, ten-stepped rotation in the Fo motor
and transient elastic deformation of the connected rotor subunits. However, the
spatial and temporal resolution of motor activities measured by smFRET were
limited by the photophysics of the FRET fluorophores. Here we evaluate the
novel FRET donor mNeonGreen as a fusion to FoF1-ATP synthase and compare it to
the previously used fluorophore EGFP. Topics of this manuscript are the
biochemical purification procedures and the activity measurements of the fully
functional mutant enzyme.


  Phylogenetic trees are pervasively used to depict evolutionary relationships.
Increasingly, researchers need to visualize large trees and compare multiple
large trees inferred for the same set of taxa (reflecting uncertainty in the
tree inference or genuine discordance among the loci analysed). Existing tree
visualization tools are however not well suited to these tasks. In particular,
side-by-side comparison of trees can prove challenging beyond a few dozen taxa.
Here, we introduce Phylo.io, a web application to visualize and compare
phylogenetic trees side-by-side. Its distinctive features are: highlighting of
similarities and differences between two trees, automatic identification of the
best matching rooting and leaf order, scalability to very large trees, high
usability, multiplatform support via standard HTML5 implementation, and
possibility to store and share visualisations. The tool can be freely accessed
at http://phylo.io. The code for the associated JavaScript library is available
at https://github.com/DessimozLab/phylo-io under an MIT open source license.


  Protein structure prediction has been a grand challenge problem in the
structure biology over the last few decades. Protein quality assessment plays a
very important role in protein structure prediction. In the paper, we propose a
new protein quality assessment method which can predict both local and global
quality of the protein 3D structural models. Our method uses both multi and
single model quality assessment method for global quality assessment, and uses
chemical, physical, geo-metrical features, and global quality score for local
quality assessment. CASP9 targets are used to generate the features for local
quality assessment. We evaluate the performance of our local quality assessment
method on CASP10, which is comparable with two stage-of-art QA methods based on
the average absolute distance between the real and predicted distance. In
addition, we blindly tested our method on CASP11, and the good performance
shows that combining single and multiple model quality assessment method could
be a good way to improve the accuracy of model quality assessment, and the
random forest technique could be used to train a good local quality assessment
model.


  For families, kinship coefficients are quantifications of the amount of
genetic sharing between a pair of individuals. These coefficients are critical
for understanding the breeding habits and genetic diversity of diploid
populations. Historically, computations of the inbreeding coefficient were used
to prohibit inbred marriages and prohibit breeding of some pairs of pedigree
animals. Such prohibitions foster genetic diversity and help prevent recessive
Mendelian disease at a population level.
  This paper gives the fastest known algorithms for computing the kinship
coefficient of a set of individuals with a known pedigree. The algorithms given
here consider the possibility that the founders of the known pedigree may
themselves be inbred, and they compute the appropriate inbreeding-adjusted
kinship coefficients. The exact kinship algorithm has running-time $O(n^2)$ for
an $n$-individual pedigree. The recursive-cut exact kinship algorithm has
running time $O(s^2m)$ where $s$ is the number of individuals in the largest
segment of the pedigree and $m$ is the number of cuts. The approximate
algorithm has running-time $O(n)$ for an $n$-individual pedigree on which to
estimate the kinship coefficients of $\sqrt{n}$ individuals from $\sqrt{n}$
founder kinship coefficients.


  Bz-423 is a promising new drug for treatment of autoimmune diseases. This
small molecule binds to subunit OSCP of the mitochondrial enzyme FoF1-ATP
synthase and modulates its catalytic activities. We investigate the binding of
Bz-423 to mitochondria in living cells and how subunit rotation in FoF1-ATP
synthase, i.e. the mechanochemical mechanism of this enzyme, is affected by
Bz-423. Therefore, the enzyme was marked selectively by genetic fusion with the
fluorescent protein EGFP to the C terminus of subunit gamma. Imaging the
threedimensional arrangement of mitochondria in living yeast cells was possible
at superresolution using structured illumination microscopy, SIM. We measured
uptake and binding of a Cy5-labeled Bz-423 derivative to mitochondrial FoF1-ATP
synthase in living yeast cells using FRET acceptor photobleaching microscopy.
Our data confirmed the binding of Cy5-labeled Bz-423 to the top of the F1
domain of the enzyme in mitochondria of living Saccharomyces cerevisiae cells.


  The reconstruction of phylogenetic trees from molecular sequence data relies
on modelling site substitutions by a Markov process, or a mixture of such
processes. In general, allowing mixed processes can result in different tree
topologies becoming indistinguishable from the data, even for infinitely long
sequences. However, when the underlying Markov process supports linear
phylogenetic invariants, then provided these are sufficiently informative, the
identifiability of the tree topology can be restored. In this paper, we
investigate a class of processes that support linear invariants once the
stationary distribution is fixed, the `equal input model'. This model
generalizes the `Felsenstein 1981' model (and thereby the Jukes--Cantor model)
from four states to an arbitrary number of states (finite or infinite), and it
can also be described by a `random cluster' process. We describe the structure
and dimension of the vector spaces of phylogenetic mixtures and of linear
invariants for any fixed phylogenetic tree (and for all trees -- the so called
`model invariants'), on any number $n$ of leaves. We also provide a precise
description of the space of mixtures and linear invariants for the special case
of $n=4$ leaves. By combining techniques from discrete random processes and
(multi-) linear algebra, our results build on a classic result that was first
established by James Lake in 1987.


  Nosocomial pathogens such as Methicillin-Resistant {\em Staphylococcus
aureus} (MRSA) and Vancomycin-resistant {\em Enterococci} (VRE) are the cause
of significant morbidity and mortality among hospital patients. It is important
to be able to assess the efficacy of control measures using data on patient
outcomes. In this paper we describe methods for analysing such data using
patient-level stochastic models which seek to describe the underlying
unobserved process of transmission. The methods are applied to detailed
longitudinal patient-level data on VRE from a study in a US hospital with eight
intensive care units (ICUs). The data comprise admission and discharge dates,
dates and results of screening tests, and dates during which precautionary
measures were in place for each patient during the study period. Results
include estimates of the efficacy of the control measures, the proportion of
unobserved patients colonized with VRE and the proportion of patients colonized
on admission.


  The specificity of molecular recognition is important to molecular
self-organization. A prominent example is the biological cell where, within a
highly crowded molecular environment, a myriad of different molecular receptor
pairs recognize their binding partner with astonishing accuracy. In thermal
equilibrium it is usually admitted that the affinity of recognizer pairs only
depends on the nature of the two binding molecules. Accordingly, Boltzmann
factors of binding energy differences relate the molecular affinities among
different target molecules that compete for the same probe. Here, we consider
the molecular recognition of short DNA oligonucleotide single strands. We show
that a better matching oligonucleotide strand can prevail against a
disproportionally more concentrated competitor that exhibits reduced affinity
due to a mismatch. The magnitude of deviation from the simple picture above may
reach several orders of magnitude. In our experiments the effective molecular
affinity of a given strand remains elevated only as long as the better matching
competitor is not present. We interpret our observations based on an
energy-barrier of entropic origin that occurs if two competing oligonucleotide
strands occupy the same probe simultaneously. In this situation the relative
binding affinities are reduced asymmetrically, which leads to an expression of
the free energy landscape that represents a formal analogue of a Landau
description of phase transitions. Our mean field description reproduces the
observations in quantitative agreement. The advantage of improved molecular
recognition comes at no energetic cost other than the design of the molecular
ensemble, and the introduction of the competitor. It will be interesting to see
if mechanisms along similar lines as exposed here, contribute to the molecular
synergy that occurs in biological systems.


  The interplay between energy efficiency and evolutionary mechanisms is
addressed. One important question is how evolutionary mechanisms can select for
the optimised usage of energy in situations where it does not lead to immediate
advantage. For example, this problem is of great importance to improve our
understanding about the major transition from unicellular to multicellular form
of life. The immediate advantage of gathering efficient individuals in an
energetic context is not clear. Although this process increases relatedness
among individuals, it also increases local competition. To address this
question, we propose a model of two competing metabolic strategies that makes
explicit reference to the resource usage. We assume the existence of an
efficient strain, which converts resource into energy at high efficiency but
displays a low rate of resource consumption, and an inefficient strain, which
consumes resource at a high rate with a low efficiency in converting it to
energy. We explore the dynamics in both well-mixed and structured populations.
The selection for optimised energy usage is measured by the likelihood of that
an efficient strain can invade a population only comprised by inefficient
strains. It is found that the region of the parameter space at which the
efficient strain can thrive in structured populations is always larger than
observed in well-mixed populations. In fact, in well-mixed populations the
efficient strain is only evolutionarily stable in the domain whereupon there is
no evolutionary dilemma. We also observe that small group sizes enhance the
chance of invasion by the efficient strain in spite of increasing the
competition among relatives. This outcome corroborates the key role played by
kin selection and shows that the group dynamics relied on group expansion,
overlapping generations and group split can balance the negative effects of
local competition.


  Overlapping genes exist in all domains of life and are much more abundant
than expected at their first discovery in the late 1970s. Assuming that the
reference gene is read in frame +0, an overlapping gene can be encoded in two
reading frames in the sense strand, denoted by +1 and +2, and in three reading
frames in the opposite strand, denoted by -0, -1 and -2. This motivated
numerous researchers to study the constraints induced by the genetic code on
the various overlapping frames, mostly based on information theory. Our focus
in this paper is on the constraints induced on two overlapping genes in terms
of amino acids, as well as polypeptides. We show that simple linear constraints
bind the amino acid composition of two proteins encoded by overlapping genes.
Novel constraints are revealed when polypeptides are considered, and not just
single amino acids. For example, in double-coding sequences with an overlapping
reading frame -2, each Tyrosine (denoted as Tyr or Y) in the overlapping frame
overlaps a Tyrosine in the reference frame +0 (and reciprocally), whereas
specific words (e.g. YY) never occur. We thus distinguish between null
constraints (YY = 0 in frame -2) and non-null constraints (Y in frame +0 <=> Y
in frame -2). Our equivalence-based constraints are symmetrical and thus enable
the characterization of the joint composition of overlapping proteins. We
describe several formal frameworks and a graph algorithm to characterize and
compute these constraints. These results yield support for understanding the
mechanisms and evolution of overlapping genes, and for developing novel
overlapping gene detection methods.


  Cells can use periodic enzyme activities to adapt to periodic environments or
existing internal rhythms and to establish metabolic cycles that schedule
biochemical processes in time. A periodically changing allocation of the
protein budget between reactions or pathways may increase the overall metabolic
efficiency. To study this hypothesis, I quantify the possible benefits of
small-amplitude enzyme rhythms in kinetic models. Starting from an
enzyme-optimised steady state, I score the effects of possible enzyme rhythms
on a metabolic objective and optimise their amplitudes and phase shifts.
Assuming small-amplitude rhythms around an optimal reference state, optimal
phases and amplitudes can be computed by solving a quadratic optimality
problem. In models without amplitude constraints, general periodic enzyme
profiles can be obtained by Fourier synthesis. The theory of optimal enzyme
rhythms combines the dynamics and economics of metabolic systems and explains
how optimal small-amplitude enzyme profiles are shaped by network structure,
kinetics, external rhythms, and the metabolic objective. The formulae show how
orchestrated enzyme rhythms can exploit synergy effects to improve metabolic
performance and that optimal enzyme profiles are not simply adapted to existing
metabolic rhythms, but that they actively shape these rhythms to improve their
own (and other enzymes') efficiency. The resulting optimal enzyme profiles
"portray" the enzymes' dynamic effects in the network: for example, enzymes
that act synergistically may be coexpressed, periodically and with some optimal
phase shifts. The theory yields optimality conditions for enzyme rhythms in
metabolic cycles, with static enzyme adaptation as a special case, and predicts
how cells should combine transcriptional and posttranslational regulation to
realise enzyme rhythms at different frequencies.


  The number of substitutions (of nucleotides, amino acids, ...) that take
place during the evolution of a sequence is a stochastic variable of
fundamental importance in the field of molecular evolution. Although the mean
number of substitutions during molecular evolution of a sequence can be
estimated for a given substitution model, no simple solution exists for the
variance of this random variable. We show in this article that the computation
of the variance is as simple as that of the mean number of substitutions for
both short and long times. Apart from its fundamental importance, this result
can be used to investigate the dispersion index R , i.e. the ratio of the
variance to the mean substitution number, which is of prime importance in the
neutral theory of molecular evolution. By investigating large classes of
substitution models, we demonstrate that although R\ge1 , to obtain R
significantly larger than unity necessitates in general additional hypotheses
on the structure of the substitution model.


  Aim: Molecular dynamics simulations and normal mode analysis are
well-established approaches to generate receptor conformational ensembles
(RCEs) for ligand docking and virtual screening. Here, we report new fast
molecular dynamics-based and normal mode analysis-based protocols combined with
conformational pocket classifications to efficiently generate RCEs. Materials
\& methods: We assessed our protocols on two well-characterized protein targets
showing local active site flexibility, dihydrofolate reductase and large
collective movements, CDK2. The performance of the RCEs was validated by
distinguishing known ligands of dihydrofolate reductase and CDK2 among a
dataset of diverse chemical decoys. Results \& discussion: Our results show
that different simulation protocols can be efficient for generation of RCEs
depending on different kind of protein flexibility.


  Current evolutionary biology models usually assume that a phenotype undergoes
gradual change. This is in stark contrast to biological intuition, which
indicates that change can also be punctuated-the phenotype can jump. Such a
jump could especially occur at speciation, i.e. dramatic change occurs that
drives the species apart. Here we derive a Central Limit Theorem for punctuated
equilibrium. We show that, if adaptation is fast, for weak convergence to
normality to hold, the variability in the occurrence of change has to disappear
with time.


  Mean-field models approximate large stochastic systems by simpler
differential equations that are supposed to approximate the mean of the larger
system. It is generally assumed that as the stochastic systems get larger
(i.e., more people or particles), they converge to the mean-field models.
Mean-field models are common in many fields, but their convergence is rarely
proved. The existing approaches rely on operator semigroups, martingales, PDEs,
or infinite systems of ODEs. We give a general proof for their convergence
using only Taylor's theorem and basic ODE results. We hope this allows applied
researchers to routinely show convergence of their mean-field models, putting
their work on a stronger foundation.


  This paper presents a compact pairwise model that describes the spread of
multi-stage epidemics on networks. The multi-stage model corresponds to a
gamma-distributed infectious period which interpolates between the classical
Markovian models with exponentially distributed infectious period and epidemics
with a constant infectious period. We show how the compact approach leads to a
system of equations whose size is independent of the range of node degrees,
thus significantly reducing the complexity of the model. Network clustering is
incorporated into the model to provide a more accurate representation of
realistic contact networks, and the accuracy of proposed closures is analysed
for different levels of clustering and number of infection stages. Our results
support recent findings that standard closure techniques are likely to perform
better when the infectious period is constant.


  Metagenomic profiling is challenging in part because of the highly uneven
sampling of the tree of life by genome sequencing projects and the limitations
imposed by performing phylogenetic inference at fixed taxonomic ranks. We
present the algorithm MetaPalette which uses long $k$-mer sizes ($k=30, 50$) to
fit a $k$-mer "palette" of a given sample to the $k$-mer palette of reference
organisms. By modeling the $k$-mer palettes of unknown organisms, the method
also gives an indication of the presence, abundance, and evolutionary
relatedness of novel organisms present in the sample. The method returns a
traditional, fixed-rank taxonomic profile which is shown on independently
simulated data to be one of the most accurate to date. Tree figures are also
returned that quantify the relatedness of novel organisms to reference
sequences and the accuracy of such figures is demonstrated on simulated
spike-ins and a metagenomic soil sample. The software implementing MetaPalette
is available at: https://github.com/dkoslicki/MetaPalette. Pre-trained
databases are included for Archaea, Bacteria, Eukaryota, and viruses.


  We analyze velocity-jump process models of persistent search for a single
target on a bounded domain. The searcher proceeds along ballistic trajectories
and is absorbed upon collision with the target boundary. When reaching the
domain boundary, the searcher chooses a random direction for its new
trajectory. For circular domains and targets, we can approximate the mean first
passage time (MFPT) using a Markov chain approximation of the search process.
Our analysis and numerical simulations reveal that the time to find the target
decreases for targets closer to the domain boundary. When there is a small
probability of direction-switching within the domain, we find the time to find
the target decreases slightly with the turning probability. We also extend our
exit time analysis to the case of partitioned domains, where there is a single
target within one of multiple disjoint subdomains. Given an average time of
transition between domains $\langle T \rangle$, we find that the optimal rate
of transition that minimizes the time to find the target obeys $\beta_{{\rm
min}} \propto 1/ \sqrt{\langle T \rangle}$.


  A rigorous understanding of how multicellular behaviors arise from the
actions of single cells requires quantitative frameworks that bridge the gap
between genetic circuits, the arrangement of cells in space, and
population-level behaviors. Here, we provide such a framework for a ubiquitous
class of multicellular systems - namely, "secrete-and-sense cells" that
communicate by secreting and sensing a signaling molecule. By using formal,
mathematical arguments and introducing the concept of a phenotype diagram, we
show how these cells tune their degrees of autonomous and collective behavior
to realize distinct single-cell and population-level phenotypes; these
phenomena have biological analogs, such as quorum sensing or paracrine
signaling. We also define the "entropy of population," a measurement of the
number of arrangements that a population of cells can assume, and demonstrate
how a decrease in the entropy of population accompanies the formation of
ordered spatial patterns. Our conceptual framework ties together diverse
systems, including tissues and microbes, with common principles.


  Specific interactions are a hallmark feature of self-assembly and
signal-processing systems in both synthetic and biological settings.
Specificity between components may arise from a wide variety of physical and
chemical mechanisms in diverse contexts, from DNA hybridization to
shape-sensitive depletion interactions. Despite this diversity, all systems
that rely on interaction specificity operate under the constraint that
increasing the number of distinct components inevitably increases off-target
binding. Here we introduce `capacity', the maximal information encodable using
specific interactions, to compare specificity across diverse experimental
systems, and to compute how specificity changes with physical parameters. Using
this framework, we find that `shape'-coding of interactions has higher capacity
than chemical (`color') coding because the strength of off-target binding is
strongly sublinear in binding site size for shapes while being linear for
colors. We also find that different specificity mechanisms, such as shape and
color, can be combined in a synergistic manner, giving a capacity greater than
the sum of the parts.


  DNA strand displacement is a key reaction in DNA homologous recombination and
DNA mismatch repair and is also heavily utilized in DNA-based computation and
locomotion. Despite its ubiquity in science and engineering, sequence-dependent
effects of displacement kinetics have not been extensively characterized. Here,
we measured toehold-mediated strand displacement kinetics using single-molecule
fluorescence in the presence of a single base pair mismatch. The apparent
displacement rate varied significantly when the mismatch was introduced in the
invading DNA strand. The rate generally decreased as the mismatch in the
invader was encountered earlier in displacement. Our data indicate that a
single base pair mismatch in the invader stalls branch migration, and
displacement occurs via direct dissociation of the destabilized incumbent
strand from the substrate strand. We combined both branch migration and direct
dissociation into a model, which we term, the concurrent displacement model,
and used the first passage time approach to quantitatively explain the salient
features of the observed relationship. We also introduce the concept of
splitting probabilities to justify that the concurrent model can be simplified
into a three-step sequential model in the presence of an invader mismatch. We
expect our model to become a powerful tool to design DNA-based reaction schemes
with broad functionality.