  We consider a model ecosystem of sessile species competing for space. In
particular, we consider the system introduced in [Mathiesen et al. Phys. Rev.
Lett. 107, 188101 (2011)] where species compete according to a fixed
interaction network with links determined by a Bernoulli process. In the limit
of a small introduction rate of new species, the model exhibits a discontinuous
transition from a high-diversity state to a low-diversity state as the
interaction probability between species, $\gamma$, is increased from zero. Here
we explore the effects of finite introduction rates and system-size on the
phase transition by utilizing efficient parallel computing. We find that the
low state appears for $\gamma>\gamma_c$. As $\gamma$ is increased further, the
high state approaches to the low state, suggesting the possibility that the two
states merge at a high $\gamma$. We find that the fraction of time spent in the
high state becomes longer with higher introduction rates, but the availability
of the two states is rather insensitive to the value of the introduction rate.
Furthermore we establish a relation between the introduction rate and the
system size, which preserves the probability for the system to remain in the
high-diversity state.


  During an epidemic control, the containment of the disease is usually
achieved through increasing devoted resource to shorten the duration of
infectiousness. However, the impact of this resource expenditure has not been
studied quantitatively. Using the well-documented cholera data, we observe
empirically that the recovery rate which is related to the duration of
infectiousness has a strong positive correlation with the average resource
devoted to the infected individuals. By incorporating this relation we build a
novel model and find that insufficient resource leads to an abrupt increase in
the infected population size, which is in marked contrast with the continuous
phase transitions believed previously. Counterintuitively, this abrupt phase
transition is more pronounced in the less contagious diseases, which usually
correspond to the most fatal ones. Furthermore, we find that even for a single
infection source, public resource needs to meet a significant amount, which is
proportional to the whole population size to ensure epidemic containment. Our
findings provide a theoretical foundation for efficient epidemic containment
strategies in the early stage.


  Nowadays, in our globalized world,the local and intercountry movements of
population have been increased. This situation makes it important for host
countries to do right predictions for the future population of their native
people as well as immigrant people. The knowledge of the attained number of
accumulated population is necessary for future planning, concerning to
education,health, job, housing, safety requirements, etc. In this work, for
updating historically well known formulas of population dynamics of a community
are revisited in the framework of compound growth and fractional calculus to
get more realistic relations. Within this context, for a time t, the population
evolution of a society which owns two different components is calculated.
Concomitant relations have been developed to provide a comparison between the
native population and the immigrant population that come into existence where
at each time interval a colonial population is joined. Eventually at time t,
the case where the native population becomes equal to immigrant population is
discussed. Moreover, the differential equation of a probability function for
the equilibrium state of accumulative population is obtained. It is seen that
the equilibrium solution can be written as a series in terms of Bernoulli
numbers. In the calculation of population dynamics, Mittag- Leffler(M-L)
function replaces the exponential function.


  Reconstruction of neuroanatomy is a fundamental problem in neuroscience.
Stochastic expression of colors in individual cells is a promising tool,
although its use in the nervous system has been limited due to various sources
of variability in expression. Moreover, the intermingled anatomy of neuronal
trees is challenging for existing segmentation algorithms. Here, we propose a
method to automate the segmentation of neurons in such (potentially
pseudo-colored) images. The method uses spatio-color relations between the
voxels, generates supervoxels to reduce the problem size by four orders of
magnitude before the final segmentation, and is parallelizable over the
supervoxels. To quantify performance and gain insight, we generate simulated
images, where the noise level and characteristics, the density of expression,
and the number of fluorophore types are variable. We also present segmentations
of real Brainbow images of the mouse hippocampus, which reveal many of the
dendritic segments.


  In multi-resolution simulations, different system components are
simultaneously modelled at different levels of resolution, these being smoothly
coupled together. In the case of enzyme systems, computationally expensive
atomistic detail is needed in the active site to capture the chemistry of
substrate binding. Global properties of the rest of the protein also play an
essential role, determining the structure and fluctuations of the binding site;
however, these can be modelled on a coarser level. Similarly, in the most
computationally efficient scheme only the solvent hydrating the active site
requires atomistic detail. We present a methodology to couple atomistic and
coarse-grained protein models, while solvating the atomistic part of the
protein in atomistic water. This allows a free choice of which protein and
solvent degrees of freedom to include atomistically, without loss of accuracy
in the atomistic description. This multi-resolution methodology can
successfully model stable ligand binding, and we further confirm its validity
via an exploration of system properties relevant to enzymatic function. In
addition to a computational speedup, such an approach can allow the
identification of the essential degrees of freedom playing a role in a given
process, potentially yielding new insights into biomolecular function.


  Structure, function and dynamics of many biomolecular systems can be
characterized by the energetic variational principle and the corresponding
systems of partial differential equations (PDEs). This principle allows us to
focus on the identification of essential energetic components, the optimal
parametrization of energies, and the efficient computational implementation of
energy variation or minimization. Given the fact that complex biomolecular
systems are structurally non-uniform and their interactions occur through
contact interfaces, their free energies are associated with various interfaces
as well, such as solute-solvent interface, molecular binding interface, lipid
domain interface, and membrane surfaces. This fact motivates the inclusion of
interface geometry, particular its curvatures, to the parametrization of free
energies. Applications of such interface geometry based energetic variational
principles are illustrated through three concrete topics: the multiscale
modeling of biomolecular electrostatics and solvation that includes the
curvature energy of the molecular surface, the formation of microdomains on
lipid membrane due to the geometric and molecular mechanics at the lipid
interface, and the mean curvature driven protein localization on membrane
surfaces. By further implicitly representing the interface using a phase field
function over the entire domain, one can simulate the dynamics of the interface
and the corresponding energy variation by evolving the phase field function,
achieving significant reduction of the number of degrees of freedom and
computational complexity. Strategies for improving the efficiency of
computational implementations and for extending applications to coarse-graining
or multiscale molecular simulations are outlined.


  Several studies highlighted the relevance of extrinsic noise in shaping cell
decision making and differentiation in molecular networks. Experimental
evidences of phenotypic differentiation are given by the presence of bimodal
distributions of gene expression levels, where the modes of the distribution
often correspond to different physiological states of the system. We
theoretically address the presence of bimodal phenotypes in the context of
microRNA (miRNA)-mediated regulation. MiRNAs are small noncoding RNA molecules
that downregulate the expression of their target mRNAs. The nature of this
interaction is titrative and induces a threshold effect: below a given target
transcription rate no mRNAs are free and available for translation. We
investigate the effect of extrinsic noise on the system by introducing a
fluctuating miRNA-transcription rate. We find that the presence of extrinsic
noise favours the presence of bimodal target distributions which can be
observed for a wider range of parameters compared to the case with intrinsic
noise only and for lower miRNA-target interaction strength. Our results suggest
that combining threshold-inducing interactions with extrinsic noise provides a
simple and robust mechanism for obtaining bimodal populations not requiring
fine tuning. We furthermore characterise the protein distributions dependence
on protein half-life.


  Membrane proteins deform the surrounding lipid bilayer, which can lead to
membrane-mediated interactions between neighboring proteins. Using the
mechanosensitive channel of large conductance (MscL) as a model system, we
demonstrate how the observed differences in protein structure can affect
membrane-mediated interactions and cooperativity among membrane proteins. We
find that distinct oligomeric states of MscL lead to distinct gateway states
for the clustering of MscL, and predict signatures of MscL structure and
spatial organization in the cooperative gating of MscL. Our modeling approach
establishes a quantitative relation between the observed shapes and cooperative
function of membrane~proteins.


  Previous kinetic models had assumed that the reaction medium was reacting at
random and without a turnover associated to thermodynamics exchanges, with a
rigid active site on the enzyme. The experimental studies show that coupling
factor 1 (CF1) from spinach chloroplasts has latent ATPase activity, which
become expressed after heat-treatment and incubation with calcium. The
sigmoidal kinetics observed on the competitive effect of glycerol on water
saturating a protein, suggests that the role of the hydration shell in the
catalytic mechanism of the CF1-ATPase, modify the number of water molecules
associated with the conformational turnover required for active site activity.
It is assume that the water associated to the hydrophilic state of the enzyme
produces a fit-in of the substrate to form (ES), follow by the catalytic action
with product formation (EP). This one induces the dissociation of water and
increases the hydrophobic attractions between R-groups. The latter, becomes the
form of the enzyme interacting with water to form the dissociated free enzyme
(E) and free product (P). Glycerol dependent suppression of the water dynamics
on two interacting sites configuration shows a change in the
H-bond-configuration. The thermodynamics modeling requires an energy
expenditure of about 4kcal/mol per each H-bond in turnover. Glycerol determined
a turnover of 14 molecules of water released from the active sites to reach the
inactive form of the enzyme. The entropy generated by turnover of the fit-in
and -out substrate and product from the protein could be dissipated-out of the
enzyme-water system itself. Coupling with the surrounding water clusters allows
to recreate H-bonds. This should involve a decrease in the number of H-bonds
present in the clusters. These changes in the mass action capability of the
water clusters could eventually become dissipated through a cooling effect.


  A system of partial differential equations is derived as a model for the
dynamics of a honey bee colony with a continuous age distribution, and the
system is then extended to include the effects of a simplified infectious
disease. In the disease-free case we analytically derive the equilibrium age
distribution within the colony and propose a novel approach for determining the
global asymptotic stability of a reduced model. Furthermore, we present a
method for determining the basic reproduction number $R_0$ of the infection;
the method can be applied to other age-structured disease models with
interacting susceptible classes. The results of asymptotic stability indicate
that a honey bee colony suffering losses will recover naturally so long as the
cause of the losses is removed before the colony collapses. Our expression for
$R_0$ has potential uses in the tracking and control of an infectious disease
within a bee colony.


  Experiments have revealed that membrane proteins can form two-dimensional
