  Timing is essential for many cellular processes, from cellular responses to
external stimuli to the cell cycle and circadian clocks. Many of these
processes are based on gene expression. For example, an activated gene may be
required to reach in a precise time a threshold level of expression that
triggers a specific downstream process. However, gene expression is subject to
stochastic fluctuations, naturally inducing an uncertainty in this
threshold-crossing time with potential consequences on biological functions and
phenotypes. Here, we consider such "timing fluctuations", and we ask how they
can be controlled. Our analytical estimates and simulations show that, for an
induced gene, timing variability is minimal if the threshold level of
expression is approximately half of the steady-state level. Timing fuctuations
can be reduced by increasing the transcription rate, while they are insensitive
to the translation rate. In presence of self-regulatory strategies, we show
that self-repression reduces timing noise for threshold levels that have to be
reached quickly, while selfactivation is optimal at long times. These results
lay a framework for understanding stochasticity of endogenous systems such as
the cell cycle, as well as for the design of synthetic trigger circuits.


  We report that protein confinement within nanoscopic vesicular compartments
corresponds to a liquid-liquid phase transition with the protein/water within
vesicle lumen interacting very differently than in bulk. We show this effect
leads to considerable structural changes on the proteins with evidence
suggesting non-alpha helical conformations. Most importantly both aspects lead
to a significant improvement on protein stability against thermal denaturation
up to 95degC at neutral pH, with little or no evidence of unfolding or reduced
enzymatic activity. The latter parameter does indeed exhibit an increase after
thermal cycling. Our results suggest that nanoscopic confinement is a promising
new avenue for the enhanced long-term storage of proteins. Moreover, our
investigations have potentially important implications for the origin of life,
since such compartmentalization may well have been critical for ensuring the
preservation of primordial functional proteins under relatively harsh
conditions, thus playing a key role in the subsequent emergence of primitive
life forms.


  Stochastic Boolean networks, or more generally, stochastic discrete networks,
are an important class of computational models for molecular interaction
networks. The stochasticity stems from the updating schedule. Standard updating
schedules include the synchronous update, where all the nodes are updated at
the same time, and the asynchronous update where a random node is updated at
each time step. The former produces a deterministic dynamics while the latter a
stochastic dynamics. A more general stochastic setting considers propensity
parameters for updating each node. Stochastic Discrete Dynamical Systems (SDDS)
is a modeling framework that considers two propensity parameters for updating
each node and uses one when the update has a positive impact on the variable,
that is, when the update causes the variable to increase its value, and uses
the other when the update has a negative impact, that is, when the update
causes it to decrease its value. This framework offers additional features for
simulations but also adds a complexity in parameter estimation of the
propensities. This paper presents a method for estimating the propensity
parameters for SDDS. The method is based on adding noise to the system using
the Google PageRank approach to make the system ergodic and thus guaranteeing
the existence of a stationary distribution. Then with the use of a genetic
algorithm, the propensity parameters are estimated. Approximation techniques
that make the search algorithms efficient are also presented and Matlab/Octave
code to test the algorithms are available at
http://www.ms.uky.edu/~dmu228/GeneticAlg/Code.html.


  Repetitive elements are important in genomic structures, functions and
regulations, yet effective methods in precisely identifying repetitive elements
in DNA sequences are not fully accessible, and the relationship between
repetitive elements and periodicities of genomes is not clearly understood. We
present an $\textit{ab initio}$ method to quantitatively detect repetitive
elements and infer the consensus repeat pattern in repetitive elements. The
method uses the measure of the distribution uniformity of nucleotides at
periodic positions in DNA sequences or genomes. It can identify periodicities,
consensus repeat patterns, copy numbers and perfect levels of repetitive
elements. The results of using the method on different DNA sequences and
genomes demonstrate efficacy and accuracy in identifying repeat patterns and
periodicities. The complexity of the method is linear with respect to the
lengths of the analyzed sequences.


  Biological pest control is increasingly used in agriculture as a an
alternative to traditional chemical pest control. In many cases, this involves
a one-off or periodic release of naturally occurring and/or genetically
modified enemies such as predators, parasitoids, or pathogens. As the
interaction between these enemies and the pest is complex and the production of
natural enemies potentially expensive, it is not surprising that both the
efficacy and economic viability of biological pest control are debated. Here,
we investigate the performance of very simple control strategies. In
particular, we show how Pareto-efficient one-off or periodic release
strategies, that optimally trade off between efficacy and economic viability,
can be devised and used to enable high efficacy at small economic costs. We
demonstrate our method on a pest-pathogen-crop model with a tunable immigration
rate of pests. By analyzing this model, we demonstrate that simple
Pareto-efficient one-off and periodic release strategies are efficacious and
simultaneously have profits that are close to the theoretical maximum obtained
by strategies optimizing only the profit. When the immigration rate of pests is
low to intermediate, one-off control strategies are sufficient, and when the
immigration of pests is high, periodic release strategies are preferable. The
methods presented here can be extended to more complex scenarios and be used to
identify promising biological pest control strategies in many circumstances.


  This paper proposes that eukaryotic cells, under severe genotoxic stress, can
commit genoautotomy (genome 'self-injury') that involves cutting and releasing
single-stranded DNA (ssDNA) fragments from double-stranded DNA and leaving
ssDNA gaps in the genome. The ssDNA gaps could be easily and precisely repaired
later. The released ssDNA fragments may play some role in the regulation of
cell cycle progression. Taken together, genoautotomy causes limited nonlethal
DNA damage, but prevents the whole genome from lethal damage, and thus should
be deemed as a eukaryotic cellular defence response to genotoxic stress.


  Cancer stem cells are controlled by developmental networks that are often
topologically indistinguishable from normal, healthy stem cells. The question
is why cancer stem cells can be both phenotypically distinct and have
morphological effects so different from normal stem cells. The difference
between cancer stem cells and normal stem cells lies not in differences their
network architecture, but rather in the spatial-temporal locality of their
activation in the genome and the resulting expression in the body. The
metastatic potential cancer stem cells is not based primarily on their network
divergence from normal stem cells, but on non-network based genetic changes
that enable the evolution of gene-based phenotypic properties of the cell that
permit its escape and travel to other parts of the body. Stem cell network
theory allows the precise prediction of stem cell behavioral dynamics and a
mathematical description of stem cell proliferation for both normal and cancer
stem cells. It indicates that the best therapeutic approach is to tackle the
highest order stem cells first, otherwise spontaneous remission of so called
cured cancers will always be a danger. Stem cell networks point to a pathway to
new methods to diagnose and cure not only stem cell cancers but cancers
generally.


  Languages and genes are both transmitted from generation to generation, with
