  Analogous to genomic sequence alignment, biological network alignment (NA)
aims to find regions of similarities between molecular networks (rather than
sequences) of different species. NA can be either local (LNA) or global (GNA).
LNA aims to identify highly conserved common subnetworks, which are typically
small, while GNA aims to identify large common subnetworks, which are typically
suboptimally conserved. We recently showed that LNA and GNA yield complementary
results: LNA has high functional but low topological alignment quality, while
GNA has high topological but low functional alignment quality. Thus, we propose
IGLOO, a new approach that integrates GNA and LNA in hope to reconcile the two.
We evaluate IGLOO against state-of-the-art LNA (NetworkBLAST, NetAligner,
AlignNemo, and AlignMCL) and GNA (GHOST, NETAL, GEDEVO, MAGNA++, WAVE, and
L-GRAAL) methods. We show that IGLOO allows for a trade-off between topological
and functional alignment quality better than the existing LNA and GNA methods
considered in our study.


  Aedes Aegypti is the vector of several deadly diseases, including Zika.
Effective and sustainable vector control measures must be deployed to keep A.
aegypti numbers under control. The distribution of A. Aegypti is subject to
spatial and climatic constraints. Using agent-based modeling, we model the
population dynamics of A. aegypti subjected to the spatial and climatic
constraints of a neighborhood in the Key West. Satellite imagery was used to
identify vegetation, houses (CO2 zones) both critical to the mosquito
lifecycle. The model replicates the seasonal fluctuation of adult population
sampled through field studies and approximates the population at a high of 986
(95% CI: [979, 993]) females and 1031 (95% CI: [1024, 1039]) males in the fall
and a low of 316 (95% CI: [313, 319]) females and 333 (95% CI: [330, 336])
males during the winter. We then simulate two biological vector control
strategies: 1) Wolbachia infection and 2) Release of Insects carrying a
Dominant Lethal gene (RIDL). Our results support the probability of sustained
Wolbachia infection within the population for two years after the year of
release. egies, our approach provides a realistic simulation environment
consisting of male and female Aedes aegypti, breeding spots, vegetation and CO2
sources.


  When choosing the timing of cross-sectional network snapshots in longitudinal
social network studies, the effect on the precision of parameter estimates
generally plays a minor role. Often the timing is opportunistic or determined
by a variety of considerations such as organizational constraints, funding, and
availability of study participants. Theory to guide the timing of network
snapshots is also missing. We use a statistical framework to relate the timing
to the precision of the parameter estimates, specifically, the sum of the
relative widths of their confidence intervals. We illustrate this
computationally using the STERGM suite of the statnet package to estimate the
parameters of the network dynamics. Analytically, we derive simple
approximations for the optimal timing when the parameters correspond to the
rates for different network events. We find that the optimal time depends the
most on the network processes with short expected durations and few expected
events such as the dissolution of ties. We apply our approximations to a simple
example of a dynamic network with formation and dissolution of iid ties.


  A major aim of evolutionary biology is to explain the respective roles of
adaptive versus non-adaptive changes in the evolution of complexity. While
selection is certainly responsible for the spread and maintenance of complex
phenotypes, this does not automatically imply that strong selection enhances
the chance for the emergence of novel traits, that is, the origination of
complexity. Population size is one parameter that alters the relative
importance of adaptive and non-adaptive processes: as population size
decreases, selection weakens and genetic drift grows in importance. Because of
this relationship, many theories invoke a role for population size in the
evolution of complexity. Such theories are difficult to test empirically
because of the time required for the evolution of complexity in biological
populations. Here, we used digital experimental evolution to test whether large
or small asexual populations tend to evolve greater complexity. We find that
both small and large---but not intermediate-sized---populations are favored to
evolve larger genomes, which provides the opportunity for subsequent increases
in phenotypic complexity. However, small and large populations followed
different evolutionary paths towards these novel traits. Small populations
evolved larger genomes by fixing slightly deleterious insertions, while large
populations fixed rare beneficial insertions that increased genome size. These
results demonstrate that genetic drift can lead to the evolution of complexity
in small populations and that purifying selection is not powerful enough to
prevent the evolution of complexity in large populations.


  Cretaceous-Cenozoic history of angiosperms led to the a certain character of
the distribution of taxa of different levels (the number of species and genera
in families, species/genera ratio in families, the number of species in the
genera). In most cases, these distributions are satisfactorily described by a
power law (Pareto distribution). In logarithmic coordinates power function is a
straight line. Empirical curves repeat this line is good enough, but on the
right side of the graph (at low volumes taxa), there is a marked deviation of
theoretical from the empirical curves. This suggests that the small volumes of
taxa should be greater for full compliance with the theoretical curves.
Modeling ratios among genera and species in families showed that only in the
case of dynamic extinction factor observed satisfactory agreement between
observed and calculated the number of species in a wide range of iterations.
This suggested that there was a differential extinction of species during the
evolution of angiosperms. This implies that the rate of extinction had to be
minimal in genera with a large number of species. On the contrary, extinction
rates may be increased by orders of magnitude with a decrease in the number of
species. As a result, large genera became getting bigger and small genera
become less. The species frequencies of distribution in the genera varied
according to a power law. The initial divergence of taxa volume, which led to
their further division into large and small, could be caused by the emergence
and expansion of herbs with their functional and adaptive capabilities.


  When predicting the fate and consequences of recurring deleterious mutations
in self-fertilising populations most models developed make the assumption that
populations have discrete non-overlapping generations. This makes them
biologically irrelevant when considering perennial species with over-lapping
generations and where mating occurs independently of the age group. Previous
models studying the effect of perennial life-histories on the genetic
properties of populations in the presence of self-fertilisation have done so
considering age-dependent selection and have found that, contrary to empirical
observations, perennial populations should exhibit lower levels of inbreeding
depression. Here we propose a simple deterministic model in continuous time
with selection at different fitness traits and feedback between population
fitness and size. We find that a perennial life-history can result in high
levels of inbreeding depression in spite of inbreeding, due to higher
frequencies of heterozygous individuals at the adult stage. We also propose
that there may be demographic advantages for self-fertilisation that are
independent of reproductive success.


  The most common assumption in evolutionary game theory is that players should
adopt a strategy that warrants the highest payoff. However, recent studies
indicate that the spatial selection for cooperation is enhanced if an
appropriate fraction of the population chooses the most common rather than the
most profitable strategy within the interaction range. Such conformity might be
due to herding instincts or crowd behavior in humans and social animals. In a
heterogeneous population where individuals differ in their degree, collective
influence, or other traits, an unanswered question remains who should conform.
Selecting conformists randomly is the simplest choice, but it is neither a
realistic nor the optimal one. We show that, regardless of the source of
heterogeneity and game parametrization, socially the most favorable outcomes
emerge if the masses conform. On the other hand, forcing leaders to conform
significantly hinders the constructive interplay between heterogeneity and
coordination, leading to evolutionary outcomes that are worse still than if
conformists were chosen randomly. We conclude that leaders must be able to
create a following for network reciprocity to be optimally augmented by
conformity. In the opposite case, when leaders are castrated and made to
follow, the failure of coordination impairs the evolution of cooperation.


  We study numerically a coagulation-fragmentation model derived by Niwa and
further elaborated by Degond et al., where a unique equilibrium distribution of
group sizes is shown to exist in both cases of continuous and discrete group
size distributions. We provide a numerical investigation of these equilibria
using three different methods to approximate the equilibrium: a recursive
algorithm based on the work of Ma et. al., a Newton method and the resolution
of the time-dependent problem. All three schemes are validated by showing that
they approximate the predicted small and large size asymptotic behaviour of the
equilibrium accurately. The recursive algorithm is used to investigate the
transition from discrete to continuous size distributions and the time
evolution scheme is exploited to show uniform convergence to equilibrium in
time and to determine convergence rates.


  In this paper results are obtained concerning the number of positive
stationary solutions in simple models of the Calvin cycle of photosynthesis and
the stability of these solutions. It is proved that there are open sets of
parameters in a model of Zhu et. al. for which there exist two positive
stationary solutions. There are never more than two isolated positive
stationary solutions but under certain explicit special conditions on the
parameters there is a whole continuum of positive stationary solutions. It is
also shown that in the set of parameter values for which two isolated positive
stationary solutions exist there is an open subset where one of the solutions
is asymptotically stable and the other is unstable. In related models , for
which it was known that more that one positive stationary solution exists, it
is proved that there are parameter values for which one of these solutions is
asymptotically stable and the other unstable. A key technical aspect of the
proofs is to exploit the fact that there is a bifurcation where the centre
manifold is one-dimensional.


  Mathematical models are important tools to study the excluded volume effects
on reaction-diffusion systems, which are known to play an important role inside
living cells. Detailed microscopic simulations with off-lattice Brownian
dynamics become computationally expensive in crowded environments. In this
paper we therefore investigate to which extent on-lattice approximations, so
called Cellular Automata models, can be used to simulate reactions and
diffusion in the presence of crowding molecules. We show that the diffusion is
most severely slowed down in the off-lattice model, since randomly distributed
obstacles effectively exclude more volume than those ordered on an artificial
grid. Crowded reaction rates can be both increased and decreased by the grid
structure and it proves important to model the molecules with realistic sizes
when excluded volume is taken into account. The grid artifacts increase with
increasing crowder density and we conclude that the computationally more
efficient on-lattice simulations are accurate approximations only for low
crowder densities.


  We present an exact, closed expression for the expected neutral Site
