  More virulent strains of influenza virus subtypes H1N1 appeared in 2007 and
H3N2 in 2011. The amino acid differences from prior less virulent strains
appear to be small when tabulated through sequence alignments and counting site
identities and similarities. Here we show how analyzing fractal hydropathic
forces responsible for globular compaction and modularity quantifies the
mutational origins of increased virulence, and also analyzes receptor sites and
N-linked glycan accretion.


  Thermodynamic tools are well suited to connecting evolution of protein
functionalities to mutations of amino acid sequences, especially for neuronal
network structures. These tools enable one to quantify changes in modular
structure and correlate them with corresponding changes in observable
properties. Here we quantify modular rodent-primate changes in amyloid
precursor protein A4 and \b{eta} amyloid fragments. These are related to
changes in cortical connectivity and to the presence (absence) of plaque
formation in primates (rodents). Two thermodynamic scales are used, descriptive
of water/air protein unfolding (old), or fractal conformational restructuring
(new). These describe complementary aspects of protein activity, at
respectively higher and lower effective temperatures.


  The enormous size and complexity of genotypic sequence space frequently
requires consideration of coarse-grained sequences in empirical models. We
develop scaling relations to quantify the effect of this coarse-graining on
properties of fitness landscapes and evolutionary paths. We first consider
evolution on a simple Mount Fuji fitness landscape, focusing on how the length
and predictability of evolutionary paths scale with the coarse-grained sequence
length and alphabet. We obtain simple scaling relations for both the weak- and
strong-selection limits, with a non-trivial crossover regime at intermediate
selection strengths. We apply these results to evolution on a biophysical
fitness landscape that describes how proteins evolve new binding interactions
while maintaining their folding stability. We combine the scaling relations
with numerical calculations for coarse-grained protein sequences to obtain
quantitative properties of the model for realistic binding interfaces and a
full amino acid alphabet.


  The ongoing effort to detect and characterize physical entanglement in
biopolymers has so far established that knots are present in many globular
proteins and also abound in viral DNA packaged inside bacteriophages. RNA
molecules, on the other hand, have not yet been systematically screened for the
occurrence of physical knots. We have accordingly undertaken the systematic
profiling of the ~6,000 RNA structures present in the protein data bank. The
search identified no more than three deeply-knotted RNA molecules. These are
ribosomal RNAs solved by cryo-em and consist of about 3,000 nucleotides.
Compared to the case of proteins and viral DNA, the observed incidence of RNA
knots is therefore practically negligible. This suggests that either
evolutionary selection, or thermodynamic and kinetic folding mechanisms act
towards minimizing the entanglement of RNA to an extent that is unparalleled by
other types of biomolecules. The properties of the three observed RNA knotting
patterns provide valuable clues for designing RNA sequences capable of
self-tying in a twist-knot fold.


  Analyses of targeted genomic sequencing data from next-generation-sequencing
(NGS) technologies typically involves mapping reads to a reference sequence or
clustering reads. For a number of species a reference genome is not available
so the analyses of targeted sequencing data, for example polymorphic structural
variation caused by mobile elements is difficult; clustering methods are
preferred for such data analysis. Clustering of reads requires a clustering
threshold parameter, which is used to compare and group reads. However,
determining the optimal clustering threshold for a read dataset is challenging
because of different sequence composition, the number of sequences present, and
also the amount of sequencing errors in the dataset. High values of the
clustering threshold parameter can falsely inflate the number of recovered
genomic regions, while low values of clustering threshold can merge reads from
distinct regions into a single cluster. Thus, an algorithm that can empirically
determine clustering threshold is needed. We propose a pipeline for clustering
genomic sequences wherein the clustering threshold is empirically determined
from the NGS data. The optimal threshold is decided based on two internal
clustering measures which assess clusters for small intra-cluster diameters and
large inter-cluster distances. We evaluate the pipeline on two simulated
datasets derived from human genome sequence simulating different genomic
regions and sequencing depth. The total number of clusters obtained from our
pipeline is closer to the actual number of reference sequences when compared to
single round of clustering. Also, the number of clusters whose consensus
sequence matches a corresponding reference sequence is higher in our pipeline.
We observe that the presence of repeat regions affects clustering accuracy.


  Deliberate deceptiveness intended to gain an advantage is commonplace in
human and animal societies. In a social dilemma, an individual may only pretend
to be a cooperator to elicit cooperation from others, while in reality he is a
defector. With this as motivation, we study a simple variant of the
evolutionary prisoner's dilemma game entailing deceitful defectors and
conditional cooperators that lifts the veil on the impact of such two-faced
behavior. Defectors are able to hide their true intentions at a personal cost,
while conditional cooperators are probabilistically successful at identifying
defectors and act accordingly. By focusing on the evolutionary outcomes in
structured populations, we observe a number of unexpected and counterintuitive
phenomena. We show that deceitful behavior may fare better if it is costly, and
that a higher success rate of identifying defectors does not necessarily favor
cooperative behavior. These results are rooted in the spontaneous emergence of
cycling dominance and spatial patterns that give rise to fascinating phase
transitions, which in turn reveal the hidden complexity behind the evolution of
deception.


  Enterococcus faecalis is the third cause of nosocomial infections. To obtain
the first comprehensive view of transcriptional organizations in this
bacterium, we used a modified RNA-seq approach enabling to discriminate primary
from processed 5'RNA ends. We also validated our approach by confirming known
features in Escherichia coli.
  We mapped 559 transcription start sites and 352 processing sites in E.
faecalis. A blind motif search retrieved canonical features of SigA- and
SigN-dependent promoters preceding TSSs mapped. We discovered 95 novel putative
regulatory RNAs, small- and antisense RNAs, and 72 transcriptional antisense
organisations.
  Presented data constitute a significant insight into bacterial RNA landscapes
and a step towards the inference of regulatory processes at transcriptional and
post-transcriptional levels in a comprehensive manner.


  Diabetes mellitus type II affects around 8 percent of the total adult
population in the world. It is the fifth leading cause of death in high income
countries and an epidemic in developing countries. We analyze protein-protein
interaction data of the pancreatic cells for normal and disease states. The
analysis exhibits overall structural similarities in the normal and disease
networks. The important differences are revealed through specific interaction
patterns and eigenvector analyses. The top contributing nodes from localized
eigenvectors as well as those being part of specific interaction patterns turn
out to be significant for the occurrence of the disease. The analysis provides
a direction for further development of novel drugs and therapies in curing the
disease by targeting specific patterns instead of a single node.


  We develop the connection between large deviation theory and more applied
approaches to stochastic hybrid systems by highlighting a common underlying
Hamiltonian structure. A stochastic hybrid system involves the coupling between
a piecewise deterministic dynamical system in $R^d$ and a time-homogeneous
Markov chain on some discrete space $\Gamma$. We assume that the Markov chain
on $\Gamma$ is ergodic, and that the discrete dynamics is much faster than the
piecewise deterministic dynamics (separation of time-scales). Using the
Perron-Frobenius theorem and the calculus-of-variations, we evaluate the rate
function of a large deviation principle in terms of a classical action, whose
Hamiltonian is given by the Perron eigenvalue of a $|\Gamma|$-dimensional
linear equation. The corresponding linear operator depends on the transition
rates of the Markov chain and the nonlinear functions of the piecewise
deterministic system. The resulting Hamiltonian is identical to one derived
using path-integrals and WKB methods. We illustrate the theory by considering
the example of stochastic ion channels.


  Reprogramming is a process of transforming differentiated cells into
pluripotent stem cells by inducing specific modifying factors in the cells.
Reprogramming is a non-equilibrium process involving a collaboration at levels
separated by orders of magnitude in time scale, namely transcription factor
binding/unbinding, protein synthesis/degradation, and epigenetic histone
modification. We propose a model of reprogramming by integrating these
temporally separated processes and show that stable states on the epigenetic
landscape should be viewed as a superposition of basin minima generated in the
fixed histone states. Slow histone modification is responsible for the narrow
valleys connecting the pluripotent and differentiated states on the epigenetic
landscape, and the pathways which largely overlap with these valleys explain
the observed heterogeneity of latencies in reprogramming. We show that histone
dynamics also creates an intermediary state observed in experiments. A change
in the mechanism of histone modification alters the pathway to bypass the
barrier, thereby accelerating the reprogramming and reducing the heterogeneity
of latencies.


  In the research scope of forest stand self-thinning, the analysis reveals a
broad picture in which the -3/2 rule takes a definite and special place. The
application of the simple geometrical model to the Douglas-fir and Scots pine
data suggests that the slope of the self-thinning curve will not remain
constant during the course of growth and self-thinning of a single forest
stand. Most probable, at the initial stages of stand growth the slope will be
less than -3/2 and at old ages of the stand the slope will be higher than -3/2.
Inevitably, a time will come when the slope is exactly equals -3/2. In other
words, the slope -3/2 is an obligatory state in the course of self-thinning of
a forest stand. At the very time of -3/2 slope two particular features coincide
with it. One is that the total bole surface area remains constant. The length
of the constancy stage would probably vary with species, initial stand
densities, their spatial arrangements, conditions of growth, and other specific
factors. Another feature of the time is that a geometric similarity in the
growth of the forest stand takes place, which is not in a contradiction with
the -3/2 rule as it had been formulated by its authors. To put it shortly, the
slope -3/2: i) is a very specific and obligatory state in the process of forest
stand growth and ii) is not an asymptote but rather a transitional point (span)
in the time of growth. These two assertions may be called a transitional status
of the -3/2 rule. The geometric model of a forest stand (Gavrikov, 2014) has
proved to be rather helpful at analyzing of real forest stand structure and
dynamics. Despite of its extreme simplicity (it uses cones as representations
of trees) the model looks like having enough similarity with real even-aged
forests since the model's predictions are often reasonably close to measured
values of power exponents.


  A ternary permutation constraint satisfaction problem (CSP) is specified by a