opportunity for differential reproduction and survivorship of forms. Here we
apply a rigorous inference framework, drawn from population genetics, to
distinguish between two broad mechanisms of language change: drift and
selection. Drift is change that results from stochasticity in transmission and
it may occur in the absence of any intrinsic difference between linguistic
forms; whereas selection is truly an evolutionary force arising from intrinsic
differences -- for example, when one form is preferred by members of the
population. Using large corpora of parsed texts spanning the 12th century to
the 21st century, we analyze three examples of grammatical changes in English:
the regularization of past-tense verbs, the rise of the periphrastic `do', and
syntactic variation in verbal negation. We show that we can reject stochastic
drift in favor of a selective force driving some of these language changes, but
not others. The strength of drift depends on a word's frequency, and so drift
provides an alternative explanation for why some words are more prone to change
than others. Our results suggest an important role for stochasticity in
language change, and they provide a null model against which selective theories
of language evolution must be compared.


  Detecting variation in the evolutionary process along chromosomes is
increasingly important as whole-genome data becomes more widely available. For
example, factors such as incomplete lineage sorting, horizontal gene transfer,
and chromosomal inversion are expected to result in changes in the underlying
gene trees along a chromosome, while changes in selective pressure and
mutational rates for different genomic regions may lead to shifts in the
underlying mutational process. We propose the split score as a general method
for quantifying support for a particular phylogenetic relationship within a
genomic data set. Because the split score is based on algebraic properties of a
matrix of site pattern frequencies, it can be rapidly computed, even for data
sets that are large in the number of taxa and/or in the length of the
alignment, providing an advantage over other methods (e.g., maximum likelihood)
that are often used to assess such support. Using simulation we explore the
properties of the split score, including its dependence on sequence length,
branch length, size of a split and its ability to detect true splits in the
underlying tree. Using a sliding window analysis, we show that split scores can
be used to detect changes in the underlying evolutionary process for
genome-scale data from primates, mosquitoes, and viruses in a computationally
efficient manner. Computation of the split score has been implemented in the
software package SplitSup.


  Background: Tumours are diverse ecosystems with persistent heterogeneity in
various cancer hallmarks like self-sufficiency of growth factor production for
angiogenesis and reprogramming of energy-metabolism for aerobic glycolysis.
This heterogeneity has consequences for diagnosis, treatment, and disease
progression.
  Methods: We introduce the double goods game to study the dynamics of these
traits using evolutionary game theory. We model glycolytic acid production as a
public good for all tumour cells and oxygen from vascularization via VEGF
production as a club good benefiting non-glycolytic tumour cells. This results
in three viable phenotypic strategies: glycolytic, angiogenic, and aerobic
non-angiogenic.
  Results: We classify the dynamics into three qualitatively distinct regimes:
(1) fully glycolytic, (2) fully angiogenic, or (3) polyclonal in all three cell
types. The third regime allows for dynamic heterogeneity even with linear
goods, something that was not possible in prior public good models that
considered glycolysis or growth-factor production in isolation.
  Conclusion: The cyclic dynamics of the polyclonal regime stress the
importance of timing for anti-glycolysis treatments like lonidamine. The
existence of qualitatively different dynamic regimes highlights the order
effects of treatments. In particular, we consider the potential of vascular
renormalization as a neoadjuvant therapy before follow up with interventions
like buffer therapy.


  We present Meraculous2, an update to the Meraculous short-read assembler that
includes (1) handling of allelic variation using "bubble" structures within the
de Bruijn graph, (2) improved gap closing, and (3) an improved scaffolding
algorithm that produces more complete assemblies without compromising
scaffolding accuracy. The speed and bandwidth efficiency of the new parallel
implementation have also been substantially improved, allowing the assembly of
a human genome to be accomplished in 24 hours on the JGI/NERSC Genepool system.
To highlight the features of Meraculous2 we present here the assembly of the
diploid human genome NA12878, and compare it with previously published
assemblies of the same data using other algorithms. The Meraculous2 assemblies
are shown to have better completeness, contiguity, and accuracy than other
published assemblies for these data. Practical considerations including
pre-assembly analyses of polymorphism and repetitiveness are described.


  Aptamers are chemically produced oligonucleotides, able to bind a variety of
targets such as drugs, proteins and pathogens with high sensitivity and
selectivity. Therefore, aptamers are largely employed for producing label-free
biosensors, with significant applications in diagnostics and drug delivery. In
particular, the anti-thrombin aptamers are biomolecules of high interest for
clinical use, because of their ability to recognize and bind the thrombin
enzyme. Among them, the DNA 15-mer thrombin-binding aptamer (TBA), has been
widely explored concerning both its structure, which was resolved with
different techniques, and its function, especially about the possibility of
using it as the active part of biosensors. This paper proposes a microscopic
model of the electrical properties of TBA and the aptamer-thrombin complex,
combining information from both structure and function. The novelty consists in
describing both the aptamer alone and the complex as an impedance network, thus
going deeper inside the issues addressed in an emerging electronics branch
known as proteotronics. The theoretical results are compared and validated with
Electrochemical Impedance Spectroscopy measurements reported in the literature.
Finally, the model suggests resistance measurements as a novel tool for testing
aptamer-target affinity.


  Water molecules inside G-protein coupled receptor have recently been
spotlighted in a series of crystal structures. To decipher the dynamics and
functional roles of internal waters in GPCR activity, we studied
A$_{\text{2A}}$ adenosine receptor using $\mu$sec-molecular dynamics
simulations. Our study finds that the amount of water flux across the
transmembrane (TM) domain varies depending on the receptor state, and that the
water molecules of the TM channel in the active state flow three times slower
than those in the inactive state. Depending on the location in solvent-protein
interface as well as the receptor state, the average residence time of water in
each residue varies from $\sim\mathcal{O}(10^2)$ psec to
$\sim\mathcal{O}(10^2)$ nsec. Especially, water molecules, exhibiting ultraslow
relaxation ($\sim\mathcal{O}(10^2)$ nsec) in the active state, are found around
the microswitch residues that are considered activity hotspots for GPCR
function. A continuous allosteric network spanning the TM domain, arising from
water-mediated contacts, is unique in the active state, underscoring the
importance of slow waters in the GPCR activation.


  This paper introduces a high-throughput software tool framework called {\it
sam2bam} that enables users to significantly speedup pre-processing for
next-generation sequencing data. The sam2bam is especially efficient on
single-node multi-core large-memory systems. It can reduce the runtime of data
pre-processing in marking duplicate reads on a single node system by 156-186x
compared with de facto standard tools. The sam2bam consists of parallel
software components that can fully utilize the multiple processors, available
memory, high-bandwidth of storage, and hardware compression accelerators if
available.
  The sam2bam provides file format conversion between well-known genome file
formats, from SAM to BAM, as a basic feature. Additional features such as
analyzing, filtering, and converting the input data are provided by {\it
plug-in} tools, e.g., duplicate marking, which can be attached to sam2bam at
runtime.
  We demonstrated that sam2bam could significantly reduce the runtime of NGS
data pre-processing from about two hours to about one minute for a whole-exome
data set on a 16-core single-node system using up to 130 GB of memory. The
sam2bam could reduce the runtime for whole-genome sequencing data from about 20
hours to about nine minutes on the same system using up to 711 GB of memory.


  Strongly correlated electrostatics of DNA systems has drawn the interest of
many groups, especially the condensation and overcharging of DNA by multivalent
counterions. By adding counterions of different valencies and shapes, one can
enhance or reduce DNA overcharging. In this papers, we focus on the effect of
multivalent co-ions, specifically divalent co-ions such as SO$_4^{2-}$. A
computational experiment of DNA condensation using Monte$-$Carlo simulation in
grand canonical ensemble is carried out where DNA system is in equilibrium with
a bulk solution containing a mixture of salt of different valency of co-ions.
Compared to system with purely monovalent co-ions, the influence of divalent
co-ions shows up in multiple aspects. Divalent co-ions lead to an increase of
monovalent salt in the DNA condensate. Because monovalent salts mostly
participate in linear screening of electrostatic interactions in the system,
more monovalent salt molecules enter the condensate leads to screening out of
short-range DNA$-$DNA like charge attraction and weaker DNA condensation free
energy. The overcharging of DNA by multivalent counterions is also reduced in
the presence of divalent co$-$ions. Strong repulsions between DNA and divalent
co-ions and among divalent co-ions themselves leads to a {\em depletion} of
negative ions near DNA surface as compared to the case without divalent
co-ions. At large distance, the DNA$-$DNA repulsive interaction is stronger in
the presence of divalent co$-$ions, suggesting that divalent co$-$ions role is
not only that of simple stronger linear screening.


  Recently, we investigated the geographical origins of Ashkenazic Jews (AJs)
and their native language Yiddish by applying a biogeographical tool, the
Geographic Population Structure (GPS), to a cohort of 367 exclusively
Yiddish-speaking and multilingual AJs genotyped on the Genochip microarray. GPS
localized most AJs along major ancient trade routes in northeastern Turkey
adjacent to primeval villages with names that may be derived from the word
"Ashkenaz." These findings were compatible with the hypothesis of an
Irano-Turko-Slavic origin for AJs and a Slavic origin for Yiddish and at odds
with the Rhineland hypothesis advocating a German origin of both. Our approach
has been recently adopted by Flegontov et al. (2016a) to trace the origin of
the Siberian Ket people and their language. Recently, Flegontov et al. (2016b)
have raised several questions concerning the accuracy of the Genochip
microarray and GPS, specifically in relation to AJs and Yiddish. Although many
of these issues have been addressed in our previous papers, we take this
opportunity to clarify the principles of the GPS approach, review the recent
biogeographical and ancient DNA findings regarding AJs, and comment on the
origin of Yiddish.


  Recently, a first step was made by the authors towards a systematic
investigation of the effect of reaction-step-size noise - uncertainty in the
step size of the reaction - on the dynamics of stochastic populations. This was
done by investigating the effect of bursty influx on the switching dynamics of
stochastic populations. Here we extend this formalism to account for bursty
reproduction processes, and improve the accuracy of the formalism to include
subleading-order corrections. Bursty reproduction appears in various contexts,
where notable examples include bursty viral production from infected cells, and
reproduction of mammals involving varying number of offspring. The main
question we quantitatively address is how bursty reproduction affects the
overall fate of the population. We consider two complementary scenarios:
population extinction and population survival; in the former a population gets
extinct after maintaining a long-lived metastable state, whereas in the latter
a population proliferates despite undergoing a deterministic drift towards
extinction. In both models reproduction occurs in bursts, sampled from an
arbitrary distribution. In the extinction problem, we show that bursty
reproduction broadens the quasi-stationary distribution of population sizes in
the metastable state, which results in an exponential decrease of the mean time
to extinction. In the survival problem, bursty reproduction yields an
exponential increase in survival probability of the population. Close to the
bifurcation limit our analytical results simplify considerably and are shown to
depend solely on the mean and variance of the burst-size distribution. Our
formalism is demonstrated on several realistic distributions which all compare
well with numerical Monte-Carlo simulations.


  Motivation: Genome rearrangement plays an important role in evolutionary
biology and has profound impacts on phenotype in organisms ranging from
microbes to humans. The mechanisms for genome rearrangement events remain
unclear. Lots of comparisons have been conducted among different species. To
reveal the mechanisms for rearrangement events, comparison of different
individuals/strains within the same species or genus (pan-genomes) is more
helpful since they are much closer to each other. Results: We study the
mechanism for inversion events via core-genome scaffold comparison of different
strains within the same species. We focus on two kinds of bacteria, Pseudomonas
aeruginosa and Escherichia coli, and investigate the inversion events among
different strains of the same specie. We find an interesting phenomenon that
long (larger than 10,000 bp) inversion regions are flanked by a pair of
Inverted Repeats (IRs) (with lengths ranging from 385 bp to 27476 bp) which are
often Insertion Sequences (ISs).This mechanism can also explain why the
breakpoint reuses for inversion events happen. We study the prevalence of the
phenomenon and find that it is a major mechanism for inversions. The other
observation is that for different rearrangement events such as transposition
and inverted block interchange, the two ends of the swapped regions are also
associated with repeats so that after the rearrangement operations the two ends
of the swapped regions remain unchanged. To our knowledge, this is the first
time such a phenomenon is reported for transposition event.


  The lengths of the telomere regions of chromosomes in a population of cells
are modelled using a chemical master equation formalism, from which the
evolution of the average number of cells of each telomere length is extracted.
In particular, the role of the telomere-elongating enzyme telomerase on these
dynamics is investigated. We show that for biologically relevant rates of cell
birth and death, one finds a critical rate, R_crit, of telomerase activity such
that the total number of cells diverges. Further, R_crit is similar in
magnitude to the rates of mitosis and cell death. The possible relationship of
this result to replicative immortality and its associated hallmark of cancer is
discussed.


  MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression
post-transcriptionally. They are involved in key biological processes and then
may play a major role in the development of human diseases including cancer, in
particular their involvement in breast cancer metastasis has been confirmed.
Recently, the authors of ref.(\cite{key1} have found that miR-429 may have a
role in the inhibition of breast cancer metastasis and have identified its
target gene CRKL as a potential candidate. In this paper, by using systems
biology tools we have shown that CRKL is involved in positive regulation of
ERK1/2 signaling pathway and contribute to the regulation of LYN through a
topological generalization of feed forward loop.


  In contrary to the common belief that only one strand of the pre-miRNA is
active (usually the 5p one that is the more abundant) while the second one
(miRNA*) is discarded, functional 5p and 3p have been observed for many miRNAs.
Among those miRNAs is miR-34a which is a target gene of the tumor suppressor
p53. In this paper we have re-examined the role of miR-34a-5p and miR-34a-3p in
the signaling pathway of p53. We found that they target overlapping sets of
genes (MDM2 and THBS1). By a GO enrichment analysis we found that THBS1 is
involved in cancer and metastasis relevant processes. We have also deduce that
p53, MDM2 and miR-34a are linked by a type 1 incoherent FFL that represent a
novel mechanism for accelerating the response of p53 to external stress
signals.


  Random heteropolymers are a minimal description of biopolymers and can
provide a theoretical framework to the investigate the formation of loops in
biophysical experiments. A two--state model provides a consistent and robust
way to study the scaling properties of loop formation in polymers of the size
of typical biological systems. Combining it with self--adjusting
simulated--tempering simulations, we can calculate numerically the looping
properties of several realizations of the random interactions within the chain.
Differently from homopolymers, random heteropolymers display at different
temperatures a continuous set of scaling exponents. The necessity of using
self--averaging quantities makes finite--size effects dominant at low
temperatures even for long polymers, shadowing the length--independent
character of looping probability expected in analogy with homopolymeric
globules. This could provide a simple explanation for the small scaling
exponents found in experiments, for example in chromosome folding.


  The $\alpha$ and $\beta$ subunits comprising the hexameric assembly of
F1-ATPase share a high degree of structural identity, though low primary
identity. Each subunit binds nucleotide in similar pockets, yet only $\beta$
subunits are catalytically active. Why? We re-examine their internal symmetry
axes and observe interesting differences. Dividing each chain into an
N-terminal head region, a C-terminal foot region, and a central torso, we
observe (1) that while the foot and head regions in all chains obtain high and
similar mobility, the torsos obtain different mobility profiles, with the
$\beta$ subunits exhibiting a higher motility compared to the $\alpha$
subunits, a trend supported by the crystallographic B-factors. The $\beta$
subunits have greater torso mobility by having fewer distributed, nonlocal
packing interactions providing a spacious and soft connectivity, and offsetting
the resultant softness with local stiffness elements, including an additional
$\beta$ sheet. (2) A loop near the nucleotide binding-domain of the $\beta$
subunits, absent in the $\alpha$ subunits, swings to create a large variation
in the occlusion of the nucleotide binding region. (3) A combination of the
softest three eigenmodes significantly reduces the RMSD between the open and
closed conformations of the $\beta$ subnits. (4) Comparisons of computed and
observed crystallographic B-factors suggest a suppression of a particular
symmetry axis in an $\alpha$ subunit. (5) Unexpectedly, the soft intra-monomer
oscillations pertain to distortions that do not create inter-monomer steric
clashes in the assembly, suggesting that structural optimization of the
assembly evolved at all levels of complexity.


  Almost all current approaches for engineering modular logic components in
synthetic biology use first-order regulators, including most CRISPR/CAS, TAL,
zinc finger, and RNA interference systems. Many practitioners understand
intuitively that second and higher order binding is necessary for scalability,
and this is easy to show for single-input single-output systems. However, no
study to date has analysed whether a more complex system, utilizing e.g.
feedback or error correction, can produce scalable computation from first-order
regulators. We prove here that first order repressor systems cannot support
bistability. In the process, we introduce a function G to measure signal
quality in molecular systems, and we show that G always decreases in dynamic
feedback systems as well as static feed-forward logic cascades of first-order
repressors. As a result, first order repressors cannot build memory or signal
buffering elements. Finally, we suggest G as a potential new property for
characterization of standard biological parts.


  Optimal sex allocation theory is one of the most intricately developed areas
of evolutionary ecology. Under a range of conditions, particularly under
population sub-division, selection favours sex being allocated to offspring
non-randomly, generating non-binomial variances of offspring group sex ratios.
Detecting non-binomial sex allocation is complicated by stochastic
developmental mortality, as offspring sex can often only be identified on
maturity with the sex of non-maturing offspring remaining unknown. We show that
current approaches for detecting non-binomiality have limited ability to detect
non-binomial sex allocation when developmental mortality has occurred. We
present a new procedure using an explicit model of sex allocation and mortality
and develop a Bayesian model selection approach (available as an R package). We
use the double and multiplicative binomial distributions to model over- and
under-dispersed sex allocation and show how to calculate Bayes factors for
comparing these alternative models to the null hypothesis of binomial sex
allocation. The ability to detect non-binomial sex allocation is greatly
increased, particularly in cases where mortality is common. The use of Bayesian
methods allows for the quantification of the evidence in favour of each
hypothesis, and our modelling approach provides an improved descriptive
capability over existing approaches. We use a simulation study to demonstrate
substantial improvements in power for detecting non-binomial sex allocation in
situations where current methods fail, and we illustrate the approach in real
scenarios using empirically obtained datasets on the sexual composition of
groups of gregarious parasitoid wasps.


  Population diversity is essential for avoiding premature convergence in
Genetic Algorithms (GAs) and for the effective use of crossover. Yet the
dynamics of how diversity emerges in populations are not well understood. We
use rigorous runtime analysis to gain insight into population dynamics and GA
performance for the ($\mu$+1) GA and the $\text{Jump}_k$ test function. We show
that the interplay of crossover and mutation may serve as a catalyst leading to
a sudden burst of diversity. This leads to improvements of the expected
optimisation time of order $\Omega(n/\log n)$ compared to mutation-only
algorithms like (1+1) EA. Moreover, increasing the mutation rate by an
arbitrarily small constant factor can facilitate the generation of diversity,
leading to speedups of order $\Omega(n)$. We also compare seven commonly used
diversity mechanisms and evaluate their impact on runtime bounds for the
($\mu$+1) GA. All previous results in this context only hold for
unrealistically low crossover probability $p_c=O(k/n)$, while we give analyses
for the setting of constant $p_c < 1$ in all but one case.
  For the typical case of constant $k > 2$ and constant $p_c$, we can compare
the resulting expected runtimes for different diversity mechanisms assuming an
optimal choice of $\mu$: $O(n^{k-1})$ for duplicate elimination/minim.,
$O(n^2\log n)$ for maximising the convex hull, $O(n\log n)$ for deterministic
crowding (assuming $p_c = k/n$), $O(n\log n)$ for maximising Hamming distance,
$O(n\log n)$ for fitness sharing, $O(n\log n)$ for single-receiver island
model.
  This proves a sizeable advantage of all variants of the ($\mu$+1) GA compared
to (1+1) EA, which requires time $\Theta(n^k)$. Experiments complement our
theoretical findings and further highlight the benefits of crossover and
diversity on $\text{Jump}_k$.


  How DNA is mapped to functional proteins is a basic question of living
matter. We introduce and study a physical model of protein evolution which
suggests a mechanical basis for this map. Many proteins rely on large-scale
motion to function. We therefore treat protein as learning amorphous matter
that evolves towards such a mechanical function: Genes are binary sequences
that encode the connectivity of the amino acid network that makes a protein.
The gene is evolved until the network forms a shear band across the protein,
which allows for long-range, soft modes required for protein function. The
evolution reduces the high-dimensional sequence space to a low-dimensional
space of mechanical modes, in accord with the observed dimensional reduction
between genotype and phenotype of proteins. Spectral analysis of the space of
$10^6$ solutions shows a strong correspondence between localization around the
shear band of both mechanical modes and the sequence structure. Specifically,
our model shows how mutations of the gene and their correlations occur at amino
acids whose interactions determine the functional mode.


  Dynamic cooperativity in monomeric enzymes is characterized in terms of a
non-Michaelis-Menten kinetic behaviour. The latter is believed to be associated
with mechanisms that include multiple reaction pathways due to enzymatic
conformational fluctuations. Recent advances in single-molecule fluorescence
spectroscopy have provided new fundamental insights on the possible mechanisms
underlying reactions catalyzed by fluctuating enzymes. Here, we present a
bottom-up approach to understand enzyme turnover kinetics at physiologically
relevant mesoscopic concentrations informed by mechanisms extracted from
single-molecule stochastic trajectories. The stochastic approach, presented
here, shows the emergence of dynamic cooperativity in terms of a slowing down
of the Michaelis-Menten (MM) kinetics resulting in negative cooperativity. For
fewer enzymes, dynamic cooperativity emerges due to the combined effects of
enzymatic conformational fluctuations and molecular discreteness. The increase
in the number of enzymes, however, suppresses the effect of enzymatic
conformational fluctuations such that dynamic cooperativity emerges solely due
to the discrete changes in the number of reacting species. These results
confirm that the turnover kinetics of fluctuating enzyme based on the
parallel-pathway MM mechanism switches over to the single-pathway MM mechanism
with the increase in the number of enzymes. For large enzyme numbers,
convergence to the exact MM equation occurs in the limit of very high substrate
concentration as the stochastic kinetics approaches the deterministic
behaviour.


  Deterministic evolutionary theory robustly predicts that populations
displaying altruistic behaviours will be driven to extinction by mutant cheats
that absorb common benefits but do not themselves contribute. Here we show that
when demographic stochasticity is accounted for, selection can in fact act in
the reverse direction to that predicted deterministically, instead favouring
cooperative behaviors that appreciably increase the carrying capacity of the
population. Populations that exist in larger numbers experience a selective
advantage by being more stochastically robust to invasions than smaller
populations, and this advantage can persist even in the presence of
reproductive costs. We investigate this general effect in the specific context
of public goods production and find conditions for stochastic selection
reversal leading to the success of public good producers. This insight,
developed here analytically, is missed by both the deterministic analysis as
well as standard game theoretic models that enforce a fixed population size.
The effect is found to be amplified by space; in this scenario we find that
selection reversal occurs within biologically reasonable parameter regimes for
microbial populations. Beyond the public good problem, we formulate a general
mathematical framework for models that may exhibit stochastic selection
reversal. In this context, we describe a stochastic analogue to r-K theory, by
which small populations can evolve to higher densities in the absence of
disturbance.


  Community ecology has traditionally relied on the competitive exclusion
principle, a piece of common wisdom in conceptual frameworks developed to
describe species assemblages. Key concepts in community ecology, such as
limiting similarity and niche partitioning, are based on competitive exclusion.
However, this classical paradigm in ecology relies on implications derived from
simple, deterministic models. Here we show how the predictions of a symmetric,
deterministic model about the way extinctions proceed can be utterly different
from the results derived from the same model when ecological drift (demographic
stochasticity) is explicitly considered. Using analytical approximations to the
steady-state conditional probabilities for assemblages with two and three
species, we demonstrate that stochastic competitive exclusion leads to a
cascade of extinctions, whereas the symmetric, deterministic model predicts a
multiple collapse of species. To test the robustness of our results, we have
studied the effect of environmental stochasticity and relaxed the species
symmetry assumption. Our conclusions highlight the crucial role of
stochasticity when deriving reliable theoretical predictions for species
community assembly.


  We analyze the paper of Nathan D. Grubaugh et al. (Nature 546, 401-405, 2017)
and find that it does not offer a convincing quantitative explanation for what
generated the temporal distribution of human Zika virus (ZIKV) cases shown in
their paper (Fig. 1d). We criticize this aspect because it is this
understanding of how human cases develop from day-today and week-to-week within
an area such as these Ground Zeros, that policymakers need in order to mitigate
future outbreaks. We present results that strongly suggest that the missing
piece is everyday human visit-revisit behavior. These results reproduce the
human outbreak data in the key areas of Miami in 2016 very well, and give
policymakers specific predictions for how changes in human flow through these
areas will affect, and hence can be used to mitigate, future ZIka outbreaks in
Miami and beyond.


  Characterizing genes with semantic information is an important process
regarding the description of gene products. In spite that complete genomes of
many organisms have been already sequenced, the biological functions of all of
their genes are still unknown. Since experimentally studying the functions of
those genes, one by one, would be unfeasible, new computational methods for
gene functions inference are needed. We present here a novel computational
approach for inferring biological function for a set of genes with previously
unknown function, given a set of genes with well-known information. This
approach is based on the premise that genes with similar behaviour should be
grouped together. This is known as the guilt-by-association principle. Thus, it
is possible to take advantage of clustering techniques to obtain groups of
unknown genes that are co-clustered with genes that have well-known semantic
information (GO annotations). Meaningful knowledge to infer unknown semantic
information can therefore be provided by these well-known genes. We provide a
method to explore the potential function of new genes according to those
currently annotated. The results obtained indicate that the proposed approach
could be a useful and effective tool when used by biologists to guide the
inference of biological functions for recently discovered genes. Our work sets
an important landmark in the field of identifying unknown gene functions
through clustering, using an external source of biological input. A simple web
interface to this proposal can be found at
http://fich.unl.edu.ar/sinc/webdemo/gamma-am/.


  The computational modelling of DNA is becoming crucial in light of new
advances in DNA nanotechnology, single-molecule experiments and in vivo DNA
tampering. Here we present a mesoscopic model for double stranded DNA (dsDNA)
at the single nucleotide level which retains the characteristic helical
structure, while being able to simulate large molecules -- up to a million base
pairs -- for time-scales which are relevant to physiological processes. This is
made possible by an efficient and highly-parallelised implementation of the
model which we discuss here. We compare the behaviour of our model with single
molecule experiments where dsDNA is manipulated by external forces or torques.
We also present some results on the kinetics of denaturation of linear DNA.


  This paper describes an an open-source software system for the automatic
conversion of NLP event representations to system biology structured data
interchange formats such as SBML and BioPAX. It is part of a larger effort to
make results of the NLP community available for system biology pathway
modelers.


  This paper evaluates the difference between human pathway curation and
current NLP systems. We propose graph analysis methods for quantifying the gap
between human curated pathway maps and the output of state-of-the-art automatic
NLP systems. Evaluation is performed on the popular mTOR pathway. Based on
analyzing where current systems perform well and where they fail, we identify
possible avenues for progress.


  We perform differential expression analysis of high-throughput sequencing
count data under a Bayesian nonparametric framework, removing sophisticated
ad-hoc pre-processing steps commonly required in existing algorithms. We
propose to use the gamma (beta) negative binomial process, which takes into
account different sequencing depths using sample-specific negative binomial
probability (dispersion) parameters, to detect differentially expressed genes
by comparing the posterior distributions of gene-specific negative binomial
dispersion (probability) parameters. These model parameters are inferred by
borrowing statistical strength across both the genes and samples. Extensive
experiments on both simulated and real-world RNA sequencing count data show
that the proposed differential expression analysis algorithms clearly
outperform previously proposed ones in terms of the areas under both the
receiver operating characteristic and precision-recall curves.


  Mathematical modelling has become an established tool for studying the
dynamics of biological systems. Current applications range from building models
that reproduce quantitative data to identifying systems with predefined
qualitative features, such as switching behaviour, bistability or oscillations.
Mathematically, the latter question amounts to identifying parameter values
associated with a given qualitative feature.
  We introduce a procedure to partition the parameter space of a parameterized
system of ordinary differential equations into regions for which the system has
a unique or multiple equilibria. The procedure is based on the computation of
the Brouwer degree, and it creates a multivariate polynomial with parameter
depending coefficients. The signs of the coefficients determine parameter
regions with and without multistationarity. A particular strength of the
procedure is the avoidance of numerical analysis and parameter sampling.
  The procedure consists of a number of steps. Each of these steps might be
addressed algorithmically using various computer programs and available
software, or manually. We demonstrate our procedure on several models of gene
transcription and cell signalling, and show that in many cases we obtain a
complete partitioning of the parameter space with respect to multistationarity.


  Studies on how to model the interplay between diseases and behavioral
responses (so-called coupled disease-behavior interaction) have attracted
increasing attention. Owing to the lack of obvious clinical evidence of
diseases, or the incomplete information related to the disease, the risks of
infection cannot be perceived and may lead to inappropriate behavioral
responses. Therefore, how to quantitatively analyze the impacts of asymptomatic
infection on the interplay between diseases and behavioral responses is of
particular importance. In this Letter, under the complex network framework, we
study the coupled disease-behavior interaction model by dividing infectious
individuals into two states: U-state (without evident clinical symptoms,
labelled as U) and I-state (with evident clinical symptoms, labelled as I). A
susceptible individual can be infected by U- or I-nodes, however, since the
U-nodes cannot be easily observed, susceptible individuals take behavioral
responses \emph{only} when they contact I-nodes. The mechanism is considered in
the improved Susceptible-Infected-Susceptible (SIS) model and the improved
Susceptible-Infected-Recovered (SIR) model, respectively. Then, one of the most
concerned problems in spreading dynamics: the epidemic thresholds for the two
models are given by two methods. The analytic results \emph{quantitatively}
describe the influence of different factors, such as asymptomatic infection,
the awareness rate, the network structure, and so forth, on the epidemic
thresholds. Moreover, because of the irreversible process of the SIR model, the
suppression effect of the improved SIR model is weaker than the improved SIS
model.


  We present the open source distributed software package Poisson-Boltzmann
Analytical Method (PB-AM), a fully analytical solution to the linearized
Poisson Boltzmann equation, for molecules represented as non-overlapping
spherical cavities. The PB-AM software package includes the generation of
outputs files appropriate for visualization using VMD, a Brownian dynamics
scheme that uses periodic boundary conditions to simulate dynamics, the ability
to specify docking criteria, and offers two different kinetics schemes to
evaluate biomolecular association rate constants. Given that PB-AM defines
mutual polarization completely and accurately, it can be refactored as a
many-body expansion to explore 2- and 3-body polarization. Additionally, the
software has been integrated into the Adaptive Poisson-Boltzmann Solver (APBS)
software package to make it more accessible to a larger group of scientists,
educators and students that are more familiar with the APBS framework.


  This paper calculates probability distributions modeling the Luria-Delbr\"uck
experiment. We show that by thinking purely in terms of generating functions,
and using a 'backwards in time' paradigm, that formulas describing various
situations can be easily obtained. This includes a generating function for
Haldane's probability distribution due to Ycart. We apply our formulas to both
simulated and real data created by looking at yeast cells acquiring an
immunization to the antibiotic canavanine.
  This paper is somewhat incomplete, having been last significantly modified in
March 29, 2014. However the first author feels that this paper has some
worthwhile ideas, and so is going to make this paper publicly available.


  We present an individual based model of cultural evolution, where interacting
agents are coded by binary strings standing for strategies for action,
blueprints for products or attitudes and beliefs. The model is patterned on an
established model of biological evolution, the Tangled Nature Model (TNM),
where a `tangle' of interactions between agents determines their reproductive
success. In addition, our agents also have the ability to copy part of each
other's strategy, a feature inspired by the Axelrod model of cultural
diversity. Unlike the latter, but similarly to the TNM, the model dynamics goes
through a series of metastable stages of increasing length, each characterized
by mutually enforcing cultural patterns. These patterns are abruptly replaced
by other patterns characteristic of the next metastable period. We analyze the
time dependence of the population and diversity in the system, show how
different cultures are formed and merge, and how their survival probability
lacks, in the model, a finite average life-time. Finally, we use historical
data on the number of car manufacturers after the introduction of the
automobile to the market, to argue that ourmodel can qualitatively reproduce
the flurry of cultural activity which follows a disruptive innovation.


  Akin to electric circuits, we construct biocircuits that are manipulated by
