  How do shared conventions emerge in complex decentralized social systems?
This question engages fields as diverse as linguistics, sociology and cognitive
science. Previous empirical attempts to solve this puzzle all presuppose that
formal or informal institutions, such as incentives for global agreement,
coordinated leadership, or aggregated information about the population, are
needed to facilitate a solution. Evolutionary theories of social conventions,
by contrast, hypothesize that such institutions are not necessary in order for
social conventions to form. However, empirical tests of this hypothesis have
been hindered by the difficulties of evaluating the real-time creation of new
collective behaviors in large decentralized populations. Here, we present
experimental results - replicated at several scales - that demonstrate the
spontaneous creation of universally adopted social conventions, and show how
simple changes in a population's network structure can direct the dynamics of
norm formation, driving human populations with no ambition for large scale
coordination to rapidly evolve shared social conventions.


  A binary phylogenetic network may or may not be obtainable from a tree by the
addition of directed edges (arcs) between tree arcs. Here, we establish a
precise and easily tested criterion (based on `2-SAT') that efficiently
determines whether or not any given network can be realized in this way.
Moreover, the proof provides a polynomial-time algorithm for finding one or
more trees (when they exist) on which the network can be based. A number of
interesting consequences are presented as corollaries; these lead to some
further relevant questions and observations, which we outline in the
conclusion.


  We study four citrate synthase homodimeric proteins within a structure-based
coarse-grained model. Two of these proteins come from thermophilic bacteria,
one from a cryophilic bacterium and one from a mesophilic organism; three are
in the closed and two in the open conformations. Even though the proteins
belong to the same fold, the model distinguishes the properties of these
proteins in a way which is consistent with experiments. For instance, the
thermophilic proteins are more stable thermodynamically than their mesophilic
and cryophilic homologues, which we observe both in the magnitude of thermal
fluctuations near the native state and in the kinetics of thermal unfolding.
The level of stability correlates with the average coordination number for
amino acids contacts and with the degree of structural compactness. The pattern
of positional fluctuations along the sequence in the closed conformation is
different than in the open conformation, including within the active site. The
modes of correlated and anticorrelated movements of pairs of amino acids
forming the active site are very different in the open and closed
conformations. Taken together, our results show that the precise location of
amino acid contacts in the native structure appears to be a critical element in
explaining the similarities and differences in the thermodynamic properties,
local flexibility and collective motions of the different forms of the enzyme.


  The study of genetic map linearization leads to a combinatorial hard problem,
called the {\em minimum breakpoint linearization} (MBL) problem. It is aimed at
finding a linearization of a partial order which attains the minimum breakpoint
distance to a reference total order. The approximation algorithms previously
developed for the MBL problem are only applicable to genetic maps in which
genes or markers are represented as signed integers. However, current genetic
mapping techniques generally do not specify gene strandedness so that genes can
only be represented as unsigned integers. In this paper, we study the MBL
problem in the latter more realistic case. An approximation algorithm is thus
developed, which achieves a ratio of $(m^2+2m-1)$ and runs in $O(n^7)$ time,
where $m$ is the number of genetic maps used to construct the input partial
order and $n$ the total number of distinct genes in these maps.


  Using coarse-grained molecular dynamics simulations, we analyze mechanically
induced dissociation and unfolding of the protein complex CD48-2B4. This
heterodimer is an indispensable component of the immunological system: 2B4 is a
receptor on natural killer cells whereas CD48 is expressed on surfaces of
various immune cells. So far, its mechanostability has not been assessed either
experimentally or theoretically. We find that the dissociation processes
strongly depend on the direction of pulling and may take place in several
pathways. Interestingly, the CD48-2B4 interface can be divided into three
distinct patches that act as units when resisting the pulling forces. At
experimentally accessible pulling speeds, the characteristic mechanostability
forces are in the range between 100 and 200 pN, depending on the pulling
direction. These characteristic forces need not be associated with tensile
forces involved in the act of separation of the complex because prior
shear-involving unraveling within individual proteins may give rise to a higher
force peak.


  This paper studies the relations among system parameters, uniqueness, and
stability of equilibria, for kinetic systems given in the form of polynomial
ODEs. Such models are commonly used to describe the dynamics of nonnegative
systems, with a wide range of application fields such as chemistry, systems
biology, process modeling or even transportation systems. Using a flux-based
description of kinetic models, a canonical representation of the set of all
possible feasible equilibria is developed. The characterization is made in
terms of strictly stable compartmental matrices to define the so-called family
of solutions. Feasibility is imposed by a set of constraints, which are linear
on a log-transformed space of complexes, and relate to the kernel of a matrix,
the columns of which span the stoichiometric subspace. One particularly
interesting representation of these constraints can be expressed in terms of a
class of monotonous decreasing functions. This allows connections to be
established with classical results in CRNT that relate to the existence and
uniqueness of equilibria along positive stoichiometric compatibility classes.
In particular, monotonicity can be employed to identify regions in the set of
possible reaction rate coefficients leading to complex balancing, and to
conclude uniqueness of equilibria for a class of positive deficiency networks.
The latter result might support constructing an alternative proof of the
well-known deficiency one theorem. The developed notions and results are
illustrated through examples.


  This paper analyzes a simplified model of viral infection and evolution using
the 'grand canonical ensemble' and formalisms from statistical mechanics and
thermodynamics to enumerate all possible viruses and to derive thermodynamic
variables for the system. We model the infection process as a series of energy
barriers determined by the genetic states of the virus and host as a function
of immune response and system temperature. We find a phase transition between a
positive temperature regime of normal replication and a negative temperature
'disordered' phase of the virus. These phases define different regimes in which
different genetic strategies are favored. Perhaps most importantly, it
demonstrates that the system has a real thermodynamic temperature. For normal
replication, this temperature is linearly related to effective temperature. The
strength of immune response rescales temperature but does not change the
observed linear relationship. For all temperatures and immunities studied, we
find a universal curve relating the order parameter to viral evolvability. Real
viruses have finite length RNA segments that encode for proteins which
determine their fitness; hence the methods put forth here could be refined to
apply to real biological systems, perhaps providing insight into immune escape,
the emergence of novel pathogens and other results of viral evolution.


  The implementation of effective protected areas is one of the central goals
of modern conservation biology. In the context of fisheries management and
marine ecosystem conservation, marine reserves often play a significant role to
achieve sustainable fisheries management. Consequently, a substantial number of
studies have been conducted to establish broad rules for the creation of MPAs,
or to test the effects of MPAs in specific regions. However, there still exist
many challenges for implementing MPAs that are effective at meeting their
goals. Deducing theoretical conditions guaranteeing that the introduction of
marine reserves will increase fisheries yields in age-structured population
dynamics is one such challenge. To derive such conditions, a simple
mathematical model is developed that follows an age-structured metapopulation
dynamics of a sedentary species. The obtained results suggest that a
sufficiently high fishing mortality rate and moderate recruitment success of an
individual's eggs is a necessary for marine reserves to increase fisheries
yields. The numerical calculations were conducted with the parameters of red
abalone (Haliotis rufescens) to visualize and to check validity of the
analytical results. They show good agreement with the analytical results, as
well as the results obtained in the previous works.


  A core topic of research in prebiotic chemistry is the search for plausible
synthetic routes that connect the building blocks of modern life such as
sugars, nucleotides, amino acids, and lipids to "molecular food sources" that
have likely been abundant on Early Earth. In a recent contribution, Albert
Eschenmoser emphasised the importance of catalytic and autocatalytic cycles in
establishing such abiotic synthesis pathways. The accumulation of intermediate
products furthermore provides additional catalysts that allow pathways to
change over time. We show here that generative models of chemical spaces based
on graph grammars make it possible to study such phenomena is a systematic
manner. In addition to repro- ducing the key steps of Eschenmoser's hypothesis
paper, we discovered previously unexplored potentially autocatalytic pathways
from HCN to glyoxylate. A cascading of autocatalytic cycles could efficiently
re-route matter, distributed over the combinatorial complex network of HCN
hydrolysation chemistry, towards a potential primordial metabolism. The
generative approach also has it intrinsic limitations: the unsupervised
expansion of the chemical space remains infeasible due to the exponential
growth of possible molecules and reactions between them. Here in particular the
combinatorial complexity of the HCN polymerisation and hydrolysation networks
forms the computational bottleneck. As a consequence, guidance of the
computational exploration by chemical experience is indispensable.


  RNA polymerase (RNAP) is molecular machine that polymerizes a RNA molecule, a
linear heteropolymer, using a single stranded DNA (ssDNA) as the corresponding
template; the sequence of monomers of the RNA is dictated by that of monomers
on the ssDNA template. While polymerizing a RNA, the RNAP walks step-by-step on
the ssDNA template in a specific direction. Thus, a RNAP can be regarded also
as a molecular motor and the sites of start and stop of its walk on the DNA
mark the two ends of the genetic message that it transcribes into RNA.
Interference of transcription of two overlapping genes is believed to regulate
the levels of their expression, i.e., the overall rate of the corresponding RNA
synthesis, through suppressive effect of one on the other. Here we model this
process as a mixed traffic of two groups of RNAP motors that are characterized
by two distinct pairs of start and stop sites. Each group polymerizes identical
copies of a RNA while the RNAs polymerized by the two groups are different.
These models, which may also be viewed as two interfering totally asymmetric
simple exclusion processes, account for all modes of transcriptional
interference in spite of their extreme simplicity. A combination of mean-field
theory and computer simulation of these models demonstrate the physical origin
of the switch-like regulation of the two interfering genes in both
co-directional and contra-directional traffic of the two groups of RNAP motors.


  Individuals might abstain from participating in an instance of an
evolutionary game for various reasons, ranging from lack of interest to risk
aversion. In order to understand the consequences of such diverse activity
patterns on the evolution of cooperation, we study a weak prisoner's dilemma
where each player's participation is probabilistic rather than certain. Players
that do not participate get a null payoff and are unable to replicate. We show
that inactivity introduces cascading failures of cooperation, which are
particularly severe on scale-free networks with frequently inactive hubs. The
drops in the fraction of cooperators are sudden, while the spatiotemporal
reorganization of compact cooperative clusters, and thus the recovery, takes
time. Nevertheless, if the activity of players is directly proportional to
their degree, or if the interaction network is not strongly heterogeneous, the
overall evolution of cooperation is not impaired. This is because inactivity
negatively affects the potency of low-degree defectors, who are hence unable to
utilize on their inherent evolutionary advantage. Between cascading failures,
the fraction of cooperators is therefore higher than usual, which lastly
balances out the asymmetric dynamic instabilities that emerge due to
intermittent blackouts of cooperative hubs.


  Extending the previous 2-gender dioecious diploid gene-mating evolution model
[arXiv:1410.3456], we attempt to answer "whether the Hardy-Weinberg global
stability and the exact analytic dynamical solutions can be found in the
generalized N-gender N-polyploid gene-mating system with an arbitrary number of
alleles?" For a 2-gender gene-mating evolution model, a pair of male and female
determines the trait of their offspring. Each of the pair contributes one
inherited character, the allele, to combine into the genotype of their
offspring. Hence, for an N-gender N-polypoid gene-mating model, each of N
different genders contributes one allele to combine into the genotype of their
offspring. We exactly solve the analytic solution of N-gender-mating
$(n+1)$-alleles governing highly-nonlinear coupled differential equations in
the genotype frequency parameter space for any positive integer N and $n$. For
an analogy, the 2-gender to N-gender gene-mating equation generalization is
analogs to the 2-body collision to the N-body collision Boltzmann equations
with continuous distribution functions of "discretized" variables instead of
"continuous" variables. We find their globally stable solution as a continuous
manifold and find no chaos. Our solution implies that the Laws of Nature, under
our assumptions, provide no obstruction and no chaos to support an N-gender
gene-mating stable system.


  We introduce Phen-Gen, a method which combines patient disease symptoms and
sequencing data with prior domain knowledge to identify the causative gene(s)
for rare disorders.


  Motivation: Proteomic mass spectrometry analysis is becoming routine in
clinical diagnostics, for example to monitor cancer biomarkers using blood
samples. However, differential proteomics and identification of peaks relevant
for class separation remains challenging.
  Results: Here, we introduce a simple yet effective approach for identifying
differentially expressed proteins using binary discriminant analysis. This
approach works by data-adaptive thresholding of protein expression values and
subsequent ranking of the dichotomized features using a relative entropy
measure. Our framework may be viewed as a generalization of the `peak
probability contrast' approach of Tibshirani et al. (2004) and can be applied
both in the two-group and the multi-group setting.
  Our approach is computationally inexpensive and shows in the analysis of a
large-scale drug discovery test data set equivalent prediction accuracy as a
random forest. Furthermore, we were able to identify in the analysis of mass
spectrometry data from a pancreas cancer study biological relevant and
statistically predictive marker peaks unrecognized in the original study.
  Availability: The methodology for binary discriminant analysis is implemented
in the R package binda, which is freely available under the GNU General Public
License (version 3 or later) from CRAN at URL
http://cran.r-project.org/web/packages/binda/ . R scripts reproducing all
described analyzes are available from the web page
http://strimmerlab.org/software/binda/ .


  C4 plants, such as maize, concentrate carbon dioxide in a specialized
compartment surrounding the veins of their leaves to improve the efficiency of
carbon dioxide assimilation. Nonlinear relationships between carbon dioxide and
oxygen levels and reaction rates are key to their physiology but cannot be
handled with standard techniques of constraint-based metabolic modeling. We
demonstrate that incorporating these relationships as constraints on reaction
rates and solving the resulting nonlinear optimization problem yields realistic
predictions of the response of C4 systems to environmental and biochemical
perturbations. Using a new genome-scale reconstruction of maize metabolism, we
build an 18000-reaction, nonlinearly constrained model describing mesophyll and
bundle sheath cells in 15 segments of the developing maize leaf, interacting
via metabolite exchange, and use RNA-seq and enzyme activity measurements to
predict spatial variation in metabolic state by a novel method that optimizes
correlation between fluxes and expression data. Though such correlations are
known to be weak in general, here the predicted fluxes achieve high correlation
with the data, successfully capture the experimentally observed base-to-tip
transition between carbon-importing tissue and carbon-exporting tissue, and
include a nonzero growth rate, in contrast to prior results from similar
methods in other systems. We suggest that developmental gradients may be
particularly suited to the inference of metabolic fluxes from expression data.


  It was reported that buffalo is a low susceptibility species resisting to
prion diseases, which are invariably fatal and highly infectious
neurodegenerative diseases that affect a wide variety of species. In molecular
structures, TSE neurodegenerative diseases are caused by the conversion from a
soluble normal cellular prion protein, predominantly with alpha-helices, into
insoluble abnormally folded infectious prions, rich in beta-sheets. This paper
studies the molecular structure and structural dynamics of buffalo prion
protein, in order to reveal the reason why buffalo are resistant to prion
diseases. We first did molecular modeling of a homology structure constructed
by one mutation at residue 143 from the Nuclear Magnetic Resonance structure of
bovine and cattle PrP(124-227); immediately we found for buffalo PrPC(124-227)
there are 5 hydrogen bonds at Asn143, but at this position bovine/cattle do not
have such hydrogen bonds. Same as that of rabbits, dogs or horses, our
molecular dynamics studies also confirmed there is a strong salt bridge
ASP178-ARG164 (O-N) keeping the beta2-alpha2 loop linked in buffalo. We also
found there is a very strong hydrogen bond SER170-TYR218 linking this loop with
the C-terminal end of alpha-helix H3. Other information such as (i) there is a
very strong salt bridge HIS187-ARG156 (N-O) linking alpha-helices H2 and H1 (if
mutation H187R is made at position 187 then the hydrophobic core of PrPC will
be exposed), (ii) at D178, there is a hydrogen bond Y169-D178 and a polar
contact R164-D178 for BufPrPC instead of a polar contact Q168-D178 for bovine
PrPC, (iii) BufPrPC owns 3-10 helices at 125-127, 152-156 and in the
beta2-alpha2 loop respectively, and (iv) in beta2-alpha2 loop there are strong
pi-contacts, etc, has been discovered.


  We consider cyclic Lotka-Volterra models with three and four strategies where
at every interaction agents play a strategy using a time-dependent probability
distribution. Agents learn from a loss by reducing the probability to play a
losing strategy at the next interaction. For that, an agent is described as an
urn containing $\beta$ balls of three respectively four types where after a
loss one of the balls corresponding to the losing strategy is replaced by a
ball representing the winning strategy. Using both mean-field rate equations
and numerical simulations, we investigate a range of quantities that allow us
to characterize the properties of these cyclic models with time-dependent
probability distributions. For the three-strategy case in a spatial setting we
observe a transition from neutrally stable to stable when changing the level of
discretization of the probability distribution. For large values of $\beta$,
yielding a good approximation to a continuous distribution, spatially
synchronized temporal oscillations dominate the system. For the four-strategy
game the system is always neutrally stable, but different regimes emerge,
depending on the size of the system and the level of discretization.


  We provide, on an extensive dataset and using several different distances,
confirmation of the hypothesis that CGR patterns are preserved along a genomic
DNA sequence, and are different for DNA sequences originating from genomes of
different species. This finding lends support to the theory that CGRs of
genomic sequences can act as graphic genomic signatures. In particular, we
compare the CGR patterns of over five hundred different 150,000 bp genomic
sequences originating from the genomes of six organisms, each belonging to one
of the kingdoms of life: H. sapiens, S. cerevisiae, A. thaliana, P. falciparum,
E. coli, and P. furiosus. We also provide preliminary evidence of this method's
applicability to closely related species by comparing H. sapiens (chromosome
21) sequences and over one hundred and fifty genomic sequences, also 150,000 bp
long, from P. troglodytes (Animalia; chromosome Y), for a total length of more
than 101 million basepairs analyzed. We compute pairwise distances between CGRs
of these genomic sequences using six different distances, and construct
Molecular Distance Maps that visualize all sequences as points in a
two-dimensional or three-dimensional space, to simultaneously display their
interrelationships. Our analysis confirms that CGR patterns of DNA sequences
from the same genome are in general quantitatively similar, while being
different for DNA sequences from genomes of different species. Our analysis of
the performance of the assessed distances uses three different quality measures
and suggests that several distances outperform the Euclidean distance, which
has so far been almost exclusively used for such studies. In particular we show
that, for this dataset, DSSIM (Structural Dissimilarity Index) and the
descriptor distance (introduced here) are best able to classify genomic
sequences.


  One of the new discoveries in quantum biology is the role of Environment
Assisted Quantum Transport (ENAQT) in excitonic transport processes. In
disordered quantum systems transport is most efficient when the environment
just destroys quantum interferences responsible for localization, but the
coupling does not drive the system to fully classical thermal diffusion yet.
This poised realm between the pure quantum and the semi-classical domains has
not been considered in other biological transport processes, such as charge
transport through organic molecules. Binding in receptor-ligand complexes is
assumed to be static as electrons are assumed to be not able to cross the
ligand molecule. We show that ENAQT makes cross ligand transport possible and
efficient between certain atoms opening the way for the reorganization of the
charge distribution on the receptor when the ligand molecule docks. This new
effect can potentially change our understanding how receptors work. We
demonstrate room temperature ENAQT on the caffeine molecule.


  We propose a model describing the synergetic feedback between word-of-mouth
(WoM) and epidemic dynamics controlled by voluntary vaccination. We combine a
game-theoretic model for the spread of WoM and a compartmental model describing
$SIR$ disease dynamics in the presence of a program of voluntary vaccination.
We evaluate and compare two scenarios, depending on what WoM disseminates: (1)
vaccine advertising, which may occur whether or not an epidemic is ongoing and
(2) epidemic status, notably disease prevalence. Understanding the synergy
between the two strategies could be particularly important for organizing
voluntary vaccination campaigns. We find that, in the initial phase of an
epidemic, vaccination uptake is determined more by vaccine advertising than the
epidemic status. As the epidemic progresses, epidemic status become
increasingly important for vaccination uptake, considerably accelerating
vaccination uptake toward a stable vaccination coverage.


  We study the dynamics of colonization of a territory by a stochastic
population at low immigration pressure. We assume a sufficiently strong Allee
effect that introduces, in deterministic theory, a large critical population
size for colonization. At low immigration rates, the average pre-colonization
population size is small thus invalidating the WKB approximation to the master
equation. We circumvent this difficulty by deriving an exact zero-flux solution
of the master equation and matching it with an approximate non-zero-flux
solution of the pertinent Fokker-Planck equation in a small region around the
critical population size. This procedure provides an accurate evaluation of the
quasi-stationary probability distribution of population sizes in the
pre-colonization state, and of the mean time to colonization, for a wide range
of immigration rates. At sufficiently high immigration rates our results agree
with WKB results obtained previously. At low immigration rates the results can
be very different.


  Within the field of phylogenetics there is growing interest in measures for
summarising the dissimilarity, or 'incongruence', of two or more phylogenetic
trees. Many of these measures are NP-hard to compute and this has stimulated a
considerable volume of research into fixed parameter tractable algorithms. In
this article we use Monadic Second Order logic (MSOL) to give alternative,
compact proofs of fixed parameter tractability for several well-known
incongruency measures. In doing so we wish to demonstrate the considerable
potential of MSOL - machinery still largely unknown outside the algorithmic
graph theory community - within phylogenetics. A crucial component of this work
is the observation that many of these measures, when bounded, imply the
existence of an 'agreement forest' of bounded size, which in turn implies that
an auxiliary graph structure, the display graph, has bounded treewidth. It is
this bound on treewidth that makes the machinery of MSOL available for proving
fixed parameter tractability. We give a variety of different MSOL formulations.
Some are based on explicitly encoding agreement forests, while some only use
them implicitly to generate the treewidth bound. Our formulations introduce a
number of "phylogenetics MSOL primitives" which will hopefully be of use to
other researchers.


  Cancer arises from successive rounds of mutations which generate tumor cells
with different genomic variation i.e. clones. For drug responsiveness and
therapeutics, it is necessary to identify the clones in tumor sample
accurately. Many methods are developed to infer tumor heterogeneity by either
computing cellular prevalence and tumor phylogeny or predicting genotype of
mutations. All methods suffer some problems e.g. inaccurate computation of
clonal frequencies, discarding clone specific genotypes etc. In the paper, we
propose a method, called- HetFHMM to infer tumor heterogeneity by predicting
clone specific genotypes and cellular prevalence. To infer clone specific
genotype, we consider the presence of multiple mutations at any genomic
location. We also tested our model on different simulated data. The results
shows that HetFHMM outperforms recent methods which infer tumor heterogeneity.
Therefore, HetFHMM is a novel approach in tumor heterogeneity research area.


  Populations of species in ecosystems are often constrained by availability of
resources within their environment. In effect this means that a growth of one
population, needs to be balanced by comparable reduction in populations of
others. In neutral models of biodiversity all populations are assumed to change
incrementally due to stochastic births and deaths of individuals. Here we
propose and model another redistribution mechanism driven by abrupt and severe
collapses of the entire population of a single species freeing up resources for
the remaining ones. This mechanism may be relevant e.g. for communities of
bacteria, with strain-specific collapses caused e.g. by invading
bacteriophages, or for other ecosystems where infectious diseases play an
important role.
  The emergent dynamics of our system is cyclic "diversity waves" triggered by
collapses of globally dominating populations. The population diversity peaks at
the beginning of each wave and exponentially decreases afterwards. Species
abundances are characterized by a bimodal time-aggregated distribution with the
lower peak formed by populations of recently collapsed or newly introduced
species, while the upper peak - species that has not yet collapsed in the
current wave. In most waves both upper and lower peaks are composed of several
smaller peaks. This self-organized hierarchical peak structure has a long-term
memory transmitted across several waves. It gives rise to a scale-free tail of
the time-aggregated population distribution with a universal exponent of 1.7.
We show that diversity wave dynamics is robust with respect to variations in
the rules of our model such as diffusion between multiple environments,
species-specific growth and extinction rates, and bet-hedging strategies.


  Enzymatic molecules that actively support many cellular processes, including
transport, cell division and cell motility, are known as motor proteins or
molecular motors. Experimental studies indicate that they interact with each
other and they frequently work together in large groups. To understand the
mechanisms of collective behavior of motor proteins we study the effect of
interactions in the transport of molecular motors along linear filaments. It is
done by analyzing a recently introduced class of totally asymmetric exclusion
processes that takes into account the intermolecular interactions via
thermodynamically consistent approach. We develop a new theoretical method that
allows us to compute analytically all dynamic properties of the system. Our
analysis shows that correlations play important role in dynamics of interacting
molecular motors. Surprisingly, we find that the correlations for repulsive
interactions are weaker and more short-range than the correlations for the
attractive interactions. In addition, it is shown that symmetry of interactions
affect dynamic properties of molecular motors. The implications of these
findings for motor proteins transport are discussed. Our theoretical
predictions are tested by extensive Monte Carlo computer simulations.


  We show that for any two values $\alpha, \beta >0 $ for which
$\alpha+\beta>1$ then there is a value $N$ so that for all $n \geq N$ the
following holds. For any binary phylogenetic tree $T$ on $n$ leaves there is a
set of $\lfloor n^\alpha \rfloor$ characters that capture $T$, and for which
each character takes at most $\lfloor n^\beta \rfloor$ distinct states. Here
`capture' means that $T$ is the unique perfect phylogeny for these characters.
Our short proof of this combinatorial result is based on the probabilistic
method.


  The functional relevance of regulating proteins is often limited to specific
binding sites such as the ends of microtubules or actin-filaments. A
localization of proteins on these functional sites is of great importance. We
present a quantitative theory for a diffusion and capture process, where
proteins diffuse on a filament and stop diffusing when reaching the filament's
end. It is found that end-association after one-dimensional diffusion is the
main source for tip-localization of such proteins. As a consequence, diffusion
and capture is highly efficient in enhancing the reaction velocity of enzymatic
reactions, where proteins and filament ends are to each other as enzyme and
substrate. We show that the reaction velocity can effectively be described
within a Michaelis-Menten framework. Together one-dimensional diffusion and
capture beats the (three-dimensional) Smoluchowski diffusion limit for the rate
of protein association to filament ends.


  Living cells use readout molecules to record the state of receptor proteins,
similar to measurements or copies in typical computational devices. But is this
analogy rigorous? Can cells be optimally efficient, and if not, why? We show
that, as in computation, a canonical biochemical readout network generates
correlations; extracting no work from these correlations sets a lower bound on
dissipation. For general input, the biochemical network cannot reach this
bound, even with arbitrarily slow reactions or weak thermodynamic driving. It
faces an accuracy-dissipation trade-off that is qualitatively distinct from and
worse than implied by the bound, and more complex steady-state copy processes
cannot perform better. Nonetheless, the cost remains close to the thermodynamic
bound unless accuracy is extremely high. Additionally, we show that
biomolecular reactions could be used in thermodynamically optimal devices under
exogenous manipulation of chemical fuels, suggesting an experimental system for
testing computational thermodynamics.


  Spatial heterogeneity plays an important role in complex ecosystem dynamics,
and therefore is also an important consideration in sustainable resource
management. However, little is known about how spatial effects can influence
management targets derived from a non-spatial harvest model. Here, we extended
the Schaefer model, a conventional non-spatial harvest model that is widely
used in resource management, to a spatially-explicit harvest model by
integrating environmental heterogeneities, as well as species exchange between
patches. By comparing the maximum sustainable yields (MSY), one of the central
management targets in resource management, obtained from the spatially extended
model with that of the conventional model, we examined the effect of spatial
heterogeneity. When spatial heterogeneity exists, we found that the Schaefer
model tends to overestimate the MSY, implying potential for causing
overharvesting. In addition, by assuming a well-mixed population in the
heterogeneous environment, we showed analytically that the Schaefer model
always overestimate the MSY, regardless of the number of patches existing. The
degree of overestimation becomes significant when spatial heterogeneity is
marked. Collectively, these results highlight the importance of integrating the
spatial structure to conduct sustainable resource management.


  Genes with similar transcriptional activation kinetics can display very
different temporal mRNA profiles due to differences in transcription time,
degradation rate and RNA processing kinetics. Recent studies have shown that a
splicing-associated RNA production delay can be significant. We introduce a
joint model of transcriptional activation and mRNA accumulation which can be
used for inference of transcription rate, RNA production delay and degradation
rate given genome-wide data from high-throughput sequencing time course
experiments. We combine a mechanistic differential equation model with a
non-parametric statistical modelling approach allowing us to capture a broad
range of activation kinetics, and use Bayesian parameter estimation to quantify
the uncertainty in the estimates of the kinetic parameters. We apply the model
to data from estrogen receptor (ER-{\alpha}) activation in the MCF-7 breast
cancer cell line. We use RNA polymerase II (pol-II) ChIP-Seq time course data
to characterise transcriptional activation and mRNA-Seq time course data to
quantify mature transcripts. We find that 11% of genes with a good signal in
the data display a delay of more than 20 minutes between completing
transcription and mature mRNA production. The genes displaying these long
delays are significantly more likely to be short. We also find a statistical
association between high delay and late intron retention in pre-mRNA data,
indicating significant splicing-associated production delays in many genes.


  In this paper we provide the derivation of a super compact pairwise model
with only 4 equations in the context of describing
susceptible-infected-susceptible (SIS) epidemic dynamics on heterogenous
networks. The super compact model is based on a new closure relation that
involves not only the average degree but also the second and third moments of
the degree distribution. Its derivation uses an a priori approximation of the
degree distribution of susceptible nodes in terms of the degree distribution of
the network. The new closure gives excellent agreement with heterogeneous
pairwise models that contain significantly more differential equations.


  This paper presents libRoadRunner, an extensible, high-performance,
cross-platform, open-source software library for the simulation and analysis of
models \ expressed using Systems Biology Markup Language (SBML). SBML is the
most widely used standard for representing dynamic networks, especially
biochemical networks. libRoadRunner supports solution of both large models and
multiple replicas of a single model on desktop, mobile and cluster computers.
libRoadRunner is a self-contained library, able to run both as a component
inside other tools via its C++ and C bindings andnteractively through its
Python interface. The Python Application Programming Interface (API) is similar
to the APIs of Matlab and SciPy, making it fast and easy to learn, even for new
users. libRoadRunner uses a custom Just-In-Time (JIT) compiler built on the
widely-used LLVM JIT compiler framework to compile SBML-specified models
directly into very fast native machine code for a variety of processors, making
it appropriate for solving very large models or multiple replicas of smaller
models. libRoadRunner is flexible, supporting the bulk of the SBML
specification (except for delay and nonlinear algebraic equations) and several
of its extensions. It offers multiple deterministic and stochastic integrators,
as well as tools for steady-state, stability analyses and flux balance
analysis. We regularly update libRoadRunner binary distributions for Mac OS X,
Linux and Windows and license them under Apache License Version 2.0.
http://www.libroadrunner.org provides online documentation, full build
instructions, binaries and a git source repository.


  We use a popular fictional disease, zombies, in order to introduce techniques
used in modern epidemiology modelling, and ideas and techniques used in the
numerical study of critical phenomena. We consider variants of zombie models,
from fully connected continuous time dynamics to a full scale exact stochastic
dynamic simulation of a zombie outbreak on the continental United States. Along
the way, we offer a closed form analytical expression for the fully connected
differential equation, and demonstrate that the single person per site two
dimensional square lattice version of zombies lies in the percolation
universality class. We end with a quantitative study of the full scale US
outbreak, including the average susceptibility of different geographical
regions.


  A wide range of studies in population genetics have employed the sample
frequency spectrum (SFS), a summary statistic which describes the distribution
of mutant alleles at a polymorphic site in a sample of DNA sequences. In
particular, recently there has been growing interest in analyzing the joint SFS
data from multiple populations to infer parameters of complex demographic
histories, including variable population sizes, population split times,
migration rates, admixture proportions, and so on. Although much methodological
progress has been made, existing SFS-based inference methods suffer from
numerical instability and high computational complexity when multiple
populations are involved and the sample size is large. In this paper, we
present new analytic formulas and algorithms that enable efficient computation
of the expected joint SFS for multiple populations related by a complex
demographic model with arbitrary population size histories (including piecewise
exponential growth). Our results are implemented in a new software package
called momi (MOran Models for Inference). Through an empirical study involving
tens of populations, we demonstrate our improvements to numerical stability and
computational complexity.


  The process of genetic recombination can be seen as a chemical reaction
network with mass-action kinetics. We review the known results on existence,
uniqueness, and global stability of an equilibrium in every compatibility class
and for all rate constants, from both the population genetics and the reaction
networks point of view.


  Any macroscopic deformation of a filamentous bundle is necessarily
accompanied by local sliding and/or stretching of the constituent filaments.
Yet the nature of the sliding friction between two aligned filaments
interacting through multiple contacts remains largely unexplored. Here, by
directly measuring the sliding forces between two bundled F-actin filaments, we
show that these frictional forces are unexpectedly large, scale logarithmically
with sliding velocity as in solid-like friction, and exhibit complex dependence
on the filaments' overlap length. We also show that a reduction of the
frictional force by orders of magnitude, associated with a transition from
solid-like friction to Stokes' drag, can be induced by coating F-actin with
polymeric brushes. Furthermore, we observe similar transitions in filamentous
microtubules and bacterial flagella. Our findings demonstrate how altering a
filament's elasticity, structure and interactions can be used to engineer
interfilament friction and thus tune the properties of fibrous composite
materials.


  In this article we use all-atom simulations to elucidate the mechanisms
underlying conformational switching and allostery within the coat protein of
the bacteriophage MS2. Assembly of most icosahedral virus capsids requires that
the capsid protein adopt different conformations at precise locations within
the capsid. It has been shown that a 19 nucleotide stem loop (TR) from the MS2
genome acts as an allosteric effector, guiding conformational switching of the
coat protein during capsid assembly. Since the principal conformational changes
occur far from the TR binding site, it is important to understand the molecular
mechanism underlying this allosteric communication. To this end, we use
all-atom simulations with explicit water combined with a path sampling
technique to sample the MS2 coat protein conformational transition, in the
presence and absence of TR-binding. The calculations find that TR binding
strongly alters the transition free energy profile, leading to a switch in the
favored conformation. We discuss changes in molecular interactions responsible
for this shift. We then identify networks of amino acids with correlated
motions to reveal the mechanism by which effects of TR binding span the
protein. The analysis predicts amino acids whose substitution by mutagenesis
could alter populations of the conformational substates or their transition
rates.


  In a diffusion-based molecular communication network, transmitters and
receivers communicate by using signalling molecules (or ligands) in a fluid
medium. This paper assumes that the transmitter uses different chemical
reactions to generate different emission patterns of signalling molecules to
represent different transmission symbols, and the receiver consists of
receptors. When the signalling molecules arrive at the receiver, they may react
with the receptors to form ligand-receptor complexes. Our goal is to study the
demodulation in this setup assuming that the transmitter and receiver are
synchronised. We derive an optimal demodulator using the continuous history of
the number of complexes at the receiver as the input to the demodulator. We do
that by first deriving a communication model which includes the chemical
reactions in the transmitter, diffusion in the transmission medium and the
ligand-receptor process in the receiver. This model, which takes the form of a
continuous-time Markov process, captures the noise in the receiver signal due
to the stochastic nature of chemical reactions and diffusion. We then adopt a
maximum a posterior framework and use Bayesian filtering to derive the optimal
demodulator. We use numerical examples to illustrate the properties of this
optimal demodulator.


  We discuss a method of approximate model reduction for networks of
biochemical reactions. This method can be applied to networks with polynomial
or rational reaction rates and whose parameters are given by their orders of
magnitude. In order to obtain reduced models we solve the problem of tropical
equilibration that is a system of equations in max-plus algebra. In the case of
networks with nonlinear fast cycles we have to solve the problem of tropical
equilibration at least twice, once for the initial system and a second time for
an extended system obtained by adding to the initial system the differential
equations satisfied by the conservation laws of the fast subsystem. The two
steps can be reiterated until the fast subsystem has no conservation laws
different from the ones of the full model. Our method can be used for formal
model reduction in computational systems biology.


  Everyday we are exposed to various chemicals via food additives, cleaning and
cosmetic products and medicines -- and some of them might be toxic. However
testing the toxicity of all existing compounds by biological experiments is
neither financially nor logistically feasible. Therefore the government
agencies NIH, EPA and FDA launched the Tox21 Data Challenge within the
"Toxicology in the 21st Century" (Tox21) initiative. The goal of this challenge
was to assess the performance of computational methods in predicting the
toxicity of chemical compounds. State of the art toxicity prediction methods
build upon specifically-designed chemical descriptors developed over decades.
Though Deep Learning is new to the field and was never applied to toxicity
prediction before, it clearly outperformed all other participating methods. In
this application paper we show that deep nets automatically learn features
resembling well-established toxicophores. In total, our Deep Learning approach
won both of the panel-challenges (nuclear receptors and stress response) as
well as the overall Grand Challenge, and thereby sets a new standard in tox
prediction.


  In cells, most of cargos are transported by motor proteins along microtubule.
Biophysically, unidirectional motion of large number of motor proteins along a
single track can be described by totally asymmetric simple exclusion process
(TASEP). From which many meaningful properties, such as the appearance of
domain wall (defined as the borderline of high density and low density of motor
protein along motion track) and boundary layers, can be obtained. However, it
is biologically obvious that a single track may be occupied by different motor
species. So previous studies based on TASEP of one particle species are not
reasonable enough to find more detailed properties of the motion of motors
along a single track. To address this problem, TASEP with two particle species
is discussed in this study. Theoretical methods to get densities of each
particle species are provided. Using these methods, phase transition related
properties of particle densities are obtained. Our analysis show that domain
wall and boundary layer of single species densities always appear
simultaneously with those of the total particle density. The height of domain
wall of total particle density is equal to the summation of those of single
species. Phase diagrams for typical model parameters are also presented. The
methods presented in this study can be generalized to analyze TASEP with more
particle species.


  Anaerobic digestion enables the water industry to treat wastewater as a
resource for generating energy and recovering valuable by-products. The
complexity of the anaerobic digestion process has motivated the development of
complex models. However, this complexity makes it intractable to pin-point
stability and emergent behaviour. Here, the widely used Anaerobic Digestion
Model No. 1 (ADM1) has been reduced to its very backbone, a syntrophic
two-tiered microbial food chain and a slightly more complex three-tiered
microbial food web, with their stability analysed as function of the inflowing
substrate concentration and dilution rate. Parameterised for phenol and
chlorophenol degradation, steady-states were always stable and non-oscillatory.
Low input concentrations of chlorophenol were sufficient to maintain
chlorophenol- and phenol-degrading populations but resulted in poor conversion
and a hydrogen flux that was too low to sustain hydrogenotrophic methanogens.
The addition of hydrogen and phenol boosted the populations of all three
organisms, resulting in the counterintuitive phenomena that (i) the phenol
degraders were stimulated by adding hydrogen, even though hydrogen inhibits
phenol degradation, and (ii) the dechlorinators indirectly benefitted from
measures that stimulated their hydrogenotrophic competitors; both phenomena
hint at emergent behaviour.


  The fall of prices of the high-throughput genome sequencing changes the
landscape of modern genomics. A number of large scale projects aimed at
sequencing many human genomes are in progress. Genome sequencing also becomes
an important aid in the personalized medicine. One of the significant side
effects of this change is a necessity of storage and transfer of huge amounts
of genomic data. In this paper we deal with the problem of compression of large
collections of complete genomic sequences. We propose an algorithm that is able
to compress the collection of 1092 human diploid genomes about 9,500 times.
This result is about 4 times better than what is offered by the other existing
compressors. Moreover, our algorithm is very fast as it processes the data with
speed 200MB/s on a modern workstation. In a consequence the proposed algorithm
allows storing the complete genomic collections at low cost, e.g., the examined
collection of 1092 human genomes needs only about 700MB when compressed, what
can be compared to about 6.7 TB of uncompressed FASTA files. The source code is
available at
http://sun.aei.polsl.pl/REFRESH/index.php?page=projects&project=gdc&subpage=about.


  The bond graph approach to modelling biochemical networks is extended to
allow hierarchical construction of complex models from simpler components. This
is made possible by representing the simpler components as thermodynamically
open systems exchanging mass and energy via ports. A key feature of this
approach is that the resultant models are robustly thermodynamically compliant:
the thermodynamic compliance is not dependent on precise numerical values of
parameters. Moreover, the models are reusable due to the well-defined interface
provided by the energy ports.
  To extract bond graph model parameters from parameters found in the
literature, general and compact formulae are developed to relate free-energy
constants and equilibrium constants. The existence and uniqueness of solutions
is considered in terms of fundamental properties of stoichiometric matrices.
The approach is illustrated by building a hierarchical bond graph model of
glycogenolysis in skeletal muscle.


  Discrete-time models are the traditional approach for capturing population
dynamics of a host-parasitoid system. Recent work has introduced a
semi-discrete framework for obtaining model update functions that connect
host-parasitoid population levels from year-to-year. In particular, this
framework uses differential equations to describe the hosts-parasitoid
interaction during the time of year where they come in contact, allowing
specific behaviors to be mechanistically incorporated into the model. We use
the semi-discrete approach to study the effects of host-feeding, which occurs
when a parasitoid consumes a potential host larva without ovipositing. Our
results show that host-feeding by itself cannot stabilize the system, and both
the host and parasitoid populations exhibit diverging oscillations similar to
the Nicholson-Bailey model. However, when combined with other stabilizing
mechanisms such as density-dependent host mortality or density-dependent
parasitoid attack rate, host-feeding expands the region of parameter space that
allows for a stable host-parasitoid equilibrium. Finally, our results show that
host-feeding causes inefficiency in the parasitoid population, which yields a
higher population of hosts per generation. This suggests that host-feeding may
have limited long-term impact in terms of suppressing host levels for
biological control applications.


  Temperature compensation is a notable property of circadian oscillators that
indicates the insensitivity of the oscillator system's period to temperature
changes; the underlying mechanism, however, is still unclear. We investigated
the influence of protein dimerization and cooperative stability in protein
degradation on the temperature compensation ability of two oscillators. Here,
cooperative stability means that high-order oligomers are more stable than
their monomeric counterparts. The period of an oscillator is affected by the
parameters of the dynamic system, which in turn are influenced by temperature.
We adopted the Repressilator and the Atkinson oscillator to analyze the
temperature sensitivity of their periods. Phase sensitivity analysis was
employed to evaluate the period variations of different models induced by
perturbations to the parameters. Furthermore, we used experimental data
provided by other studies to determine the reasonable range of parameter
temperature sensitivity. We then applied the linear programming method to the
oscillatory systems to analyze the effects of protein dimerization and
cooperative stability on the temperature sensitivity of their periods, which
reflects the ability of temperature compensation in circadian rhythms. Our
study explains the temperature compensation mechanism for circadian clocks.
Compared with the no-dimer mathematical model and linear model for protein
degradation, our theoretical results show that the nonlinear protein
degradation caused by cooperative stability is more beneficial for realizing
temperature compensation of the circadian clock.


  The origin and diversification of Sino-Tibetan populations have been a
long-standing hot debate. However, the limited genetic information of Tibetan
populations keeps this topic far from clear. In the present study, we genotyped
15 forensic autosomal STRs from 803 unrelated Tibetan individuals from Gansu
Province (635 from Gannan and 168 from Tianzhu). We combined these data with
published dataset to infer a detailed population affinities and admixture of
Sino-Tibetan populations. Our results revealed that the genetic structure of
Sino-Tibetan populations was strongly correlated with linguistic affiliations.
Although the among-population variances are relatively small, the genetic
components for Tibetan, Lolo-Burmese, and Han Chinese were quite distinctive,
especially for the Deng, Nu, and Derung of Lolo-Burmese. Southern indigenous
populations, such as Tai-Kadai and Hmong-Mien populations might have made
substantial genetic contribution to Han Chinese and Altaic populations, but not
to Tibetans. Likewise, Han Chinese but not Tibetan shared very similar genetic
makeups with Altaic populations, which did not support the North Asian origin
of Tibetan populations. The dataset generated here are also valuable for
forensic identification.


  Bovine TB is a major problem for the agricultural industry in several
countries. TB can be contracted and spread by species other than cattle and
this can cause a problem for disease control. In the UK and Ireland, badgers
are a recognised reservoir of infection and there has been substantial
discussion about potential control strategies. We present a coupling of
individual based models of bovine TB in badgers and cattle, which aims to
capture the key details of the natural history of the disease and of both
species at approximately county scale. The model is spatially explicit it
follows a very large number of cattle and badgers on a different grid size for
each species and includes also winter housing. We show that the model can
replicate the reported dynamics of both cattle and badger populations as well
as the increasing prevalence of the disease in cattle. Parameter space used as
input in simulations was swept out using Latin hypercube sampling and
sensitivity analysis to model outputs was conducted using mixed effect models.
By exploring a large and computationally intensive parameter space we show that
of the available control strategies it is the frequency of TB testing and
whether or not winter housing is practised that have the most significant
effects on the number of infected cattle, with the effect of winter housing
becoming stronger as farm size increases. Whether badgers were culled or not
explained about 5%, while the accuracy of the test employed to detect infected
cattle explained less than 3% of the variance in the number of infected cattle.


  Much of the research on the behavior of the SIS model on networks has
concerned the infinite size limit; in particular the phase transition between a
state where outbreaks can reach a finite fraction of the population, and a
state where only a finite number would be infected. For finite networks, there
is also a dynamic transition---the immortality transition---when the
per-contact transmission probability $\lambda$ reaches one. If $\lambda < 1$,
the probability that an outbreak will survive by an observation time $t$ tends
to zero as $t \rightarrow \infty$; if $\lambda = 1$, this probability is one.
We show that treating $\lambda = 1$ as a critical point predicts the
$\lambda$-dependence of the survival probability also for more moderate
$\lambda$-values. The exponent, however, depends on the underlying network.
This fact could, by measuring how a vertex' deletion changes the exponent, be
used to evaluate the role of a vertex in the outbreak. Our work also confirms
an extremely clear separation between the early die-off (from the outbreak
failing to take hold in the population) and the later extinctions
(corresponding to rare stochastic events of several consecutive transmission
events failing to occur).


  Adaptation of asexual populations is driven by beneficial mutations and
therefore the dynamics of this process, besides other factors, depend on the
distribution of beneficial fitness effects. It is known that on uncorrelated
fitness landscapes, this distribution can only be of three types: truncated,
exponential and power law. We performed extensive stochastic simulations to
study the adaptation dynamics on rugged fitness landscapes, and identified two
quantities that can be used to distinguish the underlying distribution of
beneficial fitness effects. The first quantity studied here is the fitness
difference between successive mutations that spread in the population, which is
found to decrease in the case of truncated distributions, remain nearly a
constant for exponentially decaying distributions and increase when the fitness
distribution decays as a power law. The second quantity of interest, namely,
the rate of change of fitness with time also shows quantitatively different
behaviour for different beneficial fitness distributions. The patterns
displayed by the two aforementioned quantities are found to hold for both low
and high mutation rates. We discuss how these patterns can be exploited to
determine the distribution of beneficial fitness effects in microbial
experiments.


  Bistability is considered wide-spread among bacteria and eukaryotic cells,
useful e.g. for enzyme induction, bet hedging, and epigenetic switching.
However, this phenomenon has mostly been described with deterministic dynamic
or well-mixed stochastic models. Here, we map known biological bistable systems
onto the well-characterized biochemical Schloegl model, using analytical
calculations and stochastic spatio-temporal simulations. In addition to network
architecture and strong thermodynamic driving away from equilibrium, we show
that bistability requires fine-tuning towards small cell volumes (or
compartments) and fast protein diffusion (well mixing). Bistability is thus
fragile and hence may be restricted to small bacteria and eukaryotic nuclei,
with switching triggered by volume changes during the cell cycle. For large
volumes, single cells generally loose their ability for bistable switching and
instead undergo a first-order phase transition.


  Protein-peptide interactions play a key role in cell functions. Their
structural characterization, though challenging, is important for the discovery
of new drugs. The CABS-dock web server provides an interface for modeling
protein-peptide interactions using a highly efficient protocol for the flexible
docking of peptides to proteins. While other docking algorithms require
pre-defined localization of the binding site, CABS-dock doesn't require such
knowledge. Given a protein receptor structure and a peptide sequence (and
starting from random conformations and positions of the peptide), CABS-dock
performs simulation search for the binding site allowing for full flexibility
of the peptide and small fluctuations of the receptor backbone. This protocol
was extensively tested over the largest dataset of non-redundant
protein-peptide interactions available to date (including bound and unbound
docking cases). For over 80% of bound and unbound data set cases, we obtained
models with high or medium accuracy (sufficient for practical applications).
Additionally, as optional features, CABS-dock can exclude user-selected binding
modes from docking search or to increase the level of flexibility for chosen
receptor fragments. CABS-dock is freely available as a web server at
http://biocomp.chem.uw.edu.pl/CABSdock


  The genetic structure of human populations is extraordinarily complex and of
fundamental importance to studies of anthropology, evolution, and medicine.
  As increasingly many individuals are of mixed origin, there is an unmet need
for tools that can infer multiple origins. Misclassification of such
individuals can lead to incorrect and costly misinterpretations of genomic
data, primarily in disease studies and drug trials.
  We present an advanced tool to infer ancestry that can identify the
biogeographic origins of highly mixed individuals.
  reAdmix is an online tool available at
http://chcb.saban-chla.usc.edu/reAdmix/.


  In this paper we present a generative model for protein contact networks. The
soundness of the proposed model is investigated by focusing primarily on
mesoscopic properties elaborated from the spectra of the graph Laplacian. To
complement the analysis, we study also classical topological descriptors, such
as statistics of the shortest paths and the important feature of modularity.
Our experiments show that the proposed model results in a considerable
improvement with respect to two suitably chosen generative mechanisms,
mimicking with better approximation real protein contact networks in terms of
diffusion properties elaborated from the Laplacian spectra. However, as well as
the other considered models, it does not reproduce with sufficient accuracy the
shortest paths structure. To compensate this drawback, we designed a second
step involving a targeted edge reconfiguration process. The ensemble of
reconfigured networks denotes improvements that are statistically significant.
As a byproduct of our study, we demonstrate that modularity, a well-known
property of proteins, does not entirely explain the actual network architecture
characterizing protein contact networks. In fact, we conclude that modularity,
intended as a quantification of an underlying community structure, should be
considered as an emergent property of the structural organization of proteins.
Interestingly, such a property is suitably optimized in protein contact
networks together with the feature of path efficiency.


  We are interested in prey-predator communities where the predator population
evolves much faster than the prey's (e.g. insect-tree communities). We
introduce a piecewise deterministic model for these prey-predator communities
that arises as a limit of a microscopic model when the number of predators goes
to infinity. We prove that the process has a unique invariant probability
measure and that it is exponentially ergodic. Further on, we rescale the
predator dynamics in order to model predators of smaller size. This slow-fast
system converges to a community process in which the prey dynamics is averaged
on the predator equilibria. This averaged process admits an invariant
probability measure which can be computed explicitly. We use numerical
simulations to study the convergence of the invariant probability measures of
the rescaled processes.


  It is well known that life on Earth alters its environment over evolutionary
and geological timescales. An important open question is whether this is a
result of evolutionary optimization or a universal feature of life. In the
latter case, the origin of life would be coincident with a shift in
environmental conditions. Here we present a model for the emergence of life in
which replicators are explicitly coupled to their environment through the
recycling of a finite supply of resources. The model exhibits a dynamic,
first-order phase transition from non-life to "life," where the life phase is
distinguished by selection on replicators. We show that environmental coupling
plays an important role in the dynamics of the transition. The transition
corresponds to a redistribution of matter in replicators and their environment,
driven by selection on replicators, exhibiting an explosive growth in diversity
as replicators are selected. The transition is accurately tracked by the mutual
information shared between replicators and their environment. In the absence of
successfully repartitioning system resources, the transition fails to complete,
leading to the possibility of many frustrated trials before life first emerges.
Often, the replicators that initiate the transition are not those that are
ultimately selected. The results are consistent with the view that life's
propensity to shape its environment is indeed a universal feature of
replicators, characteristic of the transition from non-life to life. We discuss
the implications of these results for understanding life's emergence and
evolutionary transitions more broadly.


  For any spatially explicit management, determining the appropriate spatial
scale of management decisions is critical to success at achieving a given
management goal. Specifically, managers must decide how much to subdivide a
given managed region: from implementing a uniform approach across the region to
considering a unique approach in each of one hundred patches and everything in
between. Spatially explicit approaches, such as the implementation of marine
spatial planning and marine reserves, are increasingly used in fishery
management. Using a spatially explicit bioeconomic model, we quantify how the
management scale affects optimal fishery profit, biomass, fishery effort, and
the fraction of habitat in marine reserves. We find that, if habitats are
randomly distributed, the fishery profit increases almost linearly with the
number of segments. However, if habitats are positively autocorrelated, then
the fishery profit increases with diminishing returns. Therefore, the true
optimum in management scale given cost to subdivision depends on the habitat
distribution pattern.


  Modeling social-ecological systems is difficult due to the complexity of
ecosystems and of individual and collective human behavior. Key components of
the social-ecological system are often over-simplified or omitted. Generalized
modeling is a dynamical systems approach that can overcome some of these
challenges. It can rigorously analyze qualitative system dynamics such as
regime shifts despite incomplete knowledge of the model's constituent
processes. Here, we review generalized modeling and use a recent study on the
Baltic Sea cod fishery's boom and collapse to demonstrate its application to
modeling the dynamics of empirical social-ecological systems. These empirical
applications demand new methods of analysis suited to larger, more complicated
generalized models. Generalized modeling is a promising tool for rapidly
developing mathematically rigorous, process-based understanding of a
social-ecological system's dynamics despite limited knowledge of the system.


  Punishment is a popular tool when governing commons in situations where free
riders would otherwise take over. It is well known that sanctioning systems,
such as the police and courts, are costly and thus can suffer from those who
free ride on other's efforts to maintain the sanctioning systems (second-order
free riders). Previous game-theory studies showed that if populations are very
large, pool punishment rarely emerges in public good games, even when
participation is optional, because of second-order free riders. Here we show
that a matching fund for rewarding cooperation leads to the emergence of pool
punishment, despite the presence of second-order free riders. We demonstrate
that reward funds can pave the way for a transition from a population of free
riders to a population of pool punishers. A key factor in promoting the
transition is also to reward those who contribute to pool punishment, yet not
abstaining from participation. Reward funds eventually vanish in raising pool
punishment, which is sustainable by punishing the second-order free riders.
This suggests that considering the interdependence of reward and punishment may
help to better understand the origins and transitions of social norms and
institutions.


  Modeling approaches of suspended, rod-like particles and recent experimental
data have shown that depletion forces display different signatures depending on
the orientation of these particles. It has been shown that axial attraction of
two rods yields contractile forces of 0.1pN that are independent of the
relative axial shift of the two rods. Here, we measured depletion-caused
interactions of actin bundles extending the phase space of single pairs of rods
to a multi-particle system. In contrast to a filament pair, we found forces up
to 3pN . Upon bundle relaxation forces decayed exponentially with a mean decay
time of 3.4s . These different dynamics are explained within the frame of a
mathematical model by taking pairwise interactions to a multi-filament scale.
The macromolecular content employed for our experiments is well below the
crowding of cells. Thus, we propose that arising forces can contribute to
biological force generation without the need to convert chemical energy into
mechanical work.


  In the rapidly evolving domain of next generation sequencing and
bioinformatics analysis, data generation is one aspect that is increasing at a
concomitant rate. The burden associated with processing large amounts of
sequencing data has emphasised the need to allocate sufficient computing
resources to complete analyses in the shortest possible time with manageable
and predictable costs. A novel method for predicting time to completion for a
popular bioinformatics software (QIIME), was developed using key variables
characteristic of the input data assumed to impact processing time. Multiple
Linear Regression models were developed to determine run time for two denoising
algorithms and a general bioinformatics pipeline. The models were able to
accurately predict clock time for denoising sequences from a naturally
assembled community dataset, but not an artificial community. Speedup and
efficiency tests for AmpliconNoise also highlighted that caution was needed
when allocating resources for parallel processing of data. Accurate modelling
of computational processing time using easily measurable predictors can assist
NGS analysts in determining resource requirements for bioinformatics software
and pipelines. Whilst demonstrated on a specific group of scripts, the
methodology can be extended to encompass other packages running on multiple
architectures, either in parallel or sequentially.


  Flagellated bacteria, such as Escherichia coli, perform directed motion in
gradients of concentration of attractants and repellents in a process called
chemotaxis. The E. coli chemotaxis signaling pathway is a model for signal
transduction, but it has unique features. We demonstrate that the need for fast
signaling necessitates high abundances of the proteins involved in this
pathway. We show that further constraints on the abundances of chemotaxis
proteins arise from the requirements of self-assembly, both of flagellar motors
and of chemoreceptor arrays. All these constraints are specific to chemotaxis,
and published data confirm that chemotaxis proteins tend to be more highly
expressed than their homologs in other pathways. Employing a chemotaxis pathway
model, we show that the gain of the pathway at the level of the response
regulator CheY increases with overall chemotaxis protein abundances. This may
explain why, at least in one E. coli strain, the abundance of all chemotaxis
proteins is higher in media with lower nutrient content. We also demonstrate
that the E. coli chemotaxis pathway is particularly robust to abundance
variations of the motor protein FliM.


  A two-valued fitness landscape is introduced for the classical Eigen's
quasispecies model. This fitness landscape can be considered as a direct
generalization of the so-called single or sharply peaked landscape. A general,
non permutation invariant quasispecies model is studied, therefore the
dimension of the problem is $2^N\times 2^N$, where $N$ is the sequence length.
It is shown that if the fitness function is equal to $w+s$ on a $G$-orbit $A$
and is equal to $w$ elsewhere, then the mean population fitness can be found as
the largest root of an algebraic equation of degree at most $N+1$. Here $G$ is
an arbitrary isometry group acting on the metric space of sequences of zeroes
and ones of the length $N$ with the Hamming distance. An explicit form of this
exact algebraic equation is given in terms of the spherical growth function of
the $G$-orbit $A$. Sufficient conditions for the so-called error threshold for
sequences of orbits are given. Motivated by the analysis of the two-valued
fitness landscapes an abstract generalization of Eigen's model is introduced
such that the sequences are identified with the points of a finite metric space
$X$ together with a group of isometries acting transitively on $X$. In
particular, a simplicial analogue of the original quasispecies model is
discussed, which can be considered as a mathematical model of the switching of
the antigenic variants for some bacteria.


  In this paper, we develop further the theory of thermostated systems along
the lines of our earlier paper. Two results are highlighted: 1) in the Markov
limit of the contracted description, a least dissipation of Helmholtz free
energy principle is established; and 2) a detailed account of the
appropriateness of this principle for nano-biology, including the evolution of
life, is presented.


  A new approach is developed to understand stability of a population and
further understanding of population momentum.


  Aedes aegypti is known as the responsible vector transmitting dengue
flavivirus. Unavailability of medication to cure the transmission of the virus
in the human blood becomes a global health issue in recent decades. World
epidemiologists are encouraged to focus on the investigation over the effective
and inexpensive way to prevent dengue transmission, i.e. mosquito control. In
this paper, we present a model depicting the dynamics of mosquito population
based on indoor-outdoor life cycle classification. The basic mosquito offspring
number was obtained and analysis of equilibria was shown. We brought along a
discussion on the application of optimal control to the model in which two
simultaneous schemes were introduced. The first scheme is done by disseminating
chemical like temephos in spots where eggs and larvae develop, meanwhile the
second scheme is done by deploying fumigation through areas where adult
mosquitoes prevalently nest, indoor as well as outdoor. A version of the
gradient-based method was presented to set up a workflow in minimizing the
objective functional with respect to some control variables. Numerical results
from the analysis of the basic mosquito offspring number with constant control
and from that with optimal control suggested that the application of fumigation
is preferable over that of temephos. It was also suggested that applying both
control schemes simultaneously gives the most significant reduction in the
population.


  We discuss the historic mortality record corresponding to the initial focus
of the yellow fever epidemic outbreak registered in Buenos Aires during the
year 1871 as compared to simulations of a stochastic population dynamics model.
This model incorporates the biology of the urban vector of yellow fever, the
mosquito Aedes aegypti, the stages of the disease in the human being as well as
the spatial extension of the epidemic outbreak. After introducing the
historical context and the restrictions it puts on initial conditions and
ecological parameters, we discuss the general features of the simulation and
the dependence on initial conditions and available sites for breeding the
vector. We discuss the sensitivity, to the free parameters, of statistical
estimators such as: final death toll, day of the year when the outbreak reached
half the total mortality and the normalized daily mortality, showing some
striking regularities. The model is precise and accurate enough to discuss the
truthfulness of the presently accepted historic discussions of the epidemic
causes, showing that there are more likely scenarios for the historic facts.


  Elementary flux modes (EFMs) are vectors defined from a metabolic reaction
network, giving the connections between substrates and products. EFMs-based
metabolic flux analysis (MFA) estimates the flux over each EFM from external
flux measurements through least-squares data fitting. In previous work we
presented an optimization method of column generation type that facilitates
EFMs-based MFA when the metabolic reaction network is so large that enumerating
all EFMs is prohibitive. In this work we extend this model by including errors
on measurements in a robust optimization framework. In the robust optimization
problem, the least-squares data fitting is minimized subject to the error on
each metabolite being as unfavourable as it can be, within a given interval. In
general, inclusion of robustness may make the optimization problem
significantly harder. However, we show that in our case the robust problem can
be stated as a convex quadratic programming problem, i.e., of the same form as
the original non-robust problem. Additionally, we demonstrate that the
column-generation technique of the non-robust problem can be extended also to
the robust problem. Furthermore, the option to indicate intervals on
metabolites that are not measured is introduced in this column generation
framework. The effect of including robustness in the model is evaluated in a
case-study, which indicated that the solutions of our non-robust problems are
in fact near-optimal also when robustness is considered. On the other hand, the
addition of intervals on unmeasured metabolites resulted in a change of optimal
solution. Implying that the inclusion of intervals on unmeasured metabolites is
more important than the direct consideration of measurement errors, this
despite up to 20% errors.


  Transcription is one of the essential processes for cells to read genetic
information encoded in genes, which is initiated by the binding of RNA
polymerase to related promoter. Experiments have found that the nucleotide
sequence of promoter has great influence on gene expression strength, or
promoter activity. In synthetic biology, one interesting question is how we can
synthesize a promoter with given activity, and which positions of promoter
sequence are important for determining its activity. In this study, based on
recent experimental data, correlations between promoter activity and its
sequence positions are analyzed by various methods. Our results show that,
except nucleotides in the two highly conserved regions, $-35$ box and $-10$
box, influences of nucleotides in other positions are also not neglectable. For
example, modifications of nucleotides around position $-19$ in spacing region
may change promoter activity in a large scale. The results of this study might
be helpful to our understanding of biophysical mechanism of gene transcription,
and may also be helpful to the design of synthetic cell factory.


  Coalescent histories are combinatorial structures that describe for a given
gene tree and species tree the possible lists of branches of the species tree
on which the gene tree coalescences take place. Properties of the number of
coalescent histories for gene trees and species trees affect a variety of
probabilistic calculations in mathematical phylogenetics. Exact and asymptotic
evaluations of the number of coalescent histories, however, are known only in a
limited number of cases. Here we introduce a particular family of species
trees, the \emph{lodgepole} species trees $(\lambda_n)_{n\geq 0}$, in which
tree $\lambda_n$ has $m=2n+1$ taxa. We determine the number of coalescent
histories for the lodgepole species trees, in the case that the gene tree
matches the species tree, showing that this number grows with $m!!$ in the
number of taxa $m$. This computation demonstrates the existence of tree
families in which the growth in the number of coalescent histories is faster
than exponential. Further, it provides a substantial improvement on the lower
bound for the ratio of the largest number of matching coalescent histories to
the smallest number of matching coalescent histories for trees with $m$ taxa,
increasing a previous bound of $(\sqrt{\pi} / 32)[(5m-12)/(4m-6)] m \sqrt{m}$
to $[ \sqrt{m-1}/(4 \sqrt{e}) ]^{m}$. We discuss the implications of our
enumerative results for phylogenetic computations.


  In contrast to the well-known destabilization of globular proteins by high
pressure, re- cent work has shown that pressure stabilizes the formation of
isolated {\alpha}-helices. However all simulations to date have obtained a
qualitatively opposite result within the experimen- tal pressure range. We show
that using a protein force field (Amber03w) parametrized in conjunction with an
accurate water model (TIP4P/2005) recovers the correct pressure- dependence and
an overall stability diagram for helix formation similar to that from experi-
ment; on the other hand, we confirm that using TIP3P water results in a very
weak pressure destabilization of helices. By carefully analyzing the
contributing factors, we show that this is not merely a consequence of
different peptide conformations sampled using TIP3P. Rather, there is a
critical role for the solvent itself in determining the dependence of total
system volume (peptide and solvent) on helix content. Helical peptide
structures exclude a smaller volume to water, relative to non-helical
structures with both the water models, but the total system volume for helical
conformations is higher than non-helical conformations with TIP3P water at low
to intermediate pressures, in contrast to TIP4P/2005 water. Our results further
emphasize the importance of using an accurate water model to study protein
folding under conditions away from standard temperature and pressure.


  Evolution in its course found a variety of solutions to the same optimisation
problem. The advent of high-throughput genomic sequencing has made available
extensive data from which, in principle, one can infer the underlying structure
on which biological functions rely. In this paper, we present a new method
aimed at extracting sites encoding structural and func- tional properties from
a set of protein primary sequences, namely a Multiple Sequence Alignment. The
method, called Critical Variable Selection, is based on the idea that subsets
of relevant sites cor- respond to subsequences that occur with a particularly
broad frequency distribution in the dataset. By applying this algorithm to in
silico sequences, to the Response Regulator Receiver and to the Voltage Sensor
Domain of Ion Channels, we show that this procedure recovers not only
information encoded in single site statistics and pairwise correlations but it
also captures dependencies going beyond pairwise correlations. The method
proposed here is complementary to Statistical Coupling Analysis, in that the
most relevant sites predicted by the two methods markedly differ. We find
robust and consistent results for datasets as small as few hundred sequences,
that reveal a hidden hierarchy of sites that is consistent with present
knowledge on biologically relevant sites and evo- lutionary dynamics. This
suggests that Critical Variable Selection is able to identify in a Multiple
Sequence Alignment a core of sites encoding functional and structural
information.


  We study the spatial patterns formed by a system of interacting particles
where the mobility of any individual is determined by the population crowding
at two different spatial scales. In this way we model the behavior of some
biological organisms (like mussels) that tend to cluster at short ranges as a
defensive strategy, and strongly disperse if there is a high population
pressure at large ranges for optimizing foraging. We perform stochastic
simulations of a particle-level model of the system, and derive and analyze a
continuous density description (a nonlinear diffusion equation). In both cases
we show that this interplay of scale-dependent-behaviors gives rise to a rich
formation of spatial patterns ranging from labyrinths to periodic cluster
arrangements. In most cases these clusters have the very peculiar appearance of
ring-like structures, i.e., organisms arranging in the perimeter of the
clusters, that we discuss in detail.


  In bacteria such as $\textit{Escherichia coli}$, DNA is compacted into a
nucleoid near the cell center, while ribosomes$-$molecular complexes that
translate messenger RNAs (mRNAs) into proteins$-$are mainly localized at the
poles. We study the impact of this spatial organization using a minimal
reaction-diffusion model for the cellular transcriptional-translational
machinery. Our model predicts that $\sim 90\%$ of mRNAs are segregated to the
poles and reveals a "circulation" of ribosomes driven by the flux of mRNAs,
from synthesis in the nucleoid to degradation at the poles. To address the
existence of non-specific, transient interactions between ribosomes and mRNAs,
we developed a novel method to efficiently incorporate such transient
interactions into reaction-diffusion equations, which allowed us to quantify
the biological implications of such non-specific interactions, e.g. for
ribosome efficiency.


  We introduce a scale-free method for testing the proportionality of branch
lengths between two phylogenetic trees that have the same topology and contain
the same set of taxa. This method scales both trees to a total length of 1 and
sums up the differences for each branch. Compared to previous methods, ours
yields a fully symmetrical score that measures proportionality without being
affected by scale. We call this score the normalized tree distance (NTD). Based
on real data, we demonstrate that NTD scores are distributed unimodally, in a
manner similar to a lognormal distribution. The NTD score can be used to, for
example, detect co-evolutionary processes and measure the accuracy of branch
length estimates.


  Societies consisting of cooperative individuals seem to require for their
continuing success that defectors be policed. The precise connection between
punishers and benefits, population structure, and division of labour, however,
remains ill-understood. Many models assume costly "peer punishment" to enforce
cooperation, but results in the economics literature suggest that this
assumption may not be generally valid. In many human and animal societies,
there is a division of labour between a purely supportive majority and a
dedicated minority of police-like enforcers. Here we present several extensions
to the Public Goods Game with punishment which allow for this possibility, and
evaluate their influence on the level of cooperative behaviour. We find that a
structure of separate subpopulations, which only interact through migration of
individuals, can have a strong effect on the evolutionary dynamics of a system
and significantly facilitate cooperation. Forcing defectors to contribute and
enabling fitness transfers to punishers both have a weak positive effect on
cooperation levels. In the presence of group competition, however, evolutionary
effects can paradoxically hinder cooperation.


  VDJ rearrangement and somatic hypermutation work together to produce
antibody-coding B cell receptor (BCR) sequences for a remarkable diversity of
antigens. It is now possible to sequence these BCRs in high throughput;
analysis of these sequences is bringing new insight into how antibodies
develop, in particular for broadly-neutralizing antibodies against HIV and
influenza. A fundamental step in such sequence analysis is to annotate each
base as coming from a specific one of the V, D, or J genes, or from an
N-addition (a.k.a. non-templated insertion). Previous work has used simple
parametric distributions to model transitions from state to state in a hidden
Markov model (HMM) of VDJ recombination, and assumed that mutations occur via
the same process across sites. However, codon frame and other effects have been
observed to violate these parametric assumptions for such coding sequences,
suggesting that a non-parametric approach to modeling the recombination process
could be useful. In our paper, we find that indeed large modern data sets
suggest a model using parameter-rich per-allele categorical distributions for
HMM transition probabilities and per-allele-per-position mutation
probabilities, and that using such a model for inference leads to significantly
improved results. We present an accurate and efficient BCR sequence annotation
software package using a novel HMM "factorization" strategy. This package,
called partis (https://github.com/psathyrella/partis/), is built on a new
general-purpose HMM compiler that can perform efficient inference given a
simple text description of an HMM.


  Knots are abundant in globular homopolymers but rare in globular proteins. To
shed new light on this long-standing conundrum, we study the influence of
sequence on the formation of knots in proteins under native conditions within
the framework of the hydrophobic-polar (HP) lattice protein model. By employing
large scale Wang-Landau simulations combined with suitable Monte Carlo trial
moves we show that, even though knots are still abundant on average, sequence
introduces large variability in the degree of self-entanglements. Moreover, we
are able to design sequences which are either almost always or almost never
knotted. Our findings serve as proof of concept that the introduction of just
one additional degree of freedom per monomer (in our case sequence) facilitates
evolution towards a protein universe in which knots are rare.


  We show that Wang-Landau sampling, combined with suitable Monte Carlo trial
moves, provides a powerful method for both the ground state search and the
determination of the density of states for the hydrophobic-polar (HP) protein
model and the interacting self-avoiding walk (ISAW) model for homopolymers. We
obtained accurate estimates of thermodynamic quantities for HP sequences with
$>100$ monomers and for ISAWs up to $>500$ monomers. Our procedure possesses an
intrinsic simplicity and overcomes the limitations inherent in more tailored
approaches making it interesting for a broad range of protein and polymer
models.


  Evolution depends on the fitness of organisms, the expected rate of
reproducing. Directly getting offspring is the most basic form of fitness, but
fitness can also be increased indirectly by helping genetically related
individuals (such as kin) to increase their fitness. The combined effect is
known as inclusive fitness. Here it is argued that a further elaboration of
fitness has evolved, particularly in humans. It is called extensive fitness and
it incorporates producing organisms that are merely similar in phenotype. The
evolvability of this mechanism is illustrated by computations on a simple model
combining heredity and behaviour. Phenotypes are driven into the direction of
high fitness through a mechanism that involves an internal estimate of fitness,
implicitly made within the organism itself. This mechanism has recently been
conjectured to be responsible for producing agency and goals. In the model,
inclusive and extensive fitness are both implemented by letting fitness
increase nonlinearly with the size of subpopulations of similar heredity (for
the indirect part of inclusive fitness) and of similar phenotype (for the
phenotypic part of extensive fitness). Populations implementing extensive
fitness outcompete populations implementing mere inclusive fitness. This occurs
because groups with similar phenotype tend to be larger than groups with
similar heredity, and fitness increases more when groups are larger. Extensive
fitness has two components, a direct component where individuals compete in
inducing others to become like them and an indirect component where individuals
cooperate and help others who are already similar to them.


  We present the mechanism of interaction of Wnt network module, which is
responsible for periodic sometogenesis, with p53 regulatory network, which is
one of the main regulators of various cellular functions, and switching of
various oscillating states by investigating p53-Wnt model. The variation in
Nutlin concentration in p53 regulating network drives the Wnt network module to
different states, stabilized, damped and sustain oscillation states, and even
to cycle arrest. Similarly, the change in Axin concentration in Wnt could able
to modulate the p53 dynamics at these states. We then solve the set of coupled
ordinary differential equations of the model using quasi steady state
approximation. We, further, demonstrate the change of p53 and Gsk3 interaction
rate, due to hypothetical catalytic reaction or external stimuli, can able to
regulate the dynamics of the two network modules, and even can control their
dynamics to protect the system from cycle arrest (apoptosis).


  The statistical estimation of phylogenies is always associated with
uncertainty, and accommodating this uncertainty is an important component of
modern phylogenetic comparative analysis. The birth-death polytomy resolver is
a method of accounting for phylogenetic uncertainty that places missing
(unsampled) taxa onto phylogenetic trees, using taxonomic information alone.
Recent studies of birds and mammals have used this approach to generate
pseudo-posterior distributions of phylogenetic trees that are complete at the
species level, even in the absence of genetic data for many species. Many
researchers have used these distributions of phylogenies for downstream
evolutionary analyses that involve inferences on phenotypic evolution,
geography, and community assembly. I demonstrate that the use of phylogenies
constructed in this fashion is inappropriate for many questions involving
traits. Because species are placed on trees at random with respect to trait
values, the birth-death polytomy resolver breaks down natural patterns of trait
phylogenetic structure. Inferences based on these trees are predictably and
often drastically biased in a direction that depends on the underlying (true)
pattern of phylogenetic structure in traits. I illustrate the severity of the
phenomenon for both continuous and discrete traits using examples from a global
bird phylogeny.


  Scoring DNA sequences against Position Weight Matrices (PWMs) is a widely
adopted method to identify putative transcription factor binding sites. While
common bioinformatics tools produce scores that can reflect the binding
strength between a specific transcription factor and the DNA, these scores are
not directly comparable between different transcription factors. Here, we
provide two different ways to find the scaling parameter $\lambda$ that allows
us to infer binding energy from a PWM score. The first approach uses a PWM and
background genomic sequence as input to estimate $\lambda$ for a specific
transcription factor, which we applied to show that $\lambda$ distributions for
different transcription factor families correspond with their DNA binding
properties. Our second method can reliably convert $\lambda$ between different
PWMs of the same transcription factor, which allows us to directly compare PWMs
that were generated by different approaches. These two approaches provide
consistent and computationally efficient ways to scale PWMs scores and estimate
transcription factor binding sites strength.