Frequency Spectrum for two neutral sites, 2-SFS, without recombination. This
spectrum is the immediate extension of the well known single site $\theta/f$
neutral SFS. Similar formulae are also provided for the case of the expected
SFS of sites that are linked to a focal neutral mutation of known frequency.
Formulae for finite samples are obtained by coalescent methods and remarkably
simple expressions are derived for the SFS of a large population, which are
also solutions of the multi-allelic Kolmogorov equations. Besides the general
interest of these new spectra, they relate to interesting biological cases such
as structural variants and introgressions. As an example, we present the
expected neutral frequency spectrum of regions with a chromosomal inversion.


  In this work, we develop first near-complete 3D models for NTD-hXPB - the
N-terminal protein domain of the human transcription factor XPB. The results
are very significant as NTD-hXPB plays a critical role in the synthesis of
proteins (specifically transcription) and DNA damage repair (specifically
nucleotide excision repair). NTD-hXPB is directly implicated in rare diseases
XP-B, XP-CS, and TTD2, whose symptoms include neurodegenerative disorders,
premature aging, and decreased fertility. NTDhXPB is also linked to anti-cancer
therapies. As a bi-product we derived 3D models of NTD-mXPB - homologue of
NTD-hXPB in mycobacterium tuberculosis aka MTB (causative agent of most cases
of tuberculosis, which surpassed HIV as #1 infectious disease killer in 2014).
These could be potential target for TB therapeutics. Our ab-initio modelling
protocol takes advantage of recent powerful advances (prediction of contact
residues) in protein structure modeling. Our protocol also includes human in
the loop, inspired by the prevailing theory that computational abilities of
human minds can be powerfully harnessed in engineering/problem-solving. Using
the developed models in this paper, we are able to propose significant insights
into (a) the role of NTD-hXPB in the process of transcription and DNA damage
repair (specifically NER), (b) the critical interactions of NTD-hXPB with
another co-labourer protein p52, (c) diseases associated with NTD-hXPB, and (d)
alterations in functionalities between NTD-hXPB (in human) and NTD-mXPB (in TB
pathogen).


  Indirect reciprocity is one of the major mechanisms of the evolution of
cooperation. Because constant monitoring and accurate evaluation in moral
assessments tend to be costly, indirect reciprocity can be exploited by cost
evaders. A recent study crucially showed that a cooperative state achieved by
indirect reciprocators is easily destabilized by cost evaders in the case with
no supportive mechanism. Here, we present a simple and widely applicable
solution that considers pre-assessment of cost evaders. In the pre-assessment,
those who fail to pay for costly assessment systems are assigned a nasty image
that leads to being rejected by discriminators. We demonstrate that considering
the pre-assessment can crucially stabilize reciprocal cooperation for a broad
range of indirect reciprocity models. In particular for the most leading social
norms we analyse the conditions under which a prosocial state becomes locally
stable.


  We present a model for evolving agents using both genetic and cultural
inheritance mechanisms. Within each agent our model maintains two distinct
information stores we call the genome and the memome. Processes of adaptation
are modeled as evolutionary processes at each level of adaptation
(phylogenetic, ontogenetic, sociogenetic). We review relevant competing models
and we show how our model improves on previous attempts to model genetic and
cultural evolutionary processes. In particular we argue our model can achieve
divergent gene-culture co-evolution.


  Sexual reproductive behavior has a necessary social coordination component as
willing and capable partners must both be in the right place at the right time.
It has recently been demonstrated that many social organizations that support
sexual reproduction can evolve in the absence of social coordination between
agents (e.g. herding, assortative mating, and natal philopatry). In this paper
we explore these results by including social transfer mechanisms to our agents
and contrasting their reproductive behavior with a control group without social
transfer mechanisms. We conclude that similar behaviors emerge in our social
learning agents as those that emerged in the non-social learning agents. Social
learners were more inclined towards natal philopatry. Social learners also
evolved a culture of eusociality including reproductive division of labor.


  Divergent cumulative cultural evolution occurs when the cultural evolutionary
trajectory diverges from the biological evolutionary trajectory. We consider
the conditions under which divergent cumulative cultural evolution can occur.
We hypothesize that two conditions are necessary. First that genetic and
cultural information are stored separately in the agent. Second cultural
information must be transferred horizontally between agents of different
generations. We implement a model with these properties and show evidence of
divergent cultural evolution under both cooperative and competitive selection
pressures.


  Protein secondary structure prediction is an important problem in
bioinformatics. Inspired by the recent successes of deep neural networks, in
this paper, we propose an end-to-end deep network that predicts protein
secondary structures from integrated local and global contextual features. Our
deep architecture leverages convolutional neural networks with different kernel
sizes to extract multiscale local contextual features. In addition, considering
long-range dependencies existing in amino acid sequences, we set up a
bidirectional neural network consisting of gated recurrent unit to capture
global contextual features. Furthermore, multi-task learning is utilized to
predict secondary structure labels and amino-acid solvent accessibility
simultaneously. Our proposed deep network demonstrates its effectiveness by
achieving state-of-the-art performance, i.e., 69.7% Q8 accuracy on the public
benchmark CB513, 76.9% Q8 accuracy on CASP10 and 73.1% Q8 accuracy on CASP11.
Our model and results are publicly available.


  One of the notable fields in studying the genetics of cancer is disease gene
identification which affects disease treatment and drug discovery. Many
researches have been done in this field. Genome-wide association studies (GWAS)
are one of them that focus on the identification of diseases-susceptible loci
on chromosomes. Recently, computational approaches, known as gene
prioritization methods have been used to identify candidate disease genes. Gene
prioritization methods integrate several data sources to discover and
prioritize the most probable candidate disease genes. In this paper, we propose
a prioritization method, called HybridRanker which is a network-based technique
and it also uses experimental data to identify candidate cancer genes. We apply
our proposed method on colorectal cancer data. It is notable to say that in
HybridRanker, for considering both local and global network information of a
protein-protein interaction network, different algorithms such as
shortest-path, random walk with restart and network propagation are exploited.
By using these algorithms, initial scores are given to genes within the
network. Furthermore, by looking through diseases with similar symptoms and
also comorbid diseases and by extracting their causing genes, the gene scores
are recalculated. We also use gene-phenotype relations for an additional
scoring of the candidate genes. Our method is validated and compared with other
prioritization methods in leave one-out cross-validation and the comparison
results show the better performance of the HybridRanker.


  The identification of disease-associated genes has recently gathered much
attention for uncovering disease complex mechanisms that could lead to new
insights into the treatment of diseases. For exploring disease-susceptible
genes, not only experimental approaches such as genome-wide association studies
(GWAS) have been used, but also computational methods. Since experimental
approaches are both time-consuming and expensive, numerous studies have
utilized computational techniques to explore disease genes. These methods use
various biological data sources and known disease genes to prioritize disease
candidate genes. In this paper, we propose a gene prioritization method
(NRSSPrioritize), which benefits from both local and global measures of a
protein-protein interaction (PPI) network and also from disease similarity
knowledge to suggest candidate genes for colorectal cancer (CRC)
susceptibility. Network Propagation, Random Walk with Restart, and Shortest
Paths are three network analysis tools that are applied to a PPI network for
the purpose of scoring candidate genes. Also, by looking through diseases with
similar symptoms to CRC and obtaining their causing genes, candidate genes are
scored in a different way. Finally, to integrate these four different scoring
schemes, Technique for Order Preference by Similarity to Ideal Solution
(TOPSIS) and Analytic Network Process (ANP) methods are applied to obtain
appropriate weights for the above four quantified measures and the weighted
summation of these measures are used to calculate the final score of each
candidate gene. NRSSPrioritize was validated by cross-validation analysis and
its results were compared with other prioritization tools, which gave the best
performance when using our proposed method.


  Renormalization Group (RG) theory provides the theoretical framework to
define Effective Theories (ETs), i.e. systematic low-resolution approximations
of arbitrary microscopic models. Markov State Models (MSMs) are shown to be
rigorous ETs for Molecular Dynamics (MD). Based on this fact, we use Real Space
RG to vary the resolution of a MSM and define an algorithm for clustering
microstates into macrostates. The result is a lower dimensional stochastic
model which, by construction, provides the optimal coarse-grained Markovian
representation of the system's relaxation kinetics. To illustrate and validate
our theory, we analyze a number of test systems of increasing complexity,
ranging from synthetic toy models to two realistic applications, built form
all-atom MD simulations. The computational cost of computing the
low-dimensional model remains affordable on a desktop computer even for
thousands of microstates.


  The flux of visitors through popular places undoubtedly influences viral
spreading -- from H1N1 and Zika viruses spreading through physical spaces such
as airports, to rumors and ideas spreading though online spaces such as
chatrooms and social media. However there is a lack of understanding of the
types of viral dynamics that can result. Here we present a minimal dynamical
model which focuses on the time-dependent interplay between the {\em mobility
through} and the {\em occupancy of} such spaces. Our generic model permits
analytic analysis while producing a rich diversity of infection profiles in
terms of their shapes, durations, and intensities. The general features of
these theoretical profiles compare well to real-world data of recent social
contagion phenomena.


  We present a new methodology for efficient and high-quality patterning of
biological reagents for surface-based biological assays. The method relies on
hydrodynamically confined nanoliter volumes of reagents to interact with the
substrate at the micrometer-length scale. We study the interplay between
diffusion, advection, and surface chemistry and present the design of a
noncontact scanning microfluidic device to efficiently present reagents on
surfaces. By leveraging convective flows, recirculation, and mixing of a
processing liquid, this device overcomes limitations of existing biopatterning
approaches, such as passive diffusion of analytes, uncontrolled wetting, and
drying artifacts. We demonstrate the deposition of analytes, showing a 2- to
5-fold increase in deposition rate together with a 10-fold reduction in analyte
consumption while ensuring less than 6% variation in pattern homogeneity on a
standard biological substrate. In addition, we demonstrate the recirculation of
a processing liquid using a microfluidic probe (MFP) in the context of a
surface assay for (i) probing 12 independent areas with a single microliter of
processing liquid and (ii) processing a 2 mm2 surface to create 170 antibody
spots of 50 x 100 {\mu}m2 area using 1.6 {\mu}L of liquid. We observe high
pattern quality, conservative usage of reagents, micrometer precision of
localization and convection-enhanced fast deposition. Such a device and method
may facilitate quantitative biological assays and spur the development of the
next generation of protein microarrays.


  In this paper, a novel West Nile Virus model looking upon the infected birds
as monitoring threshold, for the mosquitoes and birds with impulsive state
feedback control is considered. We obtain sufficient conditions of the global
asymptotical stability of the system without impulsive state feedback control
via comprehensively qualitative analysis. By using the Poincar\'e map, we
obtain that the system with impulsive state feedback control has a positive
periodic solution of order-1 or order-2 which is asymptotical stability due to
the analogue of Poincar\'e criterion, theory of differential inequalities,
differential equation geometry and so on. What's more, sufficient conditions
for existence and stability of the order one periodic solution are given by the
existence and uniqueness of the limit. Our results show that the control
measure is effective and feasible by means of numerical simulations.


  In this contribution I critically revise the alchemical reversible approach
in the context of the statistical mechanics theory of non covalent bonding in
drug receptor systems. I show that most of the pitfalls and entanglements for
the binding free energies evaluation in computer simulations are rooted in the
equilibrium assumption that is implicit in the reversible method. These
critical issues can be resolved by using a non-equilibrium variant of the
alchemical method in molecular dynamics simulations, relying on the production
of many independent trajectories with a continuous dynamical evolution of an
externally driven alchemical coordinate, completing the decoupling of the
ligand in a matter of few tens of picoseconds rather than nanoseconds. The
absolute binding free energy can be recovered from the annihilation work
distributions by applying an unbiased unidirectional free energy estimate, on
the assumption that any observed work distribution is given by a mixture of
normal distributions, whose components are identical in either direction of the
non-equilibrium process, with weights regulated by the Crooks theorem. I
finally show that the inherent reliability and accuracy of the unidirectional
estimate of the decoupling free energies, based on the production of few
hundreds of non-equilibrium independent sub-nanoseconds unrestrained alchemical
annihilation processes, is a direct consequence of the funnel-like shape of the
free energy surface in molecular recognition. An application of the technique
on a real drug-receptor system is presented in the companion paper.


  Kinesin-8 motor proteins destabilize microtubules. Their absence during cell
division is associated with disorganized mitotic chromosome movements and
chromosome loss. Despite recent work studying effects of kinesin 8s on
microtubule dynamics, it remains unclear whether the kinesin-8 mitotic
phenotypes are consequences of their effect on microtubule dynamics, their
well-established motor activity, or additional unknown functions. To better
understand the role of kinesin-8 proteins in mitosis, we have studied the
effects of deletion of the fission-yeast kinesin-8 proteins Klp5 and Klp6 on
chromosome movements and spindle length dynamics. Aberrant microtubule-driven
kinetochore pushing movements and tripolar mitotic spindles occurred in cells
lacking Klp5 but not Klp6. Kinesin-8 deletion strains showed large fluctuations
in metaphase spindle length, suggesting a disruption of spindle length
stabilization. Comparison of our results from light microscopy with a
mathematical model suggests that kinesin-8 induced effects on microtubule
dynamics, kinetochore attachment stability, and sliding force in the spindle
can explain the aberrant chromosome movements and spindle length fluctuations
seen.


  Because of the potential link between -1 programmed ribosomal frameshifting
and response of a pseudoknot (PK) RNA to force, a number of single molecule
pulling experiments have been performed on PKs to decipher the mechanism of
programmed ribosomal frameshifting. Motivated in part by these experiments, we
performed simulations using a coarse-grained model of RNA to describe the
response of a PK over a range of mechanical forces ($f$s) and monovalent salt
concentrations ($C$s). The coarse-grained simulations quantitatively reproduce
the multistep thermal melting observed in experiments, thus validating our
model. The free energy changes obtained in simulations are in excellent
agreement with experiments. By varying $f$ and $C$, we calculated the phase
diagram that shows a sequence of structural transitions, populating distinct
intermediate states. As $f$ and $C$ are changed, the stem-loop tertiary
interactions rupture first, followed by unfolding of the $3^{\prime}$-end
hairpin ($\textrm{I}\rightleftharpoons\textrm{F}$). Finally, the
$5^{\prime}$-end hairpin unravels, producing an extended state
($\textrm{E}\rightleftharpoons\textrm{I}$). A theoretical analysis of the phase
boundaries shows that the critical force for rupture scales as $\left(\log
C_{\textrm{m}}\right)^{\alpha}$ with $\alpha=1\,(0.5)$ for
$\textrm{E}\rightleftharpoons\textrm{I}$
($\textrm{I}\rightleftharpoons\textrm{F}$) transition. This relation is used to
obtain the preferential ion-RNA interaction coefficient, which can be
quantitatively measured in single-molecule experiments, as done previously for
DNA hairpins. A by-product of our work is the suggestion that the frameshift
efficiency is likely determined by the stability of the $5^{\prime}$-end
hairpin that the ribosome first encounters during translation.


  Specific protein-protein interactions are crucial in the cell, both to ensure
the formation and stability of multi-protein complexes, and to enable signal
transduction in various pathways. Functional interactions between proteins
result in coevolution between the interaction partners, causing their sequences
to be correlated. Here we exploit these correlations to accurately identify
which proteins are specific interaction partners from sequence data alone. Our
general approach, which employs a pairwise maximum entropy model to infer
couplings between residues, has been successfully used to predict the
three-dimensional structures of proteins from sequences. Thus inspired, we
introduce an iterative algorithm to predict specific interaction partners from
two protein families whose members are known to interact. We first assess the
algorithm's performance on histidine kinases and response regulators from
bacterial two-component signaling systems. We obtain a striking 0.93 true
positive fraction on our complete dataset without any a priori knowledge of
interaction partners, and we uncover the origin of this success. We then apply
the algorithm to proteins from ATP-binding cassette (ABC) transporter
complexes, and obtain accurate predictions in these systems as well. Finally,
we present two metrics that accurately distinguish interacting protein families
from non-interacting ones, using only sequence data.


  High-throughput genetic and epigenetic data are often screened for
associations with an observed phenotype. For example, one may wish to test
hundreds of thousands of genetic variants, or DNA methylation sites, for an
association with disease status. These genomic variables can naturally be
grouped by the gene they encode, among other criteria. However, standard
practice in such applications is independent screening with a universal
correction for multiplicity. We propose a Bayesian approach in which the prior
probability of an association for a given genomic variable depends on its gene,
and the gene-specific probabilities are modeled nonparametrically. This
hierarchical model allows for appropriate gene and genome-wide multiplicity
adjustments, and can be incorporated into a variety of Bayesian association
screening methodologies with negligible increase in computational complexity.
We describe an application to screening for differences in DNA methylation
between lower grade glioma and glioblastoma multiforme tumor samples from The
Cancer Genome Atlas. Software is available via the package BayesianScreening
for R at https://github.com/lockEF/BayesianScreening .


  Transient spine enlargement (3-5 min timescale) is an important event
associated with the structural plasticity of dendritic spines. Many of the
molecular mechanisms associated with transient spine enlargement have been
identified experimentally. Here, we use a systems biology approach to construct
a mathematical model of biochemical signaling and actin-mediated transient
spine expansion in response to calcium-influx due to NMDA receptor activation.
We have identified that a key feature of this signaling network is the
paradoxical signaling loop. Paradoxical components act bifunctionally in
signaling networks and their role is to control both the activation and
inhibition of a desired response function (protein activity or spine volume).
Using ordinary differential equation (ODE)-based modeling, we show that the
dynamics of different regulators of transient spine expansion including CaMKII,
RhoA, and Cdc42 and the spine volume can be described using paradoxical
signaling loops. Our model is able to capture the experimentally observed
dynamics of transient spine volume. Furthermore, we show that actin remodeling
events provide a robustness to spine volume dynamics. We also generate
experimentally testable predictions about the role of different components and
parameters of the network on spine dynamics.


  Many cell populations, exemplified by certain tumors, grow approximately
according to a Gompertzian growth model which has a slower approach to an upper
limit than that of logistic growth. Certain populations of animals and other
organisms have also recently been analyzed with the Gompertz model. This
article addresses the question of how long it takes to reduce the population
from one level to a lower one under a schedule of sudden decrements, each of
which removes a given fraction of the cell mass or population. A deterministic
periodic schedule is first examined and yields exact results for the
eradication or extinction time which is defined as that required to reduce the
number of cells to less than unity. The decrements in cell mass at each hit
could correspond to an approximation to reduction of a tumor by external beam
radiation therapy. The effects of variations in magnitude of successive
decrements, the time interval between them, the initial population size and the
intrinsic growth rate are calculated and results presented graphically.
  With a schedule governed by a Poisson process, the number of organisms or
cells satisfies a stochastic differential equation whose solution sample paths
have downward jumps as random times. The moments of the exit time then satisfy
a system of recurrent differential-difference equations. A simple
transformation results in a simpler system which has been studied both
analytically and numerically in the context of interspike intervals of a model
neuron. Results are presented for the mean eradication time for various
frequencies and magnitudes of hits and for various eventual and initial
population sizes. The standard deviation of the eradication time is also
investigated.


  We present a novel method for extracting cancer signatures by applying
statistical risk models (http://ssrn.com/abstract=2732453) from quantitative
finance to cancer genome data. Using 1389 whole genome sequenced samples from
14 cancers, we identify an "overall" mode of somatic mutational noise. We give
a prescription for factoring out this noise and source code for fixing the
number of signatures. We apply nonnegative matrix factorization (NMF) to genome
data aggregated by cancer subtype and filtered using our method. The resultant
signatures have substantially lower variability than those from unfiltered
data. Also, the computational cost of signature extraction is cut by about a
factor of 10. We find 3 novel cancer signatures, including a liver cancer
dominant signature (96% contribution) and a renal cell carcinoma signature (70%
contribution). Our method accelerates finding new cancer signatures and
improves their overall stability. Reciprocally, the methods for extracting
cancer signatures could have interesting applications in quantitative finance.


  A mystery about the origins of life is which molecular structures $-$ and
what spontaneous processes $-$ drove the autocatalytic transition from simple
chemistry to biology? Using the HP lattice model of polymer sequence spaces
leads to the prediction that random sequences of hydrophobic ($H$) and polar
($P$) monomers can collapse into relatively compact structures, exposing
hydrophobic surfaces, acting as primitive versions of today's protein
catalysts, elongating other such HP polymers, as ribosomes would now do. Such
foldamer-catalysts form an autocatalytic set, growing short chains into longer
chains that have particular sequences. The system has capacity for the
multimodality: ability to settle at multiple distinct quasi-stable states
characterized by different groups of dominating polymers. This is a testable
mechanism that we believe is relevant to the early origins of life.


  In this paper we propose an ecological resilience point of view on cancer.
This view is based on the analysis of a simple ODE model for the interactions
between cancer and normal cells. The model presents two regimes for tumor
growth. In the first, cancer arises due to three reasons: a partial corruption
of the functions that avoid the growth of mutated cells, an aggressive
phenotype of tumor cells and exposure to external carcinogenic factors. In this
case, treatments may be effective if they drive the system to the basin of
attraction of the cancer cure state. In the second regime, cancer arises
because the repair system is intrinsically corrupted. In this case, the
complete cure is not possible since the cancer cure state is no more stable,
but tumor recurrence may be delayed if treatment is prolongued. We review three
indicators of the resilience of a stable equilibrium, related with size and
shape of its basin of attraction: latitude, precariousness and resistance. A
novel method to calculate these indicators is proposed. This method is simpler
and more efficient than those currently used, and may be easily applied to
other population dynamics models. We apply this method to the model and
investigate how these indicators behave with parameters changes. Finally, we
present some simulations to illustrate how the resilience analysis can be
applied to validated models in order to obtain indicators for personalized
cancer treatments.
  Keywords: Tumor growth; Chemotherapy; Basins of Attraction; Regime shifts;
Critical transitions


  Lipid bilayers forming biological membranes are known to behave as viscous 2D
fluids on submicrometer scales; usually they contain a large number of active
protein inclusions. Recently, it has been shown [Proc. Nat. Acad. Sci. USA 112,
E3639 (2015)] that such active proteins should in- duce non-thermal fluctuating
lipid flows leading to diffusion enhancement and chemotaxis-like drift for
passive inclusions in biomembranes. Here, a detailed analytical and numerical
investigation of such effects is performed. The attention is focused on the
situations when proteins are concentrated within lipid rafts. We demonstrate
that passive particles tend to become attracted by active rafts and are
accumulated inside them.


  The ensemble of antagonistic matrices is introduced and studied. In
antagonistic matrices the entries $\mathcal A_{i,j}$ and $\mathcal A_{j,i}$ are
real and have opposite signs, or are both zero, and the diagonal is zero. This
generalization of antisymmetric matrices is suggested by the linearized
dynamics of competitive species in ecology.


  The macromolecules that encode and translate information in living systems,
DNA and RNA, exhibit distinctive structural asymmetries, including
homochirality or mirror image asymmetry and $3' - 5'$ directionality, that are
invariant across all life forms. The evolutionary advantages of these broken
symmetries remain unknown. Here we utilize a very simple model of hypothetical
self-replicating polymers to show that asymmetric autocatalytic polymers are
more successful in self-replication compared to their symmetric counterparts in
the Darwinian competition for space and common substrates. This broken-symmetry
property, called asymmetric cooperativity, arises with the maximization of a
replication potential, where the catalytic influence of inter-strand bonds on
their left and right neighbors is unequal. Asymmetric cooperativity also leads
to tentative, qualitative and simple evolution-based explanations for a number
of other properties of DNA that include four nucleotide alphabet, three
nucleotide codons, circular genomes, helicity, anti-parallel double-strand
orientation, heteromolecular base-pairing, asymmetric base compositions, and
palindromic instability, apart from the structural asymmetries mentioned above.
Our model results and tentative explanations are consistent with multiple lines
of experimental evidence, which include evidence for the presence of asymmetric
cooperativity in DNA.


  We propose a robust approach for performing automatic species-level
recognition of fossil pollen grains in microscopy images that exploits both
global shape and local texture characteristics in a patch-based matching
methodology. We introduce a novel criteria for selecting meaningful and
discriminative exemplar patches. We optimize this function during training
using a greedy submodular function optimization framework that gives a
near-optimal solution with bounded approximation error. We use these selected
exemplars as a dictionary basis and propose a spatially-aware sparse coding
method to match testing images for identification while maintaining global
shape correspondence. To accelerate the coding process for fast matching, we
introduce a relaxed form that uses spatially-aware soft-thresholding during
coding. Finally, we carry out an experimental study that demonstrates the
effectiveness and efficiency of our exemplar selection and classification
mechanisms, achieving $86.13\%$ accuracy on a difficult fine-grained species
classification task distinguishing three types of fossil spruce pollen.


  Studies of circular DNA confined to nanofluidic channels are relevant both
from a fundamental polymer-physics perspective and due to the importance of
circular DNA molecules in vivo. We here observe the unfolding of DNA from the
circular to linear configuration as a light-induced double strand break occurs,
characterize the dynamics, and compare the equilibrium conformational
statistics of linear and circular configurations. This is important because it
allows us to determine to which extent existing statistical theories describe
the extension of confined circular DNA. We find that the ratio of the
extensions of confined linear and circular DNA configurations increases as the
buffer concentration decreases. The experimental results fall between
theoretical predictions for the extended de Gennes regime at weaker confinement
and the Odijk regime at stronger confinement. We show that it is possible to
directly distinguish between circular and linear DNA molecules by measuring the
emission intensity from the DNA. Finally, we determine the rate of unfolding
and show that this rate is larger for more confined DNA, possibly reflecting
the corresponding larger difference in entropy between the circular and linear
configurations.


  Understanding and predicting outbreaks of contagious diseases are crucial to
the development of society and public health, especially for underdeveloped
countries. However, challenging problems are encountered because of complex
epidemic spreading dynamics influenced by spatial structure and human dynamics
(including both human mobility and human interaction intensity). We propose a
systematical model to depict nationwide epidemic spreading in C\^{o}te
d'Ivoire, which integrates multiple factors, such as human mobility, human
interaction intensity, and demographic features. We provide insights to aid in
modeling and predicting the epidemic spreading process by data-driven
simulation and theoretical analysis, which is otherwise beyond the scope of
local evaluation and geometrical views. We show that the requirement that the
average local basic reproductive number to be greater than unity is not
necessary for outbreaks of epidemics. The observed spreading phenomenon can be
roughly explained as a heterogeneous diffusion-reaction process by redefining
mobility distance according to the human mobility volume between nodes, which
is beyond the geometrical viewpoint. However, the heterogeneity of human
dynamics still poses challenges to precise prediction.


  We investigate disease spreading on eight empirical data sets of human
contacts (mostly proximity networks recording who is close to whom, at what
time). We compare three levels of representations of these data sets: temporal
networks, static networks and a fully connected topology. We notice that the
difference between the static and fully-connected networks -- with respect to
time to extinction and average outbreak size -- is smaller than between the
temporal and static topologies. This suggests that, for these data sets,
temporal structures influence disease spreading more than static network
structures. To explain the details in the differences between the
representations, we use 32 network measures. This study concur that long-time
temporal structures, like the turnover of nodes and links, are the most
important for the spreading dynamics.


  The computational study of conformational transitions in nucleic acids still
faces many challenges. For example, in the case of single stranded RNA
tetranucleotides, agreement between simulations and experiments is not
satisfactory due to inaccuracies in the force fields commonly used in molecular
dynamics simulations. We here use experimental data collected from
high-resolution X-ray structures to attempt an improvement of the latest
version of the AMBER force field. A modified metadynamics algorithm is used to
calculate correcting potentials designed to enforce experimental distributions
of backbone torsion angles. Replica-exchange simulations of tetranucleotides
including these correcting potentials show significantly better agreement with
independent solution experiments for the oligonucleotides containing pyrimidine
bases. Although the proposed corrections do not seem to be portable to generic
RNA systems, the simulations revealed the importance of the alpha and beta
backbone angles on the modulation of the RNA conformational ensemble. The
correction protocol presented here suggests a systematic procedure for
force-field refinement.


  Many systems in nature, from ferromagnets to flocks of birds, exhibit
ordering phenomena on the large scale. In physical systems order is
statistically robust for large enough dimensions, with relative fluctuations
due to noise vanishing with system size. Several biological systems, however,
are less stable than their physical analogues and spontaneously change their
global state on relatively short timescales. In this paper we show that there
are two crucial ingredients in these systems that enhance the effect of noise,
leading to collective changes of state: the non-symmetric nature of
interactions between individuals, and the presence of local heterogeneities in
the topology of the network. The consequences of these features can be larger
the larger the system size leading to a localization of the fluctuation modes
and a relaxation time that remains finite in the thermodynamic limit. The
system keeps changing its global state in time, being constantly driven out of
equilibrium by spontaneous fluctuations. Our results explain what is observed
in several living and social systems and are consistent with recent
experimental data on bird flocks and other animal groups.


  We focus on the long term dynamics of "killing the winner" Lotka-Volterra
models of marine communities consisting of bacteria, virus, and zooplankton.
Under suitable conditions, it is shown that there is a unique equilibrium with
all populations present which is stable, the system is permanent, and the
limiting behavior of its solutions is strongly constrained.


  It was recently observed by de Vienne et al. that a simple square root
transformation of distances between taxa on a phylogenetic tree allowed for an
embedding of the taxa into Euclidean space. While the justification for this
was based on a diffusion model of continuous character evolution along the
tree, here we give a direct and elementary explanation for it that provides
substantial additional insight. We use this embedding to reinterpret the
differences between the NJ and BIONJ tree building algorithms, providing one
illustration of how this embedding reflects tree structures in data.


  The dynamics of growth of bacterial populations has been extensively studied
for planktonic cells in well-agitated liquid culture, in which all cells have
equal access to nutrients. In the real world, bacteria are more likely to live
in physically structured habitats as colonies, within which individual cells
vary in their access to nutrients. The dynamics of bacterial growth in such
conditions is poorly understood, and, unlike that for liquid culture, there is
not a standard broadly used mathematical model for bacterial populations
growing in colonies in three dimensions (3-d). By extending the classic Monod
model of resource-limited population growth to allow for spatial heterogeneity
in the bacterial access to nutrients, we develop a 3-d model of colonies, in
which bacteria consume diffusing nutrients in their vicinity. By following the
changes in density of E.coli in liquid and embedded in glucose-limited soft
agar, we evaluate the fit of this model to experimental data. The model
accounts for the experimentally observed presence of a sub-exponential,
diffusion-limited growth regime in colonies, which is absent in liquid
cultures. The model predicts and our experiments confirm that, as a consequence
of inter-colony competition for the diffusing nutrients and of cell death,
there is a non-monotonic relationship between total number of colonies within
the habitat and the total number of individual cells in all of these colonies.
This combined theoretical-experimental study reveals that, within 3-d colonies,
E.coli cells are loosely packed, and colonies produce about 2.5 times as many
cells as the liquid culture from the same amount of nutrients. Our model
provides a baseline description of bacterial growth in 3-d, deviations from
which can be used to identify phenotypic heterogeneities and inter-cellular
interactions that further contribute to the structure of bacterial communities.


  Public Goods Games represent one of the most useful tools to study group
interactions between individuals. However, even if they could provide an
explanation for the emergence and stability of cooperation in modern societies,
they are not able to reproduce some key features observed in social and
economical interactions. The typical shape of wealth distribution - known as
Pareto Law - and the microscopic organization of wealth production are two of
them. Here, we introduce a modification to the classical formulation of Public
Goods Games that allows for the emergence of both of these features from first
principles. Unlike traditional Public Goods Games on networks, where players
contribute equally to all the games in which they participate, we allow
individuals to redistribute their contribution according to what they earned in
previous rounds. Results from numerical simulations show that not only a Pareto
distribution for the payoffs naturally emerges but also that if players don't
invest enough in one round they can act as defectors even if they are formally
cooperators. Finally, we also show that the players self-organize in a very
productive backbone that covers almost perfectly the minimum spanning tree of
the underlying interaction network. Our results not only give an explanation
for the presence of the wealth heterogeneity observed in real data but also
points to a conceptual change regarding how cooperation is defined in
collective dilemmas.


  A method is proposed to identify target states that optimize a metastability
index amongst a set of trial states and use these target states as milestones
(or core sets) to build Markov State Models (MSMs). If the optimized
metastability index is small, this automatically guarantees the accuracy of the
MSM, in the sense that the transitions between the target milestones is indeed
approximately Markovian. The method is simple to implement and use, it does not
require that the dynamics on the trial milestones be Markovian, and it also
offers the possibility to partition the system's state-space by assigning every
trial milestone to the target milestones it is most likely to visit next and to
identify transition state regions. Here the method is tested on the Gly-Ala-Gly
peptide, where it shown to correctly identify the expected metastable states in
the dihedral angle space of the molecule without \textit{a~priori} information
about these states. It is also applied to analyze the folding landscape of the
Beta3s mini-protein, where it is shown to identify the folded basin as a
connecting hub between an helix-rich region, which is entropically stabilized,
and a beta-rich region, which is energetically stabilized and acts as a kinetic
trap.


  To better understand protein-solvent interaction we have analyzed a variety
of physical and geometrical properties of the solvent-excluded surfaces (SESs)
over a large set of soluble proteins with crystal structures. We discover that
all have net negative surface charges and permanent electric dipoles. Moreover
both SES area and surface charge as well as several physical and geometrical
properties defined by them change with protein size via well-fitted power laws.
The relevance to protein-solvent interaction of these physical and geometrical
properties is supported by strong correlations between them and known
hydrophobicity scales and by their large changes upon protein unfolding. The
universal existence of negative surface charge and dipole, the characteristic
surface geometry and power laws reveal fundamental but distinct roles of
surface charge and SES in protein-solvent interaction and make it possible to
describe solvation and hydrophobic effect using theories on anion solute in
protic solvent. In particular the great significance of surface charge for
protein-solvent interaction suggests that a change of perception may be needed
since from solvation perspective folding into a native state is to optimize
surface negative charge rather than to minimize the hydrophobic surface area.


  Discrete-state stochastic models have become a well-established approach to
describe biochemical reaction networks that are influenced by the inherent
randomness of cellular events. In the last years severalmethods for accurately
approximating the statistical moments of such models have become very popular
since they allow an efficient analysis of complex networks. We propose a
generalized method of moments approach for inferring the parameters of reaction
networks based on a sophisticated matching of the statistical moments of the
corresponding stochastic model and the sample moments of population snapshot
data. The proposed parameter estimation method exploits recently developed
moment-based approximations and provides estimators with desirable statistical
properties when a large number of samples is available. We demonstrate the
usefulness and efficiency of the inference method on two case studies. The
generalized method of moments provides accurate and fast estimations of unknown
parameters of reaction networks. The accuracy increases when also moments of
order higher than two are considered. In addition, the variance of the
estimator decreases, when more samples are given or when higher order moments
are included.


  Transcription factors (TFs) exert their regulatory action by binding to DNA
with specific sequence preferences. However, different TFs can partially share
their binding sequences due to their common evolutionary origin. This
`redundancy' of binding defines a way of organizing TFs in `motif families' by
grouping TFs with similar binding preferences. Since these ultimately define
the TF target genes, the motif family organization entails information about
the structure of transcriptional regulation as it has been shaped by evolution.
Focusing on the human TF repertoire, we show that a one-parameter evolutionary
model of the Birth-Death-Innovation type can explain the TF empirical
ripartition in motif families, and allows to highlight the relevant
evolutionary forces at the origin of this organization. Moreover, the model
allows to pinpoint few deviations from the neutral scenario it assumes: three
over-expanded families (including HOX and FOX genes), a set of `singleton' TFs
for which duplication seems to be selected against, and a higher-than-average
rate of diversification of the binding preferences of TFs with a Zinc Finger
DNA binding domain. Finally, a comparison of the TF motif family organization
in different eukaryotic species suggests an increase of redundancy of binding
with organism complexity.


  We formulate an info-clustering paradigm based on a multivariate information
measure, called multivariate mutual information, that naturally extends
Shannon's mutual information between two random variables to the multivariate
case involving more than two random variables. With proper model reductions, we
show that the paradigm can be applied to study the human genome and connectome
in a more meaningful way than the conventional algorithmic approach. Not only
can info-clustering provide justifications and refinements to some existing
techniques, but it also inspires new computationally feasible solutions.


  The diffusion of a reactant to a binding target plays a key role in many
biological processes. The reaction-radius at which the reactant and target may
interact is often a small parameter relative to the diameter of the domain in
which the reactant diffuses. We develop uniform in time asymptotic expansions
in the reaction-radius of the full solution to the corresponding diffusion
equations for two separate reactant-target interaction mechanisms: the Doi or
volume reactivity model, and the Smoluchowski-Collins-Kimball partial
absorption surface reactivity model. In the former, the reactant and target
react with a fixed probability per unit time when within a specified
separation. In the latter, upon reaching a fixed separation, they
probabilistically react or the reactant reflects away from the target.
Expansions of the solution to each model are constructed by projecting out the
contribution of the first eigenvalue and eigenfunction to the solution of the
diffusion equation, and then developing matched asymptotic expansions in
Laplace-transform space. Our approach offers an equivalent, but alternative,
method to the pseudo-potential approach we previously employed in [Isaacson and
Newby 2013] for the simpler Smoluchowski pure absorption reaction mechanism. We
find that the resulting asymptotic expansions of the diffusion equation
solutions are identical with the exception of one parameter: the diffusion
limited reaction rates of the Doi and partial absorption models. This
demonstrates that for biological systems in which the reaction-radius is a
small parameter, properly calibrated Doi and partial absorption models may be
functionally equivalent.


  Although carbon nanotubes are potential candidates for DNA encapsulation and
subsequent delivery of biological payloads to living cells, the thermodynamical
spontaneity of DNA encapsulation under physiological conditions is still a
matter of debate. Using enhanced sampling techniques, we show for the first
time that, given a sufficiently large carbon nanotube, the confinement of a
double-stranded DNA segment, 5'-D(*CP*GP*CP*GP*AP*AP*TP*TP*CP*GP*CP*G)-3', is
thermodynamically favourable under physiological environments (134 mM, 310 K, 1
bar), leading to DNA-nanotube hybrids with lower free energy than the
unconfined biomolecule. A diameter threshold of 3 nm is established below which
encapsulation is inhibited. The confined DNA segment maintains its
translational mobility and exhibits the main geometrical features of the
canonical B form. To accommodate itself within the nanopore, the DNA end-to-end
length increases from 3.85 nm up to approximately 4.1 nm, due to a 0.3 nm
elastic expansion of the strand termini. The canonical Watson-Crick H-bond
network is essentially conserved throughout encapsulation, showing that the
contact between the DNA segment and the hydrophobic carbon walls results in
minor rearrangements of the nucleotides H-bonding. The results obtained here
are paramount to the usage of carbon nanotubes as encapsulation media for next
generation drug delivery technologies.


  Population structure and spatial heterogeneity are integral components of
evolutionary dynamics, in general, and of evolution of cooperation, in
particular. Structure can promote the emergence of cooperation in some
populations and suppress it in others. Here, we provide results for weak
selection to favor cooperation on regular graphs for any configuration, meaning
any arrangement of cooperators and defectors. Our results extend previous work
on fixation probabilities of single, randomly placed mutants. We find that for
any configuration cooperation is never favored for birth-death (BD) updating.
In contrast, for death-birth (DB) updating, we derive a simple, computationally
tractable formula for weak selection to favor cooperation when starting from
any configuration containing any number of cooperators and defectors. This
formula elucidates two important features: (i) the takeover of cooperation can
be enhanced by the strategic placement of cooperators and (ii) adding more
cooperators to a configuration can sometimes suppress the evolution of
cooperation. These findings give a formal account for how selection acts on all
transient states that appear in evolutionary trajectories. They also inform the
strategic design of initial states in social networks to maximally promote
cooperation. We also derive general results that characterize the interaction
of any two strategies, not only cooperation and defection.


  We present the R package SimInf which provides an efficient and very flexible
framework to conduct data-driven epidemiological modeling in realistic large
scale disease spread simulations. The framework integrates infection dynamics
in subpopulations as continuous-time Markov chains using the Gillespie
stochastic simulation algorithm and incorporates available data such as births,
deaths and movements as scheduled events at predefined time-points. Using C
code for the numerical solvers and OpenMP to divide work over multiple
processors ensures high performance when simulating a sample outcome. One of
our design goal was to make SimInf extendable and enable usage of the numerical
solvers from other R extension packages in order to facilitate complex
epidemiological research. In this paper, we provide a technical description of
the framework and demonstrate its use on some basic examples. We also discuss
how to specify and extend the framework with user-defined models.


  Co-expression network is a critical technique for the identification of
inter-gene interactions, which usually relies on all-pairs correlation (or
similar measure) computation between gene expression profiles across multiple
samples. Pearson's correlation coefficient (PCC) is one widely used technique
for gene co-expression network construction. However, all-pairs PCC computation
is computationally demanding for large numbers of gene expression profiles,
thus motivating our acceleration of its execution using high-performance
computing. In this paper, we present LightPCC, the first parallel and
distributed all-pairs PCC computation on Intel Xeon Phi (Phi) clusters. It
achieves high speed by exploring the SIMD-instruction-level and thread-level
parallelism within Phis as well as accelerator-level parallelism among multiple
Phis. To facilitate balanced workload distribution, we have proposed a general
framework for symmetric all-pairs computation by building bijective functions
between job identifier and coordinate space for the first time. We have
evaluated LightPCC and compared it to two CPU-based counterparts: a sequential
C++ implementation in ALGLIB and an implementation based on a parallel general
matrix-matrix multiplication routine in Intel Math Kernel Library (MKL) (all
use double precision), using a set of gene expression datasets. Performance
evaluation revealed that with one 5110P Phi and 16 Phis, LightPCC runs up to
$20.6\times$ and $218.2\times$ faster than ALGLIB, and up to $6.8\times$ and
$71.4\times$ faster than single-threaded MKL, respectively. In addition,
LightPCC demonstrated good parallel scalability in terms of number of Phis.
Source code of LightPCC is publicly available at
http://lightpcc.sourceforge.net.


  Statistical mechanics provides an effective framework to investigate
information processing in biochemical reactions. Within such framework
far-reaching analogies are established among (anti-) cooperative collective
behaviors in chemical kinetics, (anti-)ferromagnetic spin models in statistical
mechanics and operational amplifiers/flip-flops in cybernetics. The underlying
modeling -- based on spin systems -- has been proved to be accurate for a wide
class of systems matching classical (e.g. Michaelis--Menten, Hill, Adair)
scenarios in the infinite-size approximation. However, the current research in
biochemical information processing has been focusing on systems involving a
relatively small number of units, where this approximation is no longer valid.
Here we show that the whole statistical mechanical description of reaction
kinetics can be re-formulated via a mechanical analogy -- based on completely
integrable hydrodynamic-type systems of PDEs -- which provides explicit
finite-size solutions, matching recently investigated phenomena (e.g.
noise-induced cooperativity, stochastic bi-stability, quorum sensing). The
resulting picture, successfully tested against a broad spectrum of data,
constitutes a neat rationale for a numerically effective and theoretically
consistent description of collective behaviors in biochemical reactions.


  Cells adapt their metabolic fluxes in response to changes in the environment.
We present a framework for the systematic construction of flux-based graphs
derived from organism-wide metabolic networks. Our graphs encode the
directionality of metabolic fluxes via edges that represent the flow of
metabolites from source to target reactions. The methodology can be applied in
the absence of a specific biological context by modelling fluxes
probabilistically, or can be tailored to different environmental conditions by
incorporating flux distributions computed through constraint-based approaches
such as Flux Balance Analysis. We illustrate our approach on the central carbon
metabolism of Escherichia coli and on a metabolic model of human hepatocytes.
The flux-dependent graphs under various environmental conditions and genetic
perturbations exhibit systemic changes in their topological and community
structure, which capture the re-routing of metabolic fluxes and the varying
importance of specific reactions and pathways. By integrating constraint-based
models and tools from network science, our framework allows the study of
context-specific metabolic responses at a system level beyond standard pathway
descriptions.


  Earlier formulations of the DNA assembly problem were all in the context of
perfect assembly; i.e., given a set of reads from a long genome sequence, is it
possible to perfectly reconstruct the original sequence? In practice, however,
it is very often the case that the read data is not sufficiently rich to permit
unambiguous reconstruction of the original sequence. While a natural
generalization of the perfect assembly formulation to these cases would be to
consider a rate-distortion framework, partial assemblies are usually
represented in terms of an assembly graph, making the definition of a
distortion measure challenging. In this work, we introduce a distortion
function for assembly graphs that can be understood as the logarithm of the
number of Eulerian cycles in the assembly graph, each of which correspond to a
candidate assembly that could have generated the observed reads. We also
introduce an algorithm for the construction of an assembly graph and analyze
its performance on real genomes.


  One of the most highly debated questions in the field of animal swarming and
social behaviour, is the collective random patterns and chaotic behaviour
formed by some animal species, in particular if there is a danger. Is such a
behaviour beneficial or unfavourable for survival? Here we report on one of the
most remarkable forms of animal swarming and social behaviour - fish schooling
- from a hydrodynamic point of view. We found that some fish species do not
have preferred orientation and they swarm in a random pattern mode, despite the
excess of energy consumed. Our analyses, which includes calculations of the
hydrodynamic forces between slender bodies, show that such a behaviour enhances
the transfer of hydrodynamic information, and thus enhances the survivability
of the school. These findings support the general hypothesis that a disordered
and non-trivial collective behaviour of individuals within a nonlinear
dynamical system is essential for optimising transfer of information - an
optimisation that might be crucial for survival.


  This work is devoted to studying the dynamics of a structured population that
is subject to the combined effects of environmental stochasticity, competition
for resources, spatio-temporal heterogeneity and dispersal. The population is
spread throughout $n$ patches whose population abundances are modelled as the
solutions of a system of nonlinear stochastic differential equations living on
$[0,\infty)^n$.
  We prove that $r$, the stochastic growth rate of the total population in the
absence of competition, determines the long-term behaviour of the population.
The parameter $r$ can be expressed as the Lyapunov exponent of an associated
linearized system of stochastic differential equations. Detailed analysis shows
that if $r>0$, the population abundances converge polynomially fast to a unique
invariant probability measure on $(0,\infty)^n$, while when $r<0$, the
population abundances of the patches converge almost surely to $0$
exponentially fast. This generalizes and extends the results of Evans et al
(2014 J. Math. Biol.) and proves one of their conjectures.
  Compared to recent developments, our model incorporates very general
density-dependent growth rates and competition terms. Furthermore, we prove
that persistence is robust to small, possibly density dependent, perturbations
of the growth rates, dispersal matrix and covariance matrix of the
environmental noise. Our work allows the environmental noise driving our system
to be degenerate. This is relevant from a biological point of view since, for
example, the environments of the different patches can be perfectly correlated.
As an example we fully analyze the two-patch case, $n=2$, and show that the
stochastic growth rate is a decreasing function of the dispersion rate. In
particular, coupling two sink patches can never yield persistence, in contrast
to the results from the non-degenerate setting treated by Evans et al.


  Predicting the fate of ecologies is a daunting, albeit extremely important,
task. As part of this task one needs to develop an understanding of the
organization, hierarchies, and correlations among the species forming the
ecology. Focusing on complex food networks we present a theoretical method that
allows to achieve this understanding. Starting from the adjacency matrix the
method derives specific matrices that encode the various inter-species
relationships. The full potential of the method is achieved in a spatial
setting where one obtains detailed predictions for the emerging space-time
patterns. For a variety of cases these theoretical predictions are verified
through numerical simulations.


  Indirect reciprocity is a major mechanism in the maintenance of cooperation
among unrelated individuals. Indirect reciprocity leads to conditional
cooperation according to social norms that discriminate the good (those who
deserve to be rewarded with help) and the bad (those who should be punished by
refusal of help). Despite intensive research, however, there is no definitive
consensus on what social norms best promote cooperation through indirect
reciprocity, and it remains unclear even how those who refuse to help the bad
should be assessed. Here, we propose a new simple norm called "Staying" that
prescribes abstaining from assessment. Under the Staying norm, the image of the
person who makes the decision to give help stays the same as in the last
assessment if the person on the receiving end has a bad image. In this case,
the choice about whether or not to give help to the potential receiver does not
affect the image of the potential giver. We analyze the Staying norm in terms
of evolutionary game theory and demonstrate that Staying is most effective in
establishing cooperation compared to the prevailing social norms, which rely on
constant monitoring and unconditional assessment. The application of Staying
suggests that the strict application of moral judgment is limited.


  Expression of many genes varies as a cell transitions through different
cell-cycle stages. How coupling between stochastic expression and cell cycle
impacts cell-to-cell variability (noise) in the level of protein is not well
understood. We analyze a model, where a stable protein is synthesized in random
bursts, and the frequency with which bursts occur varies within the cell cycle.
Formulas quantifying the extent of fluctuations in the protein copy number are
derived and decomposed into components arising from the cell cycle and
stochastic processes. The latter stochastic component represents contributions
from bursty expression and errors incurred during partitioning of molecules
between daughter cells. These formulas reveal an interesting trade-off:
cell-cycle dependencies that amplify the noise contribution from bursty
expression also attenuate the contribution from partitioning errors. We
investigate existence of optimum strategies for coupling expression to the cell
cycle that minimize the stochastic component. Intriguingly, results show that a
zero production rate throughout the cell cycle, with expression only occurring
just before cell division minimizes noise from bursty expression for a fixed
mean protein level. In contrast, the optimal strategy in the case of
partitioning errors is to make the protein just after cell division. We provide
examples of regulatory proteins that are expressed only towards the end of cell
cycle, and argue that such strategies enhance robustness of cell-cycle
decisions to the intrinsic stochasticity of gene expression.


  As gene expression measurement technology is shifting from microarrays to
sequencing, the statistical tools available for their analysis must be adapted
since RNA-seq data are measured as counts. Recently, it has been proposed to
tackle the count nature of these data by modeling log-count reads per million
as continuous variables, using nonparametric regression to account for their
inherent heteroscedasticity. Adopting such a framework, we propose tcgsaseq, a
principled, model-free and efficient top-down method for detecting longitudinal
changes in RNA-seq gene sets. Considering gene sets defined a priori, tcgsaseq
identifies those whose expression vary over time, based on an original variance
component score test accounting for both covariates and heteroscedasticity
without assuming any specific parametric distribution for the transformed
counts. We demonstrate that despite the presence of a nonparametric component,
our test statistic has a simple form and limiting distribution, and both may be
computed quickly. A permutation version of the test is additionally proposed
for very small sample sizes. Applied to both simulated data and two real
datasets, the proposed method is shown to exhibit very good statistical
properties, with an increase in stability and power when compared to state of
the art methods ROAST, edgeR and DESeq2, which can fail to control the type I
error under certain realistic settings. We have made the method available for
the community in the R package tcgsaseq.


  In this paper we introduce a novel, context-free grammar, {\it
RNAFeatures$^*$}, capable of generating any RNA structure including pseudoknot
structures (pk-structure). We represent pk-structures as orientable fatgraphs,
which naturally leads to a filtration by their topological genus. Within this
framework, RNA secondary structures correspond to pk-structures of genus zero.
{\it RNAFeatures$^*$} acts on formal, arc-labeled RNA secondary structures,
called $\lambda$-structures. $\lambda$-structures correspond one-to-one to
pk-structures together with some additional information. This information
consists of the specific rearrangement of the backbone, by which a pk-structure
can be made cross-free. {\it RNAFeatures$^*$} is an extension of the grammar
for secondary structures and employs an enhancement by labelings of the symbols
as well as the production rules. We discuss how to use {\it RNAFeatures$^*$} to
obtain a stochastic context-free grammar for pk-structures, using data of RNA
sequences and structures. The induced grammar facilitates fast Boltzmann
sampling and statistical analysis. As a first application, we present an $O(n
log(n))$ runtime algorithm which samples pk-structures based on ninety tRNA
sequences and structures from the Nucleic Acid Database (NDB).


  We exhibit a fundamental relationship between measures of dynamical and
structural stability of equilibriums, arising from real dynamical systems. We
show that dynamical stability, quantified via systems local response to
external perturbations, coincides with the minimal internal perturbation able
to destabilize the equilibrium. First, by reformulating a result of control
theory, we explain that harmonic external perturbations reflect the spectral
sensitivity of the Jacobian matrix at the equilibrium, with respect to constant
changes of its coefficients. However, for this equivalence to hold, imaginary
changes of the Jacobian's coefficients have to be allowed. The connection with
dynamical stability is thus lost for real dynamical systems. We show that this
issue can be avoided, thus recovering the fundamental link between dynamical
and structural stability, by considering stochastic noise as external and
internal perturbations. More precisely, we demonstrate that a system's local
response to white-noise perturbations directly reflects the intensity of
internal white noise that it can accommodate before asymptotic mean-square
stability of the equilibrium is lost.


  We consider an SIR-type (Susceptible $\to$ Infected $\to$ Recovered)
stochastic epidemic process with multiple modes of transmission on a contact
network. The network is given by a random graph following a multilayer
configuration model where edges in different layers correspond to potentially
infectious contacts of different types. We assume that the graph structure
evolves in response to the epidemic via activation or deactivation of edges. We
derive a large graph limit theorem that gives a system of ordinary differential
equations (ODEs) describing the evolution of quantities of interest, such as
the proportions of infected and susceptible vertices, as the number of nodes
tends to infinity. Analysis of the limiting system elucidates how the coupling
of edge activation and deactivation to infection status affects disease
dynamics, as illustrated by a two-layer network example with edge types
corresponding to community and healthcare contacts. Our theorem extends some
earlier results deriving the deterministic limit of stochastic SIR processes on
static, single-layer configuration model graphs. We also describe precisely the
conditions for equivalence between our limiting ODEs and the systems obtained
via pair approximation, which are widely used in the epidemiological and
ecological literature to approximate disease dynamics on networks. Potential
applications include modeling Ebola dynamics in West Africa, which was the
motivation for this study.


  Genome-wide patterns of genetic divergence reveal mechanisms of adaptation
under gene flow. Empirical data show that divergence is mostly concentrated in
narrow genomic regions. This pattern may arise because differentiated loci
protect nearby mutations from gene flow, but recent theory suggests this
mechanism is insufficient to explain the emergence of concentrated
differentiation during biologically realistic timescales. Critically, earlier
theory neglects an inevitable consequence of genetic drift: stochastic loss of
local genomic divergence. Here we demonstrate that the rate of stochastic loss
of weak local differentiation increases with recombination distance to a
strongly diverged locus and, above a critical recombination distance, local
loss is faster than local `gain' of new differentiation. Under high migration
and weak selection this critical recombination distance is much smaller than
the total recombination distance of the genomic region under selection.
Consequently, divergence between populations increases by net gain of new
differentiation within the critical recombination distance, resulting in
tightly-linked clusters of divergence. The mechanism responsible is the balance
between stochastic loss and gain of weak local differentiation, a mechanism
acting universally throughout the genome. Our results will help to explain
empirical observations and lead to novel predictions regarding changes in
genomic architectures during adaptive divergence.


  Filopodia are ubiquitous fingerlike protrusions, spawned by many eukaryotic
cells, to probe and interact with their environments. Polymerization dynamics
of actin filaments, comprising the structural core of filopodia, largely
determine their instantaneous lengths and overall lifetimes. The polymerization
reactions at the filopodial tip require transport of G-actin, which enter the
filopodial tube from the filopodial base and diffuse toward the filament barbed
ends near the tip. Actin filaments are mechanically coupled into a tight bundle
by cross-linker proteins. Interestingly, many of these proteins are relatively
short, restricting the free diffusion of cytosolic G-actin throughout the
bundle and, in particular, its penetration into the bundle core. To investigate
the effect of steric restrictions on G-actin diffusion by the porous structure
of filopodial actin filament bundle, we used a particle-based stochastic
simulation approach. We discovered that excluded volume interactions result in
partial and then full collapse of central filaments in the bundle, leading to a
hollowed-out structure. The latter may further collapse radially due to the
activity of cross-linking proteins, hence producing conical-shaped filament
bundles. Interestingly, electron microscopy experiments on mature filopodia
indeed frequently reveal actin bundles that are narrow at the tip and wider at
the base. Overall, our work demonstrates that excluded volume effects in the
context of reaction-diffusion processes in porous networks may lead to
unexpected geometric growth patterns and complicated, history-dependent
dynamics of intermediate metastable configurations.


  Formal analysis of the emergent structural properties of dynamic networks is
largely uncharted territory. We focus here on the properties of forward
reachable sets (FRS) as a function of the underlying degree distribution and
edge duration. FRS are defined as the set of nodes that can be reached from an
initial seed via a path of temporally ordered edges; a natural extension of
connected component measures to dynamic networks. Working in a stochastic
framework, we derive closed-form expressions for the mean and variance of the
exponential growth rate of the FRS for temporal networks with both edge and
node dynamics. For networks with node dynamics, we calculate thresholds for the
growth of the FRS. The effects of finite population size are explored via
simulation and approximation. We examine how these properties vary by edge
duration and different cross-sectional degree distributions that characterize a
range of scientifically interesting normative outcomes (Poisson and Bernoulli).
The size of the forward reachable set gives an upper bound for the epidemic
size in disease transmission network models, relating this work to epidemic
modeling (Ferguson 2000, Eames 2004).


  Transcriptional regulation by transcription factors and post-transcriptional
regulation by microRNAs constitute two major modes of regulation of gene
expression. While gene expression motifs incorporating solely transcriptional
regulation are well investigated, the dynamics of motifs with dual strategies
of regulation, i.e., both transcriptional and post-transcriptional regulation,
have not been studied as extensively. In this paper, we probe the dynamics of a
four-gene motif with dual strategies of regulation of gene expression. Some of
the functional characteristics are compared with those of a two-gene motif, the
genetic toggle, employing only transcriptional regulation. Both the motifs
define positive feedback loops with the potential for bistability and
hysteresis. The four-gene motif, contrary to the genetic toggle, is found to
exhibit bistability even in the absence of cooperativity in the regulation of
gene expression. The four-gene motif further exhibits a novel dynamical feature
in which two regions of monostability with linear threshold response are
separated by a region of bistability with digital response.
  Using the linear noise approximation (LNA), we further show that the
coefficient of variation (a measure of noise), associated with the protein
levels in the steady state, has a lower magnitude in the case of the four-gene
motif as compared to the case of the genetic toggle. We next compare
transcriptional with post-transcriptional regulation from an information
theoretic perspective. We focus on two gene expression motifs, Motif 1 with
transcriptional regulation and Motif 2 with post-transcriptional regulation.
  We show that amongst the two motifs, Motif 2 has a greater capacity for
information transmission for an extended range of parameter values.


  In the first part of the paper I will present a brief review on the
Hardy-Weinberg equilibrium and its formulation in projective algebraic
geometry. In the second and last part I will discuss examples and
generalizations on the topic.


  We develop a stochastic model for the programmed frameshift of ribosomes
synthesizing a protein while moving along a mRNA template. Normally the reading
frame of a ribosome decodes successive triplets of nucleotides on the mRNA in a
step-by-step manner. We focus on the programmed shift of the ribosomal reading
frame, forward or backward, by only one nucleotide which results in a fusion
protein; it occurs when a ribosome temporarily loses its grip to its mRNA
track. Special "slippery" sequences of nucleotides and also downstream
secondary structures of the mRNA strand are believed to play key roles in
programmed frameshift. Here we explore the role of an hitherto neglected
parameter in regulating -1 programmed frameshift. Specifically, we demonstrate
that the frameshift frequency can be strongly regulated also by the density of
the ribosomes, all of which are engaged in simultaneous translation of the same
mRNA, at and around the slippery sequence. Monte Carlo simulations support the
analytical predictions obtained from a mean-field analysis of the stochastic
dynamics.


  The theory of interaction-based evolution argues that, at the most basic
level of analysis, there is a third alternative for how adaptive evolution
works besides a) accidental mutation and natural selection and b) Lamarckism,
namely, c) information provided by natural selection on the fit between the
organism and its environment is absorbed by non-accidental mutation. This
non-accidental mutation is non-Lamarckian yet useful for evolution, and is due
to evolved and continually evolving mutational mechanisms operating in the germ
cells. However, this theory has left a fundamental problem open: If mutational
mechanisms are not Lamarckian---if they are not "aware" of the environment and
the macroscale phenotype---then how could heritable novelty be due to anything
other than accidental mutation? This paper aims to address this question by
arguing the following. Mutational mechanisms can be broadly construed as
enacting local simplification operations on the DNA in germ cells, along with
gene duplication. The joint action of these mutational operations and natural
selection provides simplification under performance pressure. This joint action
creates from preexisting biological interactions new elements that have the
inherent capacity to come together into unexpected useful interactions with
other such elements, thus explaining nature's tendency for cooption. Novelty
thus arises not from a local genetic accident but from gradual network-level
evolution. Many empirical observations are explained from this perspective,
from cooption and gene fusion at the molecular level, to the evolution of
behavior and instinct at the organismal level. Finally, the nature of
mutational mechanisms and the need to study them in detail are described, and a
connection is drawn between evolution and learning.


  Cells can be considered as systems that utilize changes in thermodynamic