cutting and assembling channels through which stochastic information flows.
This diagrammatic manipulation allows us to create a method which constructs
networks by joining building blocks selected so that (a) they cover only basic
processes; (b) it is scalable to large networks; (c) the mean and
variance-covariance from the Pauli master equation form a closed system and;
(d) given the initial probability distribution, no special boundary conditions
are necessary to solve the master equation. The method aims to help with both
designing new synthetic signalling pathways and quantifying naturally existing
regulatory networks.


  Motivated by questions in biology, we investigate the stability of equilibria
of the dynamical system $\mathbf{x}^{\prime}=P(t)\nabla f(x)$ which arise as
critical points of $f$, under the assumption that $P(t)$ is positive
semi-definite. It is shown that the condition
$\int^{\infty}\lambda_{1}(P(t))~dt=\infty$, where $\lambda_{1}(P(t))$ is the
smallest eigenvalue of $P(t)$, plays a key role in guaranteeing uniform
asymptotic stability and in providing information on the basis of attraction of
those equilibria.


  We consider a system of nonlinear partial differential equations that
describes an age-structured population living in changing environment on $N$
patches. We prove existence and uniqueness of solution and analyze large time
behavior of the system in time-independent case, for periodically changing and
for irregularly varying environment. Under the assumption that every patch can
be reached from every other patch, directly or through several intermediary
patches, and that net reproductive operator has spectral radius larger than
one, we prove that population is persistent on all patches. If the spectral
radius is less or equal one, extinction on all patches is imminent.


  Recent electrophysiological experiments have shown that dopamine (D1)