clusters with regular translational and orientational protein arrangements,
which may allow cells to modulate protein function. However, the physical
mechanisms yielding supramolecular organization and collective function of
membrane proteins remain largely unknown. Here we show that bilayer-mediated
elastic interactions between membrane proteins can yield regular and
distinctive lattice architectures of protein clusters, and may provide a link
between lattice architecture and lattice function. Using the mechanosensitive
channel of large conductance (MscL) as a model system, we obtain relations
between the shape of MscL and the supramolecular architecture of MscL lattices.
We predict that the tetrameric and pentameric MscL symmetries observed in
previous structural studies yield distinct lattice architectures of MscL
clusters and that, in turn, these distinct MscL lattice architectures yield
distinct lattice activation barriers. Our results suggest general physical
mechanisms linking protein symmetry, the lattice architecture of membrane
protein clusters, and the collective function of membrane protein lattices.


  Hydrophobic thickness mismatch between integral membrane proteins and the
surrounding lipid bilayer can produce lipid bilayer thickness deformations.
Experiment and theory have shown that protein-induced lipid bilayer thickness
deformations can yield energetically favorable bilayer-mediated interactions
between integral membrane proteins, and large-scale organization of integral
membrane proteins into protein clusters in cell membranes. Within the continuum
elasticity theory of membranes, the energy cost of protein-induced bilayer
thickness deformations can be captured by considering compression and expansion
of the bilayer hydrophobic core, membrane tension, and bilayer bending,
resulting in biharmonic equilibrium equations describing the shape of lipid
bilayers for a given set of bilayer-protein boundary conditions. Here we
develop a combined analytic and numerical methodology for the solution of the
equilibrium elastic equations associated with protein-induced lipid bilayer
deformations. Our methodology allows accurate prediction of thickness-mediated
protein interactions for arbitrary protein symmetries at arbitrary protein
separations and relative orientations. We provide exact analytic solutions for
cylindrical integral membrane proteins with constant and varying hydrophobic
thickness, and develop perturbative analytic solutions for non-cylindrical
protein shapes. We complement these analytic solutions, and assess their
accuracy, by developing both finite element and finite difference numerical
solution schemes. Taken together, the work presented here puts into place an
analytic and numerical framework which allows calculation of bilayer-mediated
elastic interactions between integral membrane proteins for the complicated
protein shapes suggested by structural biology and at the small protein
separations most relevant for the crowded membrane environments provided by
living cells.


  In recent experiments [T. Basta et al., Proc. Natl. Acad. Sci. U.S.A. 111,
670 (2014)] lipids and membrane proteins were observed to self-assemble into
membrane protein polyhedral nanoparticles (MPPNs) with a well-defined
polyhedral protein arrangement and characteristic size. We develop a model of
MPPN self-assembly in which the preferred symmetry and size of MPPNs emerge
from the interplay of protein-induced lipid bilayer deformations, topological
defects in protein packing, and thermal effects. With all model parameters
determined directly from experiments, our model correctly predicts the observed
symmetry and size of MPPNs. Our model suggests how key lipid and protein
properties can be modified to produce a range of MPPN symmetries and sizes in
experiments.


  An active participation of players in evolutionary games depends on several
factors, ranging from personal stakes to the properties of the interaction
network. Diverse activity patterns thus have to be taken into account when
studying the evolution of cooperation in social dilemmas. Here we study the
weak prisoner's dilemma game, where the activity of each player is determined
in a probabilistic manner either by its degree or by its payoff. While
degree-correlated activity introduces cascading failures of cooperation that
are particularly severe on scale-free networks with frequently inactive hubs,
payoff-correlated activity provides a more nuanced activity profile, which
ultimately hinders systemic breakdowns of cooperation. To determine optimal
conditions for the evolution of cooperation, we introduce an exponential decay
to payoff-correlated activity that determines how fast the activity of a player
returns to its default state. We show that there exists an intermediate decay
rate, at which the resolution of the social dilemma is optimal. This can be
explained by the emerging activity patterns of players, where the inactivity of
hubs is compensated effectively by the increased activity of average-degree
players, who through their collective influence in the network sustain a higher
level of cooperation. The sudden drops in the fraction of cooperators observed
with degree-correlated activity therefore vanish, and so does the need for the
lengthy spatiotemporal reorganization of compact cooperative clusters. The
absence of such asymmetric dynamic instabilities thus leads to an optimal
resolution of social dilemmas, especially when the conditions for the evolution
of cooperation are strongly adverse.


  Mapping gene expression as a quantitative trait using whole genome-sequencing
and transcriptome analysis allows to discover the functional consequences of
genetic variation. We developed a novel method and ultra-fast software Findr
for higly accurate causal inference between gene expression traits using
cis-regulatory DNA variations as causal anchors, which improves current methods
by taking into account hidden confounders and weak regulations. Findr
outperformed existing methods on the DREAM5 Systems Genetics challenge and on
the prediction of microRNA and transcription factor targets in human
lymphoblastoid cells, while being nearly a million times faster. Findr is
publicly available at https://github.com/lingfeiwang/findr


  Consensus methods are widely used for combining phylogenetic trees into a