entropy as information. Therefore, they serve as useful models for
investigating the relationships between entropy production and information
transmission, i.e., signal transduction. Based on the hypothesis that cells
apply a chemical reaction cascade for the most efficient transduction of
information, we adopted a coding design that minimizes the number of bits per
concentration of molecules that are employed for information transduction. As a
result, the average rate of entropy production is uniform across all cycles in
a cascade reaction. Thus, the entropy production rate can be a valuable measure
for the quantification of intracellular signal transduction.


  We consider a very general stochastic model for an SIR epidemic on a network
which allows an individual's infectious period, and the time it takes to
contact each of its neighbours after becoming infected, to be correlated. We
write down the message passing system of equations for this model and prove,
for the first time, that it has a unique feasible solution. We also generalise
an earlier result by proving that this solution provides a rigorous upper bound
for the expected epidemic size (cumulative number of infection events) at any
fixed time $t>0$.
  We specialise these results to a homogeneous special case where the graph
(network) is symmetric. The message passing system here reduces to just four
equations. We prove that cycles in the network inhibit the spread of infection,
and derive important epidemiological results concerning the final epidemic size
and threshold behaviour for a major outbreak. For Poisson contact processes,
this message passing system is equivalent to a non-Markovian pair approximation
model, which we show has well-known pairwise models as special cases.
  We show further that a sequence of message passing systems, starting with the