modulation of pyramidal cells in prefrontal cortex reduces spike frequency
adaptation and enhances NMDA transmission. Using four models, from
multicompartmental to integrate and fire, we examine the effects of these
modulations on sustained (delay) activity in a reverberatory network. We find
that D1 modulation may enable robust network bistability yielding selective
reverberation among cells that code for a particular item or location. We
further show that the tuning curve of such cells is sharpened, and that
signal-to-noise ratio is increased. We postulate that D1 modulation affects the
tuning of "memory fields" and yield efficient distributed dynamic
representations.


  Using an analytically solvable stochastic model, we study the properties of a
simple genetic circuit consisting of multiple copies of an self-regulating
gene. We analyse how the variation in gene copy number and the mutations
changing the auto-regulation strength affect the steady-state distribution of
protein concentration.
  We predict that one-reporter assay, an experimental method where the
extrinsic noise level is inferred from the comparison of expression variance of
a single and duplicated reporter gene, may give an incorrect estimation of the
extrinsic noise contribution when applied to self-regulating genes.
  We also show that an imperfect duplication of an auto-activated gene,
changing the regulation strength of one of the copies, may lead to a hybrid,
binary+graded response of these genes to external signal.
  The analysis of relative changes in mean gene expression before and after
duplication suggests that evolutionary accumulation of gene duplications may
non-trivially depend on the inherent noisiness of a given gene, quantified by
maximal mean frequency of bursts.
  Moreover, we find that the dependence of gene expression noise on gene copy