subset Pi of the symmetric group S_3. An instance of such a problem consists of
a set of variables V and a set of constraints C, where each constraint is an
ordered triple of distinct elements from V. The goal is to construct a linear
order alpha on V such that, for each constraint (a,b,c) in C, the ordering of
a,b,c induced by alpha is in Pi. Excluding symmetries and trivial cases there
are 11 such problems, and their complexity is well known. Here we consider the
variant of the problem, denoted 2-Pi, where we are allowed to construct two
linear orders alpha and beta and each constraint needs to be satisfied by at
least one of the two. We give a full complexity classification of all 11 2-Pi
problems, observing that in the switch from one to two linear orders the
complexity landscape changes quite abruptly and that hardness proofs become
rather intricate. We then focus on one of the 11 problems in particular, which
is closely related to the '2-Caterpillar Compatibility' problem in the
phylogenetics literature. We show that this particular CSP remains hard on
three linear orders, and also in the biologically relevant case when we swap
three linear orders for three phylogenetic trees, yielding the '3-Tree
Compatibility' problem. Due to the biological relevance of this problem we also
give extremal results concerning the minimum number of trees required, in the
worst case, to satisfy a set of rooted triplet constraints on n leaf labels.


  Here we generalize our previous model of molecular motors trafficking
subdiffusing cargos in viscoelastic cytosol by (i) including mechanochemical
coupling between cyclic conformational fluctuations of the motor protein driven
by the reaction of ATP hydrolysis and its translational motion within the
simplest two-state model of hand-over-hand motion of kinesin, and also (ii) by
taking into account the anharmonicity of the tether between the motor and cargo
(its maximally possible extension length). It is shown that the major earlier
results such as occurrence of normal versus anomalous transport depending on
the amplitude of binding potential, cargo size and the motor turnover frequency
not only survive in this more realistic model, but the results also look very
similar for the correspondingly adjusted parameters. However, this more
realistic model displays a substantially larger thermodynamic efficiency due to
a bidirectional mechanochemical coupling. For realistic parameters, the maximal
thermodynamic efficiency can be transiently about 50% as observed for kinesins,
and even larger, surprisingly also in a novel strongly anomalous (sub)transport
regime, where the motor enzymatic turnovers become also anomalously slow and
cannot be characterized by a turnover rate. Here anomalously slow dynamics of
the cargo enforces anomalously slow cyclic kinetics of the motor protein.


  The magnitude and variability of Earth's biodiversity have puzzled scientists
ever since paleontologic fossil databases became available. We identify and
study a model of interdependent species where both endogenous and exogenous
impacts determine the nonstationary extinction dynamics. The framework provides
an explanation for the qualitative difference of marine and continental
biodiversity growth. In particular, the stagnation of marine biodiversity may
result from a global transition from an imbalanced to a balanced state of the
species dependency network. The predictions of our framework are in agreement
with paleontologic databases.


  Here we present ComPPI, a cellular compartment specific database of proteins
