  When designing genetic circuits, the typical primitives used in major
existing modelling formalisms are gene interaction graphs, where edges between
genes denote either an activation or inhibition relation. However, when
designing experiments, it is important to be precise about the low-level
mechanistic details as to how each such relation is implemented. The rule-based
modelling language Kappa allows to unambiguously specify mechanistic details
such as DNA binding sites, dimerisation of transcription factors, or
co-operative interactions. However, such a detailed description comes with
complexity and computationally costly execution. We propose a general method
for automatically transforming a rule-based program, by eliminating
intermediate species and adjusting the rate constants accordingly. Our method
consists of searching for those interaction patterns known to be amenable to
equilibrium approximations (e.g. Michaelis-Menten scheme). The reduced model is
efficiently obtained by static inspection over the rule-set, and it represents
a particular theoretical limit of the original model. The Bhattacharyya
distance is proposed as a metric to estimate the reduction error for a given
observable. The tool is tested on a detailed rule-based model of a
$\lambda$-phage switch, which lists $96$ rules and $16$ agents. The reduced
model has $11$ rules and $5$ agents, and provides a dramatic reduction in
simulation time of several orders of magnitude.


  A quantitative understanding of organism-level behavior requires predictive
models that can capture the richness of behavioral phenotypes, yet are simple
enough to connect with underlying mechanistic processes. Here we investigate
the motile behavior of nematodes at the level of their translational motion on
surfaces driven by undulatory propulsion. We broadly sample the nematode
behavioral repertoire by measuring motile trajectories of the canonical lab
strain $C. elegans$ N2 as well as wild strains and distant species. We focus on
trajectory dynamics over timescales spanning the transition from ballistic
(straight) to diffusive (random) movement and find that salient features of the
motility statistics are captured by a random walk model with independent
dynamics in the speed, bearing and reversal events. We show that the model
parameters vary among species in a correlated, low-dimensional manner
suggestive of a common mode of behavioral control and a trade-off between
exploration and exploitation. The distribution of phenotypes along this primary
mode of variation reveals that not only the mean but also the variance varies
considerably across strains, suggesting that these nematode lineages employ
contrasting ``bet-hedging'' strategies for foraging.


  The human T cell repertoire is generated by the rearrangement of variable
(V), diversity (D) and joining (J) segments on the T cell receptor (TCR) loci.
To determine whether the structural ordering of these gene segments on the TCR
loci contributes to the observed clonal frequencies, the TCR loci were examined
for self-similarity and periodicity in terms of gene segment organization.
Logarithmic transformation of numeric sequence order demonstrated that the V
and J gene segments for both T cell receptor alpha (TRA) and beta (TRB) loci
were arranged in a self-similar manner when the spacing between adjacent
segments was considered as a function of the size of the neighboring gene
segment. The ratio of genomic distance between either the J (in TRA) or D (in
TRB) segments and successive V segments on these loci declined logarithmically.
Accounting for the gene segments occurring on helical DNA molecules, in a
logarithmic distribution, sine and cosine functions of the log transformed
angular coordinates of the start and stop nucleotides of successive TCR gene
segments showed an ordered progression across the locus, supporting a
log-periodic organization. T cell clonal frequencies, based on V and J segment
usage, from three normal stem cell donors plotted against the respective
segment locations on TRB locus demonstrated a periodic variation. We
hypothesize that this quasi-periodic variation in T cell clonal repertoire may
be influenced by the location of the gene segments on the logarithmically
scaled TCR loci. Interactions between the two strands of DNA in the double
helix may influence the probability of gene segment usage by means of either
constructive or destructive interference resulting from the superposition of
the two helices, impacting probability of DNA recombination.


  Intracellular cargos which are transported by molecular motors move
stochastically along cytoskeleton filaments. In particular for bidirectionally
transported cargos it is an open question whether the characteristics of their
motion can result from pure stochastic fluctuations or whether some
coordination of the motors is needed. The results of a mean-field model of
cargo-motors dynamics, which was proposed by M\"uller et al.[1] suggest the
existence of high motility states which would result from a stochastic
tug-of-war. Here we analyze a non-mean field extension of their model, that
takes explicitly the position of each motor into account. We find that high
motility states then disappear. We consider also a mutual motor-motor
activation, as an explicit mechanism of motor coordination. We show that the
results of the mean-field model are recovered only in case of a strong
motor-motor activation in the limit of a high number of motors.


  Under models of isolation-by-distance, population structure is determined by
the probability of identity-by-descent between pairs of genes according to the
geographic distance between them. Well established analytical results indicate
that the relationship between geographical and genetic distance depends mostly
on the neighborhood size of the population, $N_b = 4{\pi}{\sigma}^2 D_e$, which
represents a standardized measure of dispersal. To test this prediction, we
model local dispersal of haploid individuals on a two-dimensional torus using
four dispersal kernels: Rayleigh, exponential, half-normal and triangular. When
neighborhood size is held constant, the distributions produce similar patterns
of isolation-by-distance, confirming predictions. Considering this, we propose
that the triangular distribution is the appropriate null distribution for
isolation-by-distance studies. Under the triangular distribution, dispersal is
uniform within an area of $4{\pi}{\sigma}^2$ (i.e. the neighborhood area),
which suggests that the common description of neighborhood size as a measure of
a local panmictic population is valid for popular families of dispersal
distributions. We further show how to draw from the triangular distribution
efficiently and argue that it should be utilized in other studies in which
computational efficiency is important


  We model the evolution of eukaryotic protein-protein interaction (PPI)
networks. In our model, PPI networks evolve by two known biological mechanisms:
(1) Gene duplication, which is followed by rapid diversification of duplicate
interactions. (2) Neofunctionalization, in which a mutation leads to a new
interaction with some other protein. Since many interactions are due to simple
surface compatibility, we hypothesize there is an increased likelihood of
interacting with other proteins in the target protein's neighborhood. We find
good agreement of the model on 10 different network properties compared to
high-confidence experimental PPI networks in yeast, fruit flies, and humans.
Key findings are: (1) PPI networks evolve modular structures, with no need to
invoke particular selection pressures. (2) Proteins in cells have on average
about 6 degrees of separation, similar to some social networks, such as
human-communication and actor networks. (3) Unlike social networks, which have
a shrinking diameter (degree of maximum separation) over time, PPI networks are
predicted to grow in diameter. (4) The model indicates that evolutionarily old
proteins should have higher connectivities and be more centrally embedded in
their networks. This suggests a way in which present-day proteomics data could
provide insights into biological evolution.


  Superoxide reductases (SORs) are superoxide (O2-)-detoxifying enzymes that
catalyse the reduction of O2- into hydrogen peroxide. Three different classes
of SOR have been reported on the basis of the presence or not of an additional
N-terminal domain. They all share a similar active site, with an unusual
non-heme Fe atom coordinated by four equatorial histidines and one axial
cysteine residues. Crucial catalytic reaction intermediates of SOR are
purported to be Fe(3+)-(hydro)peroxo species. Using resonance Raman
spectroscopy, we compared the vibrational properties of the Fe3+ active site of
two different classes of SOR, from Desulfoarculus baarsii and Treponema
pallidum, along with their ferrocyanide and their peroxo complexes. In both
species, rapid treatment with H2O2 results in the stabilization of a side-on
high spin Fe(3+)-(eta(2)-OO) peroxo species. Comparison of these two peroxo
species reveals significant differences in vibrational frequencies and bond
strengths of the Fe-O2 (weaker) and O-O (stronger) bonds for the T. pallidum
enzyme. Thus, the two peroxo adducts in these two SORs have different
stabilities which are also seen to be correlated with differences in the Fe-S
coordination strengths as gauged by the Fe-S vibrational frequencies. This was
interpreted from structural variations in the two active sites, resulting in
differences in the electron donating properties of the trans cysteine ligand.
Our results suggest that the structural differences observed in the active site
of different classes of SORs should be a determining factor for the rate of
release of the iron-peroxo intermediate during enzymatic turnover.


  The superoxide radical O(2)(-.) is a toxic by-product of oxygen metabolism.
Two O(2)(-.) detoxifying enzymes have been described so far, superoxide
dismutase and superoxide reductase (SOR), both forming H2O2 as a reaction
product. Recently, the SOR active site, a ferrous iron in a [Fe(2+) (N-His)(4)
(S-Cys)] pentacoordination, was shown to have the ability to form a complex
with the organometallic compound ferrocyanide. Here, we have investigated in
detail the reactivity of the SOR-ferrocyanide complex with O(2)(-.) by pulse
and gamma-ray radiolysis, infrared, and UV-visible spectroscopies. The complex
reacts very efficiently with O(2)(-.). However, the presence of the
ferrocyanide adduct markedly modifies the reaction mechanism of SOR, with the
formation of transient intermediates different from those observed for SOR
alone. A one-electron redox chemistry appears to be carried out by the
ferrocyanide moiety of the complex, whereas the SOR iron site remains in the
reduced state. Surprisingly, the toxic H2O2 species is no longer the reaction
product. Accordingly, in vivo experiments showed that formation of the
SOR-ferrocyanide complex increased the antioxidant capabilities of SOR
expressed in an Escherichia coli sodA sodB recA mutant strain. Altogether,
these data describe an unprecedented O(2)(-.) detoxification activity,
catalyzed by the SOR-ferrocyanide complex, which does not conduct to the
production of the toxic H2O2 species.


  Superoxide reductase SOR is a non-heme iron metalloenzyme that detoxifies