number and auto-regulation strength may qualitatively differ, e.g. in
monotonicity, depending on whether the noise is measured by Fano factor or
coefficient of variation. Thus, experimentally-based hypotheses linking gene
expression noise and evolutionary optimisation may be ambiguous as they are
dependent on the particular function chosen to quantify noise.


  We discuss fitness landscapes and how they can be modified to account for
co-evolution. We are interested in using the landscape as a way to model
rational decision making in a toy economic system. We develop a model very
similar to the Tangled Nature Model of Christensen et. al. that we call the
Tangled Decision Model. This is a natural setting for our discussion of
co-evolutionary fitness landscapes. We use a Monte Carlo step to simulate
decision making and investigate two different decision making procedures.


  Recently there has been renewed interest in phylogenetic inference methods
based on phylogenetic invariants, alongside the related Markov invariants.
Broadly speaking, both these approaches give rise to polynomial functions of
sequence site patterns that, in expectation value, either vanish for particular
evolutionary trees (in the case of phylogenetic invariants) or have well
understood transformation properties (in the case of Markov invariants). While
both approaches have been valued for their intrinsic mathematical interest, it
is not clear how they relate to each other, and to what extent they can be used
as practical tools for inference of phylogenetic trees.
  In this paper, by focusing on the special case of binary sequence data and
quartets of taxa, we are able to view these two different polynomial-based
approaches within a common framework. To motivate the discussion, we present
three desirable statistical properties that we argue any phylogenetic method
should satisfy: (1) sensible behaviour under reordering of input sequences; (2)
stability as the taxa evolve independently according to a Markov process; and
(3) ability to detect if the conditions of a continuous-time process are
violated. Motivated by these statistical properties, we develop and explore
several new phylogenetic inference methods. In particular, we develop a
statistical bias-corrected version of the Markov invariants approach which
satisfies all three properties. We also extend previous work by showing that
the phylogenetic invariants can be implemented in such a way as to satisfy
property (3). A simulation study shows that, in comparison to other methods,
our new proposed approach based on bias-corrected Markov invariants is
extremely powerful for phylogenetic inference.


  Controlling pest insects is a challenge of main importance to preserve crop
production. In the context of Integrated Pest Management (IPM) programs, we
develop a generic model to study the impact of mating disruption control using
an artificial female pheromone to confuse males and adversely affect their
mating opportunities. Consequently the reproduction rate is diminished leading
to a decline in the population size. For more efficient control, trapping is
used to capture the males attracted to the artificial pheromone. The model,
derived from biological and ecological assumptions, is governed by a system of
ODEs. A theoretical analysis of the model without control is first carried out
to establish the properties of the endemic equilibrium. Then, control is added
and the theoretical analysis of the model enables to identify threshold values
of pheromone which are practically interesting for field applications. In
particular, we show that there is a threshold above which the global asymptotic
stability of the trivial equilibrium is ensured, i.e. the population goes to
extinction. Finally we illustrate the theoretical results via numerical
experiments.


  Information theory provides a useful tool to understand the evolution of
complex nonlinear systems and their sustainability. In particular, Fisher
Information (FI) has been evoked as a useful measure of sustainability and the
variability of dynamical systems including self-organising systems. By
utilising FI, we investigate the sustainability of the logistic model for
different perturbations in the positive and/or negative feedback. Specifically,
we consider different oscillatory modulations in the parameters for positive
and negative feedbacks and investigate their effect on the evolution of the
system and Probability Density Functions (PDFs). Depending on the relative time
scale of the perturbation to the response time of the system (the linear growth
rate), we demonstrate the maintenance of the initial condition for a long time,
manifested by a broad bimodal PDF. We present the analysis of FI in different
cases and elucidate its implications for the sustainability of population
dynamics. We also show that a purely oscillatory growth rate can lead to a
finite amplitude solution while self-organisation of these systems can break
down with an exponentially growing solution due to the fluctuation in negative
feedback.


  We study equilibria in domains with boundaries for a first-order aggregation
model that includes social interactions and exogenous forces. Such equilibrium
solutions can be connected or disconnected, the latter consisting in a delta
concentration on the boundary and a free swarm component in the interior of the
domain. Equilibria are stationary points of an energy functional, and stable
configurations are local minimizers of this functional. We find a one-parameter
family of disconnected equilibrium configurations which are not energy
minimizers; the only stable equilibria are the connected states. Nevertheless,
we demonstrate that in certain cases the dynamical evolution, along the
gradient flow of the energy functional, tends to overwhelmingly favour the
formation of (unstable) disconnected equilibria.


  Empirical observations in marine ecosystems have suggested a balance of
biological and advection time scales as a possible explanation of species
coexistence. To characterise this scenario, we measure the time to fixation in
neutrally evolving populations in chaotic flows. Contrary to intuition the
variation of time scales does not interpolate straightforwardly between the
no-flow and well-mixed limits; instead we find that fixation is the slowest at
intermediate Damk\"ohler numbers, indicating long-lasting coexistence of
species. Our analysis shows that this slowdown is due to spatial organisation
on an increasingly modularised network. We also find that diffusion can either
slow down or speed up fixation, depending on the relative time scales of flow
and evolution.


  The response of lipid bilayers to osmotic stress is an important part of
cellular function. Previously, in [Oglecka et al. 2014], we reported that
cell-sized giant unilamellar vesicles (GUVs) exposed to hypotonic media,
respond to the osmotic assault by undergoing a cyclical sequence of swelling
and bursting events, coupled to the membrane's compositional degrees of
freedom. Here, we seek to deepen our quantitative understanding of the
essential pulsatile behavior of GUVs under hypotonic conditions, by advancing a
comprehensive theoretical model for vesicle dynamics. The model quantitatively
captures our experimentally measured swell-burst parameters for
single-component GUVs, and reveals that thermal fluctuations enable rate
dependent pore nucleation, driving the dynamics of the swell-burst cycles. We
further identify new scaling relationships between the pulsatile dynamics and
GUV properties. Our findings provide a fundamental framework that has the
potential to guide future investigations on the non-equilibrium dynamics of
vesicles under osmotic stress.


  The living-together of distinct organisms in a single termite nest along with