single estimate of the evolutionary tree for a group of species. As more taxa
are added, the new source trees may begin to tell a different evolutionary
story when restricted to the original set of taxa. However, if the new trees,
restricted to the original set of taxa, were to agree exactly with the earlier
trees, then we might hope that their consensus would either agree with or
resolve the original consensus tree. In this paper, we ask under what
conditions consensus methods exist that are 'future proof' in this sense. While
we show that some methods (e.g. Adams consensus) have this property for
specific types of input, we also establish a rather surprising `no-go' theorem:
there is no 'reasonable' consensus method that satisfies the future-proofing
property in general. We then investigate a second notion of 'future proofing'
for consensus methods, in which trees (rather than taxa) are added, and
establish some positive and negative results. We end with some questions for
future work.


  Summary: We introduce RBPBind, a web-based tool for the quantitative
prediction of RNA-protein interactions. Given a user-specified RNA and a
protein selected from a set of several common RNA-binding proteins, RBPBind
computes the binding curve and effective binding constant of the reaction in
question. The server also computes the probability that, at a given
concentration of protein, a protein molecule will bind to any particular
nucleotide along the RNA. The software for RBPBind is adapted from the Vienna
RNA package, whose source code has been modified to accommodate the effects of
single stranded RNA binding proteins. RBPBind thus fully incorporates the
effect of RNA secondary structure on protein-RNA interactions. Availability:
Our web server is available at http://bioserv.mps.ohio-state.edu/RBPBind


  We consider a continuous population whose dynamics is described by the
standard stationary Fleming-Viot process, so that the genealogy of $n$
uniformly sampled individuals is distributed as the Kingman $n$-coalescent. In
this note, we study some genealogical properties of this population when the
sample is conditioned to fall entirely into a subpopulation with most recent
common ancestor (MRCA) shorter than $\varepsilon$. First, using the comb
representation of the total genealogy (Lambert & Uribe Bravo 2016), we show
that the genealogy of the descendance of the MRCA of the sample on the
timescale $\varepsilon$ converges as $\varepsilon\to 0$. The limit is the
so-called Brownian coalescent point process (CPP) stopped at an independent
Gamma random variable with parameter $n$, which can be seen as the genealogy at
a large time of the total population of a rescaled critical birth-death
process, biased by the $n$-th power of its size. Secondly, we show that in this
limit the coalescence times of the $n$ sampled individuals are i.i.d. uniform
random variables in $(0,1)$. These results provide a coupling between two
standard models for the genealogy of a random exchangeable population: the
Kingman coalescent and the Brownian CPP.


  Recently developed deep learning techniques have significantly improved the
accuracy of various speech and image recognition systems. In this paper we
adapt some of these techniques for protein secondary structure prediction. We
first train a series of deep neural networks to predict eight-class secondary
structure labels given a protein's amino acid sequence information and find
that using recent methods for regularization, such as dropout and weight-norm
constraining, leads to measurable gains in accuracy. We then adapt recent
convolutional neural network architectures--Inception, ReSNet, and DenseNet
with Batch Normalization--to the problem of protein structure prediction. These
convolutional architectures make heavy use of multi-scale filter layers that
simultaneously compute features on several scales, and use residual connections
to prevent underfitting. Using a carefully modified version of these
architectures, we achieve state-of-the-art performance of 70.0% per amino acid
accuracy on the public CB513 benchmark dataset. Finally, we explore additions
from sequence-to-sequence learning, altering the model to make its predictions
conditioned on both the protein's amino acid sequence and its past secondary
structure labels. We introduce a new method of ensembling such a conditional
model with our convolutional model, an approach which reaches 70.6% Q8 accuracy
on CB513. We argue that these results can be further refined for larger boosts
in prediction accuracy through more sophisticated attempts to control
overfitting of conditional models. We aim to release the code for these
experiments as part of the TensorFlow repository.


  We consider the classic Moran process modeling the spread of genetic
mutations, as extended to structured populations by Lieberman et al.\ (Nature,
2005). In this process, individuals are the vertices of a connected graph $G$.
Initially, there is a single mutant vertex, chosen uniformly at random. In each
step, a random vertex is selected for reproduction with a probability
proportional to its fitness: mutants have fitness $r>1$, while non-mutants have
fitness 1. The vertex chosen to reproduce places a copy of itself to a
uniformly random neighbor in $G$, replacing the individual that was there. The
process ends when the mutation either reaches fixation (i.e., all vertices are
mutants), or gets extinct. The principal quantity of interest is the
probability with which each of the two outcomes occurs.
  A problem that has received significant attention recently concerns the
existence of families of graphs, called strong amplifiers of selection, for
which the fixation probability tends to 1 as the order $n$ of the graph
increases, and the existence of strong suppressors of selection, for which this
probability tends to 0. For the case of directed graphs, it is known that both
strong amplifiers and suppressors exist. For the case of undirected graphs,
however, the problem has remained open, and the general belief has been that
neither strong amplifiers nor suppressors exist. In this paper we disprove this
belief, by providing the first examples of such graphs. The strong amplifier we
present has fixation probability $1-\tilde O(n^{-1/3})$, and the strong
suppressor has fixation probability $\tilde O(n^{-1/4})$. Both graph
constructions are surprisingly simple. We also prove a general upper bound of
$1-\tilde \Omega(n^{-1/3})$ on the fixation probability of any undirected
graph. Hence, our strong amplifier is existentially optimal.


  Aging is associated with the accumulation of damage throughout a persons
life. Individual health can be assessed by the Frailty Index (FI). The FI is
calculated simply as the proportion $f$ of accumulated age related deficits
relative to the total, leading to a theoretical maximum of $f \leq 1$.
Observational studies have generally reported a much more stringent bound, with
$f \leq f_{max} <1$. The value of $f_{max}$ in observational studies appears to
be non-universal, but $f_{max} \approx 0.7$ is often reported. A previously
developed network model of individual aging was unable to recover $f_{max}<1$
while retaining the other observed phenomenology of increasing $f$ and
mortality rates with age. We have developed a computationally accelerated
network model that also allows us to tune the scale-free network exponent
$\alpha$. The network exponent $\alpha$ significantly affects the growth of
mortality rates with age. However, we are only able to recover $f_{max}$ by
also introducing a deficit sensitivity parameter $1-q$, which is equivalent to
a false-negative rate $q$. Our value of $q=0.3$ is comparable to finite
sensitivities of age-related deficits with respect to mortality that are often
reported in the literature. In light of non-zero $q$, we use mutual information
$I$ to provide a non-parametric measure of the predictive value of the FI with
respect to individual mortality. We find that $I$ is only modestly degraded by
$q<1$, and this degradation is mitigated when increasing number of deficits are
included in the FI. We also find that the information spectrum, i.e. the mutual
information of individual deficits vs connectivity, has an approximately
power-law dependence that depends on the network exponent $\alpha$. Mutual
information $I$ is therefore a useful tool for characterizing the network
topology of aging populations.


  A timely immunization can be effective against certain diseases and can save
thousands of lives. However, for some diseases it has been difficult, so far,
to develop an efficient vaccine. Malaria, a tropical disease caused by a
parasite of the genus Plasmodium, is one example. Bioinformatics has opened the
way to new lines of experimental investigation One example is reverse
vaccinology that aims to identify antigens that are capable of generating an
immune response in a given organism using in silico studies. In this study we
applied a reverse vaccinology methodology using a bioinformatics pipeline. We
obtained 45 potential linear B cells consensus epitopes from the whole genome
of P. falciparum 3D7 that can be used as candidates for malaria vaccines. The
direct implication of the results obtained is to open the way to experimentally
validate more epitopes to increase the efficiency of the available treatments
against malaria and to explore the methodology in other diseases.


  Selective sweeps affect neutral genetic diversity through hitchhiking. While
this effect is limited to the local genomic region of the sweep in panmictic
populations, we find that in spatially-extended populations the combined
effects of many unlinked sweeps can affect patterns of ancestry (and therefore
neutral genetic diversity) across the whole genome. Even low rates of sweeps
can be enough to skew the spatial locations of ancestors such that neutral
mutations that occur in an individual living outside a small region in the
center of the range have virtually no chance of fixing in the population.


  Samples of multiple complete genome sequences contain vast amounts of
information about the evolutionary history of populations, much of it in the
associations among polymorphisms at different loci. Current methods that take
advantage of this linkage information rely on models of recombination and
coalescence, limiting the sample sizes and populations that they can analyze.
We introduce a method, Minimal-Assumption Genomic Inference of Coalescence
(MAGIC), that reconstructs key features of the evolutionary history, including
the distribution of coalescence times, by integrating information across
genomic length scales without using an explicit model of recombination,
demography or selection. Using simulated data, we show that MAGIC's performance
is comparable to PSMC' on single diploid samples generated with standard
coalescent and recombination models. More importantly, MAGIC can also analyze
arbitrarily large samples and is robust to changes in the coalescent and
recombination processes. Using MAGIC, we show that the inferred coalescence
time histories of samples of multiple human genomes exhibit inconsistencies
with a description in terms of an effective population size based on
single-genome data.


  Intrinsically Disordered Proteins (IDPs) perform a broad range of biological
functions. Their relevance has motivated intense research activity seeking to
characterize their sequence/structure/function relationships. However, the
conformational plasticity of these molecules hampers the application of
traditional structural approaches, and new tools and concepts are being
developed to address the challenges they pose. Small-Angle Scattering (SAS) is
a structural biology technique that probes the size and shape of disordered
proteins and their complexes with other biomolecules. The low-resolution nature
of SAS can be compensated with specially designed computational tools and its
combined interpretation with complementary structural information. In this
review, we describe recent advances in the application of SAS to disordered
proteins and highly flexible complexes and discuss current challenges.


  We demonstrate that with two small modifications, the popular dielectric
continuum model is capable of predicting, with high accuracy, ion solvation
thermodynamics in numerous polar solvents, and ion solvation free energies in
water--co-solvent mixtures. The first modification involves perturbing the
macroscopic dielectric-flux interface condition at the solute--solvent
interface with a nonlinear function of the local electric field, giving what we
have called a solvation-layer interface condition (SLIC). The second
modification is a simple treatment of the microscopic interface potential
(static potential). We show that the resulting model exhibits high accuracy
without the need for fitting solute atom radii in a state-dependent fashion.
Compared to experimental results in nine water--co-solvent mixtures, SLIC
predicts transfer free energies to within 2.5 kJ/mol. The co-solvents include
both protic and aprotic species, as well as biologically relevant denaturants
such as urea and dimethylformamide. Furthermore, our results indicate that the
interface potential is essential to reproduce entropies and heat capacities.
The present work, together with previous studies of SLIC illustrating its
accuracy for biomolecules in water, indicates it as a promising dielectric
continuum model for accurate predictions of molecular solvation in a wide range
of conditions.


  We propose a mathematical model for cholera with treatment through
quarantine. The model is shown to be both epidemiologically and mathematically
well posed. In particular, we prove that all solutions of the model are
positive and bounded; and that every solution with initial conditions in a
certain meaningful set remains in that set for all time. The existence of
unique disease-free and endemic equilibrium points is proved and the basic
reproduction number is computed. Then, we study the local asymptotic stability
of these equilibrium points. An optimal control problem is proposed and
analyzed, whose goal is to obtain a successful treatment through quarantine. We
provide the optimal quarantine strategy for the minimization of the number of
infectious individuals and bacteria concentration, as well as the costs
associated with the quarantine. Finally, a numerical simulation of the cholera
outbreak in the Department of Artibonite (Haiti), in 2010, is carried out,
illustrating the usefulness of the model and its analysis.


  Decomposition of biomolecular reaction networks into pathways is a powerful
approach to the analysis of metabolic and signalling networks. Current
approaches based on analysis of the stoichiometric matrix reveal information
about steady-state mass flows (reaction rates) through the network. In this
work we show how pathway analysis of biomolecular networks can be extended
using an energy-based approach to provide information about energy flows
through the network. This energy-based approach is developed using the
engineering-inspired bond graph methodology to represent biomolecular reaction
networks.
  The approach is introduced using glycolysis as an exemplar; and is then
applied to analyse the efficiency of free energy transduction in a biomolecular
cycle model of a transporter protein (Sodium-Glucose Transport Protein 1,
SGLT1).
  The overall aim of our work is to present a framework for modelling and
analysis of biomolecular reactions and processes which considers energy flows
and losses as well as mass transport.


  The long-term evolution of epidemic processes depends crucially on the
structure of contact networks. As empirical evidence indicates that human
populations exhibit strong community organization, we investigate here how such
mesoscopic configurations affect the likelihood of epidemic recurrence. Through
numerical simulations on real social networks and theoretical arguments using
spectral methods, we demonstrate that highly contagious diseases that would
have otherwise died out rapidly can persist indefinitely for an optimal range
of modularity in contact networks.


  The construction of reciprocal relationships requires cooperative
interactions during the initial meetings. However, cooperative behavior with
strangers is risky because the strangers may be exploiters. In this study, we
show that people increase the likelihood of cooperativeness of strangers by
using lightweight non-risky interactions in risky situations based on the
analysis of a social network game (SNG). They can construct reciprocal
relationships in this manner. The interactions involve low-cost signaling
because they are not generated at any cost to the senders and recipients.
Theoretical studies show that low-cost signals are not guaranteed to be
reliable because the low-cost signals from senders can lie at any time.
However, people used low-cost signals to construct reciprocal relationships in
an SNG, which suggests the existence of mechanisms for generating reliable,
low-cost signals in human evolution.


  A system level view of cellular processes for human and several organisms can
be cap- tured by analyzing molecular interaction networks. A molecular
interaction network formed of differentially expressed genes and their
interactions helps to understand key players behind disease development. So, if
the functions of these genes are blocked by altering their interactions, it
would have a great impact in controlling the disease. Due to this promising
consequence, the problem of inferring disease causing genes and their pathways
has attained a crucial position in computational biology research. However,
considering the huge size of interaction networks, executing computations can
be costly. Review of literatures shows that the methods proposed for finding
the set of disease causing genes could be assessed in terms of their accuracy
which a perfect algorithm would find. Along with accuracy, the time complexity
of the method is also important, as high time complexities would limit the
number of pathways that could be found within a pragmatic time interval.


  The simplest epidemiologic model composed by mutually exclusive compartments
SIR (susceptible-infected-susceptible) is presented to describe a reality. From
health concerns to situations related with marketing, informatics or even
sociology, several are the fields that are using this epidemiological model as
a first approach to better understand a situation. In this paper, the basic
transmission model is analyzed, as well as simple tools that allows us to
extract a great deal of information about possible solutions. A set of
applications - traditional and new ones - is described to show the importance
of this model.


  The lactose operon in Escherichia coli was the first known gene regulatory
network, and it is frequently used as a prototype for new modeling paradigms.
Historically, many of these modeling frameworks use differential equations.
More recently, Stigler and Veliz-Cuba proposed a Boolean network model that
captures the bistability of the system and all of the biological steady states.
In this paper, we model the well-known arabinose operon in E. coli with a
Boolean network. This has several complex features not found in the lac operon,
such as a protein that is both an activator and repressor, a DNA looping
mechanism for gene repression, and the lack of inducer exclusion by glucose.
For 11 out of 12 choices of initial conditions, we use computational algebra
and Sage to verify that the state space contains a single fixed point that
correctly matches the biology. The final initial condition, medium levels of
arabinose and no glucose, successfully predicts the system's bistability.
Finally, we compare the state space under synchronous and asynchronous update,
and see that the former has several artificial cycles that go away under a
general asynchronous update.


  Rapidly mutating pathogens may be able to persist in the population and reach
an endemic equilibrium by escaping hosts' acquired immunity. For such diseases,
multiple biological, environmental and population-level mechanisms determine
the dynamics of the outbreak, including pathogen's epidemiological traits (e.g.
transmissibility, infectious period and duration of immunity), seasonality,
interaction with other circulating strains and hosts' mixing and spatial
fragmentation. Here, we study a susceptible-infected-recovered-susceptible
model on a metapopulation where individuals are distributed in subpopulations
connected via a network of mobility flows. Through extensive numerical
simulations, we explore the phase space of pathogen's persistence and map the
dynamical regimes of the pathogen following emergence. Our results show that
spatial fragmentation and mobility play a key role in the persistence of the
disease whose maximum is reached at intermediate mobility values. We describe
the occurrence of different phenomena including local extinction and emergence
of epidemic waves, and assess the conditions for large scale spreading.
Findings are highlighted in reference to previous works and to real scenarios.
Our work uncovers the crucial role of hosts' mobility on the ecological
dynamics of rapidly mutating pathogens, opening the path for further studies on
disease ecology in the presence of a complex and heterogeneous environment.


  The transitional and well-developed regimes of turbulent shear flows exhibit
a variety of remarkable scaling laws that are only now beginning to be
systematically studied and understood. In the first part of this article, we
summarize recent progress in understanding the friction factor of turbulent
flows in rough pipes and quasi-two-dimensional soap films, showing how the data
obey a two-parameter scaling law known as roughness-induced criticality, and
exhibit power-law scaling of friction factor with Reynolds number that depends
on the precise form of the nature of the turbulent cascade. These results hint
at a non-equilibrium fluctuation-dissipation relation that applies to turbulent
flows. The second part of this article concerns the lifetime statistics in
smooth pipes around the transition, showing how the remarkable
super-exponential scaling with Reynolds number reflects deep connections
between large deviation theory, extreme value statistics, directed percolation
and the onset of coexistence in predator-prey ecosystems. Both these phenomena
reflect the way in which turbulence can be fruitfully approached as a problem
in non-equilibrium statistical mechanics.


  Based on a detailed crossbridge model for individual myosin II motors, we
systematically study the influence of mechanical load and adenosine
triphosphate (ATP) concentration on small myosin II ensembles made from
different isoforms. For skeletal and smooth muscle myosin II, which are often
used in actomyosin gels that reconstitute cell contractility, fast forward
movement is restricted to a small region of phase space with low mechanical
load and high ATP concentration, which is also characterized by frequent
ensemble detachment. At high load, these ensembles are stalled or move
backwards, but forward motion can be restored by decreasing ATP concentration.
In contrast, small ensembles of nonmuscle myosin II isoforms, which are found
in the cytoskeleton of nonmuscle cells, are hardly affected by ATP
concentration due to the slow kinetics of the bound states. For all isoforms,
the thermodynamic efficiency of ensemble movement increases with decreasing ATP
concentration, but this effect is weaker for the nonmuscle myosin II isoforms.


  The ability of the adaptive immune system to respond to arbitrary pathogens
stems from the broad diversity of immune cell surface receptors (TCRs). This
diversity originates in a stochastic DNA editing process (VDJ recombination)
that acts each time a new immune cell is created from a stem cell. By analyzing
T cell sequence repertoires taken from the blood and thymus of mice of
different ages, we quantify the significant changes in this process that occur
in development from embryo to young adult. We find a rapid increase with age in
the number of random insertions in the VDJ recombination process, leading to a
dramatic increase in diversity. Since the blood accumulates thymic output over
time, blood repertoires are mixtures of different statistical recombination
processes and, by unraveling the mixture statistics, we can obtain a clear
picture of the time evolution of the early immune system. Sequence repertoire
analysis also allows us to detect the effect of selection on the output of the
VDJ recombination process. The effects we find are nearly identical between
thymus and blood, suggesting that they mainly reflect selection for proper
folding of the TCR receptor protein.


  In a previous article, an algorithm for identifying therapeutic targets in
Boolean networks modeling pathological mechanisms was introduced. In the
present article, the improvements made on this algorithm, named kali, are
described. These improvements are i) the possibility to work on asynchronous
Boolean networks, ii) a finer assessment of therapeutic targets and iii) the
possibility to use multivalued logic. kali assumes that the attractors of a
dynamical system, such as a Boolean network, are associated with the phenotypes
of the modeled biological system. Given a logic-based model of pathological
mechanisms, kali searches for therapeutic targets able to reduce the
reachability of the attractors associated with pathological phenotypes, thus
reducing their likeliness. kali is illustrated on an example network and used
on a biological case study. The case study is a published logic-based model of
bladder tumorigenesis from which kali returns consistent results. However, like
any computational tool, kali can predict but can not replace human expertise:
it is a supporting tool for coping with the complexity of biological systems in
the field of drug discovery.


  Cardiovascular diseases (CVD) represent a major health issue in patients with
schizophrneia (SCZ) and bipolar disorder (BD), but the exact nature of
cardiometabolic (CM) abnormalities involved and the underlying mechanisms
remain unclear. Using polygenic risk scores (PRS) and LD score regression, we
investigated the shared genetic bases of SCZ and BD with a panel of 28
cardiometabolic traits. We performed Mendelian randomization (MR) to elucidate
casual relationships between the two groups of disorders. The analysis was
based on large-scale meta-analyses of genome-wide association studies (GWAS).
We also identified the potential shared genetic variants by a statistical
approach based on local true discovery rates, and inferred the pathways
involved. We found polygenic associations of SCZ with glucose metabolism
abnormalities, adverse adipokine profiles, increased wait-hip ratio and raised
visceral adiposity. However, BMI showed inverse genetic correlation and
polygenic link with SCZ. On the other hand, we observed polygenic associations
with an overall favorable CM profile in BD. MR analysis showed that SCZ may be
causally linked to raised triglyceride and that lower fasting glucose may be
linked to BD; otherwise MR did not reveal other significant causal
relationships in general. We also identified numerous SNPs and pathways shared
between SCZ/BD with cardiometabolic traits, some of which are related to
inflammation or the immune system. In conclusion, SCZ patients may be
genetically associated with several CM abnormalities independent of medication
side-effects, and proper surveillance and management of CV risk factors may be
required from the onset of the disease. On the other hand, CM abnormalities in
BD are more likely to be secondary.


  Single particle reconstruction (SPR) from cryo-electron microscopy (EM) is a
technique in which the 3D structure of a molecule needs to be determined from
its contrast transfer function (CTF) affected, noisy 2D projection images taken
at unknown viewing directions. One of the main challenges in cryo-EM is the
typically low signal to noise ratio (SNR) of the acquired images. 2D
classification of images, followed by class averaging, improves the SNR of the
resulting averages, and is used for selecting particles from micrographs and
for inspecting the particle images. We introduce a new affinity measure, akin
to the Mahalanobis distance, to compare cryo-EM images belonging to different
defocus groups. The new similarity measure is employed to detect similar
images, thereby leading to an improved algorithm for class averaging. We
evaluate the performance of the proposed class averaging procedure on synthetic
datasets, obtaining state of the art classification.


  Using a series of computer simulations we have demonstrated a scenario of the
early evolution of the bird-type primitive language. We do not assume wise
agents who can use a grammar and manage an evolution without a grammar.
Duplication and mutation of phrases is our strategy. Such a strategy is seen in
wide classes of living phenomena.


  Estimating phylogenetic trees is an important problem in evolutionary
biology, environmental policy and medicine. Although trees are estimated, their
uncertainties are discarded by mathematicians working in tree space. Here we
explicitly model the multivariate uncertainty of tree estimates. We consider
both the cases where uncertainty information arises extrinsically (through
covariate information) and intrinsically (through the tree estimates
themselves). The importance of accounting for tree uncertainty in tree space is
demonstrated in two case studies. In the first instance, differences between
gene trees are small relative to their uncertainties, while in the second, the
differences are relatively large. Our main goal is visualization of tree
uncertainty, and we demonstrate advantages of our method with respect to
reproducibility, speed and preservation of topological differences compared to
visualization based on multidimensional scaling. The proposal highlights that
phylogenetic trees are estimated in an extremely high-dimensional space,
resulting in uncertainty information that cannot be discarded. Most
importantly, it is a method that allows biologists to diagnose whether
differences between gene trees are biologically meaningful, or due to
uncertainty in estimation.


  This chapter takes as its departure point a neural level theory of insight
that arose from studies of the sparse, distributed, content-addressable
architecture of associative memory. It is argued that convergent thought is
most fruitfully characterized in terms of, not the generation of a single
correct solution (as it is conventionally construed), but using concepts in
their most compact form by activating neural cell assemblies that respond to
their most typical properties. This allows them to be deployed in a
conventional manner such that effort is reserved for exploring causal
relationships. Conversely, it is argued that divergent thought is most
fruitfully characterized in terms of, not the generation of multiple solutions
(as it is conventionally construed), but using concepts in a form that is,
albeit expanded, constrained by the situation, by activating neural cell
assemblies that respond to context-specific atypical properties. This allows
them to be deployed in a manner that is conducive to exploring unconventional
yet potentially relevant associations, and unearthing potentially useful
relationships of correlation. Thus, divergent thought can involve as few as one
idea. This proposal is compatible with widespread beliefs that (1) most
creative tasks require not many solutions but one, yet entail both divergent
and convergent thinking, and (2) not all problems with multiple solutions
require creative thinking, and conversely, some problems with single solution
do require creative thought. The chapter discusses how the ability to shift
between convergent and divergent modes of thought may have evolved, and it
concludes with educational and vocational implications.


  In a comment on our manuscript "Strong selection significantly increases
epistatic interactions in the long-term evolution of a protein", Dr. Crona
challenges our assertion that shared entropy (that is, information) between two
residues implies epistasis between those residues, by constructing an explicit
example of three loci (say A, B, and C), where A and B are epistatically linked
(leading to shared entropy between A and B), and A and C also depend
epistatically (leading to shared entropy between A and C), so that loci B and C
are correlated (share entropy).


  Bacterial mobility is powered by rotation of helical flagellar filaments
driven by rotary motors. Flagellin isolated from {\it Salmonella Typhimurium}
SJW1660 strain, which differs by a point mutation from the wild-type strain,
assembles into straight filaments in which flagellin monomers are arranged into
left-handed helix. Using small-angle X-ray scattering (SAXS) and osmotic stress
methods, we investigated the high-resolution structure of SJW1660 flagellar
filaments as well as intermolecular forces that govern their assembly into
dense hexagonal bundle. The scattering data were fitted to high-resolution
models, which took into account the atomic structure of the flagellin subunits.
The analysis revealed the exact helical arrangement and the super-helical twist
of the flagellin subunits within the filaments. Under osmotic stress the
filaments formed $2D$ hexagonal bundles. Monte-Carlo simulations and continuum
theories were used to analyze the scattering data from hexagonal arrays,
revealing how bulk modulus, as well as how the deflection length depends on the
applied osmotic stress. Scattering data from aligned flagellar bundles
confirmed the predicated structure-factor scattering peak line-shape.
Quantitative analysis of the measured equation of state of the bundles revealed
the contributions of the electrostatic, hydration, and elastic interactions to
the intermolecular interactions associated with bundling of straight
semi-flexible flagellar filaments.}%1 {Insert Received for publication Date and
in final form Date.


  This paper presents a novel extension of the edge-based compartmental model
for epidemics with arbitrary distributions of transmission and recovery times.
Using the message passing approach we also derive a new pairwise-like model for
epidemics with Markovian transmission and an arbitrary recovery period. The new
pairwise-like model allows one to formally prove that the message passing and
edge-based compartmental models are equivalent in the case of Markovian
transmission and arbitrary recovery processes. The edge-based and message
passing models are conjectured to also be equivalent for arbitrary transmission
processes; we show the first step of a full proof of this. The new
pairwise-like model encompasses many existing well-known models that can be
obtained by appropriate reductions. It is also amenable to a relatively
straightforward numerical implementation. We test the theoretical results by
comparing the numerical solutions of the various pairwise-like models to
results based on explicit stochastic network simulations.


  This is a Reply to comments published in Physics of Life Reviews, on our
article "Evolutionary game theory using agent-based methods" (Physics of Life
Reviews, 2016, arXiv:1404.0994).


  Intrinsically disordered proteins (IDPs) and proteins with intrinsically
disordered regions (IDRs) govern a daunting number of physiological processes.
For such proteins, molecular mechanisms governing their interactions with
proteins involved in signal transduction pathways remain unclear. Using the
folded, calcium-loaded calmodulin (CaM) interaction with the calcineurin
regulatory IDP as a prototype for IDP-mediated signal transduction events, we
uncover the interplay of IDP structure and electrostatic interactions in
determining the kinetics of protein-protein association. Using an array of
biophysical approaches including stopped-flow and computational simulation, we
quantify the relative contributions of electrostatic interactions and
conformational ensemble characteristics in determining association kinetics of
calmodulin (CaM) and the calcineurin regulatory domain (CaN RD). Our chief
findings are that CaM/CN RD association rates are strongly dependent on ionic
strength and that observed rates are largely determined by the
electrostatically-driven interaction strength between CaM and the narrow CaN RD
calmodulin recognition domain. These studies highlight a molecular mechanism of
controlling signal transduction kinetics that may be utilized in wide-ranging
signaling cascades that involve IDPs.


  Dynamic instability of microtubules is considered using frameworks of
non-linear thermodynamics and non-equilibrium reaction-diffusion systems.
Stochastic assembly/disassembly phases in the polymerization dynamics of
microtubules are treated as a result of collective clusterization of
microdefects (holes in structure). The model explains experimentally observed
power law dependence of catastrophe frequency from the microtubule growth rate.
Additional reaction-diffusion-precipitation model is developed to account for
kinetic limitations in microtubule dynamics. It is shown that large scale
periodic microtubules length fluctuations are accompanied by concentration
autowaves. We built corresponding parametric diagram mapping areas of
stationary, non-stationary and metastable solutions. The loss of stability for
the stationary solutions happens through bifurcation of Andronov-Hopf. Using
parametric diagram we classify cytostatic effect of microtubule stabilizing
drugs in four major classes and analyze their compatibility and possible
synergistic effect for cancer treatment therapy.


  Genetic stability is a key factor in maintaining, survival and reproduction
of biological cells. It relies on many processes, but one of the most important
is a {\it homologous recombination}, in which the repair of breaks in
double-stranded DNA molecules is taking place with a help of several specific
proteins. In bacteria this task is accomplished by RecA proteins that are
active as nucleoprotein filaments formed on single-stranded segments of DNA. A
critical step in the homologous recombination is a search for a corresponding
homologous region on DNA, which is called a {\it homology search}. Recent
single-molecule experiments clarified some aspects of this process, but its
molecular mechanisms remain not well understood. We developed a quantitative
theoretical approach to analyze the homology search. It is based on a
discrete-state stochastic model that takes into account the most relevant
physical-chemical processes in the system. Using a method of first-passage
processes, a full dynamic description of the homology search is presented. It
is found that the search dynamics depends on the degree of extension of DNA
molecules and on the size of RecA nucleoprotein filaments, in agreement with
experimental single-molecule measurements of DNA pairing by RecA proteins. Our
theoretical calculations, supported by extensive Monte Carlo computer
simulations, provide a molecular description of the mechanisms of the homology
search.


  We introduce and analyze coupled, multi-strain epidemic models designed to
simulate the emergence and dissemination of mutant (e.g. drug-resistant)
pathogen strains. In particular, we investigate the mathematical and biological
properties of a general class of multi-strain epidemic models in which the
infectious compartments of each strain are coupled together in a general
manner. We derive explicit expressions for the basic reproduction number of
each strain and highlight their importance in regulating the system dynamics
(e.g. the potential for an epidemic outbreak) and the existence of nonnegative
endemic solutions. Importantly, we find that the basic reproduction number of
each strain is independent of the mutation rates between the strains --- even
under quite general assumptions for the form of the infectious compartment
coupling. Moreover, we verify that the coupling term promotes strain
coexistence (as an extension of the competitive exclusion principle) and
demonstrate that the strain with the greatest reproductive capacity is not
necessarily the most prevalent. Finally, we briefly discuss the implications of
our results for public health policy and planning.


  We study the Moran process as adapted by Lieberman, Hauert and Nowak. This is
a model of an evolving population on a graph or digraph where certain
individuals, called "mutants" have fitness r and other individuals, called
non-mutants have fitness 1. We focus on the situation where the mutation is
advantageous, in the sense that r>1. A family of digraphs is said to be
strongly amplifying if the extinction probability tends to 0 when the Moran
process is run on digraphs in this family. The most-amplifying known family of
digraphs is the family of megastars of Galanis et al. We show that this family
is optimal, up to logarithmic factors, since every strongly-connected n-vertex
digraph has extinction probability Omega(n^(-1/2)). Next, we show that there is
an infinite family of undirected graphs, called dense incubators, whose
extinction probability is O(n^(-1/3)). We show that this is optimal, up to
constant factors. Finally, we introduce sparse incubators, for varying edge
density, and show that the extinction probability of these graphs is O(n/m),
where m is the number of edges. Again, we show that this is optimal, up to
constant factors.


  Engineering systems modelling and analysis based on the bond graph approach
has been applied to biomolecular systems. In this context, the notion of a
Faraday-equivalent chemical potential is introduced which allows chemical
potential to be expressed in an analogous manner to electrical volts thus
allowing engineering intuition to be applied to biomolecular systems. Redox
reactions, and their representation by half-reactions, are key components of
biological systems which involve both electrical and chemical domains. A bond
graph interpretation of redox reactions is given which combines bond graphs
with the Faraday- equivalent chemical potential. This approach is particularly
relevant when the biomolecular system implements chemoelectrical transduction -
for example chemiosmosis within the key metabolic pathway of mitochondria:
oxidative phosphorylation. An alternative way of implementing computational
modularity using bond graphs is introduced and used to give a physically based
model of the mitochondrial electron transport chain. To illustrate the overall
approach, this model is analysed using the Faraday-equivalent chemical
potential approach and engineering intuition is used to guide affinity
equalisation: a energy based analysis of the mitochondrial electron transport
chain.


  Single type and multitype branching process have been used to study the
dynamics of a variety of stochastic birth-death type phenomena in biology and
physics. Their use in epidemiology goes back to Whittle's study of a
Susceptible--Infected--Recovered (SIR) model in the 1950s. In the case of an
SIR model, the presence of only one infectious class allows for the use of
single type branching processes. Multitype branching processes allow for
multiple infectious classes and have latterly been used to study metapopulation
models of disease. In this article, we develop a Continuous Time Markov Chain
(CTMC) model of Infectious Salmon Anemia virus in two patches, two CTMC models
in one patch and companion multitype branching process (MTBP) models. The CTMC
models are related to deterministic models which inform the choice of
parameters. The probability of extinction is computed for the CTMC via
numerical methods and approximated by the MTBP in the supercritical regime. The
stochastic models are treated as toy models and the parameter choices are made
to highlight regions of the parameter space where CTMC and MTBP agree or
disagree, without regard to biological significance. Partial extinction events
are defined and their relevance discussed. A case is made for calculating the
probability of such events, noting that MTBP's are not suitable for making
these calculations.


  The molecular mechanism of ion channel gating and substrate modulation is
elusive for many voltage gated ion channels, such as eukaryotic sodium ones.
The understanding of channel functions is a pressing issue in molecular
biophysics and biology. Mathematical modeling, computation and analysis of
membrane channel charge transport have become an emergent field and give rise
to significant contributions to our understanding of ion channel gating and
function. This review summarizes recent progresses and outlines remaining
challenges in mathematical modeling, simulation and analysis of ion channel
charge transport. One of our focuses is the Poisson-Nernst-Planck (PNP) model
and its generalizations. Specifically, the basic framework of the PNP system
and some of its extensions, including size effects, ion-water interactions,
coupling with density functional theory and relation to fluid flow models. A
reduced theory, the Poisson- Boltzmann-Nernst-Planck (PBNP) model, and a
differential geometry based ion transport model are also discussed. For proton
channel, a multiscale and multiphysics Poisson-Boltzmann-Kohn-Sham (PBKS) model
is presented. We show that all of these ion channel models can be cast into a
unified variational multiscale framework with a macroscopic continuum domain of
the solvent and a microscopic discrete domain of the solute. The main strategy
is to construct a total energy functional of a charge transport system to
encompass the polar and nonpolar free energies of solvation and chemical
potential related energies. Current computational algorithms and tools for
numerical simulations and results from mathematical analysis of ion channel
systems are also surveyed.


  Asymptomatic individuals in the context of malarial disease refers to
subjects who carry a parasite load but do not show clinical symptoms. A correct
understanding of the influence of asymptomatic individuals on transmission
dynamics will provide a comprehensive description of the complex interplay
between the definitive host (female \textit{Anopheles} mosquito), intermediate
host (human) and agent (\textit{Plasmodium} parasite). The goal of this article
is to conduct a rigorous mathematical analysis of a new compartmentalized
malaria model accounting for asymptomatic human hosts for the purpose of
calculating the basic reproductive number ($\mathcal{R}_0$), and determining
the bifurcations that might occur at the onset of disease free equilibrium. A
point of departure of this model from others appearing in literature is that
the asymptomatic compartment is decomposed into two mutually disjoint
sub-compartments by making use of the naturally acquired immunity (NAI) of the
population under consideration. After deriving the model, a qualitative
analysis is carried out to classify the stability of the equilibria of the
system. Our results show that the dynamical system is locally asymptotically
stable provided that $\mathcal{R}_0<1$. However this stability is not global,
owning to the occurrence of a sub-critical bifurcation in which additional
non-trivial sub-threshold equilibrium solutions appear in response to a
specified parameter being perturbed. To ensure that the model does not undergo
a backward bifurcation, we demand that an auxiliary parameter denoted
$\Lambda<1$ in addition to the threshold constraint $\mathcal{R}_0<1$. The
authors hope that this qualitative analysis will fill in the gaps of what is
currently known about asymptomatic malaria and aid in designing strategies that
assist the further development of malaria control and eradication efforts.


  Several modern genomic technologies, such as DNA-Methylation arrays, measure
spatially registered probes that number in the hundreds of thousands across
multiplechromosomes. The measured probes are by themselves less interesting
scientifically; instead scientists seek to discover biologically interpretable
genomic regions comprised of contiguous groups of probes which may act as
biomarkers of disease or serve as a dimension-reducing pre-processing step for
downstream analyses. In this paper, we introduce an unsupervised feature
learning technique which maps technological units (probes) to biological units
(genomic regions) that are common across all subjects. We use ideas from fusion
penalties and convex clustering to introduce a method for Spatial Convex
Clustering, or SpaCC. Our method is specifically tailored to detecting
multi-subject regions of methylation, but we also test our approach on the
well-studied problem of detecting segments of copy number variation. We
formulate our method as a convex optimization problem, develop a massively
parallelizable algorithm to find its solution, and introduce automated
approaches for handling missing values and determining tuning parameters.
Through simulation studies based on real methylation and copy number variation
data, we show that SpaCC exhibits significant performance gains relative to
existing methods. Finally, we illustrate SpaCC's advantages as a pre-processing
technique that reduces large-scale genomics data into a smaller number of
genomic regions through several cancer epigenetics case studies on subtype
discovery, network estimation, and epigenetic-wide association.


  Reproduction and natural selection are the key elements of life. In order to
reproduce, the genetic material must be doubled, separated and placed into two
new daughter cells, each containing a complete set of chromosomes and
organelles. In mitosis, transition from one process to the next is guided by
intricate surveillance mechanisms, known as the mitotic checkpoints.
Dis-regulation of cell division through checkpoint malfunction can lead to
developmental defects and contribute to the development or progression of
tumors. This review approaches two important mitotic checkpoints, the spindle
assembly checkpoint (SAC) and the spindle position checkpoint (SPOC). The
highly conserved spindle assembly checkpoint (SAC) controls the onset of
anaphase by preventing premature segregation of the sister chromatids of the
duplicated genome, to the spindle poles. In contrast, the spindle position
checkpoint (SPOC), in the budding yeast S. cerevisiae, ensures that during
asymmetric cell division mitotic exit does not occur until the spindle is
properly aligned with the cell polarity axis. Although there are no known
homologs, there is indication that functionally similar checkpoints exist also
in animal cells.


  The highly conserved spindle assembly checkpoint (SAC) ensures that the
sister chromatids of the duplicated genome are not separated and distributed to
the spindle poles before all chromosomes have been properly linked to the
microtubules of the mitotic spindle. Biochemically, the SAC delays cell cycle
progression by preventing activation of the anaphase-promoting complex (APC/C)
or cyclosome; whose activation by Cdc20 is required for sister-chromatid
separation, which marks the transition into anaphase. In response to activation
of the checkpoint, various species control the activity of both APC/C and
Cdc20. However, the underlying regulatory pathways remain largely elusive. In
this study, five possible model variants of APC/C regulation were constructed,
namely BubR1, Mad2, MCC, MCF2 and an all-pathways model variant. These models
are validated with experimental data from the literature. A wide range of
parameter values have been tested to find critical values of the APC binding
rate. The results show that all variants are able to capture the wild type
behaviour of the APC. However, only one model variant, which included both MCC
as well as BubR1 as potent inhibitors of the APC, was able to reproduce both
wild type and mutant type behaviour of APC regulation. The presented work has
successfully distinguished between five competing dynamical models of the same
biological system using a systems biology approach. Furthermore, the results
suggest that systems-level approach is vital for molecular biology and could
also be used for compare the pathways of relevance with the objective to
generate hypotheses and improve our understanding.


  Sexual partnerships that overlap in time (concurrent relationships) may play
a significant role in the HIV epidemic, but the precise effect is unclear. We
derive edge-based compartmental models of disease spread in idealized dynamic
populations with and without concurrency to allow for an investigation of its
effects. Our models assume that partnerships change in time and individuals
enter and leave the at-risk population. Infected individuals transmit at a
constant per-partnership rate to their susceptible partners. In our idealized
populations we find regions of parameter space where the existence of
concurrent partnerships leads to substantially faster growth and higher
equilibrium levels, but also regions in which the existence of concurrent
partnerships has very little impact on the growth or the equilibrium.
Additionally we find mixed regimes in which concurrency significantly increases
the early growth, but has little effect on the ultimate equilibrium level.
Guided by model predictions, we discuss general conditions under which
concurrent relationships would be expected to have large or small effects in
real-world settings. Our observation that the impact of concurrency saturates
suggests that concurrency-reducing interventions may be most effective in
populations with low to moderate concurrency.


  In this article, we investigate the evolution of sexual diploid populations
which are hosts for active TE families. Our purpose is to explore the
relationship between the environmental change, that influences such population
and activity of those TEs that are present in genomes.


  This paper is concerned with the analysis of vaccination strategies in a
stochastic SIR (susceptible $\to$ infected $\to$ removed) model for the spread
of an epidemic amongst a population of individuals with a random network of
social contacts that is also partitioned into households. Under various vaccine
action models, we consider both household-based vaccination schemes, in which
the way in which individuals are chosen for vaccination depends on the size of
the households in which they reside, and acquaintance vaccination, which
targets individuals of high degree in the social network. For both types of
vaccination scheme, assuming a large population with few initial infectives, we
derive a threshold parameter which determines whether or not a large outbreak
can occur and also the probability and fraction of the population infected by
such an outbreak. The performance of these schemes is studied numerically,
focusing on the influence of the household size distribution and the degree
distribution of the social network. We find that acquaintance vaccination can
significantly outperform the best household-based scheme if the degree
distribution of the social network is heavy-tailed. For household-based
schemes, when the vaccine coverage is insufficient to prevent a major outbreak
and the vaccine is imperfect, we find situations in which both the probability
and size of a major outbreak under the scheme which minimises the threshold
parameter are \emph{larger} than in the scheme which maximises the threshold
parameter.


  We recall the perturbation expansion for Michaelis-Menten kinetics, beyond
the standard quasi-steady-state approximation (sQSSA). Against this background,
we are able to appropriately apply the alternative approach to the study of
singularly perturbed differential equations that is based on the
renormalization group (SPDERG), by clarifying similarities and differences. In
the present demanding situation, we directly renormalize the bare initial
condition value for the substrate. Our main results are: i) the 2nd order
SPDERG uniform approximations to the correct solutions contain, up to 1st
order, the same outer components as the known perturbation expansion ones; ii)
the differential equation to be solved for the derivation of the 1st order
outer substrate component is simpler within the SPDERG approach; iii) the
approximations better reproduce the numerical solutions of the original problem
in a region encompassing the matching one, because of the 2nd order terms in
the inner components, calculated here for the first time to our knowledge: iv)
the refined SPDERG uniform approximations, that we propose, give the correct
asymptotically vanishing solutions, too, and allow to obtain results nearly
indistinguishable from the solutions of the original problem in a large part of
the whole relevant time window, even in the studied unfavourable kinetic
constant case, for an expansion parameter value as large as {\epsilon} = 0.5.


  Inverse statistical approaches to determine protein structure and function
from Multiple Sequence Alignments (MSA) are emerging as powerful tools in
computational biology. However the underlying assumptions of the relationship
between the inferred effective Potts Hamiltonian and real protein structure and
energetics remain untested so far. Here we use lattice protein model (LP) to
benchmark those inverse statistical approaches. We build MSA of highly stable
sequences in target LP structures, and infer the effective pairwise Potts
Hamiltonians from those MSA. We find that inferred Potts Hamiltonians reproduce
many important aspects of 'true' LP structures and energetics. Careful analysis
reveals that effective pairwise couplings in inferred Potts Hamiltonians depend
not only on the energetics of the native structure but also on competing folds;
in particular, the coupling values reflect both positive design (stabilization
of native conformation) and negative design (destabilization of competing
folds). In addition to providing detailed structural information, the inferred
Potts models used as protein Hamiltonian for design of new sequences are able
to generate with high probability completely new sequences with the desired
folds, which is not possible using independent-site models. Those are
remarkable results as the effective LP Hamiltonians used to generate MSA are
not simple pairwise models due to the competition between the folds. Our
findings elucidate the reasons for the success of inverse approaches to the
modelling of proteins from sequence data, and their limitations.


  Here we study ODE epidemic models with spread of awareness, assuming that a
certain proportion of the hosts will become aware of the ongoing outbreak upon
recovery. This study builds on W. Just and J. Saldana's work in [1], and is
conducted under the same framework, while addressing the influence of the
awareness gained from direct experience of the disease.
  In [1], the authors investigated the question whether preventive behavioral
response triggered by awareness of the infection is sufficient to prevent
future flare-ups from low endemic levels if awareness decays over time. They
showed that if all the hosts experienced infection return directly to the
susceptible compartment upon recovery, such oscillations are ruled out in
Susceptible-Aware-Infectious Susceptible models with a single compartment of
aware hosts, but can occur if two distinct compartments of aware hosts who
differ in their willingness to alert other susceptible hosts are considered.
Qualitatively, the models studied here produce the same results when we assume
that recovery from the disease may or even will convey awareness from direct
experience.


  A pattern Recognition of a probability distribution of amino acids is
obtained for selected families of proteins. The mathematical model is derived
from a theory of protein families formation which is derived from application
of a Pauli's master equation method.


  We study the tandem duplication distance between binary sequences and their
roots. In other words, the quantity of interest is the number of tandem
duplication operations of the form $\seq x = \seq a \seq b \seq c \to \seq y =
\seq a \seq b \seq b \seq c$, where $\seq x$ and $\seq y$ are sequences and
$\seq a$, $\seq b$, and $\seq c$ are their substrings, needed to generate a
binary sequence of length $n$ starting from a square-free sequence from the set
$\{0,1,01,10,010,101\}$. This problem is a restricted case of finding the
duplication/deduplication distance between two sequences, defined as the
minimum number of duplication and deduplication operations required to
transform one sequence to the other. We consider both exact and approximate
tandem duplications. For exact duplication, denoting the maximum distance to
the root of a sequence of length $n$ by $f(n)$, we prove that $f(n)=\Theta(n)$.
For the case of approximate duplication, where a $\beta$-fraction of symbols
may be duplicated incorrectly, we show that the maximum distance has a sharp
transition from linear in $n$ to logarithmic at $\beta=1/2$. We also study the
duplication distance to the root for sequences with a given root and for
special classes of sequences, namely, the de Bruijn sequences, the Thue-Morse
sequence, and the Fibbonaci words. The problem is motivated by genomic tandem
duplication mutations and the smallest number of tandem duplication events
required to generate a given biological sequence.


  In this work we analyze the post monsoon Dengue outbreaks by analyzing the
transient and long term dynamics of Dengue incidences and its environmental
correlates in Ahmedabad city in western India from 2005-2012. We calculate the
reproduction number $R_p$ using the growth rate of post monsoon Dengue
outbreaks and biological parameters like host and vector incubation periods and
vector mortality rate, and its uncertainties are estimated through Monte-Carlo
simulations by sampling parameters from their respective probability
distributions. Reduction in Female Aedes mosquito density required for an
effective prevention of Dengue outbreaks is also calculated. The non stationary
pattern of Dengue incidences and its climatic correlates like rainfall
temperature is analyzed through Wavelet based methods. We find that the mean
time lag between peak of monsoon and Dengue is 9 weeks. Monsoon and Dengue
cases are phase locked from 2008-2012 in the 16-32 weeks band. The duration of
post monsoon outbreak has been increasing every year, especially post 2008,
even though the intensity and duration of monsoon has been decreasing.
Temperature and Dengue incidences show correlations in the same band, but phase
lock is not stationary.


  Chemical reaction networks with generalized mass-action kinetics lead to
power-law dynamical systems. As a simple example, we consider the Lotka
reactions with two chemical species and arbitrary power-law kinetics. We study
existence, uniqueness, and stability of the positive equilibrium, in
particular, we characterize its global asymptotic stability in terms of the
kinetic orders.


  We introduce the package PhylogeneticTrees for Macaulay2 which allows users
to compute phylogenetic invariants for group-based tree models. We provide some
background information on phylogenetic algebraic geometry and show how the
package PhylogeneticTrees can be used to calculate a generating set for a
phylogenetic ideal as well as a lower bound for its dimension. Finally, we show
how methods within the package can be used to compute a generating set for the
join of any two ideals.


  Specific interactions between receptors and their target ligands in the
presence of non-target ligands are crucial for biological processes such as T
cell ligand discrimination. To discriminate between the target and non-target
ligands, cells have to increase specificity by amplifying the small differences
in affinity among ligands. In addition, sensitivity to ligand concentration and
quick discrimination are also important to detect low amounts of target ligands
and facilitate fast cellular decision making after ligand recognition. In this
work, we find that ultraspecificity is naturally derived from a well-known
mechanism for zero-order ultrasensitivity to concentration. We also show that
this mechanism can produce an optimal balance of specificity, sensitivity, and
quick discrimination. Furthermore, we show that a model for insensitivity to
large number of non-terget ligands can be naturally derived from the
ultraspecificity model. Zero-order ultraspecificity may provide alternative way
to understand ligand discrimination from the viewpoint of the nonlinear
properties of biochemical reactions.


  Global, population-wide oscillations in models of cyclic dominance may result
in the collapse of biodiversity due to the accidental extinction of one species
in the loop. Previous research has shown that such oscillations can emerge if
the interaction network has small-world properties, and more generally, because
of long-range interactions among individuals or because of mobility. But
although these features are all common in nature, global oscillations are
rarely observed in actual biological systems. This begets the question what is
the missing ingredient that would prevent local oscillations to synchronize
across the population to form global oscillations. Here we show that, although
heterogeneous species-specific invasion rates fail to have a noticeable impact
on species coexistence, randomness in site-specific invasion rates successfully
hinders the emergence of global oscillations and thus preserves biodiversity.
Our model takes into account that the environment is often not uniform but
rather spatially heterogeneous, which may influence the success of microscopic
dynamics locally. This prevents the synchronization of locally emerging
oscillations, and ultimately results in a phenomenon where one type of
randomness is used to mitigate the adverse effects of other types of randomness
in the system.


  In human relations individuals' gender and age play a key role in the
structures and dynamics of their social arrangements. In order to analyze the
gender preferences of individuals in interaction with others at different
stages of their lives we study a large mobile phone dataset. To do this we
consider four fundamental gender-related caller and callee combinations of
human interactions, namely male to male, male to female, female to male, and
female to female, which together with age, kinship, and different levels of
friendship give rise to a wide scope of human sociality. Here we analyse the
relative strength of these four types of interaction using a large dataset of
mobile phone communication records. Our analysis suggests strong age dependence
for an ego of one gender choosing to call an individual of either gender. We
observe a strong opposite sex bonding across most of their reproductive age.
However, older women show a strong tendency to connect to another female that
is one generation younger in a way that is suggestive of the
\emph{grandmothering effect}. We also find that the relative strength among the
four possible interactions depends on phone call duration. For calls of medium
and long duration, opposite gender interactions are significantly more probable
than same gender interactions during the reproductive years, suggesting
potential emotional exchange between spouses. By measuring the fraction of
calls to other generations we find that mothers tend to make calls more to
their daughters than to their sons, whereas fathers make calls more to their
sons than to their daughters. For younger people, most of their calls go to
same generation alters, while older people call the younger people more
frequently, which supports the suggestion that \emph{affection flows downward}.


  Sensors are the first element of the pathways that control the response of
cells to their environment. After chemical, the next most important cue is
mechanical, and protein complexes that produce or enable a chemical signal in
response to a mechanical stimulus are called mechanosensors. There is a sharp
distinction between sensing an external force or pressure/tension applied to
the cell, and sensing the mechanical stiffness of the environment. We call the
first mechanosensitivity of the 1st kind, and the latter mechanosensitivity of
the 2nd kind. There are two variants of protein complexes that act as
mechanosensors of the 2nd kind: producing either a one-off or a reversible
action. The latent complex of TGF-$\beta$ is an example of the one-off action:
on the release of active TGF-$\beta$ signal, the complex is discarded and needs
to be replaced. In contrast, focal adhesion kinase (FAK) in a complex with
integrin is a reversible mechanosensor, which initiates the chemical signal in
its active phosphorylated conformation, but can spontaneously return to its
closed folded conformation. Here we study the physical mechanism of the
reversible mechanosensor of the 2nd kind, using FAK as a practical example. We
find how the rates of conformation changes depend on the substrate stiffness
and the pulling force applied from the cell cytoskeleton. The results compare
well with the phenotype observations of cells on different substrates.


  In this paper we propose a straightforward operational definition of variants
of disordered proteins, taking the human proteome as a case study. The focus is
on a distinction between mostly unstructured proteins and proteins which
contain long unstructured regions accommodated in an overall folded structure.
In particular we distinguish: i) Not disordered proteins (NDPs), that either
have all their residues ordered or do not have disordered segments longer than
30 residues nor more than 30% of disordered residues; ii) Proteins with
intrinsically disordered regions (IDRPs), that have at least one disordered
domain longer than 30 residues, but disordered in less than 30% of their
residues; iii) Proteins that are intrinsically disordered (IDPs), that have
both at least one disordered segment longer than 30 residues and that are
disordered on more than 30% of their residues; iv) Proteins with fragmented
disorder (FRAG_IDPs), that do not have a disordered fragment longer than 30
residues but that, nevertheless, have at least 30% or more of their residues
predicted as disordered. The potential use of these variants is checked over
several groups of disease-related proteins. Our main conclusions point out that
IDRPs are more similar to NDPs than to IDPs. IDRPs and NDPs have a similar
functional repertoire and probably share a lock-and-key mechanism of
interaction with substrates. IDRPs and IDPs are differently present among human
disease-related proteins. IDRPs probably do not play a specific role in the
development of complex diseases, since their frequency is similar in
disease-related proteins and in the entire human proteome. IDPs can play a role
in the emergence of cancer, neurodegenerative, thyroid and liver disease.


  This paper focuses on pattern matching in the DNA sequence. It was inspired
by a previously reported method that proposes encoding both pattern and
sequence using prime numbers. Although fast, the method is limited to rather
small pattern lengths, due to computing precision problem. Our approach
successfully deals with large patterns, due to our implementation that uses
modular arithmetic. In order to get the results very fast, the code was adapted
for multithreading and parallel implementations. The method is reduced to
assembly language level instructions, thus the final result shows significant
time and memory savings compared to the reference algorithm.


  The Wright-Fisher model is the most popular population model for describing
the behaviour of evolutionary systems with a finite population size.
Approximations to the model have commonly been used for the analysis of
time-resolved genome sequence data, but the model itself has rarely been tested
against genomic data. Here, we evaluate the extent to which it can be inferred
as the correct model given experimental data. Given genome-wide data from an
evolutionary experiment, we validate the Wright-Fisher model as the better
model for variance in a finite population in contrast to a Gaussian model of
allele frequency propagation. However, we note a range of circumstances under
which the Wright-Fisher model cannot be correctly identified. We discuss the
potential for more rapid approximations to the Wright-Fisher model.


  Recent computational efforts have shown that the current potential energy
models used in molecular dynamics are not accurate enough to describe the
conformational ensemble of RNA oligomers and suggest that molecular dynamics
should be complemented with experimental data. We here propose a scheme based
on the maximum entropy principle to combine simulations with bulk experiments.
In the proposed scheme the noise arising from both the measurements and the
forward models used to back calculate the experimental observables is
explicitly taken into account. The method is tested on RNA nucleosides and is
then used to construct chemically consistent corrections to the Amber RNA force
field that allow a large set of experimental data on nucleosides and
dinucleosides to be correctly reproduced. The transferability of these
corrections is assessed against independent data on tetranucleotides and
displays a previously unreported agreement with experiments. This procedure can
be applied to enforce multiple experimental data on multiple systems in a
self-consistent framework thus suggesting a new paradigm for force field
refinement.


  Long, flexible physical filaments are naturally tangled and knotted, from
macroscopic string down to long-chain molecules. The existence of knotting in a
filament naturally affects its configuration and properties, and may be very
stable or disappear rapidly under manipulation and interaction. Knotting has
been previously identified in protein backbone chains, for which these
mechanical constraints are of fundamental importance to their molecular
functionality, despite their being open curves in which the knots are not
mathematically well defined; knotting can only be identified by closing the
termini of the chain somehow. We introduce a new method for resolving knotting
in open curves using virtual knots, a wider class of topological objects that
do not require a classical closure and so naturally capture the topological
ambiguity inherent in open curves. We describe the results of analysing
proteins in the Protein Data Bank by this new scheme, recovering and extending
previous knotting results, and identifying topological interest in some new
cases. The statistics of virtual knots in protein chains are compared with
those of open random walks and Hamiltonian subchains on cubic lattices,
identifying a regime of open curves in which the virtual knotting description
is likely to be important.


  In 3-D reconstruction problems, the image data obtained from cryo electron
microscopy is the projection of many heterogeneous instances of the object
under study (e.g., a virus). When the object is heterogeneous but has an
overall symmetry, it is natural to describe the object as stochastic with
symmetrical statistics. This paper presents a maximum likelihood reconstruction
approach which allows each object to lack symmetry while constraining the {\it
statistics} of the ensemble of objects to have symmetry. This algorithm is
demonstrated on bacteriophage HK97 and is contrasted with an existing algorithm
in which each object, while still heterogeneous, has the symmetry.
Reconstruction results show that the proposed algorithm eliminates
long-standing distortions in previous heterogeneity calculations associated
with symmetry axes, and provides estimates that make more biologically sense
than the estimates of existing algorithms.


  This paper is about numerical control of HIV propagation. The contribution of
the paper is threefold: first, a novel model of HIV propagation is proposed;
second, the methods from numerical optimal control are successfully applied to
the developed model to compute optimal control profiles; finally, the computed
results are applied to the real problem yielding important and practically
relevant results.


  Methicillin-resistant Staphylococcus aureus (MRSA) is a difficult-to-treat
infection that only in the European Union affects about 150,000 patients and
causes extra costs of 380 million Euros annually to the health-care systems.
Increasing efforts have been taken to mitigate the epidemics and to avoid
potential outbreaks in low endemic settings. Understanding the population
dynamics of MRSA through modeling is essential to identify the causal
mechanisms driving the epidemics and to generalize conclusions to different
contexts. We develop an innovative high-resolution spatiotemporal contact
network model of interactions between patients to reproduce the hospital
population in the context of the Stockholm County in Sweden and simulate the
spread of MRSA within this population. Our model captures the spatial and
temporal heterogeneity caused by human behavior and by the dynamics of mobility
within wards and hospitals. We estimate that in this population the epidemic
threshold is at about 0.008. We also identify that these heterogeneous contact
patterns cause the emergence of super-spreader patients and a polynomial growth
of the epidemic curve. We finally study the effect of standard intervention
control strategies and identify that screening is more effective than improved
hygienic in order to cause smaller or null outbreaks.


  The ability of a honeybee swarm to select the best nest site plays a
fundamental role in determining the future colony's fitness. To date, the
nest-site selection process has mostly been modelled and theoretically analysed
for the case of binary decisions. However, when the number of alternative nests
is larger than two, the decision process dynamics qualitatively change. In this
work, we extend previous analyses of a value-sensitive decision-making
mechanism to a decision process among N nests. First, we present the
decision-making dynamics in the symmetric case of N equal-quality nests. Then,
we generalise our findings to a best-of-N decision scenario with one superior
nest and N-1 inferior nests, previously studied empirically in bees and ants.
Whereas previous binary models highlighted the crucial role of inhibitory
stop-signalling, the key parameter in our new analysis is the relative time
invested by swarm members in individual discovery and in signalling behaviours.
Our new analysis reveals conflicting pressures on this ratio in symmetric and
best-of-N decisions, which could be solved through a time-dependent signalling
strategy. Additionally, our analysis suggests how ecological factors
determining the density of suitable nest sites may have led to selective
pressures for an optimal stable signalling ratio.