and their interactions enabling an extensive, compartmentalized protein-protein
interaction network analysis (http://ComPPI.LinkGroup.hu). ComPPI enables the
user to filter biologically unlikely interactions, where the two interacting
proteins have no common subcellular localizations and to predict novel
properties, such as compartment-specific biological functions. ComPPI is an
integrated database covering four species (S. cerevisiae, C. elegans, D.
melanogaster and H. sapiens). The compilation of nine protein-protein
interaction and eight subcellular localization data sets had four curation
steps including a manually built, comprehensive hierarchical structure of more
than 1600 subcellular localizations. ComPPI provides confidence scores for
protein subcellular localizations and protein-protein interactions. ComPPI has
user-friendly search options for individual proteins giving their subcellular
localization, their interactions and the likelihood of their interactions
considering the subcellular localization of their interacting partners.
Download options of search results, whole proteomes, organelle-specific
interactomes and subcellular localization data are available on its website.
Due to its novel features, ComPPI is useful for the analysis of experimental
results in biochemistry and molecular biology, as well as for proteome-wide
studies in bioinformatics and network science helping cellular biology,
medicine and drug design.


  The class of epistatic fitness landscapes is much more diverse than the class
of non-epistatic landscapes, and so it stands to reason that there exist
dynamical phenomena that can only be realized in the presence of epistasis.
Here, we compare evolutionary dynamics on all finite epistatic landscapes
versus all finite non-epistatic landscapes, under weak mutation. We first
analyze the mean fitness trajectory - that is, the time course of the expected
fitness of a population. We show that for any epistatic fitness landscape and
starting genotype, there always exists a non-epistatic fitness landscape and
starting genotype that produces the exact same mean fitness trajectory. Thus,
surprisingly, the space of mean fitness trajectories that can be realized by
epistatic landscapes is no more diverse than the space of mean fitness
trajectories that can be realized by non-epistatic landscapes. On the other
hand, we show that epistatic fitness landscapes can produce dynamics in the
time-evolution of the variance in fitness across replicate populations and in
the time-evolution of the expected number of substitutions that cannot be
produced by any non-epistatic landscape. These results on identifiability have
implications for efforts to infer epistasis from the types of data often
measured in experimental populations.


  The problem addressed here is that of simultaneous treatment of several gene
expression datasets, possibly collected under different experimental conditions
and/or platforms. Using robust statistics, a large scale statistical analysis
has been conducted over $20$ datasets downloaded from the Gene Expression
Omnibus repository. The differences between datasets are compared to the
variability inside a given dataset. Evidence that meaningful biological
information can be extracted by merging different sources is provided.


  For practical construction of complex synthetic genetic networks able to
perform elaborate functions it is important to have a pool of relatively simple
"bio-bricks" with different functionality which can be compounded together. To
complement engineering of very different existing synthetic genetic devices
such as switches, oscillators or logical gates, we propose and develop here a
design of synthetic multiple input distributed classifier with learning
ability. Proposed classifier will be able to separate multi-input data, which
are inseparable for single input classifiers. Additionally, the data classes
could potentially occupy the area of any shape in the space of inputs. We study
two approaches to classification, including hard and soft classification and
confirm the schemes of genetic networks by analytical and numerical results.


  Motivated by the fact that the same social dilemma can be perceived
differently by different players, we here study evolutionary multigames in
structured populations. While the core game is the weak prisoner's dilemma, a
fraction of the population adopts either a positive or a negative value of the
sucker's payoff, thus playing either the traditional prisoner's dilemma or the
snowdrift game. We show that the higher the fraction of the population adopting
a different payoff matrix, the more the evolution of cooperation is promoted.
The microscopic mechanism responsible for this outcome is unique to structured
populations, and it is due to the payoff heterogeneity, which spontaneously
introduces strong cooperative leaders that give rise to an asymmetric strategy
imitation flow in favor of cooperation. We demonstrate that the reported
evolutionary outcomes are robust against variations of the interaction network,
and they also remain valid if players are allowed to vary which game they play
over time. These results corroborate existing evidence in favor of
heterogeneity-enhanced network reciprocity, and they reveal how different
perceptions of social dilemmas may contribute to their resolution.


  Protein machines often exhibit long range interplay between different sites
in order to achieve their biological tasks. We investigate and characterize the
non--linear energy localization and the basic mechanisms of energy transfer in
protein devices. By studying two different model protein machines, with
different biological functions, we show that genuinely non--linear phenomena
are responsible for energy transport between the different machine sites
involved in the biological functions. The energy transfer turns out to be
extremely efficient from an energetic point of view: by changing the energy
initially provided to the model device, we identify a well defined range of
energies where the time for the energy transport to occur is minimal and the
amount of transferred energy is maximum. Furthermore, by introducing an
implicit solvent, we show that the energy is localized on the internal residues
of the protein structure, thus minimizing the dissipation.


  Elementary flux modes (EFMs) are pathways through a metabolic reaction
network that connect external substrates to products. Using EFMs, a metabolic
network can be transformed into its macroscopic counterpart, in which the
internal metabolites have been eliminated and only external metabolites remain.
In EFMs-based metabolic flux analysis (MFA) experimentally determined external
fluxes are used to estimate the flux of each EFM. It is in general prohibitive
to enumerate all EFMs for complex networks, since the number of EFMs increases
rapidly with network complexity. In this work we present an optimization-based
method that dynamically generates a subset of EFMs and solves the EFMs-based
MFA problem simultaneously. The obtained subset contains EFMs that contribute
to the optimal solution of the EFMs-based MFA problem. The usefulness of our
method was examined in a case-study using data from a Chinese hamster ovary
cell culture and two networks of varied complexity. It was demonstrated that
the EFMs-based MFA problem could be solved at a low computational cost, even
for the more complex network. Additionally, only a fraction of the total number
of EFMs was needed to compute the optimal solution.


  This work introduces novel unconditionally stable operator splitting methods
for solving the time dependent nonlinear Poisson-Boltzmann (NPB) equation for
the electrostatic analysis of solvated biomolecules. In a pseudo-transient
continuation solution of the NPB equation, a long time integration is needed to
reach the steady state. This calls for time stepping schemes that are stable
and accurate for large time increments. The existing alternating direction
implicit (ADI) methods for the NPB equation are known to be conditionally
stable, although being fully implicit. To overcome this difficulty, we propose
several new operator splitting schemes, in both multiplicative and additive
styles, including locally one-dimensional (LOD) schemes and additive operator
splitting (AOS) schemes. The proposed schemes become much more stable than the
ADI methods, and some of them are indeed unconditionally stable in dealing with
solvated proteins with source singularities and non-smooth solutions.
Numerically, the orders of convergence in both space and time are found to be
one. Nevertheless, the precision in calculating the electrostatic free energy
is low, unless a small time increment is used. Further accuracy improvements
are thus considered. After acceleration, the optimized LOD method can produce a
reliable energy estimate by integrating for a small and fixed number of time
steps. Since one only needs to solve a tridiagonal linear system in each
independent one dimensional process, the overall computation is very efficient.
The unconditionally stable LOD method scales linearly with respect to the
number of atoms in the protein studies, and is over 20 times faster than the
conditionally stable ADI methods.


  (An updated version of this manuscript has been accepted to Scientific
Reports in 2016, please refer to http://www.nature.com/articles/srep31900)
  The highly anticipated transition from next generation sequencing (NGS) to
third generation sequencing (3GS) has been difficult primarily due to high
error rates and excessive sequencing cost. The high error rates make the
assembly of long erroneous reads of large genomes challenging because existing
software solutions are often overwhelmed by error correction tasks. Here we
report a hybrid assembly approach that simultaneously utilizes NGS and 3GS data
to address both issues. We gain advantages from three general and basic design
principles: (i) Compact representation of the long reads lead to efficient
alignments. (ii) Base-level errors can be skipped; structural errors need to be
detected and corrected. (iii) Structurally correct 3GS reads are assembled and
polished. In our implementation, preassembled NGS contigs are used to derive
the compact representation of the long reads, which established an algorithmic
conversion from a de Bruijn graph to an overlap graph, the two major assembly
paradigms. Moreover, since NGS and 3GS data can compensate each other, our
hybrid assembly approach reduces both of their sequencing requirements.
Experiments show that our software is able to assemble mammalian-sized genomes
orders of magnitude more efficiently in time than existing methods, while
saving about half of the sequencing cost.


  Biogeographical regions (geographically distinct assemblages of species and
communities) constitute a cornerstone for ecology, biogeography, evolution and
conservation biology. Species turnover measures are often used to quantify
biodiversity patterns, but algorithms based on similarity and clustering are
highly sensitive to common biases and intricacies of species distribution data.
Here we apply a community detection approach from network theory that
incorporates complex, higher order presence-absence patterns. We demonstrate
the performance of the method by applying it to all amphibian species in the
world (c. 6,100 species), all vascular plant species of the USA (c. 17,600),
and a hypothetical dataset containing a zone of biotic transition. In
comparison with current methods, our approach tackles the challenges posed by
transition zones and succeeds in identifying a larger number of commonly
recognised biogeographical regions. This method constitutes an important
advance towards objective, data derived identification and delimitation of the
world's biogeographical regions.


  Robust control theory has been successfully applied to numerous real-world
problems using a small set of devices called {\it controllers}. However, the
real systems represented by networks contain unreliable components and modern
robust control engineering has not addressed the problem of structural changes
on a large network. Here, we introduce the concept of structurally robust
control of complex networks and provide a concrete example using an algorithmic
framework that is widely applied in engineering. The developed analytical
tools, computer simulations and real network analyses lead herein to the
discovery that robust control can be achieved in scale-free networks with
exactly the same order of controllers required in a standard non-robust
configuration by adjusting only the minimum degree. The presented methodology
also addresses the probabilistic failure of links in real systems, such as
neural synaptic unreliability in {\it C. elegans}, and suggests a new direction
to pursue in studies of complex networks in which control theory has a role.


  Comprehensive models of stochastic, clonally reproducing populations are
defined in terms of general branching processes, allowing birth during maternal
life, as for higher organisms, or by splitting, as in cell division. The
populations are assumed to start small, by mutation or immigration, reproduce
supercritically while smaller than the habitat carrying capacity but
subcritically above it. Such populations establish themselves with a
probability wellknown from branching process theory. Once established, they
grow up to a band around the carrying capacity in a time that is logarithmic in
the latter, assumed large. There they prevail during a time period whose
duration is exponential in the carrying capacity. Even populations whose life
style is sustainble in the sense that the habitat carrying capacity is not
eroded but remains the same, ultimately enter an extinction phase, which again
lasts for a time logarithmic in the carrying capacity. However, if the habitat
can carry a population which is large, say millions of individuals, and it
manages to avoid early extinction, time in generations to extinction will be
exorbitantly long, and during it, population composition over ages, types,
lineage etc. will have time to stabilise. This paper aims at an exhaustive
description of the life cycle of such populations, from inception to
extinction, extending and overviewing earlier results. We shall also say some
words on persistence times of populations with smaller carrying capacities and
short life cycles, where the population may indeed be in danger in spite of not
eroding its environment.


  We study the regulating mechanism of p53 on the properties of cell cycle
dynamics in the light of the proposed model of interacting p53 and cell cycle
networks via p53. Irradiation (IR) introduce to p53 compel p53 dynamics to
suffer different phases, namely oscillating and oscillation death (stabilized)
phases. The IR induced p53 dynamics undergo collapse of oscillation with
collapse time \Delta t which depends on IR strength. The stress p53 via IR
drive cell cycle molecular species MPF and cyclin dynamics to different states,
namely, oscillation death, oscillations of periods, chaotic and sustain
oscillation in their bifurcation diagram. We predict that there could be a
critical \Delta t_c induced by p53 via IR_c, where, if \Delta t < \Delta t_c
the cell cycle may come back to normal state, otherwise it will go to cell
cycle arrest (apoptosis).


  We present a new method for exploring cancer gene expression data based on
tools from algebraic topology. Our method selects a small relevant subset from
tens of thousands of genes while simultaneously identifying nontrivial higher
order topological features, i.e., holes, in the data. We first circumvent the
problem of high dimensionality by dualizing the data, i.e., by studying genes
as points in the sample space. Then we select a small subset of the genes as
landmarks to construct topological structures that capture persistent, i.e.,
topologically significant, features of the data set in its first homology
group. Furthermore, we demonstrate that many members of these loops have been
implicated for cancer biogenesis in scientific literature. We illustrate our
method on five different data sets belonging to brain, breast, leukemia, and
ovarian cancers.


  Network alignment can be used to transfer functional knowledge between
conserved regions of different networks. Typically, existing methods use a node
cost function (NCF) to compute similarity between nodes in different networks
and an alignment strategy (AS) to find high-scoring alignments with respect to
the total NCF over all aligned nodes (or node conservation). But, they then
evaluate quality of their alignments via some other measure that is different
than the node conservation measure used to guide the alignment construction
process. Typically, one measures the amount of conserved edges, but only after
alignments are produced. Hence, a recent attempt aimed to directly maximize the
amount of conserved edges while constructing alignments, which improved
alignment accuracy. Here, we aim to directly maximize both node and edge
conservation during alignment construction to further improve alignment
accuracy. For this, we design a novel measure of edge conservation that (unlike
existing measures that treat each conserved edge the same) weighs each
conserved edge so that edges with highly NCF-similar end nodes are favored. As
a result, we introduce a novel AS, Weighted Alignment VotEr (WAVE), which can
optimize any measures of node and edge conservation, and which can be used with
any NCF or combination of multiple NCFs. Using WAVE on top of established
state-of-the-art NCFs leads to superior alignments compared to the existing
methods that optimize only node conservation or only edge conservation or that
treat each conserved edge the same. And while we evaluate WAVE in the
computational biology domain, it is easily applicable in any domain.


  We present an explicit solution to a classic model of cell-population growth
introduced by Luria and Delbrueck 70 years ago to study the emergence of
mutations in bacterial populations. In this model a wild-type population is
assumed to grow exponentially in a deterministic fashion. Proportional to the
wild-type population size, mutants arrive randomly and initiate new
sub-populations of mutants that grows stochastically according to a
supercritical birth and death process. We give an exact expression for the
generating function of the total number of mutants at a given wild type
population size. We present a simple expression for the probability of finding
no mutants, and a recursion formula for the probability of finding a given
number of mutants. In the "large population-small mutation"-limit we recover
recent results of Kessler and Levin for a fully stochastic version of the
process.


  One of the most impressive features of moving animal groups is their ability
to perform sudden coherent changes in travel direction. While this collective
decision can be a response to an external perturbation, such as the presence of
a predator, recent studies show that such directional switching can also emerge
from the intrinsic fluctuations in the individual behaviour. However, the cause
and the mechanism by which such collective changes of direction occur are not
fully understood yet. Here, we present an experimental study of spontaneous
collective turns in natural flocks of starlings. We employ a recently developed
tracking algorithm to reconstruct three-dimensional trajectories of each
individual bird in the flock for the whole duration of a turning event. Our
approach enables us to analyze changes in the individual behavior of every
group member and reveal the emergent dynamics of turning. We show that
spontaneous turns start from individuals located at the elongated edges of the
flocks, and then propagate through the group. We find that birds on the edges
deviate from the mean direction of motion much more frequently than other
individuals, indicating that persistent localized fluctuations are the crucial
ingredient for triggering a collective directional change. Finally, we
quantitatively show that birds follow equal radius paths during turning
allowing the flock to change orientation and redistribute risky locations among
group members. The whole process of turning is a remarkable example of how a
self-organized system can sustain collective changes and reorganize, while
retaining coherence.


  The strand symmetric model is a phylogenetic model designed to reflect the
symmetry inherent in the double-stranded structure of DNA. We show that the set
of known phylogenetic invariants for the general strand symmetric model of the
three leaf claw tree entirely defines the ideal. This knowledge allows one to
determine the vanishing ideal of the general strand symmetric model of any
trivalent tree. Our proof of the main result is computational. We use the fact
that the Zariski closure of the strand symmetric model is the secant variety of
a toric variety to compute the dimension of the variety. We then show that the
known equations generate a prime ideal of the correct dimension using
elimination theory.


  In the last few decades, numerous experiments have shown that humans do not
always behave so as to maximize their material payoff. Cooperative behavior
when non-cooperation is a dominant strategy (with respect to the material
payoffs) is particularly puzzling. Here we propose a novel approach to explain
cooperation, assuming what Halpern and Pass (2013) call "translucent players".
Typically, players are assumed to be "opaque", in the sense that a deviation by
one player does not affect the strategies used by other players. But a player
may believe that if he switches from one strategy to another, the fact that he
chooses to switch may be visible to the other players. For example, if he
chooses to defect in Prisoner's Dilemma, the other player may sense his guilt.
We show that by assuming translucent players, we can recover many of the
regularities observed in human behavior in well-studied games such as
Prisoner's Dilemma, Traveler's Dilemma, Bertrand Competition, and the Public
Goods game.


  This paper compares simulation results and experimental published data
concerning the interaction between $\alpha_2\beta_1$ integrin and collagen.This
phenomenon is characterized by a competition among steric, reciprocal
orientation constraints of the two structures and long range electrostatic
force. The published experimental results shown that the
$\textrm{k}_{\textrm{on}}$ rate, for the 42-mer collagen fragment and the
$\alpha_2\textrm{I}$ integrin domain complex, ranges in order of magnitude from
$10^4-10^5$ M$^{-1}$s$^{-1}$. This is the lower bound of the interval defined
by the diffusion limited regime and the orientational constraint. The
electrostatic affinity between 21-mer collagen fragment and the
$\alpha_2\textrm{I}$ integrin domain due to the divalent cations
($\textrm{Mg}^{2+}$, $\textrm{Co}^{2+}$, $\textrm{Mn}^{2+}$ and
$\textrm{Ca}^{2+}$) has been expressed in terms of the Pearson correlation
coefficient between the electrostatic desolvation potential of the ligand and
the electrostatic interaction potential of the receptor. The obtained values of
the Pearson coefficient clearly reveals the complementarity between the
$\alpha_2\textrm{I}$ integrin domain and 21-mer collagen in the metal ion
dependent site. Atomistic Brownian dynamics simulations applied in the
diffusion limited regime, correctly reproduce the binding kinetics, of the
$\alpha_2\textrm{I}$ integrin domain and the 21-mer collagen in presence of
$\textrm{Mg}^{2+}$, $\textrm{Co}^{2+}$ and $\textrm{Mn}^{2+}$, according to
published experimental data.


  Fundamental properties of macroscopic gene-mating dynamic evolutionary
systems are investigated. We focus on a single locus, any number of alleles in
a two-gender dioecious population, for a large class of systems within
population genetics. Our governing equations are time-dependent differential
equations labeled by a set of genotype frequencies. Our equations are uniquely
derived from 4 assumptions within any population: (1) a closed system; (2)
average-and-random mating process (mean-field behavior); (3) Mendelian
inheritance; (4) exponential growth/death. Although our equations are nonlinear
with time-evolutionary dynamics, we have obtained an exact analytic
time-dependent solution and an exactly solvable model. From the
phenomenological viewpoint, any initial parameter of genotype frequencies of a
closed system will eventually approach a stable fixed point. Under time
evolution, we show (1) the monotonic behavior of genotype frequencies, (2) any
genotype or allele that appears in the population will never become extinct,
(3) the Hardy-Weinberg law, and (4) the global stability without chaos in the
parameter space. To demonstrate the experimental evidence for our theory, as an
example, we show a mapping from the data of blood type genotype frequencies of
world ethnic groups to our stable fixed-point solutions. This fixed-point
Hardy-Weinberg manifold, attracting any initial point in any Euclidean fiber
bounded within the genotype frequency space to the fixed point where this fiber
is attached. The stable base manifold and its attached fibers form a fiber
bundle, which fills in the whole genotype frequency space completely. We can
define the genetic distance of two populations as their geodesic distance on
the equilibrium manifold. The modification of our theory under the process of
natural selection and mutation is addressed.


  Traditional mathematical models of photosynthesis are based on mass action
kinetics of light reactions. This approach requires the modeller to enumerate
all the possible state combinations of the modelled chemical species. This
leads to combinatorial explosion in the number of reactions although the
structure of the model could be expressed more compactly. We explore the use of
rule-based modelling, in particular, a simplified variant of Kappa, to
compactly capture and automatically reduce existing mathematical models of
photosynthesis. Finally, the reduction procedure is implemented in BioNetGen
language and demonstrated on several ODE models of photosynthesis processes.
This is an extended version of the paper published in proceedings of 5th
International Workshop on Static Analysis and Systems Biology (SASB) 2014.


  A differential equations model is developed for the 2014 Ebola epidemics in
Sierra Leone, Liberia, and Guinea. The model describes the dynamic interactions
of the susceptible and infected populations of these countries. The model
incorporates the principle features of contact tracing, namely, the number of
contacts per identified infectious case, the likelihood that a traced contact
is either incubating or infectious, and the efficiency of the contact tracing
process.The model is first fitted to current cumulative reported case data in
each country. The data fitted simulations are then projected forward in time,
with varying parameter regimes corresponding to contact tracing efficiencies.
These projections quantify the importance of the identification, isolation, and
contact tracing processes for containment of the epidemics.


  Earlier this year, we published a scathing critique of a paper by Mendez et
al. (2013) in which the claim was made that a Y chromosome was 237,000-581,000
years old. Elhaik et al. (2014) also attacked a popular article in Scientific
American by the senior author of Mendez et al. (2013), whose title was "Sex
with other human species might have been the secret of Homo sapiens's [sic]
success" (Hammer 2013). Five of the 11 authors of Mendez et al. (2013) have now
written a "rebuttal," and we were allowed to reply. Unfortunately, our reply
was censored for being "too sarcastic and inflamed." References were removed,
meanings were castrated, and a dedication in the Acknowledgments was deleted.
Now, that the so-called rebuttal by 45% of the authors of Mendez et al. (2013)
has been published together with our vasectomized reply, we decided to make
public our entire reply to the so called "rebuttal." In fact, we go one step
further, and publish a version of the reply that has not even been
self-censored.


  We study an open-chain totally asymmetric exclusion process (TASEP) with
stochastic gates present at the two boundaries. The gating dynamics has been
modeled keeping the physical system of ion-channel gating in mind. These gates
can randomly switch between an open state and a closed state. In the open
state, the gates are highly permeable such that any particle arriving at the
gate immediately passes through. In the closed state, a particle gets trapped
at the gate and cannot pass through until the gate switches open again. We
calculate the phase-diagram of the system and find important and non-trivial
differences with the phase-diagram of a regular open-chain TASEP. In
particular, depending on switching rates of the two gates, the system may or
may not admit a maximal current phase. Our analytic calculation within
mean-field theory captures the main qualitative features of our Monte Carlo
simulation results. We also perform a refined mean-field calculation where the
correlations at the boundaries are taken into account. This theory shows
significantly better quantitative agreement with our simulation results.


  Neutral models which assume ecological equivalence between species provide
null models for community assembly. In Hubbell's Unified Neutral Theory of
Biodiversity (UNTB), many local communities are connected to a single
metacommunity through differing immigration rates. Our ability to fit the full
multi-site UNTB has hitherto been limited by the lack of a computationally
tractable and accurate algorithm. We show that a large class of neutral models
with this mainland-island structure but differing local community dynamics
converge in the large population limit to the hierarchical Dirichlet process.
Using this approximation we developed an efficient Bayesian fitting strategy
for the multi-site UNTB. We can also use this approach to distinguish between
neutral local community assembly given a non-neutral metacommunity distribution
and the full UNTB where the metacommunity too assembles neutrally. We applied
this fitting strategy to both tropical trees and a data set comprising 570,851
sequences from 278 human gut microbiomes. The tropical tree data set was
consistent with the UNTB but for the human gut neutrality was rejected at the
whole community level. However, when we applied the algorithm to gut microbial
species within the same taxon at different levels of taxonomic resolution, we
found that species abundances within some genera were almost consistent with
local community assembly. This was not true at higher taxonomic ranks. This
suggests that the gut microbiota is more strongly niche constrained than
macroscopic organisms, with different groups adopting different functional
roles, but within those groups diversity may at least partially be maintained
by neutrality.We also observed a negative correlation between body mass index
and immigration rates within the family Ruminococcaceae.


  Translation is a crucial step in gene expression. During translation,
macromolecules called ribosomes "read" the mRNA strand in a sequential manner
and produce a corresponding protein. Translation is known to consume most of
the cell's energy. Maximizing the protein production rate in mRNA translation,
subject to the bounded biomolecullar budget, is thus an important problem in
both biology and biotechnology. We consider this problem using a mathematical
model for mRNA translation called the ribosome flow model (RFM). For an mRNA
strand with $n$ sites the RFM includes $n$ state-variables that encode the
normalized ribosomal density at each site, and $n+1$ positive parameters: the
initiation rate and elongation rates along the chain. An affine constraint on
these rates is used to model the bounded cellular budget. We show that for a
homogeneous constraint the rates that maximize the steady-state protein
production rate have a special structure. They are symmetric with respect to
the middle of the chain, and monotonically increase as we move towards the
center of the chain. The ribosomal densities corresponding to the optimal rates
monotonically decrease along the chain. We discuss some of the biological
implications of these results.


  Recent advances in imaging technology now provide us with 3D images of
developing organs. These can be used to extract 3D geometries for simulations
of organ development. To solve models on growing domains, the displacement
fields between consecutive image frames need to be determined. Here we develop
and evaluate different landmark-free algorithms for the determination of such
displacement fields from image data. In particular, we examine minimal
distance, normal distance, diffusion-based and uniform mapping algorithms and
test these algorithms with both synthetic and real data in 2D and 3D. We
conclude that in most cases the normal distance algorithm is the method of
choice and wherever it fails, diffusion-based mapping provides a good
alternative.


  Intrinsically disordered proteins are fascinating the community of protein
science since the last decade, at least. There is a well-established line of
research that intends to reveal the crucial role played by intrinsically
disordered proteins (IDPs) in the development of human diseases. The main
argument is that IDPs are differentially more present in groups of
disease-related proteins. In this note we compare the frequency of disorder in
human proteins, both disease-related and not. The frequency of disorder is
comparable in the two sub-groups of proteins. Disorder is widespread in human
proteins, but it is not a specific pre-requisite of proteins involved in the
development of cancer, cardiovascular diseases, diabetes and neurodegenerative
diseases. A tendency of cancer-related proteins to be statistically more
disordered than the rest of human proteins is confirmed.


  In this paper we consider epidemic models of directly transmissible SIR
(susceptible $\to$ infective $\to$ recovered) and SEIR (with an additional
latent class) infections in fully-susceptible populations with a social
structure, consisting either of households or of households and workplaces. We
review most reproduction numbers defined in the literature for these models,
including the basic reproduction number $R_0$ introduced in the companion paper
of this, for which we provide a simpler, more elegant derivation. Extending
previous work, we provide a complete overview of the inequalities among these
reproduction numbers and resolve some open questions. Special focus is put on
the exponential-growth-associated reproduction number $R_r$, which is loosely
defined as the estimate of $R_0$ based on the observed exponential growth of an
emerging epidemic obtained when the social structure is ignored. We show that
for the vast majority of the models considered in the literature $R_r \geq R_0$
when $R_0 \ge 1$ and $R_r \leq R_0$ when $R_0 \le 1$. We show that, in contrast
to models without social structure, vaccination of a fraction $1-1/R_0$ of the
population, chosen uniformly at random, with a perfect vaccine is usually
insufficient to prevent large epidemics. In addition, we provide significantly
sharper bounds than the existing ones for bracketing the critical vaccination
coverage between two analytically tractable quantities, which we illustrate by
means of extensive numerical examples.


  Cooperation is usually represented as a Prisoner's Dilemma game. Although
individual self-interest may not favour cooperation, cooperation can evolve if,
for example, players interact multiple times adjusting their behaviour
accordingly to opponent's previous action. To analyze population dynamics,
replicator equation has been widely used under several versions. Although it is
usually stated that a strategy called Generous-tit-for-tat is the winner within
the reactive strategies set, here we show that this result depends on
replicator's version and on the number of available strategies, stemming from
the fact that a dynamics system is also defined by the number of available
strategies and not only by the model version. Using computer simulations and
analytical arguments, we show that Generous-tit-for-tat victory is found only
if the number of strategies available is not too large, with defection winning
otherwise.


  Wealthy individuals may be less tempted to defect than those with
comparatively low payoffs. To take this into consideration, we introduce
coevolutionary success-driven multigames in structured populations. While the
core game is always the weak prisoner's dilemma, players whose payoffs from the
previous round exceed a threshold adopt only a minimally low temptation to
defect in the next round. Along with the strategies, the perceived strength of
the social dilemma thus coevolves with the success of each individual player.
We show that the lower the threshold for using the small temptation to defect,
the more the evolution of cooperation is promoted. Importantly, the promotion
of cooperation is not simply due to a lower average temptation to defect, but
rather due to a dynamically reversed direction of invasion along the interfaces
that separate cooperators and defectors on regular networks. Conversely, on
irregular networks, in the absence of clear invasion fronts, the promotion of
cooperation is due to intermediate-degree players. At sufficiently low
threshold values, these players accelerate the erosion of defectors and
significantly shorten the fixation time towards more cooperative stationary
states. Coevolutionary multigames could thus be the new frontier for the swift
resolution of social dilemmas.


  To maintain a constant cell size, dividing cells have to coordinate cell
cycle events with cell growth. This coordination has for long been supposed to
rely on the existence of size thresholds determining cell cycle progression
[1]. In budding yeast, size is controlled at the G1/S transition [11]. In
agreement with this hypothesis, the size at birth influences the time spent in
G1: smaller cells have a longer G1 period [3]. Nevertheless, even though cells
born smaller have a longer G1, the compensation is imperfect and they still bud
at smaller cell sizes. In bacteria, several recent studies have shown that the
incremental model of size control, in which size is controlled by addition of a
constant volume (in contrast to a size threshold), is able to quantitatively
explain the experimental data on 4 different bacterial species [6, 5, 6, 7].
Here, we report on experimental results for the budding yeast Saccharomyces
cerevisiae, finding, surprisingly, that cell size control in this organism is
very well described by the incremental model, suggesting a common strategy for
cell size control with bacteria. Additionally, we argue that for S. cerevisiae
the volume increment is not added from birth to division, but rather between
two budding events.


  PLINK 1 is a widely used open-source C/C++ toolset for genome-wide
association studies (GWAS) and research in population genetics. However, the
steady accumulation of data from imputation and whole-genome sequencing studies
has exposed a strong need for even faster and more scalable implementations of
key functions. In addition, GWAS and population-genetic data now frequently
contain probabilistic calls, phase information, and/or multiallelic variants,
none of which can be represented by PLINK 1's primary data format.
  To address these issues, we are developing a second-generation codebase for
PLINK. The first major release from this codebase, PLINK 1.9, introduces
extensive use of bit-level parallelism, O(sqrt(n))-time/constant-space
Hardy-Weinberg equilibrium and Fisher's exact tests, and many other algorithmic
improvements. In combination, these changes accelerate most operations by 1-4
orders of magnitude, and allow the program to handle datasets too large to fit
in RAM. This will be followed by PLINK 2.0, which will introduce (a) a new data
format capable of efficiently representing probabilities, phase, and
multiallelic variants, and (b) extensions of many functions to account for the
new types of information.
  The second-generation versions of PLINK will offer dramatic improvements in
performance and compatibility. For the first time, users without access to
high-end computing resources can perform several essential analyses of the
feature-rich and very large genetic datasets coming into use.


  We consider the relationship between stationary distributions for stochastic
models of reaction systems and Lyapunov functions for their deterministic
counterparts. Specifically, we derive the well known Lyapunov function of
reaction network theory as a scaling limit of the non-equilibrium potential of
the stationary distribution of stochastically modeled complex balanced systems.
We extend this result to general birth-death models and demonstrate via example
that similar scaling limits can yield Lyapunov functions even for models that
are not complex or detailed balanced, and may even have multiple equilibria.


  We investigate in vitro fibrillation kinetics of the hormone peptide glucagon
at various concentrations using confocal microscopy and determine the glucagon
fibril persistence length $60 \mu\textrm{m}$. At all concentrations we observe
that periods of individual fibril growth are interrupted by periods of stasis.
The growth probability is large at high and low concentrations and is reduced
for intermediate glucagon concentrations. To explain this behavior we propose a
simple model, where fibrils come in two forms, one built entirely from glucagon
monomers and one entirely from glucagon trimers. The opposite building blocks
act as fibril growth blockers, and this generic model reproduces experimental
behavior well.


  An adaptive network model using SIS epidemic propagation with link-type
dependent link activation and deletion is considered. Bifurcation analysis of
the pairwise ODE approximation and the network-based stochastic simulation is
carried out, showing that three typical behaviours may occur; namely,
oscillations can be observed besides disease-free or endemic steady states. The
oscillatory behaviour in the stochastic simulations is studied using Fourier
analysis, as well as through analysing the exact master equations of the
stochastic model. A compact pairwise approximation for the dynamic network case
is also developed and, for the case of link-type independent rewiring, the
outcome of epidemics and changes in network structure are concurrently
presented in a single bifurcation diagram. By going beyond simply comparing
simulation results to mean-field models, our approach yields deeper insights
into the observed phenomena and help better understand and map out the
limitations of mean-field models.


  Network alignment aims to find regions of topological or functional
similarities between networks. In computational biology, it can be used to
transfer biological knowledge from a well-studied species to a poorly-studied
species between aligned network regions. Typically, existing network aligners
first compute similarities between nodes in different networks (via a node cost
function) and then aim to find a high-scoring alignment (node mapping between
the networks) with respect to "node conservation", typically the total node
cost function over all aligned nodes. Only after an alignment is constructed,
the existing methods evaluate its quality with respect to an alternative
measure, such as "edge conservation". Thus, we recently aimed to directly
optimize edge conservation while constructing an alignment, which improved
alignment quality. Here, we approach a novel idea of maximizing both node and
edge conservation, and we also approach this idea from a novel perspective, by
aligning optimally edges between networks first in order to improve node cost
function needed to then align well nodes between the networks. In the process,
unlike the existing measures of edge conservation that treat each conserved
edge the same, we favor conserved edges that are topologically similar over
conserved edges that are topologically dissimilar. We show that our novel
method, which we call GRaphlet Edge AlignmenT (GREAT), improves upon
state-of-the-art methods that aim to optimize node conservation only or edge
conservation only.


  In this paper, we focus on susceptible-infected-susceptible dynamics on
metapopulation networks, where nodes represent subpopulations, and where agents
diffuse and interact. Recent studies suggest that heterogeneous network
structure between elements plays an important role in determining the threshold
of infection rate at the onset of epidemics, a fundamental quantity governing
the epidemic dynamics. We consider the general case in which the infection rate
at each node depends on its population size, as shown in recent empirical
observations. We first prove that a sufficient condition for the endemic
threshold (i.e., its upper bound), previously derived based on a mean-field
approximation of network structure, also holds true for arbitrary networks. We
also derive an improved condition showing that networks with the rich-club
property (i.e., high connectivity between nodes with a large number of links)
are more prone to disease spreading. The dependency of infection rate on
population size introduces a considerable difference between this upper bound
and estimates based on mean-field approximations, even when degree-degree
correlations are considered. We verify the theoretical results with numerical
simulations.


  The study of configurational parameters of deformed DNA is a relevant problem
in research of such important biological process as double helix compactization
in cell. The deformations accompanied with local disruptions of the regular
macromolecule structure cause significant bending of the double helix, or
kinks. In this paper an approach for Kratky-Porod model to calculate
persistence length of DNA macromolecule with kinks is developed. The presented
approach considers kinks of arbitrary configuration, including two basic types
of kinks, type 1 - sharp kink caused by unstacking a single base pair step, and
type 2 - intrinsic-induced kink that involves several base pairs. Within
developed approach analytical expressions for persistence length, coil size and
gyration radius of kinky double helix were obtained.


  Identity-by-descent (IBD) is a fundamental concept in genetics with many
applications. In a common definition, two haplotypes are said to contain an IBD
segment if they share a segment that is inherited from a recent shared common
ancestor without intervening recombination. Long IBD segments (> 1cM) can be
efficiently detected by a number of algorithms using high-density SNP array
data from a population sample. However, these approaches detect IBD based on
contiguous segments of identity-by-state, and such segments may exist due to
the conflation of smaller, nearby IBD segments. We quantified this effect using
coalescent simulations, finding that nearly 40% of inferred segments 1-2cM long
are results of conflations of two or more shorter segments, under demographic
scenarios typical for modern humans. This biases the inferred IBD segment
length distribution, and so can affect downstream inferences. We observed this
conflation effect universally across different IBD detection programs and human
demographic histories, and found inference of segments longer than 2cM to be
much more reliable (less than 5% conflation rate). As an example of how this
can negatively affect downstream analyses, we present and analyze a novel
estimator of the de novo mutation rate using IBD segments, and demonstrate that
the biased length distribution of the IBD segments due to conflation can lead
to inflated estimates if the conflation is not modeled. Understanding the
conflation effect in detail will make its correction in future methods more
tractable.


  In post-transcriptional regulation, an mRNA molecule is bound by many
proteins and/or miRNAs to modulate its function. To enable combinatorial gene
regulation, these binding partners of an RNA must communicate with each other,
exhibiting cooperativity. Even in the absence of direct physical interactions
between the binding partners, such cooperativity can be mediated through RNA
secondary structures, since they affect the accessibility of the binding sites.
Here we propose a quantitative measure of this structure-mediated cooperativity
that can be numerically calculated for an arbitrary RNA sequence. Focusing on
an RNA with two binding sites, we derive a characteristic difference of free
energy differences, i.e. \Delta\Delta G, as a measure of the effect of the
occupancy of one binding site on the binding strength of another. We apply this
measure to a large number of human and C. elegans mRNAs, and find that
structure-mediated cooperativity is a generic feature. Interestingly, this
cooperativity not only affects binding sites in close proximity along the
sequence but also configurations in which one binding site is located in the
5'UTR and the other is located in the 3'UTR of the mRNA. Furthermore, we find
that this end-to-end cooperativity is determined by the UTR sequences while the
sequences of the coding regions are irrelevant.


  Recent experimental results on the effect of miRNA on the decay of its target
mRNA have been analyzed against a previously hypothesized single molecule
degradation pathway. According to that hypothesis, the silencing complex
(miRISC) first interacts with its target mRNA and then recruits the protein
complexes associated with NOT1 and PAN3 to trigger deadenylation (and
subsequent degradation) of the target mRNA. Our analysis of the experimental
decay patterns allowed us to refine the structure of the degradation pathways
at the single molecule level. Surprisingly, we found that if the previously
hypothesized network was correct, only about 7% of the target mRNA would be
regulated by the miRNA mechanism, which is inconsistent with the available
knowledge. Based on systematic data analysis, we propose the alternative
hypothesis that NOT1 interacts with miRISC before binding to the target mRNA.
Moreover, we show that when miRISC binds alone to the target mRNA, the mRNA is
degraded more slowly, probably through a deadenylation-independent pathway. The
new biochemical pathway we propose both fits the data and paves the way for new
experimental work to identify new interactions.


  We present a high-throughput microfluidic platform capable of quantitating up
to 384 biomarkers in 4 distinct samples by immunoassay. The microfluidic device
contains 384 unit cells, which can be individually programmed with pairs of
capture and detection antibody. Samples are quantitated in each unit cell by
four independent MITOMI detection areas, allowing four samples to be analyzed
in parallel for a total of 1,536 assays per device. We show that the device can
be pre-assembled and stored for weeks at elevated temperature and we performed
proof-of-concept experiments simultaneously quantitating IL-6, IL-1\b{eta},
TNF-{\alpha}, PSA, and GFP. Finally, we show that the platform can be used to
identify functional antibody combinations by screening 64 antibody combinations
requiring up to 384 unique assays per device.


  ${\bf Background:}$ The early growth dynamics of the West African Ebola virus
epidemic has been qualitatively different for Guinea, Sierra Leone and Liberia.
However, it is important to understand these disparate dynamics as trends of a
single epidemic spread over regions with similar geographic and cultural
aspects, with likely common parameters for transmission rates and the
reproduction number $R_0$. ${\bf Methods:}$ We combine a discrete, stochastic
SEIR model with a two-scale community network model to demonstrate that the
different regional trends may be explained by different community mixing rates.
Heuristically, the effect of different community mixing rates may be understood
as the observation that two individuals infected by the same chain of
transmission are more likely to know one another in a less-mixed community.
Saturation effects occur as the contacts of an infected individual are more
likely to already be exposed by the same chain of transmission. ${\bf
Results:}$ The effects of community mixing, together with the effects of
stochasticity, can explain the qualitative difference in the growth of Ebola
virus cases in each country, and why the probability of large outbreaks may
have recently increased. An increase in the rate of Ebola cases in Guinea in
late August, and a local fitting of the transient dynamics of the Ebola cases
in Liberia, suggests that the epidemic in Liberia has been more severe, and the
epidemic in Guinea is worsening, due to discrete seeding events as the epidemic
spreads into new communities. ${\bf Conclusions:}$ A relatively simple network
model provides insight on the role of local effects such as saturation that
would be difficult to otherwise quantify. Our results predict that exponential
growth of an epidemic is driven by the exposure of new communities,
underscoring the importance of control measures that limit this spread.


  We consider the cell population dynamics with $n$ different phenotypes. Both
the Markovian branching process model (stochastic model) and the ordinary
differential equation (ODE) system model (deterministic model) are presented,
and exploited to investigate the dynamics of the phenotypic proportions. We
will prove that in both models, these proportions will tend to constants
regardless of initial population states ("phenotypic equilibrium") under weak
conditions, which explains the experimental phenomenon in Gupta et al.'s paper.
We also prove that Gupta et al.'s explanation is the ODE model under a special
assumption. As an application, we will give sufficient and necessary conditions
under which the proportion of one phenotype tends to $0$ (die out) or $1$
(dominate). We also extend our results to non-Markovian cases.


  Arterial blood pressure is a key vital sign for the health of the human body.
As such, accurate and reproducible measurement techniques are necessary for
successful diagnosis. Blood pressure measurement is an example of molecular
communication in regulated biological systems. In general, communication in
regulated biological systems is difficult because the act of encoding
information about the state of the system can corrupt the message itself. In
this paper, we propose three strategies to cope with this problem to facilitate
reliable molecular communication links: communicate from the outskirts; build
it in; and leave a small footprint. Our strategies---inspired by communication
in natural biological systems---provide a classification to guide the design of
molecular communication mechanisms in synthetic biological systems. We
illustrate our classification using examples of the first two strategies in
natural systems. We then consider a molecular link within a model based on the
Michaelis-Menten kinetics. In particular, we compute the capacity of the link,
which reveals the potential of communicating using our leave a small footprint
strategy. This provides a way of identifying whether the molecular link can be
improved without affecting the function, and a guide to the design of synthetic
biological systems.


  Epigenome modulation in response to the environment potentially provides a
mechanism for organisms to adapt, both within and between generations. However,
neither the extent to which this occurs, nor the molecular mechanisms involved
are known. Here we investigate DNA methylation variation in Swedish Arabidopsis
thaliana accessions grown at two different temperatures. Environmental effects
on DNA methylation were limited to transposons, where CHH methylation was found
to increase with temperature. Genome-wide association mapping revealed that the
extensive CHH methylation variation was strongly associated with genetic
variants in both cis and trans, including a major trans-association close to
the DNA methyltransferase CMT2. Unlike CHH methylation, CpG gene body
methylation (GBM) on the coding region of genes was not affected by growth
temperature, but was instead strongly correlated with the latitude of origin.
Accessions from colder regions had higher levels of GBM for a significant
fraction of the genome, and this was correlated with elevated transcription
levels for the genes affected. Genome-wide association mapping revealed that
this effect was largely due to trans-acting loci, a significant fraction of
which showed evidence of local adaptation. These findings constitute the first
direct link between DNA methylation and adaptation to the environment, and
provide a basis for further dissecting how environmentally driven and
genetically determined epigenetic variation interact and influence organismal
fitness.


  This paper is concerned with Crump-Mode-Jagers branching processes,
describing spread of an epidemic depending on the proportion of the population
that is vaccinated. Births in the branching process are aborted independently
with a time-dependent probability given by the fraction of the population
vaccinated. Stochastic monotonicity and continuity results for a wide class of
functions (e.g., extinction time and total number of births over all time)
defined on such a branching process are proved using coupling arguments,
leading to optimal vaccination schemes to control corresponding functions
(e.g., duration and final size) of epidemic outbreaks. The theory is
illustrated by applications to the control of the duration of mumps outbreaks
in Bulgaria.


  Homeostasis is a running theme in biology. Often achieved through feedback
regulation strategies, homeostasis allows living cells to control their
internal environment as a means for surviving changing and unfavourable
environments. While many endogenous homeostatic motifs have been studied in
living cells, some other motifs may remain under-explored or even undiscovered.
At the same time, known regulatory motifs have been mostly analyzed at the
deterministic level, and the effect of noise on their regulatory function has
received low attention. Here we lay the foundation for a regulation theory at
the molecular level that explicitly takes into account the noisy nature of
biochemical reactions and provides novel tools for the analysis and design of
robust homeostatic circuits. Using these ideas, we propose a new regulation
motif, which we refer to as {\em antithetic integral feedback, and demonstrate
its effectiveness as a strategy for generically regulating a wide class of
reaction networks. By combining tools from probability and control theory, we
show that the proposed motif preserves the stability of the overall network,
steers the population of any regulated species to a desired set point, and
achieves robust perfect adaptation -- all with low prior knowledge of reaction
rates. Moreover, our proposed regulatory motif can be implemented using a very
small number of molecules and hence has a negligible metabolic load.
Strikingly, the regulatory motif exploits stochastic noise, leading to enhanced
regulation in scenarios where noise-free implementations result in
dysregulation. Finally, we discuss the possible manifestation of the proposed
antithetic integral feedback motif in endogenous biological circuits and its
realization in synthetic circuits.


  Ebola virus disease is a lethal human and primate disease that currently
requires a particular attention from the international health authorities due
to important outbreaks in some Western African countries and isolated cases in
the United Kingdom, the USA and Spain. Regarding the emergency of this
situation, there is a need of development of decision tools, such as
mathematical models, to assist the authorities to focus their efforts in
important factors to eradicate Ebola. In this work, we propose a novel
deterministic spatial-temporal model, called Be-CoDiS (Between-Countries
Disease Spread), to study the evolution of human diseases within and between
countries. The main interesting characteristics of Be-CoDiS are the
consideration of the movement of people between countries, the control measure
effects and the use of time dependent coefficients adapted to each country.
First, we focus on the mathematical formulation of each component of the model
and explain how its parameters and inputs are obtained. Then, in order to
validate our approach, we consider two numerical experiments regarding the
2014-15 Ebola epidemic. The first one studies the ability of the model in
predicting the EVD evolution between countries starting from the index cases in
Guinea in December 2013. The second one consists of forecasting the evolution
of the epidemic by using some recent data. The results obtained with Be-CoDiS
are compared to real data and other models outputs found in the literature.
Finally, a brief parameter sensitivity analysis is done. A free Matlab version
of Be-CoDiS is available at: http://www.mat.ucm.es/momat/software.htm


  Quantitative plant disease resistance is believed to be more durable than
qualitative resistance, since it exerts less selective pressure on the
pathogens. However, the process of progressive pathogen adaptation to
quantitative resistance is poorly understood, which makes it difficult to
predict its durability or to derive principles for its sustainable deployment.
Here, we study the dynamics of pathogen adaptation in response to quantitative
plant resistance affecting pathogen reproduction rate and its carrying
capacity. We developed a stochastic model for the continuous evolution of a
pathogen population within a quantitatively resistant host. We assumed that
pathogen can adapt to a host by the progressive restoration of reproduction
rate or of carrying capacity, or of both. Our model suggests that a combination
of QTLs affecting distinct pathogen traits was more durable if the evolution of
repressed traits was antagonistic. Otherwise, quantitative resistance that
depressed only pathogen reproduction was more durable. In order to decelerate
the progressive pathogen adaptation, QTLs that decrease the pathogen's ability
to extend must be combined with QTLs that decrease the spore production per
lesion or the infection efficiency or that increase the latent period. Our
theoretical framework can help breeders to develop principles for sustainable
deployment of quantitative trait loci..


  In systems biology effective models are widely used due to the complexity of
biological system. They result from a coarse-graining process which employs
specific assumptions. Frequently one does not start with a model taking all
details into account and then performs a coarse-graining process, but rather
one starts right away with the effective equations and often the underlying
assumptions remain hidden or unclear. We exemplify the analysis of an effective
model by analyzing a time delay equation for the JAK/STAT5 signaling pathway
and show how one can avoid wrong conclusions and obtain a deeper understanding
of the biological system . By analyzing the assumptions leading to a
coarse-grained model one might be able to gain new insight into the involved
biological processes. Further, the compliance of the model with experimental
data can be considered as a validation of the assumptions made in the
derivation of the mathematical equations.


  Networks of biochemical reactions, like cellular metabolic networks, are kept
in non-equilibrium steady states by the exchange fluxes connecting them to the
environment. In most cases, feasible flux configurations can be derived from
minimal mass-balance assumptions upon prescribing in- and out-take fluxes. Here
we consider the problem of inferring intracellular flux patterns from
extracellular metabolite levels. Resorting to a thermodynamic out of
equilibrium variational principle to describe the network at steady state, we
show that the switch from fermentative to oxidative phenotypes in cells can be
characterized in terms of the glucose, lactate, oxygen and carbon dioxide
concentrations. Results obtained for an exactly solvable toy model are fully
recovered for a large scale reconstruction of human catabolism. Finally we
argue that, in spite of the many approximations involved in the theory,
available data for several human cell types are well described by the predicted
phenotypic map of the problem.


  We consider a stochastic model for an evolving population. We show that in
the presence of genotype extinctions the population dies out for a low mutation
probability but may survive for a high mutation probability. This turns upside
down the widely held belief that above a certain mutation threshold a
population cannot survive.


  Btrim is a fast and lightweight software to trim adapters and low quality
regions in reads from ultra high-throughput next-generation sequencing
machines. It also can reliably identify barcodes and assign the reads to the
original samples. Based on a modified Myers's bit-vector dynamic programming
algorithm, Btrim can handle indels in adapters and barcodes. It removes low
quality regions and trims off adapters at both or either end of the reads. A
typical trimming of 30M reads with two sets of adapter pairs can be done in
about a minute with a small memory footprint. Btrim is a versatile stand-alone
tool that can be used as the first step in virtually all next-generation
sequence analysis pipelines. The program is available at
\url{http://graphics.med.yale.edu/trim/}.


  Similar to intelligent multicellular neural networks controlling human
brains, even single cells surprisingly are able to make intelligent decisions
to classify several external stimuli or to associate them. This happens because
of the fact that gene regulatory networks can perform as perceptrons, simple
intelligent schemes known from studies on Artificial Intelligence. We study the
role of genetic noise in intelligent decision making at the genetic level and
show that noise can play a constructive role helping cells to make a proper
decision. We show this using the example of a simple genetic classifier able to
classify two external stimuli.


  We present a method to obtain numerically accurate values of configurational
free energies of semiflexible macromolecular systems, based on the technique of
thermodynamic integration combined with normal-mode analysis of a reference
system subject to harmonic constraints. Compared with previous free-energy
calculations that depend on a reference state, our approach introduces two
innovations, namely the use of internal coordinates to constrain the reference
states and the ability to freely select these reference states. As a
consequence, it is possible to explore systems that undergo substantially
larger fluctuations than those considered in previous calculations, including
semiflexible biopolymers having arbitrary ratios of contour length L to
persistence length P. To validate the method, high accuracy is demonstrated for
free energies of prime DNA knots with L/P=20 and L/P=40, corresponding to DNA
lengths of 3000 and 6000 base pairs, respectively. We then apply the method to
study the free-energy landscape for a model of a synaptic nucleoprotein complex
containing a pair of looped domains, revealing a bifurcation in the location of
optimal synapse (crossover) sites. This transition is relevant to target-site
selection by DNA-binding proteins that occupy multiple DNA sites separated by
large linear distances along the genome, a problem that arises naturally in
gene regulation, DNA recombination, and the action of type-II topoisomerases.


  DNA methylation is an intensely studied epigenetic mark implicated in many
biological processes of direct clinical relevance. While sequencing based
technologies are increasingly allowing high resolution measurements of DNA
methylation, statistical modelling of such data is still challenging. In
particular, statistical identification of differentially methylated regions
across different conditions poses unresolved challenges in accounting for
spatial correlations within the statistical testing procedure.
  We propose a non-parametric, kernel-based method, M3D, to detect higher-order
changes in methylation profiles, such as shape, across pre-defined regions. The
test statistic explicitly accounts for differences in coverage levels between
samples, thus handling in a principled way a major confounder in the analysis
of methylation data. Empirical tests on real and simulated data sets show an
increased power compared to established methods, as well as considerable
robustness with respect to coverage and replication levels.
  Code is implemented in R, and is freely available via Bioconductor package
M3D.


  Deoxycitidine in solution exists as two tautomers one of which is an
uncanonical imino one. The latter can dominate with such derivatives as
5-methyl, 5-hydroxymethyl- and 5-formylcytosine. The imino tautomer potentially
is able to form a hoosteen GC base pair. To detect such pair, it is suggested
to use 1H15N NMR. Formation of GC-Hoogsteen base pair with imino tautomer of
cytosine can be a reason for epigenetic effects of 5-methyl- and
5-hydroxymethylcytosine.


  Protein fibril accumulation at interfaces is an important step in many
physiological processes and neurodegenerative diseases as well as in designing
materials. Here we show, using $\beta$-lactoglobulin fibrils as a model, that
semiflexible fibrils exposed to a surface do not possess the Gaussian
distribution of curvatures characteristic for wormlike chains, but instead
exhibit a spontaneous curvature, which can even lead to ring-like
conformations. The long-lived presence of such rings is confirmed by atomic
force microscopy, cryogenic scanning electron microscopy and passive probe
particle tracking at air- and oil-water interfaces. We reason that this
spontaneous curvature is governed by structural characteristics on the
molecular level and is to be expected when a chiral and polar fibril is placed
in an inhomogeneous environment such as an interface. By testing
$\beta$-lactoglobulin fibrils with varying average thicknesses, we conclude
that fibril thickness plays a determining role in the propensity to form rings.


  It is well-known that gene activation/deactivation dynamics may be a major
source of randomness in genetic networks, also in the case of large
concentrations of the transcription factors. In this work, we investigate the
effect of realistic extrinsic noises acting on gene deactivation in a common
network motif - the positive feedback of a transcription factor on its own
synthesis - under a variety of settings, i.e., distinct cellular types,
distribution of proteins and properties of the external perturbations. At
variance with standard models where the perturbations are Gaussian unbounded,
we focus here on bounded extrinsic noises to better mimic biological reality.
Our results suggest that the gene switching velocity is a key parameter to
modulate the response of the network. Simulations suggest that, if the gene
switching is fast and many proteins are present, an irreversible noise-induced
first order transition is observed as a function of the noise intensity. If
noise intensity is further increased a second order transition is also
observed. When gene switching is fast, a similar scenario is observed even when
few proteins are present, provided that larger cells are considered, which is
mainly influenced by noise autocorrelation. On the contrary, if the gene
switching is slow, no fist order transitions are observed. In the concluding
remarks possible implications of the irreversible transitions for cellular
differentiation are briefly discussed.


  Reconciliation methods aim at recovering macro evolutionary events and at
localizing them in the species history, by observing discrepancies between gene
family trees and species trees. In this article we introduce an Integer Linear
Programming (ILP) approach for the NP-hard problem of computing a most
parsimonious time-consistent reconciliation of a gene tree with a species tree
when dating information on speciations is not available. The ILP formulation,
which builds upon the DTL model, returns a most parsimonious reconciliation
ranging over all possible datings of the nodes of the species tree. By studying
its performance on plausible simulated data we conclude that the ILP approach
is significantly faster than a brute force search through the space of all
possible species tree datings. Although the ILP formulation is currently
limited to small trees, we believe that it is an important proof-of-concept
which opens the door to the possibility of developing an exact, parsimony based
approach to dating species trees. The software (ILPEACE) is freely available
for download.


  The formation of a collectively moving group benefits individuals within a
population in a variety of ways such as ultra-sensitivity to perturbation,
collective modes of feeding, and protection from environmental stress. While
some collective groups use a single organizing principle, others can
dynamically shift the behavior of the group by modifying the interaction rules
at the individual level. The surface-dwelling bacterium Myxococcus xanthus
forms dynamic collective groups both to feed on prey and to aggregate during
times of starvation. The latter behavior, termed fruiting-body formation,
involves a complex, coordinated series of density changes that ultimately lead
to three-dimensional aggregates comprising hundreds of thousands of cells and
spores. This multi-step developmental process most likely involves several
different single-celled behaviors as the population condenses from a loose,
two-dimensional sheet to a three-dimensional mound. Here, we use
high-resolution microscopy and computer vision software to spatiotemporally
track the motion of thousands of individuals during the initial stages of
fruiting body formation. We find that a combination of cell-contact-mediated
alignment and internal timing mechanisms drive a phase transition from
exploratory flocking, in which cell groups move rapidly and coherently over
long distances, to a reversal-mediated localization into streams, which act as
slow-spreading, quasi-one-dimensional nematic fluids. These observations lead
us to an active liquid crystal description of the myxobacterial development
cycle.


  Taylors Law (TL) describes the scaling relationship between the mean and
variance of populations as a power-law. TL is widely observed in ecological
systems across space and time with exponents varying largely between 1 and 2.
Many ecological explanations have been proposed for TL but it is also commonly
observed outside ecology. We propose that TL arises from the constraining
influence of two primary variables: the number of individuals and the number of
censuses or sites. We show that most possible configurations of individuals
among censuses or sites produce the power-law form of TL with exponents between
1 and 2. This feasible set approach suggests that TL is a statistical pattern
driven by two constraints, providing an a priori explanation for this
ubiquitous pattern. However, the exact form of any specific mean-variance
relationship cannot be predicted in this way, i.e., this approach does a poor
job of predicting variation in the exponent, suggesting that TL may still
contain ecological information.


  Seasonal influenza epidemics cause consistent, considerable, widespread loss
annually in terms of economic burden, morbidity, and mortality. With access to
accurate and reliable forecasts of a current or upcoming influenza epidemic's
behavior, policy makers can design and implement more effective
countermeasures. We developed a framework for in-season forecasts of epidemics
using a semiparametric Empirical Bayes framework, and applied it to predict the
weekly percentage of outpatient doctors visits for influenza-like illness, as
well as the season onset, duration, peak time, and peak height, with and
without additional data from Google Flu Trends, as part of the CDC's 2013--2014
"Predict the Influenza Season Challenge". Previous work on epidemic modeling
has focused on developing mechanistic models of disease behavior and applying
time series tools to explain historical data. However, these models may not
accurately capture the range of possible behaviors that we may see in the
future. Our approach instead produces possibilities for the epidemic curve of
the season of interest using modified versions of data from previous seasons,
allowing for reasonable variations in the timing, pace, and intensity of the
seasonal epidemics, as well as noise in observations. Since the framework does
not make strict domain-specific assumptions, it can easily be applied to other
diseases as well. Another important advantage of this method is that it
produces a complete posterior distribution for any desired forecasting target,
rather than mere point predictions. We report prospective
influenza-like-illness forecasts that were made for the 2013--2014 U.S.
influenza season, and compare the framework's cross-validated prediction error
on historical data to that of a variety of simpler baseline predictors.


  Advances of modern sensing and sequencing technologies generate a deluge of
high dimensional space-temporal physiological and next-generation sequencing
(NGS) data. Physiological traits are observed either as continuous random
functions, or on a dense grid and referred to as function-valued traits. Both
physiological and NGS data are highly correlated data with their inherent
order, spacing, and functional nature which are ignored by traditional
summary-based univariate and multivariate regression methods designed for
quantitative genetic analysis of scalar trait and common variants. To capture
morphological and dynamic features of the data and utilize their dependent
structure, we propose a functional linear model (FLM) in which a trait curve is
modeled as a response function, the genetic variation in a genomic region or
gene is modeled as a functional predictor, and the genetic effects are modeled
as a function of both time and genomic position (FLMF) for genetic analysis of
function-valued trait with both GWAS and NGS data. By extensive simulations, we
demonstrate that the FLMF has the correct type 1 error rates and much higher
power to detect association than the existing methods. The FLMF is applied to
sleep data from Starr County health studies where oxygen saturation were
measured in 22,670 seconds on average for 833 individuals. We found 65 genes
that were significantly associated with oxygen saturation functional trait with
P-values ranging from 2.40E-06 to 2.53E-21. The results clearly demonstrate
that the FLMF substantially outperforms the traditional genetic models with
scalar trait.


  Discovering patterns in networks of protein-protein interactions (PPIs) is a
central problem in systems biology. Alignments between these networks aid
functional understanding as they uncover important information, such as
evolutionary conserved pathways, protein complexes and functional orthologs.
The objective of a multiple network alignment is to create clusters of nodes
that are evolutionarily conserved and functionally consistent across all
networks. Unfortunately, the alignment methods proposed thus far do not fully
meet this objective, as they are guided by pairwise scores that do not utilize
the entire functional and topological information across all networks. To
overcome this weakness, we propose FUSE, a multiple network aligner that
utilizes all functional and topological information in all PPI networks. It
works in two steps. First, it computes novel similarity scores of proteins
across the PPI networks by fusing from all aligned networks both the protein
wiring patterns and their sequence similarities. It does this by using
Non-negative Matrix Tri-Factorization (NMTF). When we apply NMTF on the five
largest and most complete PPI networks from BioGRID, we show that NMTF finds a
larger number of protein pairs across the PPI networks that are functionally
conserved than can be found by using protein sequence similarities alone. This
demonstrates complementarity of protein sequence and their wiring patterns in
the PPI networks. In the second step, FUSE uses a novel maximum weight
k-partite matching approximation algorithm to find an alignment between
multiple networks. We compare FUSE with the state of the art multiple network
aligners and show that it produces the largest number of functionally
consistent clusters that cover all aligned PPI networks. Also, FUSE is more
computationally efficient than other multiple network aligners.


  Bistability is ubiquitous in biological systems. For example, bistability is
found in many reaction networks that involve the control and execution of
important biological functions, such as signalling processes. Positive feedback
loops, composed of species and reactions, are necessary for bistability, and
generally for multi-stationarity, to occur. These loops are therefore often
used to illustrate and pinpoint the parts of a multi-stationary network that
are relevant (`responsible') for the observed multi-stationarity. However
positive feedback loops are generally abundant in reaction networks but not all
of them are important for subsequent interpretation of the network's dynamics.
  We present an automated procedure to determine the relevant positive feedback
loops of a multi-stationary reaction network. The procedure only reports the
loops that are relevant for multi-stationarity (that is, when broken
multi-stationarity disappears) and not all positive feedback loops of the
network. We show that the relevant positive feedback loops must be understood
in the context of the network (one loop might be relevant for one network, but
cannot create multi-stationarity in another). Finally, we demonstrate the
procedure by applying it to several examples of signaling processes, including
a ubiquitination and an apoptosis network, and to models extracted from the
Biomodels database.
  We have developed and implemented an automated procedure to find relevant
positive feedback loops in reaction networks. The results of the procedure are
useful for interpretation and summary of the network's dynamics.