the termite builder colony, is emblematic in its ecological and evolutionary
significance. On top of preserving biodiversity, these interspecific and
intraspecific symbioses provide useful examples of interindividual associations
thought to underly transitions in organic evolution. Being interindividual in
nature, such processes may involve emergent phenomena and hence call for
analytical solutions provided by computing tools and modelling, as opposed to
classical biological methods of analysis. Here we provide selected examples of
such solutions, showing that termite studies may profit from a symbiotic-like
link with computing science to open up wide and new research avenues in ecology
and evolution.


  The majority of mammalian genomic transcripts do not directly code for
proteins and it is currently believed that most of these are not under
evolutionary constraint. However given the abundance non-coding RNA (ncRNA) and
its strong affinity for inter-RNA binding, these molecules have the potential
to regulate proteins in a highly distributed way, similar to artificial neural
networks. We explore this analogy by devising a simple architecture for a
biochemical network that can function as an associative memory. We show that
the steady state solution for this chemical network has the same structure as
an associative memory neural network model. By allowing the choice of
equilibrium constants between different ncRNA species, the concentration of
unbound ncRNA can be made to follow any pattern and many patterns can be stored
simultaneously. The model is studied numerically and within certain parameter
regimes it functions as predicted. Even if the starting concentration pattern
is quite different, it is shown to converge to the original pattern most of the
time. The network is also robust to mutations in equilibrium constants. This
calls into question the criteria for deciding if a sequence is under
evolutionary constraint.


  In this paper, we address the formal characterization of targets triggering
cellular trans-differentiation in the scope of Boolean networks with
asynchronous dynamics. Given two fixed points of a Boolean network, we are
interested in all the combinations of mutations which allow to switch from one
fixed point to the other, either possibly, or inevitably. In the case of
existential reachability, we prove that the set of nodes to (permanently) flip
are only and necessarily in certain connected components of the interaction
graph. In the case of inevitable reachability, we provide an algorithm to
identify a subset of possible solutions.


  We consider (N,r) games of predation with N species and r < N prey and
predators, acting in a cyclic way. Further basic reactions include
reproduction, decay and diffusion over a regular grid, without a hard
constraint on the occupation number per site, so in a "bosonic" implementation.
For special combinations of N and r and appropriate parameter choices we
observe games within games, that is, different coexisting games, depending on
the spatial resolution. As concrete and simplest example we analyze the (6,3)
game. Once the players segregate from a random initial distribution, domains
emerge, which effectively play a (2,1)-game on the coarse scale of domain
diameters, while agents inside the domains play (3,1) (rock-paper-scissors),
leading to spiral formation with species chasing each other. As the (2,1)-game
has a winner in the end, the coexistence of domains is transient, while agents
inside the remaining domain coexist, until demographic fluctuations lead to
extinction of all but one species in the very end. This means that we observe a
dynamical generation of multiple space and time scales with emerging
re-organization of players upon segregation, starting from a simple set of
rules on the smallest scale and changed rules from the coarser perspective.
These observations are based on Gillespie simulations. We discuss the
deterministic limit derived from a van Kampen expansion. The linear stability
analysis predicts the number of forming domains, their composition in terms of
species; it also explains the instability of interfaces between domains, which
drives their extinction; spiral patterns are identified as motion along
heteroclinic cycles. The numerical solutions reproduce the observed patterns of
the Gillespie simulations including even extinction events, so that the
mean-field analysis here is very conclusive, which is due to the specific
implementation of rules.


  This paper uses the recent idea that the fundamental haploid-diploid
lifecycle of eukaryotic organisms implements a rudimentary form of learning
within evolution. A general approach for evolutionary computation is here
derived that differs from all previous known work using diploid
representations. The primary role of recombination is also changed from that
previously considered in both natural and artificial evolution under the new
view. Using well-known abstract tuneable models it is shown that varying
fitness landscape ruggedness varies the benefit of the new approach.


  We are at a stage in our evolution where we do not yet know if we will ever
communicate with intelligent beings that have evolved on other planets, yet we
are intelligent and curious enough to wonder about this. We find ourselves
wondering about this at the very beginning of a long era in which stellar
luminosity warms many planets, and by our best models, continues to provide
equally good opportunities for intelligent life to evolve. By simple Bayesian
reasoning, if, as we believe, intelligent life forms have the same propensity
to evolve later on other planets as we had to evolve on ours, it follows that
they will likely not pass through a similar wondering stage in their evolution.
This suggests that the future holds some kind of interstellar communication
that will serve to inform newly evolved intelligent life forms that they are
not alone before they become curious.


  In this paper, we explain why the chaotic model (CM) of Bahi and Michel
(2008) accurately simulates gene mutations over time. First, we demonstrate
that the CM model is a truly chaotic one, as defined by Devaney. Then, we show
that mutations occurring in gene mutations have the same chaotic dynamic, thus
making the use of chaotic models relevant for genome evolution.


  Among the unsolved problems in computational biology, protein folding is one
of the most interesting challenges. To study this folding, tools like neural
networks and genetic algorithms have received a lot of attention, mainly due to
the NP-completeness of the folding process. The background idea that has given
rise to the use of these algorithms is obviously that the folding process is
predictable. However, this important assumption is disputable as chaotic
properties of such a process have been recently highlighted. In this paper,
which is an extension of a former work accepted to the 2011 International Joint
Conference on Neural Networks (IJCNN11), the topological behavior of a
well-known dynamical system used for protein folding prediction is evaluated.
It is mathematically established that the folding dynamics in the 2D
hydrophobic-hydrophilic (HP) square lattice model, simply called "the 2D model"
in this document, is indeed a chaotic dynamical system as defined by Devaney.
Furthermore, the chaotic behavior of this model is qualitatively and
quantitatively deepened, by studying other mathematical properties of disorder,
namely: the indecomposability, instability, strong transitivity, and constants
of expansivity and sensitivity. Some consequences for both biological paradigms
and structure prediction using this model are then discussed. In particular, it
is shown that some neural networks seems to be unable to predict the evolution
of this model with accuracy, due to its complex behavior.


  Properties of cooperation's probability function in Prisoner`s Dilemma have
impact on evolution of game. Basic model defines that probability of
cooperation depends linearly, both on the player's altruism and the co-player's
reputation. I propose modification of activation function to smooth one
(hyperbolic tangent with scaling parameter a, which corresponds to its shape)
and observe three phases for different range of a. (1) For small a, strategies
seem to randomly change in time and situation of mixed choices (one cooperates
and second defects) dominate. (2) For medium a, players choose only one
strategy for given period of time (the common state can switch to opposite one
with some probability). (3) For large a, mixed strategy (once defect, once
cooperate) is coexisting with common strategies and no change is allowed. I
believe that proposed function characterizes better socio-economical phenomena
and especially phase 1 and 2 contain most of human behavior.


  It is often presumed, that life evolves relatively fast on planets with