superoxide radical in microorganisms. Its active site consists of an unusual
non-heme Fe2+ center in a [His4 Cys1] square pyramidal pentacoordination, with
the axial cysteine ligand proposed to be an essential feature in catalysis. Two
NH peptide groups from isoleucine 118 and histidine 119 establish H-bondings
with the sulfur ligand (Desulfoarculus baarsii SOR numbering). In order to
investigate the catalytic role of these H-bonds, the isoleucine 118 residue of
the SOR from Desulfoarculus baarsii was mutated into alanine, aspartate or
serine residues. Resonance Raman spectroscopy showed that the mutations
specifically induced an increase of the strength of the Fe3+-S(Cys) and
S-C$\beta$(Cys) bonds as well as a change in conformation of the cysteinyl side
chain, which was associated with the alteration of the NH hydrogen bonding to
the sulfur ligand. The effects of the isoleucine mutations on the reactivity of
SOR with $O2\bullet$- were investigated by pulse radiolysis. These studies
showed that the mutations induced a specific increase of the pKa of the first
reaction intermediate, recently proposed to be an $Fe2+-O2\bullet-$ species.
These data were supported by DFT calculations carried out on three models of
the $Fe2+-O2\bullet-$ intermediate, with one, two or no H-bonds on the sulfur
ligand. Our results demonstrated that the hydrogen bonds between the NH
(peptide) and the cysteine ligand tightly control the rate of protonation of
the $Fe2+-O2\bullet-$ reaction intermediate to form an Fe3+-OOH species.


  The NAD(P)H:flavin oxidoreductase from Escherichia coli, named Fre, is a
monomer of 26.2 kDa that catalyzes the reduction of free flavins using NADPH or
NADH as electron donor. The enzyme does not contain any prosthetic group but
accommodates both the reduced pyridine nucleotide and the flavin in a ternary
complex prior to oxidoreduction. The specificity of the flavin reductase for
the pyridine nucleotide was studied by steady-state kinetics using a variety of
NADP analogs. Both the nicotinamide ring and the adenosine part of the
substrate molecule have been found to be important for binding to the
polypeptide chain. However, in the case of NADPH, the 2'-phosphate group
destabilized almost completely the interaction with the adenosine moiety.
Moreover, NADPH and NMNH are very good substrates for the flavin reductase, and
we have shown that both these molecules bind to the enzyme almost exclusively
by the nicotinamide ring. This provides evidence that the flavin reductase
exhibits a unique mode for recognition of the reduced pyridine nucleotide. In
addition, we have shown that the flavin reductase selectively transfers the
pro-R hydrogen from the C-4 position of the nicotinamide ring and is therefore
classified as an A-side-specific enzyme.


  Superoxide reductase (SOR) is a small metalloenzyme that catalyzes reduction