homogeneous one just described, converges to the deterministic
Kermack-McKendrick equations for this stochastic model. For Poisson contact and
recovery, we show that this convergence is monotone, from which it follows that
the message passing system (and hence also the pairwise model) here provides a
better approximation to the expected epidemic size at time $t>0$ than the
Kermack-McKendrick model.


  Nosocomial outbreaks of bacteria are well-documented. Based on these
incidents, and the heavy usage of antibiotics in hospitals, it has been assumed
that antibiotic resistance evolves in hospital environments. To test this
assumption, we studied resistance phenotypes of bacteria collected from patient
isolates at a community hospital over a 2.5-year period. A graphical model
analysis shows no association between resistance and patient information other
than time of arrival. This allows us to focus on time course data.
  We introduce a Hospital Transmission Model, based on negative binomial delay.
Our main contribution is a statistical hypothesis test called the Nosocomial
Evolution of Resistance Detector (NERD). It calculates the significance of
resistance trends occurring in a hospital. It can inform hospital staff about
the effects of various practices and interventions, can help detect clonal
outbreaks, and is available as an R-package.
  We applied the NERD method to each of the 16 antibiotics in the study via 16
hypothesis tests. For 13 of the antibiotics, we found that the hospital
environment had no significant effect upon the evolution of resistance; the
hospital is merely a piece of the larger picture. The p-values obtained for the
other three antibiotics (Cefepime, Ceftazidime and Gentamicin) indicate that
particular care should be taken in hospital practices with these antibiotics.
One of the three, Ceftazidime, was significant after accounting for multiple
hypotheses, indicating a trend of decreased resistance for this drug.


  Understanding protein-protein interactions is central to our understanding of