clement conditions, at least in its basic forms, and that extended periods of
habitability are subsequently needed for the evolution of higher life forms.
Many planets are however expected to be only transiently habitable. On a large
set of otherwise suitable planets life will therefore just not have the time to
develop on its own to a complexity level as it did arise on earth with the
cambrian explosion. The equivalent of a cambrian explosion may however have the
chance to unfold on transiently habitable planets if it would be possible to
fast forward evolution by 3-4 billion years (with respect to terrestrial
timescales). We argue here, that this is indeed possible when seeding the
candidate planet with the microbial lifeforms, bacteria and unicellular
eukaryotes alike, characterizing earth before the cambrian explosion. An
interstellar mission of this kind, denoted the `Genesis project', could be
carried out by a relatively low-cost robotic microcraft equipped with a
on-board gene laboratory for the in situ synthesis of the microbes.
  We review here our current understanding of the processes determining the
timescales shaping the geo-evolution of an earth-like planet, the prospect of
finding Genesis candidate planets and selected issues regarding the mission
layout. Discussing the ethical aspects connected with a Genesis mission, which
would be expressively not for human benefit, we will also touch the risk that a
biosphere incompatibility may arise in the wake of an eventual manned
exploration of a second earth.


  Various genome evolutionary models have been proposed these last decades to
predict the evolution of a DNA sequence over time, essentially described using
a mutation matrix. By essence, all of these models relate the evolution of DNA
sequences to the computation of the successive powers of the mutation matrix.
To make this computation possible, hypotheses are assumed for the matrix, such
as symmetry and time-reversibility, which are not compatible with mutation
rates that have been recently obtained experimentally on genes ura3 and can1 of
the Yeast Saccharomyces cerevisiae. In this work, authors investigate
systematically the possibility to relax either the symmetry or the
time-reversibility hypothesis of the mutation matrix, by investigating all the
possible matrices of size 2*2 and 3*3. As an application example, the
experimental study on the Yeast Saccharomyces cerevisiae has been used in order
to deduce a simple mutation matrix, and to compute the future evolution of the
rate purine/pyrimidine for $ura3$ on the one hand, and of the particular
behavior of cytosines and thymines compared to purines on the other hand.


  The Notch pathway plays a vital role in determining whether cells in the
intestinal epithelium adopt a secretory or an absorptive phenotype. Cell fate
specification is coordinated via Notch's interaction with the canonical Wnt
pathway. Here, we propose a new mathematical model of the Notch and Wnt
pathways, in which the Hes1 promoter acts as a hub for pathway crosstalk.
Computational simulations of the model can assist in understanding how healthy
intestinal tissue is maintained, and predict the likely consequences of
biochemical knockouts upon cell fate selection processes. Chemical reaction
network theory (CRNT) is a powerful, generalised framework which assesses the
capacity of our model for monostability or multistability, by analysing
properties of the underlying network structure without recourse to specific
parameter values or functional forms for reaction rates. CRNT highlights the
role of beta-catenin in stabilising the Notch pathway and damping oscillations,
demonstrating that Wnt-mediated actions on the Hes1 promoter can induce
dynamical transitions in the Notch system, from multistability to
monostability. Time-dependent model simulations of cell pairs reveal the
stabilising influence of Wnt upon the Notch pathway, in which beta-catenin- and
Dsh-mediated action on the Hes1 promoter are key in shaping the subcellular
dynamics. Where Notch-mediated transcription of Hes1 dominates, there is Notch
oscillation and maintenance of fate flexibility; Wnt-mediated transcription of
Hes1 favours bistability akin to cell fate selection. Cells could therefore
regulate the proportion of Wnt- and Notch-mediated control of the Hes1 promoter
to coordinate the timing of cell fate selection as they migrate through the
intestinal epithelium and are subject to reduced Wnt stimuli.


  We show that the recently introduced logarithmic metrics used to predict
disease arrival times on complex networks are approximations of more general
network-based measures derived from random walks theory. Using the daily
air-traffic transportation data we perform numerical experiments to compare the
infection arrival time with this alternative metric that is obtained by
accounting for multiple walks instead of only the most probable path. The
comparison with direct simulations of arrival times reveals a higher
correlation compared to the shortest path approach used previously. In addition
our method allows to connect fundamental observables in epidemic spreading with
the cumulant generating function of the hitting time for a Markov chain. Our
results provides a general and computationally efficient approach to the
problem using only algebraic methods.


  Malaria remains endemic in tropical areas, especially in Africa. For the
evaluation of new tools and to further ourunderstanding of host-parasite
interactions, knowing the environmental risk of transmission-even at a very
local scale-isessential. The aim of this study was to assess how malaria
transmission is influenced and can be predicted by local climaticand
environmental factors. As the entomological part of a cohort study of 650
newborn babies in nine villages in the ToriBossito district of Southern Benin
between June 2007 and February 2010, human landing catches were performed to
assessthe density of malaria vectors and transmission intensity. Climatic
factors as well as household characteristics were recordedthroughout the study.
Statistical correlations between Anopheles density and environmental and
climatic factors weretested using a three-level Poisson mixed regression model.
The results showed both temporal variations in vector density(related to season
and rainfall), and spatial variations at the level of both village and house.
These spatial variations could belargely explained by factors associated with
the house's immediate surroundings, namely soil type, vegetation index andthe
proximity of a watercourse. Based on these results, a predictive regression
model was developed using a leave-one-outmethod, to predict the spatiotemporal
variability of malaria transmission in the nine villages. This study points up
theimportance of local environmental factors in malaria transmission and
describes a model to predict the transmission risk ofindividual children, based
on environmental and behavioral characteristics.


  Vector-transmitted diseases such as dengue fever and chikungunya have been
spreading rapidly in many parts of the world. The Zika virus has been known
since 1947 and invaded South America in 2013. It can be transmitted not only by
(mosquito) vectors but also directly through sexual contact. Zika has developed
into a serious global health problem because, while most cases are asymptomatic
or very light, babies born to Zika - infected mothers may develop microcephaly
and other very serious birth defects.
  We formulate and analyze two epidemic models for vector-transmitted diseases,
one appropriate for dengue and chikungunya fever outbreaks and one that
includes direct transmission appropriate for Zika virus outbreaks. This is
especially important because the Zika virus is the first example of a disease
that can be spread both indirectly through a vector and directly (through
sexual contact). In both cases, we obtain expressions for the basic
reproduction number and show how to use the initial exponential growth rate to
estimate the basic reproduction number. However, for the model that includes
direct transmission some additional data would be needed to identify the
fraction of cases transmitted directly. Data for the 2015 Zika virus outbreak
in Barranquilla, Colombia has been used to fit parameters to the model
developed here and to estimate the basic reproduction number.


  Herbert Simon's classic rich-get-richer model is one of the simplest
empirically supported mechanisms capable of generating heavy-tail size
distributions for complex systems. Simon argued analytically that a population
of flavored elements growing by either adding a novel element or randomly
replicating an existing one would afford a distribution of group sizes with a
power-law tail. Here, we show that, in fact, Simon's model does not produce a
simple power law size distribution as the initial element has a dominant
first-mover advantage, and will be overrepresented by a factor proportional to
the inverse of the innovation probability. The first group's size discrepancy
cannot be explained away as a transient of the model, and may therefore be many
orders of magnitude greater than expected. We demonstrate how Simon's analysis
was correct but incomplete, and expand our alternate analysis to quantify the
variability of long term rankings for all groups. We find that the expected
time for a first replication is infinite, and show how an incipient group must
break the mechanism to improve their odds of success. We present an example of
citation counts for a specific field that demonstrates a first-mover advantage
consistent with our revised view of the rich-get-richer mechanism. Our findings
call for a reexamination of preceding work invoking Simon's model and provide
an expanded understanding going forward.


  Microbial populations adapt to their environment by acquiring advantageous
mutations, but in the early twentieth century, questions about how these
organisms acquire mutations arose. The experiment of Salvador Luria and Max
Delbr\"uck that won them a Nobel Prize in 1969 confirmed that mutations don't
occur out of necessity, but instead can occur many generations before there is
a selective advantage, and thus organisms follow Darwinian evolution instead of
Lamarckian. Since then, new areas of research involving microbial evolution
have spawned as a result of their experiment. Determining the mutation rate of
a cell is one such area. Probability distributions that determine the number of
mutants in a large population have been derived by D. E. Lea, C. A. Coulson,
and J. B. S. Haldane. However, not much work has been done when time of cell
division is dependent on the cell age, and even less so when cell division is
asymmetric, which is the case in most microbial populations. Using probability
generating function methods, we rigorously construct a probability distribution
for the cell population size given a life-span distribution for both mother and
daughter cells, and then determine its asymptotic growth rate. We use this to
construct a probability distribution for the number of mutants in a large cell
population, which can be used with likelihood methods to estimate the cell
mutation rate.


  The present study, investigated some reproductive and fertility parameters in
some wild populations, originating from the North Tunisia (4 populations of
Medicago ciliaris and 3 populations of Medicago intertexta). Previous finding
revealed that these species are genetically distinct and easily recognized by
the number of flowers per inflorescence and pod dimensions. However,
biometrical traits and isozyme patterns intermediacy between these two species
had detected the existence of a potential spontaneous interspecific hybrid
originating from Sedjnane locality in Tunisia. Indeed, the present work has
shown significant decrease of pollen fertility and seed production for this
population when compared to the others(pollen viability 75%, pollen
germinability 8% and pod production=9%). These results suggested a possible
natural interspecific hybrid and confirming introgressive hybridization
possibility between M. intertexta and M. ciliaris.


  Multiple biological processes are driven by oscillatory gene expression at
different time scales. Pulsatile dynamics are thought to be widespread, and
single-cell live imaging of gene expression has lead to a surge of dynamic,
possibly oscillatory, data for different gene networks. However, the regulation
of gene expression at the level of an individual cell involves reactions
between finite numbers of molecules, and this can result in inherent randomness
in expression dynamics, which blurs the boundaries between aperiodic
fluctuations and noisy oscillators. Thus, there is an acute need for an
objective statistical method for classifying whether an experimentally derived
noisy time series is periodic. Here we present a new data analysis method that
combines mechanistic stochastic modelling with the powerful methods of
non-parametric regression with Gaussian processes. Our method can distinguish
oscillatory gene expression from random fluctuations of non-oscillatory
expression in single-cell time series, despite peak-to-peak variability in
period and amplitude of single-cell oscillations. We show that our method
outperforms the Lomb-Scargle periodogram in successfully classifying cells as
oscillatory or non-oscillatory in data simulated from a simple genetic
oscillator model and in experimental data. Analysis of bioluminescent live cell
imaging shows a significantly greater number of oscillatory cells when
luciferase is driven by a {\it Hes1} promoter (10/19), which has previously
been reported to oscillate, than the constitutive MoMuLV 5' LTR (MMLV) promoter
(0/25). The method can be applied to data from any gene network to both
quantify the proportion of oscillating cells within a population and to measure
the period and quality of oscillations. It is publicly available as a MATLAB
package.


  Stochastic fluctuations of molecule numbers are ubiquitous in biological
systems. Important examples include gene expression and enzymatic processes in
living cells. Such systems are typically modelled as chemical reaction networks
whose dynamics are governed by the Chemical Master Equation. Despite its simple
structure, no analytic solutions to the Chemical Master Equation are known for
most systems. Moreover, stochastic simulations are computationally expensive,
making systematic analysis and statistical inference a challenging task.
Consequently, significant effort has been spent in recent decades on the
development of efficient approximation and inference methods. This article
gives an introduction to basic modelling concepts as well as an overview of
state of the art methods. First, we motivate and introduce deterministic and
stochastic methods for modelling chemical networks, and give an overview of
simulation and exact solution methods. Next, we discuss several approximation
methods, including the chemical Langevin equation, the system size expansion,
moment closure approximations, time-scale separation approximations and hybrid
methods. We discuss their various properties and review recent advances and
remaining challenges for these methods. We present a comparison of several of
these methods by means of a numerical case study and highlight some of their
respective advantages and disadvantages. Finally, we discuss the problem of
inference from experimental data in the Bayesian framework and review recent
methods developed the literature. In summary, this review gives a
self-contained introduction to modelling, approximations and inference methods
for stochastic chemical kinetics.


  Computational docking is the core process of computer-aided drug design; it
aims at predicting the best orientation and conformation of a small drug
molecule when bound to a target large protein receptor. The docking quality is
typically measured by a scoring function: a mathematical predictive model that
produces a score representing the binding free energy and hence the stability
of the resulting complex molecule. We analyze the performance of both learning
techniques on the scoring power, the ranking power, docking power, and
screening power using the PDBbind 2013 database. For the scoring and ranking
powers, the proposed learning scoring functions depend on a wide range of
features (energy terms, pharmacophore, intermolecular) that entirely
characterize the protein-ligand complexes. For the docking and screening
powers, the proposed learning scoring functions depend on the intermolecular
features of the RF-Score to utilize a larger number of training complexes. For
the scoring power, the DL\_RF scoring function achieves Pearson's correlation
coefficient between the predicted and experimentally measured binding
affinities of 0.799 versus 0.758 of the RF scoring function. For the ranking
power, the DL scoring function ranks the ligands bound to fixed target protein
with accuracy 54% for the high-level ranking and with accuracy 78% for the
low-level ranking while the RF scoring function achieves (46% and 62%)
respectively. For the docking power, the DL\_RF scoring function has a success
rate when the three best-scored ligand binding poses are considered within 2
\AA\ root-mean-square-deviation from the native pose of 36.0% versus 30.2% of
the RF scoring function. For the screening power, the DL scoring function has
an average enrichment factor and success rate at the top 1% level of (2.69 and
6.45%) respectively versus (1.61 and 4.84%) respectively of the RF scoring
function.


  Stochastic compartmental models are important tools for understanding the
course of infectious diseases epidemics in populations and in prospective
evaluation of intervention policies. However, calculating the likelihood for
discretely observed data from even simple models -- such as the ubiquitous
susceptible-infectious-removed (SIR) model -- has been considered
computationally intractable, since its formulation almost a century ago.
Recently researchers have proposed methods to circumvent this limitation
through data augmentation or approximation, but these approaches often suffer
from high computational cost or loss of accuracy. We develop the mathematical
foundation and an efficient algorithm to compute the likelihood for discretely
observed data from a broad class of stochastic compartmental models. We also
give expressions for the derivatives of the transition probabilities using the
same technique, making possible inference via Hamiltonian Monte Carlo (HMC). We
use the 17th century plague in Eyam, a classic example of the SIR model, to
compare our recursion method to sequential Monte Carlo, analyze using HMC, and
assess the model assumptions. We also apply our direct likelihood evaluation to
perform Bayesian inference for the 2014-2015 Ebola outbreak in Guinea. The
results suggest that the epidemic infectious rates have decreased since October
2014 in the Southeast region of Guinea, while rates remain the same in other
regions, facilitating understanding of the outbreak and the effectiveness of
Ebola control interventions.


  Natural protein sequences that self-assemble to form globular structures are
compact with high packing densities in the folded states. It is known that
proteins unfold upon addition of denaturants, adopting random coil structures.
The dependence of the radii of gyration on protein size in the folded and
unfolded states obeys the same scaling laws as synthetic polymers. Thus, one
might surmise that the mechanism of collapse in proteins and polymers ought to
be similar. However, because the number of amino acids in single domain
proteins is not significantly greater than about two hundred, it has not been
resolved if the unfolded states of proteins are compact under conditions that
favor the folded states - a problem at the heart of how proteins fold. By
adopting a theory used to derive polymer-scaling laws, we find that the
propensity for the unfolded state of a protein to be compact is universal and
is encoded in the contact map of the folded state. Remarkably, analysis of over
2000 proteins shows that proteins rich in $\beta$-sheets have greater tendency
to be compact than $\alpha$-helical proteins. The theory provides insights into
the reasons for the small size of single domain proteins and the physical basis
for the origin of multi-domain proteins. Application to non-coding RNA
molecules show that they have evolved to collapse sharing similarities to
$\beta$-sheet proteins. An implication of our theory is that the evolution of
natural foldable sequences is guided by the requirement that for efficient
folding they should populate minimum energy compact states under folding
conditions. This concept also supports the compaction selection hypothesis used
to rationalize the unusually condensed states of viral RNA molecules.


  In order to interpret recent experimental studies of the dependence of
conductance of ds-DNA as the DNA is pulled from the 3'end1-3'end2 ends, which
find a sharp conductance jump for a very short (4.5 %) stretching length, we
carried out multiscale modeling, to predict the conductance of dsDNA as it is
mechanically stretched to promote various structural polymorphisms. We
calculate the current along the stretched DNA using a combination of molecular
dynamics simulations, non-equilibrium pulling simulations, quantum mechanics
calculations, and kinetic Monte Carlo simulations. For 5'end1-5'end2
attachments we find an abrupt jump in the current within a very short
stretching length (6 $ \AA $ or 17 %) leading to a melted DNA state. In
contrast, for 3'end1-3'end2 pulling it takes almost 32$ \AA $ (84 %) of
stretching to cause a similar jump in the current. Thus, we demonstrate that
charge transport in DNA can occur over stretching lengths of several
nanometers. We find that this unexpected behaviour in the B to S conformational
DNA transition arises from highly inclined base pair geometries that result
from this pulling protocol. We find that the dramatically different conductance
behaviors for two different pulling protocols arise from the nature of how
hydrogen bonds of DNA base pairs break.


  The ability to reason beyond established knowledge allows Organic Chemists to
solve synthetic problems and to invent novel transformations. Here, we propose
a model which mimics chemical reasoning and formalises reaction prediction as
finding missing links in a knowledge graph. We have constructed a knowledge
graph containing 14.4 million molecules and 8.2 million binary reactions, which
represents the bulk of all chemical reactions ever published in the scientific
literature. Our model outperforms a rule-based expert system in the reaction
prediction task for 180,000 randomly selected binary reactions. We show that
our data-driven model generalises even beyond known reaction types, and is thus
capable of effectively (re-) discovering novel transformations (even including
transition-metal catalysed reactions). Our model enables computers to infer
hypotheses about reactivity and reactions by only considering the intrinsic
local structure of the graph, and because each single reaction prediction is
typically achieved in a sub-second time frame, our model can be used as a
high-throughput generator of reaction hypotheses for reaction discovery.


  We study a stochastic model of gene expression, in which protein production
has a form of random bursts whose size distribution is arbitrary, whereas
protein decay is a first-order reaction. We find exact analytical expressions
for the time evolution of the cumulant-generating function for the most general
case when both the burst size probability distribution and the model parameters
depend on time in an arbitrary (e.g. oscillatory) manner, and for arbitrary
initial conditions. We show that in the case of periodic external activation
and constant protein degradation rate, the response of the gene is analogous to
the RC low-pass filter, where slow oscillations of the external driving have a
greater effect on gene expression than the fast ones. We also demonstrate that
the $n$-th cumulant of the protein number distribution depends on the $n$-th
moment of the burst size distribution. We use these results to show that
different measures of noise (coefficient of variation, Fano factor, fractional
change of variance) may vary in time in a different manner. Therefore, any
biological hypothesis of evolutionary optimization based on the nonmonotonicity
of a chosen measure of noise must justify why it assumes that biological
evolution quantifies noise in that particular way. Finally, we show that not
only for exponentially distributed burst sizes but also for a wider class of
burst size distributions (e.g. Dirac delta and gamma) the control of gene
expression level by burst frequency modulation gives rise to proportional
scaling of variance of the protein number distribution to its mean, whereas the
control by amplitude modulation implies proportionality of protein number
variance to the mean squared.


  We propose a model solely based on actin treadmilling and polymerisation
which describes many characteristic states of actin wave formation: spots,
spirals and travelling waves. In our model, as in experiments on cell
recovering motility following actin depolymerisation, we choose an isotropic
low density initial condition; polymerisation of actin filaments then raises
the density towards the Onsager threshold where they align. We show that this
alignment, in turn, destabilizes the isotropic phase and generically induces
transient actin spots or spirals as part of the dynamical pathway towards a
polarized phase which can either be uniform or consist of a series of
actin-wave trains (flocks). Our results uncover a universal route to actin wave
formation in the absence of any system specific nonlinear biochemistry, and it
may help understand the mechanism underlying the observation of actin spots and
waves in vivo. They also suggest a minimal setup to design similar patterns in
vitro.


  Phosphorylation, the enzyme-mediated addition of a phosphate group to a
molecule, is a ubiquitous chemical mechanism in biology. Multisite
phosphorylation, the addition of phosphate groups to multiple sites of a single
molecule, may be distributive or processive. Distributive systems can be
bistable, while processive systems were recently shown to be globally stable.
However, this global convergence result was proven only for a specific
mechanism of processive phosphorylation/dephosphorylation (namely, all
catalytic reactions are reversible). Accordingly, we generalize this result to
allow for processive phosphorylation networks in which each reaction may be
irreversible, and also to account for possible product inhibition. We
accomplish this by defining an all-encompassing processive network that
encapsulates all of these schemes, and then appealing to recent results of
Marcondes de Freitas, Wiuf, and Feliu that assert global convergence by way of
monontone systems theory and network/graph reductions (which correspond to
removal of intermediate complexes). Our results form a case study into the
question of when global convergence is preserved when reactions and/or
intermediate complexes are added to or removed from a network.


  In nematodes, environmental or physiological perturbations alter death's
scaling of time. In human cancer, genetic perturbations alter death's curvature
of time. Those changes in scale and curvature follow the constraining contours
of death's invariant geometry. I show that the constraints arise from a
fundamental extension to the theories of randomness, invariance and scale. A
generalized Gompertz law follows. The constraints imposed by the invariant
Gompertz geometry explain the tendency of perturbations to stretch or bend
death's scaling of time. Variability in death rate arises from a combination of
constraining universal laws and particular biological processes.