of O(2)(*)(-) to H(2)O(2) and thus provides an antioxidant mechanism against
superoxide radicals. Its active site contains an unusual mononuclear ferrous
center, which is very efficient during electron transfer to O(2)(*)(-)
[Lombard, M., Fontecave, M., Touati, D., and Nivi{\`e}re, V. (2000) J. Biol.
Chem. 275, 115-121]. The reaction of the enzyme from Desulfoarculus baarsii
with superoxide was studied by pulse radiolysis methods. The first step is an
extremely fast bimolecular reaction of superoxide reductase with superoxide,
with a rate constant of (1.1 +/- 0.3) x 10(9) M(-1) s(-1). A first intermediate
is formed which is converted to a second one at a much slower rate constant of
500 +/- 50 s(-1). Decay of the second intermediate occurs with a rate constant
of 25 +/- 5 s(-1). These intermediates are suggested to be iron-superoxide and
iron-peroxide species. Furthermore, the role of glutamate 47 and lysine 48,
which are the closest charged residues to the vacant sixth iron coordination
site, has been investigated by site-directed mutagenesis. Mutation of glutamate
47 into alanine has no effect on the rates of the reaction. On the contrary,
mutation of lysine 48 into an isoleucine led to a 20-30-fold decrease of the
rate constant of the bimolecular reaction, suggesting that lysine 48 plays an
important role during guiding and binding of superoxide to the iron center II.
In addition, we report that expression of the lysine 48 sor mutant gene hardly
restored to a superoxide dismutase-deficient Escherichia coli mutant the
ability to grow under aerobic conditions.


  Superoxide reductases (SORs) are small metalloenzymes, which catalyze
reduction of O2*- to H2O2. The reaction of the enzyme from Treponema pallidum
with superoxide was studied by pulse radiolysis methods. The first step is an
extremely fast bi-molecular reaction of the ferrous center with O2, with a rate
constant of 6 x 10 (8) M(-1) s(-1). A first intermediate is formed which is
converted to a second one with a slower rate constant of 4800 s(-1). This
latter value is 10 times higher than the corresponding one previously reported
in the case of SOR from Desulfoarculus baarsii. The reconstituted spectra for
the two intermediates are consistent with formation of transient iron-peroxide
species.


  Background: Protein-protein interaction (PPI) network analyses are highly
valuable in deciphering and understanding the intricate organisation of
cellular functions. Nevertheless, the majority of available protein-protein
interaction networks are context-less, i.e. without any reference to the
spatial, temporal or physiological conditions in which the interactions may
occur. In this work, we are proposing a protocol to infer the most likely
protein-protein interaction (PPI) network in human macrophages. Results: We
integrated the PPI dataset from the Agile Protein Interaction DataAnalyzer
(APID) with different meta-data to infer a contextualized macrophage-specific
interactome using a combination of statistical methods. The obtained
interactome is enriched in experimentally verified interactions and in proteins
involved in macrophage-related biological processes (i.e. immune response
activation, regulation of apoptosis). As a case study, we used the
contextualized interactome to highlight the cellular processes induced upon
Mycobacterium tuberculosis infection. Conclusion: Our work confirms that
contextualizing interactomes improves the biological significance of
bioinformatic analyses. More specifically, studying such inferred network
rather than focusing at the gene expression level only, is informative on the
processes involved in the host response. Indeed, important immune features such
as apoptosis are solely highlighted when the spotlight is on the protein
interaction level.


  The flavin reductase ActVB is involved in the last step of actinorhodin
biosynthesis in Streptomyces coelicolor. Although ActVB can be isolated with
some FMN bound, this form was not involved in the flavin reductase activity. By
studying the ferric reductase activity of ActVB, we show that its FMN-bound
form exhibits a proper enzymatic activity of reduction of iron complexes by
NADH. This shows that ActVB active site exhibits a dual property with regard to
the FMN. It can use it as a substrate that goes in and off the active site or
as a cofactor to provide an electron transferase activity to the polypeptide.


  The two-component flavin-dependent monooxygenases belong to an emerging class
of enzymes involved in oxidation reactions in a number of metabolic and
biosynthetic pathways in microorganisms. One component is a NAD(P)H:flavin
oxidoreductase, which provides a reduced flavin to the second component, the
proper monooxygenase. There, the reduced flavin activates molecular oxygen for
substrate oxidation. Here, we study the flavin reductase ActVB and ActVA-ORF5
gene product, both reported to be involved in the last step of biosynthesis of
the natural antibiotic actinorhodin in Streptomyces coelicolor. For the first
time we show that ActVA-ORF5 is a FMN-dependent monooxygenase that together
with the help of the flavin reductase ActVB catalyzes the oxidation reaction.
The mechanism of the transfer of reduced FMN between ActVB and ActVA-ORF5 has
been investigated. Dissociation constant values for oxidized and reduced flavin
(FMNox and FMNred) with regard to ActVB and ActVA-ORF5 have been determined.
The data clearly demonstrate a thermodynamic transfer of FMNred from ActVB to
ActVA-ORF5 without involving a particular interaction between the two protein
components. In full agreement with these data, we propose a reaction mechanism
in which FMNox binds to ActVB, where it is reduced, and the resulting FMNred
moves to ActVA-ORF5, where it reacts with O2 to generate a flavinperoxide
intermediate. A direct spectroscopic evidence for the formation of such species
within ActVA-ORF5 is reported.


  The active site of superoxide reductase SOR consists of an Fe2+ center in an
unusual [His4 Cys1] square-pyramidal geometry. It specifically reduces
superoxide to produce H2O2. Here, we have reacted the SOR from Desulfoarculus
baarsii directly with H2O2. We have found that its active site can transiently
stabilize an Fe3+-peroxo species that we have spectroscopically characterized
by resonance Raman. The mutation of the strictly conserved Glu47 into alanine
results in a stabilization of this Fe3+-peroxo species, when compared to the
wild-type form. These data support the hypothesis that the reaction of SOR
proceeds through such a Fe3+-peroxo intermediate. This also suggests that Glu47
might serve to help H2O2 release during the reaction with superoxide.


  Some sulfate-reducing and microaerophilic bacteria rely on the enzyme
superoxide reductase (SOR) to eliminate the toxic superoxide anion radical
(O2*-). SOR catalyses the one-electron reduction of O2*- to hydrogen peroxide
at a nonheme ferrous iron center. The structures of Desulfoarculus baarsii SOR
(mutant E47A) alone and in complex with ferrocyanide were solved to 1.15 and
1.7 A resolution, respectively. The latter structure, the first ever reported
of a complex between ferrocyanide and a protein, reveals that this
organo-metallic compound entirely plugs the SOR active site, coordinating the
active iron through a bent cyano bridge. The subtle structural differences
between the mixed-valence and the fully reduced SOR-ferrocyanide adducts were
investigated by taking advantage of the photoelectrons induced by X-rays. The
results reveal that photo-reduction from Fe(III) to Fe(II) of the iron center,
a very rapid process under a powerful synchrotron beam, induces an expansion of
the SOR active site.


  A complex contagion is an infectious process in which individuals may require
multiple transmissions before changing state. These are used to model behaviors
if an individual only adopts a particular behavior after perceiving a consensus
among others. We may think of individuals as beginning inactive and becoming
active after contact with a sufficient number of active partners. These have
been studied in a number of cases, but analytic models for the dynamic spread
of complex contagions are typically complex. Here we study the dynamics of the
Watts Threshold Model (WTM) assuming transmission occurs in continuous time as
a Poisson process, or in discrete time where individuals transmit to all
partners in the time step following their activation. We adapt techniques
developed for infectious disease modeling to develop and analyze analytic
models for the dynamics of the WTM in Configuration Model networks and a class
of random clustered (triangle-based) networks. The resulting model is
relatively simple and compact. We use it to gain insights into the contagion
dynamics. In the infinite population limit, we derive conditions under which
cascades happen with an arbitrarily small initial proportion active, confirming
a hypothesis of Watts for this case. We also observe hybrid phase transitions
when cascades are not possible for small initial conditions, but occur for
large enough initial conditions. We derive sufficient conditions for this
hybrid phase transition to occur. We show that in many cases, if the hybrid
phase transition occurs, all individuals eventually become active. Finally, we
discuss the role clustering plays in facilitating or impeding the spread and
find that the hypothesis of Watts that was confirmed in Configuration Model
networks does not hold in general. This approach allows us to unify many
existing disparate observations and derive new results.


  We study inhomogeneous host-pathogen dynamics to model the global amphibian
population extinction in a lake basin system. The lake basin system is modeled
as quenched disorder. In this model we show that once the pathogen arrives at
the lake basin it spreads from one lake to another, eventually spreading to the
entire lake basin system in a wave like pattern. The extinction time has been
found to depend on the steady state host population and pathogen growth rate.
Linear estimate of the extinction time is computed. The steady state host
population shows a threshold behavior in the interaction strength for a given
growth rate.


  Cells must coordinate adjustments in genome expression to accommodate changes
in their environment. We hypothesized that the amount of transcriptome change
is proportional to the amount of environmental change. To capture the effects
of environmental changes on the transcriptome, we compared transcriptome
diversities (defined as the Shannon entropy of frequency distribution) of
silkworm fat-body tissues cultured with several concentrations of
phenobarbital. Although there was no proportional relationship, we did identify
a drug concentration tipping point between 0.25 and 1.0 mM. Cells cultured in
media containing lower drug concentrations than the tipping point showed
uniformly high transcriptome diversities, while those cultured at higher drug
concentrations than the tipping point showed uniformly low transcriptome
diversities. The plasticity of transcriptome diversity was corroborated by
cultivations of fat bodies in MGM-450 insect medium without phenobarbital and
in 0.25 mM phenobarbital-supplemented MGM-450 insect medium after previous
cultivation (cultivation for 80 hours in MGM-450 insect medium without
phenobarbital, followed by cultivation for 10 hours in 1.0 mM
phenobarbital-supplemented MGM-450 insect medium). Interestingly, the
transcriptome diversities of cells cultured in media containing 0.25 mM
phenobarbital after previous cultivation (cultivation for 80 hours in MGM-450
insect medium without phenobarbital, followed by cultivation for 10 hours in
1.0 mM phenobarbital-supplemented MGM-450 insect medium) were different from
cells cultured in media containing 0.25 mM phenobarbital after previous
cultivation (cultivation for 80 hours in MGM-450 insect medium without
phenobarbital). This hysteretic phenomenon of transcriptome diversities
indicates multi-stability of the genome expression system.


  We propose a metric which can be used to compute the amount of heritable
variation enabled by a given dynamical system. A distribution of selection
pressures is used such that each pressure selects a particular fixed point via
competitive exclusion in order to determine the corresponding distribution of
potential fixed points in the population dynamics. This metric accurately
detects the number of species present in artificially prepared test systems,
and furthermore can correctly determine the number of heritable sets in
clustered transition matrix models in which there are no clearly defined
genomes. Finally, we apply our metric to the GARD model and show that it
accurately reproduces prior measurements of the model's heritability.


  We consider an infinitely large population under stabilising selection and
mutation in which the allelic effects determining a polygenic trait vary
between loci. We obtain analytical expressions for the stationary genetic
variance as a function of the distribution of effects, mutation rate and
selection coefficient. We also study the dynamics of the allele frequencies,
focussing on short-term evolution of the phenotypic mean as it approaches the
optimum after an environmental change. We find that when most effects are
small, the genetic variance does not change appreciably during adaptation, and
the time until the phenotypic mean reaches the optimum is short if the number
of loci is large. However, when most effects are large, the change of the
variance during the adaptive process cannot be neglected. In this case, the
short-term dynamics may be described by that of a single locus of large effect.
Our results may be used to understand polygenic selection driving rapid
adaptation.


  A large number of studies have shown the existence of metabolic covalent
modifications in different molecular structures, able to store biochemical
information that is not encoded by the DNA. Some of these covalent mark
patterns can be transmitted across generations (epigenetic changes). Recently,
the emergence of Hopfield-like attractor dynamics has been observed in the
self-organized enzymatic networks, which have the capacity to store functional
catalytic patterns that can be correctly recovered by the specific input
stimuli. The Hopfield-like metabolic dynamics are stable and can be maintained
as a long-term biochemical memory. In addition, specific molecular information
can be transferred from the functional dynamics of the metabolic networks to
the enzymatic activity involved in the covalent post-translational modulation
so that determined functional memory can be embedded in multiple stable
molecular marks. Both the metabolic dynamics governed by Hopfield-type
attractors (functional processes) and the enzymatic covalent modifications of
determined molecules (structural dynamic processes) seem to represent the two
stages of the dynamical memory of cellular metabolism (metabolic memory).
Epigenetic processes appear to be the structural manifestation of this cellular
metabolic memory. Here, a new framework for molecular information storage in
the cell is presented, which is characterized by two functionally and
molecularly interrelated systems: a dynamic, flexible and adaptive system
(metabolic memory) and an essentially conservative system (genetic memory). The
molecular information of both systems seems to coordinate the physiological
development of the whole cell.


  Inspired by biological evolution, we consider the following so-called
accessibility percolation problem: The vertices of the unoriented
$n$-dimensional binary hypercube are assigned independent $U(0, 1)$ weights,
referred to as fitnesses. A path is considered accessible if fitnesses are
strictly increasing along it. We prove that the probability that the global
fitness maximum is accessible from the all zeroes vertex converges to
$1-\frac{1}{2}\ln\left(2+\sqrt{5}\right)$ as $n\rightarrow\infty$. Moreover, we
prove that if one conditions on the location of the fitness maximum being $v$,
then provided $v$ is not too close to the all zeroes vertex in Hamming
distance, the probability that $v$ is accessible converges to a function of
this distance divided by $n$ as $n\rightarrow\infty$. This resolves a
conjecture by Berestycki, Brunet and Shi in almost full generality.
  As a second result we show that, for any graph, accessibility percolation can
equivalently be formulated in terms of first-passage site percolation. This
connection is of particular importance for the study of accessibility
percolation on trees.


  In order to conduct a statistical analysis on a given set of phylogenetic
gene trees, we often use a distance measure between two trees. In a statistical
distance-based method to analyze discordance between gene trees, it is a key to
decide "biological meaningful" and "statistically well-distributed" distance
between trees. Thus, in this paper, we study the distributions of the three
tree distance metrics: the edge difference, the path difference, and the
precise $K$ interval cospeciation distance, between two trees: first, we focus
on distributions of the three tree distances between two random unrooted trees
with $n$ leaves ($n \geq 4$); and then we focus on the distributions the three
tree distances between a fixed rooted species tree with $n$ leaves and a random
gene tree with $n$ leaves generated under the coalescent process with given the
species tree. We show some theoretical results as well as simulation study on
these distributions.


  The antifungal effects of orthorhombic (~10 nm; spherical) and monoclinic
(~50 nm; tetrapod) sulfur nanoparticles (SNPs) were studied against the
NADPH-dependent glutathione reductase (GR) mediated xenobiotic detoxification
system (GSH-GSSG) in filamentous fungi (Aspergillus niger as a model organism).
Both the SNPs induced significant reduction in fungal growth and spore
formation, and also introduced marked deformities on the surface of
conidiophores at their sub-inhibitory concentrations. A genome wide
transcriptome profile then revealed manifold reduction in the expression of
GSH-GSSG transcripts among SNPs treated fungal isolates, which is unusual for
the micron sized elemental sulfur but probably effective in terms of antifungal
efficacy.


  Enzyme costs play a major role in the choice of metabolic routes, both in
evolution and bioengineering. Given desired fluxes, necessary enzyme levels can
be estimated based on known rate laws and on a principle of minimal enzyme
cost. With logarithmic metabolite levels as free variables, enzyme cost
functions and constraints in optimality and sampling problems can be handled
easily. The set of feasible metabolite profiles forms a polytope in
log-concentration space, whose points represent all possible steady states of a
kinetic model. We show that enzyme cost is a convex function on this polytope.
This makes enzyme cost minimization - finding optimal enzyme profiles and
corresponding metabolite profiles that realize a desired flux at a minimal cost
- a convex optimization problem.


  Binding of a ligand on a protein changes the flexibility of certain parts of
the protein, which directly affects its function. These changes are not the
same at each point, some parts become more flexible and some others become
stiffer. Here, an equation is derived that gives the stiffness map for
proteins. The model is based on correlations of fluctuations of pairs of points
that need to be evaluated by molecular dynamics simulations. The model is also
cast in terms of the Gaussian Network Model and changes of stiffness upon
dimerization of AKT1 are evaluated as an example.


  Although, several factors have been attributed to thermostability, the
stabilization strategies used by proteins are still enigmatic. Studies on
recombinant xylanase which has the ubiquitous (\b{eta}/{\alpha})8 TIM
(Triosephosphate isomerase) barrel fold showed that, just a single extreme
N-terminus mutation (V1L) markedly enhanced the thermostability by 5 {\deg}C
without loss of catalytic activity whereas another mutation, V1A at the same
position decreased the stability by 2 {\deg}C. Based on computational analysis
of their crystal structures including residue interaction network, we
established a link between N- to C-terminal contacts and protein stability. We
demonstrate that augmenting of N- to C-terminal non-covalent interactions is
associated with the enhancement of protein stability. We propose that the
strategy of mutations at the termini could be exploited with a view to modulate
stability without compromising on enzymatic activity, or in general, protein
function, in diverse folds where N- and C-termini are in close proximity.
Finally, we discuss the implications of our results for the development of
therapeutics involving proteins and for designing effective protein engineering
strategies.


  Long, elastic filaments cross-linked and deformed by active molecular motors
occur in various natural settings. The overall macroscopic mechanical response
of such a composite network depends on the coupling between the active and the
passive properties of the underlying constituents and nonlocal interactions
between different parts of the composite. In a simple one dimensional system,
using a mean field model, it has been shown that the combination of motor
activity and finite filament extensibility yields a persistence length scale
over which strain decays. Here we study a similar system, in the complementary
limit of strong noise and moderate extensibility, using Brownian multi-particle
collision dynamics-based numerical simulations that includes the coupling
between motor kinetics and local filament extensibility. While the numerical
model shows deviations from the mean field predictions due to the presence of
strong active noise caused by the variations in individual motor activity,
several qualitative features are still retained. Specifically, for fixed motor
attachment and detachment rates, the decay is length is set by the ratio of the
passive elasticity to the active shear resistance generated by attached motors.
Our study generalizes the notion of persistence in passive thermal systems to
actively driven systems with testable predictions.


  Crosslinked semi-flexible and flexible filaments that are actively deformed
by molecular motors occur in various natural settings, such as the ordered
eukaryotic flagellum, and the disordered cytoskeleton. The deformation of these
composite systems is driven by active motor forces and resisted by passive
filament elasticity, and structural constraints due to permanent cross-links.
Using a mean field theory for a one-dimensional ordered system, we show that
the combination of motor activity and finite filament extensibility yields a
characteristic persistence length scale over which active strain decays. This
decay length is set by the ability of motors to respond to combination of the
weak extensional elasticity, passive shear resistance and the viscoelastic
properties of the motor assembly, and generalizes the notion of persistence in
purely thermal filaments to active systems.


  Two sequentially Markov coalescent models (SMC and SMC') are available as
tractable approximations to the ancestral recombination graph (ARG). We present
a Markov process describing coalescence at two fixed points along a pair of
sequences evolving under the SMC'. Using our Markov process, we derive a number
of new quantities related to the pairwise SMC', thereby analytically
quantifying for the first time the similarity between the SMC' and ARG. We use
our process to show that the joint distribution of pairwise coalescence times
at recombination sites under the SMC' is the same as it is marginally under the
ARG, which demonstrates that the SMC' is, in a particular well-defined,
intuitive sense, the most appropriate first-order sequentially Markov
approximation to the ARG. Finally, we use these results to show that population
size estimates under the pairwise SMC are asymptotically biased, while under
the pairwise SMC' they are approximately asymptotically unbiased.


  The bacterial flagellar motor (BFM) is responsible for driving bacterial
locomotion and chemotaxis, fundamental processes in pathogenesis and biofilm
formation. In the BFM, torque is generated at the interface between
transmembrane proteins (stators) and a rotor. It is well-established that the
passage of ions down a transmembrane gradient through the stator complex
provides the energy needed for torque generation. However, the physics involved
in this energy conversion remain poorly understood. Here we propose a
mechanically specific model for torque generation in the BFM. In particular, we
identify two fundamental forces involved in torque generation: electrostatic
and steric. We propose that electrostatic forces serve to position the stator,
while steric forces comprise the actual 'power stroke'. Specifically, we
predict that ion-induced conformational changes about a proline 'hinge' residue
in an $\alpha$-helix of the stator are directly responsible for generating the
power stroke. Our model predictions fit well with recent experiments on a
single-stator motor. Furthermore, we propose several experiments to elucidate
the torque-speed relationship in motors where the number of stators may not be
constant. The proposed model provides a mechanical explanation for several
fundamental features of the flagellar motor, including: torque-speed and
speed-ion motive force relationships, backstepping, variation in step sizes,
and the puzzle of swarming experiments.


  In biological cells and novel diagnostic devices biochemical receptors need
to be sensitive to extremely small concentration changes of signaling
molecules. The accuracy of such molecular signaling is ultimately limited by
the counting noise imposed by the thermal diffusion of molecules. Many
macromolecules and organelles transiently bind to molecular motors and are then
actively transported. We here show that a random albeit directed delivery of
signaling molecules to within a typical diffusion distance to the receptor
reduces the correlation time of the counting noise, effecting an improved
sensing precision. The conditions for this active focusing are indeed
compatible with observations in living cells. Our results are relevant for a
better understanding of molecular cellular signaling and the design of novel
diagnostic devices.


  Network dynamics are typically presented as a time series of network
properties captured at each period. The current approach examines the dynamical
properties of transmission via novel measures on an integrated, temporally
extended network representation of interaction data across time. Because it
encodes time and interactions as network connections, static network measures
can be applied to this "temporal web" to reveal features of the dynamics
themselves. Here we provide the technical details and apply it to agent-based
implementations of the well-known SEIR and SEIS epidemiological models.


  Fisheries management agencies around the world collect age data for the
purpose of assessing the status of natural resources in their jurisdiction.
Estimates of mortality rates represent a key information to assess the
sustainability of fish stocks exploitation. Contrary to medical research or
manufacturing where survival analysis is routinely applied to estimate failure
rates, survival analysis has seldom been applied in fisheries stock assessment
despite similar purposes between these fields of applied statistics. In this
paper, we developed hazard functions to model the dynamic of an exploited fish
population. These functions were used to estimate all parameters necessary for
stock assessment (including natural and fishing mortality rates as well as gear
selectivity) by maximum likelihood using age data from a sample of catch. This
novel application of survival analysis to fisheries stock assessment was tested
by Monte Carlo simulations to assert that it provided un-biased estimations of
relevant quantities. The method was applied to data from the Queensland
(Australia) sea mullet (Mugil cephalus) commercial fishery collected between
2007 and 2014. It provided, for the first time, an estimate of natural
mortality affecting this stock: 0.22 $\pm$ 0.08 year$^{-1}$.


  Normal anxiety is considered an adaptive response to the possible presence of
danger, but is susceptible to dysregulation. Anxiety disorders are prevalent at
high frequency in contemporary human societies, yet impose substantial
disability upon their sufferers. This raises a puzzle: why has evolution left
us vulnerable to anxiety disorders? We develop a signal detection model in
which individuals must learn how to calibrate their anxiety responses: they
need to learn which cues indicate danger in the environment. We derive the
optimal strategy for doing so, and find that individuals face an inevitable
exploration-exploitation tradeoff between obtaining a better estimate of the
level of risk on one hand, and maximizing current payoffs on the other. Because
of this tradeoff, a subset of the population can become trapped in a state of
self-perpetuating over-sensitivity to threatening stimuli, even when
individuals learn optimally. This phenomenon arises because when individuals
become too cautious, they stop sampling the environment and fail to correct
their misperceptions, whereas when individuals become too careless they
continue to sample the environment and soon discover their mistakes. Thus,
over-sensitivity to threats becomes common whereas under-sensitivity becomes
rare. We suggest that this process may be involved in the development of
excessive anxiety in humans.


  The combinatorics of RNA plays a central role in biology. Mathematical
biologists have several commonly-used models for RNA: words in a fixed alphabet
(representing the primary sequence of nucleotides) and plane trees
(representing the secondary structure, or folding of the RNA sequence). This
paper considers an augmented version of the standard model of plane trees, one
that incorporates some observed constraints on how the folding can occur. In
particular we assume the alphabet consists of complementary pairs, for instance
the Watson-Crick pairs A-U and C-G of RNA.
  Given a word in the alphabet, a valid plane tree is a tree for which, when
the word is folded around the tree, each edge matches two complementary
letters. Consider the graph whose vertices are valid plane trees for a fixed
word and whose edges are given by Condon, Heitsch, and Hoos's local moves. We
prove this graph is connected.
  We give an explicit algorithm to construct a valid plane tree from a primary
sequence, assuming that at least one valid plane tree exists. The tree produced
by our algorithm has other useful characterizations, including a uniqueness
condition defined by local moves. We also study enumerative properties of valid
plane trees, analyzing how the number of valid plane trees depends on the
choice of sequence length and alphabet size. Finally we show that the
proportion of words with at least one valid plane tree goes to zero as the word
size increases. We also give some open questions.


  It is common to use a compartmental, fluid model described by a system of
ordinary differential equations (ODEs) to model disease spread. In addition to
their simplicity, these models are also the mean-field approximations of more
accurate stochastic models of disease spread on contact networks. For the
simplest case of a stochastic susceptible-infected-susceptible (SIS) process
(infection with recovery) on a complete network, it has been shown that the
fraction of infected nodes converges to the mean-field ODE as the number of
nodes increases. However the proofs are not simple, requiring sophisticated
probability, partial differential equations (PDE), or infinite systems of ODEs.
We provide a short proof in this case for convergence in mean-square on finite
time-intervals using a system of two ODEs and a moment inequality and also
obtain the first lower bound on the expected fraction of infected nodes.


  Caspases are a family of cysteine proteases that play essential roles during
apoptosis, and we presume some of them may also protect the cell from oxidative
stress. We found that the absence of yeast caspase 1(Yca1)in Saccharomyces
cerevisiae leads to a more intense burst of mitochondrial reactive oxygen
species (ROS) In addition, compared to wild type yeast cells, the ability of
yca1 mutant cells to maintain mitochondrial activity is significantly reduced
after either oxidative stress treatment or aging. During mitochondrial ROS
burst, deletion of the yca1 gene delayed structural damage of a green
fluorescent protein (GFP) reporter bound in the inner mitochondrial membrane.
This work implies that yeast caspase 1 is closely connected to the oxidative
stress response. We speculate that Yca1 can discriminate proteins damaged by
oxidation and accelerate their hydrolysis to attenuate the ROS burst.


  We propose a population model for TB-HIV/AIDS coinfection transmission
dynamics, which considers antiretroviral therapy for HIV infection and
treatments for latent and active tuberculosis. The HIV-only and TB-only
sub-models are analyzed separately, as well as the TB-HIV/AIDS full model. The
respective basic reproduction numbers are computed, equilibria and stability
are studied. Optimal control theory is applied to the TB-HIV/AIDS model and
optimal treatment strategies for co-infected individuals with HIV and TB are
derived. Numerical simulations to the optimal control problem show that non
intuitive measures can lead to the reduction of the number of individuals with
active TB and AIDS.


  Stochastic kinetic models of genetic expression are able to describe protein
fluctuations. A comparative study of the canonical and a feedback model is
given here by using stochastic simulation methods. The feedback model is
skeleton model implementation of the circular gene hypothesis, which suggests
the interaction between the synthesis and degradation of mRNA. Qualitative and
quantitative changes in the shape and in the numerical characteristics of the
stationary distributions suggest that more combined experimental and
theoretical studies should be done to uncover the details of the kinetic
mechanism of gene expression.


  In recent years single molecule force spectroscopy has opened a new avenue to
provide profiles of the complex energy landscape of biomolecules. In this
field, quantitative analyses of the data employing sound theoretical models,
have played a major role in interpreting data and anticipating outcomes of
experiments. Here, we explain how by using temperature as a variable in
mechanical unfolding of biomolecules in force spectroscopy, the roughness of
the energy landscape can be measured without making any assumptions about the
underlying reaction coordinate. Estimates of other aspects of the energy
landscape such as free energy barriers or the transition state (TS) locations
could depend on the precise model used to analyze the experimental data. We
illustrate the inherent difficulties in obtaining the transition state location
from loading rate or force-dependent unfolding rates. Because the transition
state moves as the force or the loading rate is varied, which is reminiscent of
the Hammond effect, it is in general difficult to invert the experimental data.
The independence of the TS location on force holds good only for brittle or
hard biomolecules whereas the TS location changes considerably if the molecule
is soft or plastic. Finally, we discuss the goodness of the end-to-end distance
(or pulling) coordinate of the molecule as a surrogate reaction coordinate in a
situation like force-induced ligand-unbinding from a gated molecular receptor
as well as force-quench refolding of an RNA hairpin.


  When the body gets infected by a pathogen or receives a vaccine dose, the
immune system develops pathogen-specific immunity. Induced immunity decays in
time and years after recovery/vaccination the host might become susceptible
again. Exposure to the pathogen in the environment boosts the immune system
thus prolonging the duration of the protection. Such an interplay of within
host and population level dynamics poses significant challenges in rigorous
mathematical modeling of immuno-epidemiology. The aim of this paper is twofold.
First, we provide an overview of existing models for waning of
disease/vaccine-induced immunity and immune system boosting. Then a new
modeling approach is proposed for SIRVS dynamics, monitoring the immune status
of individuals and including both waning immunity and immune system boosting.
We show that some previous models can be considered as special cases or
approximations of our framework.


  In this paper, decision theory was used to derive Bayes and minimax decision
rules to estimate allelic frequencies and to explore their admissibility.
Decision rules with uniformly smallest risk usually do not exist and one
approach to solve this problem is to use the Bayes principle and the minimax
principle to find decision rules satisfying some general optimality criterion
based on their risk functions. Two cases were considered, the simpler case of
biallelic loci and the more complex case of multiallelic loci. For each locus,
the sampling model was a multinomial distribution and the prior was a Beta
(biallelic case) or a Dirichlet (multiallelic case) distribution. Three loss
functions were considered: squared error loss (SEL), Kulback-Leibler loss (KLL)
and quadratic error loss (QEL). Bayes estimators were derived under these three
loss functions and were subsequently used to find minimax estimators using
results from decision theory. The Bayes estimators obtained from SEL and KLL
turned out to be the same. Under certain conditions, the Bayes estimator
derived from QEL led to an admissible minimax estimator (which was also equal
to the maximum likelihood estimator). The SEL also allowed finding admissible
minimax estimators. Some estimators had uniformly smaller variance than the MLE
and under suitable conditions the remaining estimators also satisfied this
property. In addition to their statistical properties, the estimators derived
here allow variation in allelic frequencies, which is closer to the reality of
finite populations exposed to evolutionary forces.


  As researchers collect increasingly large molecular data sets to reconstruct
the Tree of Life, the heterogeneity of signals in the genomes of diverse
organisms poses challenges for traditional phylogenetic analysis. A class of
phylogenetic methods known as "species tree methods" have been proposed to
directly address one important source of gene tree heterogeneity, namely the
incomplete lineage sorting or deep coalescence that occurs when evolving
lineages radiate rapidly, resulting in a diversity of gene trees from a single
underlying species tree. Although such methods are gaining in popularity, they
are being adopted with caution in some quarters, in part because of an
increasing number of examples of strong phylogenetic conflict between
concatenation or supermatrix methods and species tree methods. Here we review
theory and empirical examples that help clarify these conflicts. Thinking of
concatenation as a special case of the more general model provided by the
multispecies coalescent can help explain a number of differences in the
behavior of the two methods on phylogenomic data sets. Recent work suggests
that species tree methods are more robust than concatenation approaches to some
of the classic challenges of phylogenetic analysis, including rapidly evolving
sites in DNA sequences, base compositional heterogeneity and long branch
attraction. We show that approaches such as binning, designed to augment the
signal in species tree analyses, can distort the distribution of gene trees and
are inconsistent. Computationally efficient species tree methods that
incorporate biological realism are a key to phylogenetic analysis of whole
genome data.


  Although mutations drive the evolutionary process, the rates at which the
mutations occur are themselves subject to evolutionary forces. Our purpose here
is to understand the role of selection and random genetic drift in the
evolution of mutation rates, and we address this question in asexual
populations at mutation-selection equilibrium neglecting selective sweeps.
Using a multitype branching process, we calculate the fixation probability of a
rare nonmutator in a large asexual population of mutators, and find that a
nonmutator is more likely to fix when the deleterious mutation rate of the
mutator population is high. Compensatory mutations in the mutator population
are found to decrease the fixation probability of a nonmutator when the
selection coefficient is large. But, surprisingly, the fixation probability
changes nonmonotonically with increasing compensatory mutation rate when the
selection is mild. Using these results for the fixation probability and a
drift-barrier argument, we find a novel relationship between the mutation rates
and the population size. We also discuss the time to fix the nonmutator in an
adapted population of asexual mutators, and compare our results with
experiments.


  Salamanders (urodela) have among the largest vertebrate genomes, ranging in
size from 10 to over 80 pg. The urodela are divided into ten extant families
each with a characteristic range in genome size. Although changes in genome
size often occur randomly and in the absence of selection pressure, non-random
patterns of genome size variation are evident among specific vertebrate
lineages. Here we report that genome size in salamander families varies
inversely with species richness and other ecological factors: clades that began
radiating earlier (older crown age) tend to have smaller genomes, higher levels
of diversity and larger geographical ranges. These observations support the
hypothesis that urodel families with larger genomes either have a lower
propensity to diversify or are more vulnerable to extinction than families with
smaller genomes.


  The analysis of the three-dimensional structure of proteins is an important
topic in molecular biochemistry. Structure plays a critical role in defining
the function of proteins and is more strongly conserved than amino acid
sequence over evolutionary timescales. A key challenge is the identification
and evaluation of structural similarity between proteins; such analysis can aid
in understanding the role of newly discovered proteins and help elucidate
evolutionary relationships between organisms. Computational biologists have
developed many clever algorithmic techniques for comparing protein structures,
however, all are based on heuristic optimization criteria, making statistical
interpretation somewhat difficult. Here we present a fully probabilistic
framework for pairwise structural alignment of proteins. Our approach has
several advantages, including the ability to capture alignment uncertainty and
to estimate key "gap" parameters which critically affect the quality of the
alignment. We show that several existing alignment methods arise as maximum a
posteriori estimates under specific choices of prior distributions and error
models. Our probabilistic framework is also easily extended to incorporate
additional information, which we demonstrate by including primary sequence
information to generate simultaneous sequence-structure alignments that can
resolve ambiguities obtained using structure alone. This combined model also
provides a natural approach for the difficult task of estimating evolutionary
distance based on structural alignments. The model is illustrated by comparison
with well-established methods on several challenging protein alignment
examples.


  Intransitivity is supposed to be a main reason for deficits in coevolutionary
progress and inheritable superiority. Besides, coevolutionary dynamics is
characterized by interactions yielding subjective fitness, but aiming at
solutions that are superior with respect to an objective measurement. Such an
approximation of objective fitness may be, for instance, generalization
performance. In the paper a link between rating-- and ranking--based measures
of intransitivity and fitness landscapes that can address the dichotomy between
subjective and objective fitness is explored. The approach is illustrated by
numerical experiments involving a simple random game with continuously tunable
degree of randomness.


  We typically think of cells as responding to external signals independently
by regulating their gene expression levels, yet they often locally exchange
information and coordinate. Can such spatial coupling be of benefit for
conveying signals subject to gene regulatory noise? Here we extend our
information-theoretic framework for gene regulation to spatially extended
systems. As an example, we consider a lattice of nuclei responding to a
concentration field of a transcriptional regulator (the "input") by expressing
a single diffusible target gene. When input concentrations are low, diffusive
coupling markedly improves information transmission; optimal gene activation
functions also systematically change. A qualitatively new regulatory strategy
emerges where individual cells respond to the input in a nearly step-like
fashion that is subsequently averaged out by strong diffusion. While motivated
by early patterning events in the Drosophila embryo, our framework is
generically applicable to spatially coupled stochastic gene expression models.


  We study the effects of multiple binding sites in the promoter of a genetic
oscillator. We evaluate the regulatory function of a promoter with multiple
binding sites in the absence of cooperative binding, and consider different
hypotheses for how the number of bound repressors affects transcription rate.
Effective Hill exponents of the resulting regulatory functions reveal an
increase in the nonlinearity of the feedback with the number of binding sites.
We identify optimal configurations that maximize the nonlinearity of the
feedback. We use a generic model of a biochemical oscillator to show that this
increased nonlinearity is reflected in enhanced oscillations, with larger
amplitudes over wider oscillatory ranges. Although the study is motivated by
genetic oscillations in the zebrafish segmentation clock, our findings may
reveal a general principle for gene regulation.


  Spherical viral shells with icosahedral symmetry have been considered as
quasicrystalline tilings. Similarly to known Caspar-Klug quasi-equivalence
theory, the presented approach also minimizes the number of conformations
necessary for the protein molecule bonding with its neighbors in the shell, but
is based on different geometrical principles. It is assumed that protein
molecule centers are located at vertices of tiles with identical edges, and the
number of different tile types is minimal. Idealized coordinates of
nonequivalent by symmetry protein positions in six various capsid types are
obtained. The approach describes in a uniform way both the structures
satisfying the well-known Caspar-Klug geometrical model and the structures
contradicting this model.


  Motivation: Sequence mapping is the cornerstone of modern genomics. However,
most existing sequence mapping algorithms are insufficiently general.
  Results: We introduce context schemes: a method that allows the unambiguous
recognition of a reference base in a query sequence by testing the query for
substrings from an algorithmically defined set. Context schemes only map when
there is a unique best mapping, and define this criterion uniformly for all
reference bases. Mappings under context schemes can also be made stable, so
that extension of the query string (e.g. by increasing read length) will not
alter the mapping of previously mapped positions. Context schemes are general
in several senses. They natively support the detection of arbitrary complex,
novel rearrangements relative to the reference. They can scale over orders of
magnitude in query sequence length. Finally, they are trivially extensible to
more complex reference structures, such as graphs, that incorporate additional
variation. We demonstrate empirically the existence of high performance context
schemes, and present efficient context scheme mapping algorithms.
  Availability and Implementation: The software test framework created for this
work is available from
https://registry.hub.docker.com/u/adamnovak/sequence-graphs/.
  Contact: benedict@soe.ucsc.edu
  Supplementary Information: Six supplementary figures and one supplementary
section are available with the online version of this article.


  Background. The calculation of diffusion-controlled ligand binding rates is
important for understanding enzyme mechanisms as well as designing enzyme
inhibitors. We demonstrate the accuracy and effectiveness of a Lagrangian
particle-based method, smoothed particle hydrodynamics (SPH), to study
diffusion in biomolecular systems by numerically solving the time-dependent
Smoluchowski equation for continuum diffusion.
  Results. The numerical method is first verified in simple systems and then
applied to the calculation of ligand binding to an acetylcholinesterase
monomer. Unlike previous studies, a reactive Robin boundary condition (BC),
rather than the absolute absorbing (Dirichlet) boundary condition, is
considered on the reactive boundaries. This new boundary condition treatment
allows for the analysis of enzymes with "imperfect" reaction rates. Rates for
inhibitor binding to mAChE are calculated at various ionic strengths and
compared with experiment and other numerical methods. We find that imposition
of the Robin BC improves agreement between calculated and experimental reaction
rates.
  Conclusions. Although this initial application focuses on a single monomer
system, our new method provides a framework to explore broader applications of
SPH in larger-scale biomolecular complexes by taking advantage of its
Lagrangian particle-based nature.


  In this paper we study disease spread over a randomly switched network, which
is modeled by a stochastic switched differential equation based on the so
called $N$-intertwined model for disease spread over static networks. Assuming
that all the edges of the network are independently switched, we present
sufficient conditions for the convergence of infection probability to zero.
Though the stability theory for switched linear systems can naively derive a
necessary and sufficient condition for the convergence, the condition cannot be
used for large-scale networks because, for a network with $n$ agents, it
requires computing the maximum real eigenvalue of a matrix of size exponential
in $n$. On the other hand, our conditions that are based also on the spectral
theory of random matrices can be checked by computing the maximum real
eigenvalue of a matrix of size exactly $n$.


  The Cancer Genome Atlas project was initiated by the National Cancer
Institute in order to characterize the genomes of hundreds of tumors of various
cancer types. While much effort has been put into detecting somatic genomic
variation in these data, somatic structural variation induced by the activity
of transposable element insertions has not been reported. Transposable elements
(TEs) are particularly relevant in cancer in part because of several known
cases in which a TE insertion is directly linked to cancer formation and
studies linking the epigenetic status of retrotransposons to carcinogenesis and
patient outcome. Additionally, evidence for somatic retrotransposition in
eukaryotic genomes suggests that some tissues and therefore some cancer types
may be disposed to increased retrotransposition. We built upon previous work to
develop a highly efficient computational pipeline for the detection of
non-reference mobile ele- ment insertions from high-throughput paired-end whole
genome sequencing data that is capable of detecting breakpoints through a local
assembly strategy. Using this, we analyzed 33 whole genome tumor datasets with
paired normal samples from TCGA across 3 different cancer types: glioblastoma
multiforme (GBM), ovarian serous cystoadenocarcinoma (OV) and colorectal ade-
nocarcinoma (COAD). We detected 72 insertions in colon samples, almost all of
them LINE-1 elements, and none in GBM or OV. The amount of somatic
retrotransposition varies widely between samples with 61 insertions present in
one case. The lack of somatic retrotransposon insertions in GBM and OV samples
suggests that TE activity in cancer is restricted to certain cancer types.


  The aim of this paper is to give an introductory review of the cryoelectron
microscopy as a complex data source for the most of the system biology
branches, including the most perspective non-local approaches known as
"localomics" and "dynamomics". A brief summary of various cryoelectron
mi-croscopy methods and corresponding system biological ap-proaches is given in
the text. The above classification can be considered as a useful framework for
the primary comprehen-sions about cryoelectron microscopy aims and instrumental
tools. We do not discuss any of these concepts in details, but merely point out
that their methodological complexity follows only from the structure-functional
complexity of biological systems which are investigated in this manner. We also
postu-late that one can employ some of the cryoelectron microscopic techniques
not only for observation, but also for modification and structural
refunctionalization of some biological and similar soft matter objects and
microscopic samples. In other worlds, we start with the cryoelectron microscopy
as a tool for the sys-tem biology and progress to its applying as an instrument
for system biology and functional biomimetics; i.e. "system cryobi-ology" goes
over into "synthetic cryobiology" or "cryogenic biomimetics". All these
conclusions can be deduced from the most recent works of the latest years,
including just submitted foreign papers. This article provides an up-to-date
description of the conceptual basis for the novel view on the computational
cryoelectron microscopy (in silico) approaches and the data mining principles
which lie at the very foundation of modern structural analysis and
reconstruction.


  This article explains how natural selection works and how it has been
inappropriately applied to the description of cultural change. It proposes an
alternative evolutionary explanation for cultural evolution that describes it
in terms of communal exchange.


  In food webs, many interacting species coexist despite the restrictions
imposed by the competitive exclusion principle and apparent competition. For
the generalized Lotka-Volterra equations, sustainable coexistence necessitates
nonzero determinant of the interaction matrix. Here we show that this
requirement is equivalent to demanding that each species be part of a
non-overlapping pairing, which substantially constrains the food web structure.
We demonstrate that a stable food web can always be obtained if a
non-overlapping pairing exists. If it does not, the matrix rank can be used to
quantify the lack of niches, corresponding to unpaired species. For the species
richness at each trophic level, we derive the food web assembly rules, which
specify sustainable combinations. In neighboring levels, these rules allow the
higher level to avert competitive exclusion at the lower, thereby incorporating
apparent competition. In agreement with data, the assembly rules predict high
species numbers at intermediate levels and thinning at the top and bottom.
Using comprehensive food web data, we demonstrate how omnivores or parasites
with hosts at multiple trophic levels can loosen the constraints and help
obtain coexistence in food webs. Hence, omnivory may be the glue that keeps
communities intact even under extinction or ecological release of species.


  In this paper, we introduce the software suite, Hermes, which provides fast,
novel algorithms for RNA secondary structure kinetics. Using the fast Fourier
transform to efficiently compute the Boltzmann probability that a secondary
structure S of a given RNA sequence has base pair distance x [resp. y] from
reference structure A [resp. B], Hermes computes the exact kinetics of folding
from A to B in this coarse-grained model. In particular, Hermes computes the
mean first passage time from the transition probability matrix by using matrix
inversion, and also computes the equilibrium time from the rate matrix by using
spectral decomposition. Due to the model granularity and the speed of Hermes,
it is capable of determining secondary structure refolding kinetics for large
RNA sequences, beyond the range of other methods. Comparative benchmarking of
Hermes with other methods indicates that Hermes provides refolding kinetics of
accuracy suitable for use in computational design of RNA, an important area of
synthetic biology. Source code and documentation for Hermes are available at
http://bioinformatics. bc.edu/clotelab/Hermes/.


  Darwin claims in the {\em Origin} that similarity is evidence for common
ancestry, but that adaptive similarities are "almost valueless" as evidence.
This claim seems reasonable for some adaptive similarities but not for others.
Here we clarify and evaluate these and related matters by using the law of
likelihood as an analytic tool and by considering mathematical models of three
evolutionary processes -- directional selection, stabilizing selection, and
drift. Our results apply both to Darwin's theory of evolution and to modern
evolutionary biology.


  Biological systems are modular, and this modularity affects the evolution of
biological systems over time and in different environments. We here develop a
theory for the dynamics of evolution in a rugged, modular fitness landscape. We
show analytically how horizontal gene transfer couples to the modularity in the
system and leads to more rapid rates of evolution at short times. The model, in
general, analytically demonstrates a selective pressure for the prevalence of
modularity in biology. We use this model to show how the evolution of the
influenza virus is affected by the modularity of the proteins that are
recognized by the human immune system. Approximately 25\% of the observed rate
of fitness increase of the virus could be ascribed to a modular viral
landscape.


  Biological functions are carried out by groups of interacting molecules,
cells or tissues, known as communities. Membership in these communities may
overlap when biological components are involved in multiple functions. However,
traditional clustering methods detect non-overlapping communities. These
detected communities may also be unstable and difficult to replicate, because
traditional methods are sensitive to noise and parameter settings. These
aspects of traditional clustering methods limit our ability to detect
biological communities, and therefore our ability to understand biological
functions.
  To address these limitations and detect robust overlapping biological
communities, we propose an unorthodox clustering method called SpeakEasy which
identifies communities using top-down and bottom-up approaches simultaneously.
Specifically, nodes join communities based on their local connections, as well
as global information about the network structure. This method can quantify the
stability of each community, automatically identify the number of communities,
and quickly cluster networks with hundreds of thousands of nodes.
  SpeakEasy shows top performance on synthetic clustering benchmarks and
accurately identifies meaningful biological communities in a range of datasets,
including: gene microarrays, protein interactions, sorted cell populations,
electrophysiology and fMRI brain imaging.


  Networks have become an important tool for infectious disease epidemiology.
Most previous theoretical studies of transmission network models have either
considered simple Markovian dynamics at the individual level, or have focused
on the invasion threshold and final outcome of the epidemic. Here, we provide a
general theory for early real-time behaviour of epidemics on large
configuration model networks (i.e. static and locally unclustered), in
particular focusing on the computation of the Malthusian parameter that
describes the early exponential epidemic growth. Analytical, numerical and
Monte-Carlo methods under a wide variety of Markovian and non-Markovian
assumptions about the infectivity profile are presented. Numerous examples
provide explicit quantification of the impact of the network structure on the
temporal dynamics of the spread of infection and provide a benchmark for
validating results of large scale simulations.


  This manuscript has been written to address questions related to our recent
publication (Science 347:78-81, 2015). We appreciate the many reactions to this
paper that have been communicated to us, either privately or publicly. The
following addresses several of the most important statistical and technical
issues related to our analysis and conclusions. Our responses to non-technical
questions are available at
http://www.hopkinsmedicine.org/news/media/releases/bad_luck_of_random_mutations_plays_predominant_role_in_cancer_study_shows


  It has been hypothesized that in the era just before the last universal
common ancestor emerged, life on earth was fundamentally collective. Ancient
life forms shared their genetic material freely through massive horizontal gene
transfer (HGT). At a certain point, however, life made a transition to the
modern era of individuality and vertical descent. Here we present a minimal
model for this hypothesized "Darwinian transition." The model suggests that
HGT-dominated dynamics may have been intermittently interrupted by
selection-driven processes during which genotypes became fitter and decreased
their inclination toward HGT. Stochastic switching in the population dynamics
with three-point (hypernetwork) interactions may have destabilized the
HGT-dominated collective state and led to the emergence of vertical descent and
the first well-defined species in early evolution. A nonlinear analysis of a
stochastic model dynamics covering key features of evolutionary processes (such
as selection, mutation, drift and HGT) supports this view. Our findings thus
suggest a viable route from early collective evolution to the start of
individuality and vertical Darwinian evolution, enabling the emergence of the
first species.


  The phase structure of a homopolymer chain is investigated in terms of a
universal theoretical model, designed to describe the infrared limit of slow
spatial variations. The effects of chirality are studied and compared with the
influence of a short-range attractive interaction between monomers, at various
ambient temperature values. In the high temperature limit the homopolymer chain
is in the self-avoiding random walk phase. But at low temperatures, two
different phases are possible: When short-range attractive interactions
dominate over chirality, the chain collapses into a space- filling
conformation. But when the attractive interactions become weaker, there is a
low temperature unfolding transition and the chain becomes like a straight rod.
Between the high temperature and low temperature limits, several intermediate
states are observed. For sufficiently high values of short-range attraction,
the conventional {\theta}-regime is observed between the self-avoiding random
walk phase and the space filling collapsed phase. But when chirality increases,
there is a trasition from the {\theta}-regime to a pseudogap state. Moreover, a
regime akin the {\theta}-regime is identified between the pseudogap state, and
the low temperature phase where the chain is like a straight rod. Applications
to polymers and proteins, in particular collagen, are suggested.


  This is a review of William Feller's important contributions to mathematical
biology. The seminal paper [Feller1951] "Diffusion processes in genetics" was
particularly influential on the development of stochastic processes at the
interface to evolutionary biology, and interesting ideas in this direction
(including a first characterization of what is nowadays known as "Feller's
branching diffusion") already shaped up in the paper [Feller 1939] (written in
German) "The foundations of a probabistic treatment of Volterra's theory of the
struggle for life". Feller's article "On fitness and the cost of natural
selection" [Feller 1967] contains a critical analysis of the concept of
"genetic load".


  The decreasing costs and increasing speed and accuracy of DNA sample
collection, preparation, and sequencing has rapidly produced an enormous volume
of genetic data. However, fast and accurate analysis of the samples remains a
bottleneck. Here we present D$^{4}$RAGenS, a genetic sequence identification
algorithm that exhibits the Big Data handling and computational power of the
Dynamic Distributed Dimensional Data Model (D4M). The method leverages linear
algebra and statistical properties to increase computational performance while
retaining accuracy by subsampling the data. Two run modes, Fast and Wise, yield
speed and precision tradeoffs, with applications in biodefense and medical
diagnostics. The D$^{4}$RAGenS analysis algorithm is tested over several
datasets, including three utilized for the Defense Threat Reduction Agency
(DTRA) metagenomic algorithm contest.


  The study of games and their equilibria is central to developing insights for
understanding many socio-economic phenomena. Here we present a dynamical
systems view of the equilibria of two-person, payoff-symmetric games. In
particular, using this perspective, we discuss the differences between two
solution concepts for such games - namely, those of Nash equilibrium and
co-action equilibrium. For the Nash equilibrium, we show that the dynamical
view can provide an equilibrium refinement, selecting one equilibrium among
several possibilities, thereby solving the issue of multiple equilibria that
appear in some games. We illustrate in detail this dynamical perspective by
considering three well known 2-person games namely the Prisoner's Dilemma, game
of Chicken and the Stag-Hunt. We find that in all of these cases, co-action
equilibria tends to correspond to `nicer' strategies than those corresponding
to Nash equilibria.


  1. Predicting space use patterns of animals from their interactions with the
environment is fundamental for understanding the effect of habitat changes on
ecosystem functioning. Recent attempts to address this problem have sought to
unify resource selection analysis, where animal space use is derived from
available habitat quality, and mechanistic movement models, where detailed
movement processes of an animal are used to predict its emergent utilization
distribution. Such models bias the animal's movement towards patches that are
easily available and resource-rich, and the result is a predicted probability
density at a given position being a function of the habitat quality at that
position. However, in reality, the probability that an animal will use a patch
of the terrain tends to be a function of the resource quality in both that
patch and the surrounding habitat. 2. We propose a mechanistic model where this
non-local effect of resources naturally emerges from the local movement
processes, by taking into account the relative utility of both the habitat
where the animal currently resides and that of where it is moving. We give
statistical techniques to parametrize the model from location data and
demonstrate application of these techniques to GPS data of caribou in
Newfoundland. 3. Steady-state animal probability distributions arising from the
model have complex patterns that cannot be expressed simply as a function of
the local quality of the habitat. In particular, large areas of good habitat
are used more intensively than smaller patches of equal quality habitat,
whereas isolated patches are used less frequently. 4. Whilst we focus on
habitats in this study, our modelling framework can be readily used with any
environmental covariates and therefore represents a unification of mechanistic
modelling and step selection approaches to understanding animal space use.


  The Ebola epidemic in West Africa is the largest ever recorded, with over
27,000 cases and 11,000 deaths as of June 2015. The public health response was
challenged by difficulties with disease surveillance, which impacted subsequent
analysis and decision-making regarding optimal interventions. We developed a
stage-structured model of Ebola virus disease (EVD). A key feature of the model
is that it includes a generalized correction term accounting for factors such
as the fraction of cases reported and fraction of the population at risk (e.g.
due to contact patterns, interventions, spatiotemporal spread, pre-existing
immunity, asymptomatic cases, etc.). We generated a range of short-term
forecasts for Guinea, Liberia, and Sierra Leone, which we then validated using
subsequent data. We also used the model to examine the uncertainty in the
relative contributions to transmission by the different stages of infection
(early, late, and funeral). We found that a wide range of forecasted
trajectories fit approximately equally well to the early data. However, by
including the correction factor term the best-fit models correctly forecasted
EVD activity for all three countries, both individually and for all countries
combined. In particular, the model correctly forecasted the slow-down in
Liberia, as well as the continued exponential growth in Sierra Leone through
November 2014. Parameter unidentifiability issues hindered estimation of the
relative contributions of each stage of transmission from incidence and deaths
data alone, which poses a challenge in determining optimal intervention
strategies, and underscores the need for additional data collection. Even with
these limited data, however, it is still possible to accurately capture and
predict the epidemic dynamics by using a simplified correction term that
approximately accounts for the complex underlying factors driving disease
spread.


  It has been suggested that quantum coherence in the selectivity filter of ion
channel may play a key role in fast conduction and selectivity of ions.
However, it has not been clearly elucidated yet why classical coherence is not
sufficient for this purpose. In this paper, we investigate the classical
vibrational coherence between carbonyl groups oscillations in the selectivity
filter of KcsA ion channels based on the data obtained from molecular dynamics
simulations. Our results show that classical coherence plays no effective role
in fast ionic conduction.


  We present a Bayesian method for characterizing the mating system of
populations reproducing through a mixture of self-fertilization and random
outcrossing. Our method uses patterns of genetic variation across the genome as
a basis for inference about pure hermaphroditism, androdioecy, and gynodioecy.
We extend the standard coalescence model to accommodate these mating systems,
accounting explicitly for multilocus identity disequilibrium, inbreeding
depression, and variation in fertility among mating types. We incorporate the
Ewens Sampling Formula (ESF) under the infinite-alleles model of mutation to
obtain a novel expression for the likelihood of mating system parameters. Our
Markov chain Monte Carlo (MCMC) algorithm assigns locus-specific mutation
rates, drawn from a common mutation rate distribution that is itself estimated
from the data using a Dirichlet Process Prior (DPP) model. Among the parameters
jointly inferred are the population-wide rate of self-fertilization,
locus-specific mutation rates, and the number of generations since the most
recent outcrossing event for each sampled individual.


  Given any finite and closed chemical reaction system, it is possible to
efficiently determine whether or not it contains a `self-sustaining and
collectively autocatalytic' subset of reactions, and to find such subsets when
they exist. However, for systems that are potentially open-ended (for example,
when no prescribed upper bound is placed on the complexity or size/length of
molecules types), the theory developed for the finite case breaks down. We
investigate a number of subtleties that arise in such systems that are absent
in the finite setting, and present several new results.


  In this article we formulate new models for coupled systems of bulk-surface
reaction-diffusion equations on stationary volumes. The bulk reaction-diffusion
equations are coupled to the surface reaction-diffusion equations through
linear Robin-type boundary conditions. We then state and prove the necessary
conditions for diffusion-driven instability for the coupled system. Due to the
nature of the coupling between bulk and surface dynamics, we are able to
decouple the stability analysis of the bulk and surface dynamics. Under a
suitable choice of model parameter values, the bulk reaction-diffusion system
can induce patterning on the surface independent of whether the surface
reaction-diffusion system produces or not, patterning. On the other hand, the
surface reaction-diffusion system can not generate patterns everywhere in the
bulk in the absence of patterning from the bulk reaction-diffusion system. For
this case, patterns can only be induced in regions close to the surface
membrane. Various numerical experiments are presented to support our
theoretical findings. Our most revealing numerical result is that, Robin-type
boundary conditions seem to introduce a boundary layer coupling the bulk and
surface dynamics.


  Astrophysical ionizing radiation events have been recognized as a potential
threat to life on Earth, primarily through depletion of stratospheric ozone and
subsequent increase in surface-level solar ultraviolet radiation. Simulations
of the atmospheric effects of a variety of events (such as supernovae,
gamma-ray bursts, and solar proton events) have been previously published,
along with estimates of biological damage at Earth's surface. In this work, we
employed the TUV radiative transfer model to expand and improve calculations of
surface-level irradiance and biological impacts following an ionizing radiation
event. We considered changes in surface-level UVB, UVA, and photosynthetically
active radiation (visible light) for clear-sky conditions and fixed aerosol
parameter values. We also considered a wide range of biological effects on
organisms ranging from humans to phytoplankton. We found that past work
overestimated UVB irradiance, but that relative estimates for increase in
exposure to DNA damaging radiation are still similar to our improved
calculations. We also found that the intensity of biologically damaging
radiation varies widely with organism and specific impact considered; these
results have implications for biosphere-level damage following astrophysical
ionizing radiation events. When considering changes in surface-level visible
light irradiance, we found that, contrary to previous assumptions, a decrease
in irradiance is only present for a short time in very limited geographical
areas; instead we found a net increase for most of the modeled time-space
region. This result has implications for proposed climate changes associated
with ionizing radiation events.


  Single particle electron cryomicroscopy (cryo-EM) allows for structures of
proteins and protein complexes to be determined from images of non-crystalline
specimens. Cryo-EM data analysis requires electron microscope images of
randomly oriented ice-embedded protein particles to be rotated and translated
to allow for coherent averaging when calculating three-dimensional (3D)
structures. Rotation of 2D images is usually done with the assumption that the
magnification of the electron microscope is the same in all directions.
However, due to electron optical aberrations, this condition is not met with
some electron microscopes when used with the settings necessary for cryo-EM
with a direct detector device (DDD) camera. Correction of images by linear
interpolation in real space has allowed high-resolution structures to be
calculated from cryo-EM images for symmetric particles. Here we describe and
compare a simple real space method, a simple Fourier space method, and a
somewhat more sophisticated Fourier space method to correct images for a
measured anisotropy in magnification. Further, anisotropic magnification causes
contrast transfer function (CTF) parameters estimated from image power spectra
to have an apparent systematic astigmatism. To address this problem we develop
an approach to adjust CTF parameters measured from distorted images so that
they can be used with corrected images. The effect of anisotropic magnification
on CTF parameters provides a simple way of detecting magnification anisotropy
in cryo-EM datasets.


  Visualization and analysis of molecular profiling data together with
biological networks are able to provide new mechanistical insights into
biological functions. Currently, high-throughput data are usually visualized on
top of predefined network layouts which are not always adapted to a given data
analysis task. We developed a Cytoscape app which allows to construct
biological network layouts based on the data from molecular profiles imported
as values of nodes attributes. DeDaL is a Cytoscape 3.0 app which uses linear
and non-linear algorithms of dimension reduction to produce data-driven network
layouts based on multidimensional data (typically gene expression). DeDaL
implements several data pre-processing and layout post-processing steps such as
continuous morphing between two arbitrary network layouts and aligning one
network layout with respect to another one by rotating and mirroring. Combining
these possibilities facilitates creating insightful network layouts
representing both structural network features and the correlation patterns in
multivariate data. DeDaL is the first method allowing to construct biological
network layouts from high-throughput data. DeDaL is freely available for
downloading together with step-by-step tutorial at
http://bioinfo-out.curie.fr/projects/dedal/.


  Marine protected areas (MPAs) have attracted much attention as a tool for
sustainable fisheries management, restoring depleted fisheries stocks and
maintaining ecosystems. However, even with total exclusion of fishing effort,
depleted stocks sometimes show little or no recovery over a long time period.
Here, using a mathematical model, we show that multiple stable states may hold
the key to understanding the tendency for fisheries stocks to recover because
of MPAs. We find that MPAs can have either a positive effect or almost no
effect on the recovery of depleted fishing stocks, depending on the fish
migration patterns and the fishing policies. MPAs also reinforce ecological
resilience, particularly for migratory species. In contrast to previous
reports, our results show that MPAs have small or sometimes negative effects on
the recovery of sedentary species. Unsuitable MPA planning might result in low
effectiveness or even deterioration of the existing condition.


  While most current high-throughput DNA sequencing technologies generate short
reads with low error rates, emerging sequencing technologies generate long
reads with high error rates. A basic question of interest is the tradeoff
between read length and error rate in terms of the information needed for the
perfect assembly of the genome. Using an adversarial erasure error model, we
make progress on this problem by establishing a critical read length, as a
function of the genome and the error rate, above which perfect assembly is
guaranteed. For several real genomes, including those from the GAGE dataset, we
verify that this critical read length is not significantly greater than the
read length required for perfect assembly from reads without errors.


  Evolutionary dynamics has been classically studied for homogeneous
populations, but now there is a growing interest in the non-homogenous case.
One of the most important models has been proposed by Lieberman, Hauert and
Nowak, adapting to a weighted directed graph the classical process described by
Moran. The Markov chain associated with the graph can be modified by erasing
all non-trivial loops in its state space, obtaining the so-called Embedded
Markov chain (EMC). The fixation probability remains unchanged, but the
expected time to absorption (fixation or extinction) is reduced. In this paper,
we shall use this idea to compute asymptotically the average fixation
probability for complete bipartite graphs. To this end, we firstly review some
recent results on evolutionary dynamics on graphs trying to clarify some
points. We also revisit the 'Star Theorem' proved by Lieberman, Hauert and
Nowak for the star graphs. Theoretically, EMC techniques allow fast computation
of the fixation probability, but in practice this is not always true. Thus, in
the last part of the paper, we compare this algorithm with the standard Monte
Carlo method for some kind of complex networks.


  We first generalise ideas discussed by Kiss et al. (2015) to prove a theorem
for generating exact closures (here expressing joint probabilities in terms of
their constituent marginal probabilities) for susceptible-infectious-removed
(SIR) dynamics on arbitrary graphs (networks). For Poisson transmission and
removal processes, this enables us to obtain a systematic reduction in the
number of differential equations needed for an exact `moment closure'
representation of the underlying stochastic model. We define `transmission
blocks' as a possible extension of the block concept in graph theory and show
that the order at which the exact moment closure representation is curtailed is
the size of the largest transmission block. More generally, approximate
closures of the hierarchy of moment equations for these dynamics are typically
defined for the first and second order yielding mean-field and pairwise models
respectively. It is frequently implied that, in principle, closed models can be
written down at arbitrary order if only we had the time and patience to do
this. However, for epidemic dynamics on networks, these higher-order models
have not been defined explicitly. Here we unambiguously define hierarchies of
approximate closed models that can utilise subsystem states of any order, and
show how well-known models are special cases of these hierarchies.


  We consider an environment where players are involved in a public goods game
and must decide repeatedly whether to make an individual contribution or not.
However, players lack strategically relevant information about the game and
about the other players in the population. The resulting behavior of players is
completely uncoupled from such information, and the individual strategy
adjustment dynamics are driven only by reinforcement feedbacks from each
player's own past. We show that the resulting "directional learning" is
sufficient to explain cooperative deviations away from the Nash equilibrium. We
introduce the concept of k-strong equilibria, which nest both the Nash
equilibrium and the Aumann-strong equilibrium as two special cases, and we show
that, together with the parameters of the learning model, the maximal
k-strength of equilibrium determines the stationary distribution. The
provisioning of public goods can be secured even under adverse conditions, as
long as players are sufficiently responsive to the changes in their own payoffs
and adjust their actions accordingly. Substantial levels of public cooperation
can thus be explained without arguments involving selflessness or social
preferences, solely on the basis of uncoordinated directional (mis)learning.


  This short note provides a simple formal proof of a folklore result in
statistical phylogenetics concerning the convergence of bootstrap support for a
tree and its edges.


  Mathematical modeling is now used commonly in the analysis of signaling
networks. With advances in high resolution microscopy, the spatial location of
different signaling molecules and the spatio-temporal dynamics of signaling
microdomains are now widely acknowledged as key features of biochemical signal
transduction. Reaction-diffusion mechanisms are commonly used to model such
features, often with a heavy reliance on numerical simulations to obtain
results. However, simulations are parameter dependent and may not be able to
provide an understanding of the full range of the system responses. Analytical
approaches on the other hand provide a framework to study the entire phase
space. In this tutorial, we provide a largely analytical method for studying
reaction-diffusion models and analyzing their stability properties. Using two
representative biological examples, we demonstrate how this approach can guide
experimental design.