almost all complex biological processes. Computational tools exploiting rapidly
growing genomic databases to characterize protein-protein interactions are
urgently needed. Such methods should connect multiple scales from evolutionary
conserved interactions between families of homologous proteins, over the
identification of specifically interacting proteins in the case of multiple
paralogs inside a species, down to the prediction of residues being in physical
contact across interaction interfaces. Statistical inference methods detecting
residue-residue coevolution have recently triggered considerable progress in
using sequence data for quaternary protein structure prediction; they require,
however, large joint alignments of homologous protein pairs known to interact.
The generation of such alignments is a complex computational task on its own;
application of coevolutionary modeling has in turn been restricted to proteins
without paralogs, or to bacterial systems with the corresponding coding genes
being co-localized in operons. Here we show that the Direct-Coupling Analysis
of residue coevolution can be extended to connect the different scales, and
simultaneously to match interacting paralogs, to identify inter-protein
residue-residue contacts and to discriminate interacting from noninteracting
families in a multiprotein system. Our results extend the potential
applications of coevolutionary analysis far beyond cases treatable so far.


  In this work we consider a model of epidemic spreading coupled with an
opinion dynamics in a fully-connected population. Regarding the opinion
dynamics, the individuals may be in two distinct states, namely in favor or
against a vaccination campaign. Individuals against the vaccination follow a
standard SIS model, whereas the pro-vaccine individuals can also be in a third
compartment, namely Vaccinated. In addition, the opinions change according to
the majority-rule dynamics in groups with three individuals. We also consider
that the vaccine can give permanent or temporary immunization to the
individuals. By means of analytical calculations and computer simulations, we
show that the opinion dynamics can drastically affect the disease propagation,
and that the engagement of the pro-vaccine individuals can be crucial for
stopping the epidemic spreading. The full numerical code for simulate the model
is available from the authors' webpage.


  All living systems -- from the origin of life to modern cells -- rely on a
set of biochemical reactions that are simultaneously self-sustaining and
autocatalytic. This notion of an autocatalytic set has been formalized
graph-theoretically (as `RAF'), leading to mathematical results and
polynomial-time algorithms that have been applied to simulated and real
chemical reaction systems. In this paper, we investigate the emergence of
autocatalytic sets in polymer models when the catalysis rate of each molecule
type is drawn from some probability distribution. We show that although the
average catalysis rate $f$ for RAFs to arise depends on this distribution, a
universal linear upper and lower bound for $f$ (with increasing system size)
still applies. However, the probability of the appearance (and size) of
autocatalytic sets can vary widely, depending on the particular catalysis
distribution. We use simulations to explore how tight the mathematical bounds
are, and the reasons for the observed variations. We also investigate the
impact of inhibition (where molecules can also inhibit reactions) on the
emergence of autocatalytic sets, deriving new mathematical and algorithmic
results.