  The Poland-Scheraga model describes the denaturation transition of two
complementary - in particular, equally long - strands of DNA, and it has
enjoyed a remarkable success both for quantitative modeling purposes and at a
more theoretical level. The solvable character of the homogeneous version of
the model is one of features to which its success is due. In the bio-physical
literature a generalization of the model, allowing different length and non
complementarity of the strands, has been considered and the solvable character
extends to this substantial generalization. We present a mathematical analysis
of the homogeneous generalized Poland-Scheraga model. Our approach is based on
the fact that such a model is a homogeneous pinning model based on a bivariate
renewal process, much like the basic Poland-Scheraga model is a pinning model
based on a univariate, i.e. standard, renewal. We present a complete analysis
of the free energy singularities, which include the localization-delocalization
critical point and (in general) other critical points that have been only
partially captured in the physical literature. We obtain also precise estimates
on the path properties of the model.


  The ability to estimate the evolutionary distance between extant genomes
plays a crucial role in many phylogenomic studies. Often such estimation is
based on the parsimony assumption, implying that the distance between two
genomes can be estimated as the rearrangement distance equal the minimal number
of genome rearrangements required to transform one genome into the other.
However, in reality the parsimony assumption may not always hold, emphasizing
the need for estimation that does not rely on the rearrangement distance. The
distance that accounts for the actual (rather than minimal) number of
rearrangements between two genomes is often referred to as the true
evolutionary distance. While there exists a method for the true evolutionary
distance estimation, it however assumes that genomes can be broken by
rearrangements equally likely at any position in the course of evolution. This
assumption, known as the random breakage model, has recently been refuted in
favor of the more rigorous fragile breakage model postulating that only certain
"fragile" genomic regions are prone to rearrangements.
  We propose a new method for estimating the true evolutionary distance between
two genomes under the fragile breakage model. We evaluate the proposed method
on simulated genomes, which show its high accuracy. We further apply the
proposed method for estimation of evolutionary distances within a set of five
yeast genomes and a set of two fish genomes.


  Motivation: Prediction of phenotypes from high-dimensional data is a crucial
task in precision biology and medicine. Many technologies employ genomic
biomarkers to characterize phenotypes. However, such elements are not
sufficient to explain the underlying biology. To improve this, pathway analysis
techniques have been proposed. Nevertheless, such methods have shown lack of
accuracy in phenotypes classification. Results: Here we propose a novel
methodology called MITHrIL (Mirna enrIched paTHway Impact anaLysis) for the
analysis of signaling pathways, which has built on top of the work of Tarca et
al., 2009. MITHrIL extends pathways by adding missing regulatory elements, such
as microRNAs, and their interactions with genes. The method takes as input the
expression values of genes and/or microRNAs and returns a list of pathways
sorted according to their deregulation degree, together with the corresponding
statistical significance (p-values). Our analysis shows that MITHrIL
outperforms its competitors even in the worst case. In addition, our method is
able to correctly classify sets of tumor samples drawn from TCGA. Availability:
MITHrIL is freely available at the following URL:
http://alpha.dmi.unict.it/mithril/


  Active diffusion of intracellular components is emerging as an important
process in cell biology. This process is mediated by complex assemblies of
molecular motors and cytoskeletal filaments that drive force generation in the
cytoplasm and facilitate enhanced motion. The kinetics of molecular motors have
been precisely characterized in-vitro by single molecule approaches, however,
their in-vivo behavior remains elusive. Here, we study the active diffusion of
vesicles in mouse oocytes, where this process plays a key role in nuclear
positioning during development, and combine an experimental and theoretical
framework to extract molecular-scale force kinetics (force, power-stroke, and
velocity) of the in-vivo active process. Assuming a single dominant process, we
find that the nonequilibrium activity induces rapid kicks of duration $\tau
\sim$ 300 $\mu$s resulting in an average force of $F \sim$ 0.4 pN on vesicles
in in-vivo oocytes, remarkably similar to the kinetics of in-vitro myosin-V.
Our results reveal that measuring in-vivo active fluctuations allows extraction
of the molecular-scale activity in agreement with single-molecule studies and
demonstrates a mesoscopic framework to access force kinetics.


  The structure and quantitative composition of the cerebral cortex are
interrelated with its computational capacity. Empirical data analyzed here
indicate a certain hierarchy in local cortical composition. Specifically,
neural wire, i.e., axons and dendrites take each about 1/3 of cortical space,
spines and glia/astrocytes occupy each about $(1/3)^{2}$, and capillaries
around $(1/3)^{4}$. Moreover, data analysis across species reveals that these
fractions are roughly brain size independent, which suggests that they could be
in some sense optimal and thus important for brain function. Is there any
principle that sets them in this invariant way? This study first builds a model
of local circuit in which neural wire, spines, astrocytes, and capillaries are
mutually coupled elements and are treated within a single mathematical
framework. Next, various forms of wire minimization rule (wire length, surface
area, volume, or conduction delays) are analyzed, of which, only minimization
of wire volume provides realistic results that are very close to the empirical
cortical fractions. As an alternative, a new principle called "spine economy
maximization" is proposed and investigated, which is associated with
maximization of spine proportion in the cortex per spine size that yields
equally good but more robust results. Additionally, a combination of wire cost
and spine economy notions is considered as a meta-principle, and it is found
that this proposition gives only marginally better results than either pure
wire volume minimization or pure spine economy maximization, but only if spine
economy component dominates. In sum, these results suggest that the new spine
economy principle may be important for brain evolutionary design in a broader
context.


  An increasing number of complex systems are now modeled as networks of
coupled dynamical entities. Nonlinearity and high-dimensionality are hallmarks
of the dynamics of such networks but have generally been regarded as obstacles
to control. Here I discuss recent advances on mathematical and computational
approaches to control high-dimensional nonlinear network dynamics under general
constraints on the admissible interventions. I also discuss the potential of
network control to address pressing scientific problems in various disciplines.


  Selection, mutation and random drift affect the dynamics of allele
frequencies and consequently of quantitative traits. While the macroscopic
dynamics of quantitative traits can be measured, the underlying allele
frequencies are typically unobserved. Can we understand how the macroscopic
observables evolve without following these microscopic processes? The problem
has previously been studied by analogy with statistical mechanics: the allele
frequency distribution at each time is approximated by the stationary form,
which maximizes entropy. We explore the limitations of this method when
mutation is small ($4N\!\mu<1$) so that populations are typically close to
fixation and we extend the theory in this regime to account for changes in
mutation strength. We consider a single diallelic locus under either
directional selection, or with over-dominance, and then generalise to multiple
unlinked biallelic loci with unequal effects. We find that the maximum entropy
approximation is remarkably accurate, even when mutation and selection change
rapidly.


  In this review we discuss different mathematical models of gene regulatory
networks as relevant to the onset and development of cancer. After discussion
of alternative modelling approaches, we use a paradigmatic two-gene network to
focus on the role played by time delays in the dynamics of gene regulatory
networks. We contrast the dynamics of the reduced model arising in the limit of
fast mRNA dynamics with that of the full model. The review concludes with the
discussion of some open problems.


  We introduce and investigate an SIS-type model for the spread of an
infectious disease, where the infected population is structured with respect to
the different strain of the virus/bacteria they are carrying. Our aim is to
capture the interesting scenario when individuals infected with different
strains cause secondary (new) infections at different rates. Therefore, we
consider a nonlinear infection process, which generalises the bilinear process
arising from the classic mass-action assumption. Our main motivation is to
study competition between different strains of a virus/bacteria. From the
mathematical point of view, we are interested whether the nonlinear infection
process leads to a well-posed model. We use a semilinear formulation to show
global existence and positivity of solutions up to a critical value of the
exponent in the nonlinearity. Furthermore, we establish the existence of the
endemic steady state for particular classes of nonlinearities.


  Networks are ubiquitous throughout science and engineering. A number of
methods, including some from our own group, have explored how one goes about
computing or predicting the dynamics of networks given information about
internal models of individual nodes and network connectivity, possibly with
additional information provided by statistical or descriptive metrics that
characterize the network. But what can be inferred about network dynamics when
there is no knowledge or information about the internal model or dynamics of
participating nodes? Here, we explore how connected subsets of nodes
competitively interact in order to activate a common downstream node they
connect into. We achieve this by assuming a simple set of rules borrowed from
neurophysiology. The model we develop reflects a local process from which
global network dynamics emerges. We call this model a competitive refractory
dynanics model. It is derived from a consideration of spatial and temporal
summation in biological neurons, whereby summating post synaptic potentials
(PSPs) along the dendritic tree contribute towards the membrane potential at
the initial segment reaching a threshold potential. We first show how the
'winning node' or set of 'winning' nodes that achieve activation of a
downstream node is computable by the model. We then derive a formal definition
of optimized network signaling within our framework. We define a ratio between
the signaling latencies on the edges of the network and the internal time it
takes individual nodes to process incoming signals. We show that an optimal
ratio is one where the speed of information propagation between connected nodes
does not exceed the internal dynamic time scale of the nodes. We then show how
we can use these results to arrive at a unique interpretation for the
prevalence of the small world network topology in natural and engineered
systems.


  The transition from petascale to exascale computers is characterized by
substantial changes in the computer architectures and technologies. The
research community relying on computational simulations is being forced to
revisit the algorithms for data generation and analysis due to various
concerns, such as higher degrees of concurrency, deeper memory hierarchies,
substantial I/O and communication constraints. Simulations today typically save
all data to analyze later. Simulations at the exascale will require us to
analyze data as it is generated and save only what is really needed for
analysis, which must be performed predominately in-situ, i.e., executed
sufficiently fast locally, limiting memory and disk usage, and avoiding the
need to move large data across nodes.
  In this paper, we present a distributed method that enables in-situ data
analysis for large protein folding trajectory datasets. Traditional trajectory
analysis methods currently follow a centralized approach that moves the
trajectory datasets to a centralized node and processes the data only after
simulations have been completed. Our method, on the other hand, captures
conformational information in-situ using local data only while reducing the
storage space needed for the part of the trajectory under consideration. This
method processes the input trajectory data in one pass, breaks from the
centralized approach of traditional analysis, avoids the movement of trajectory
data, and still builds the global knowledge on the formation of individual
$\alpha$-helices or $\beta$-strands as trajectory frames are generated.


  We study the geometry of metrics and convexity structures on the space of
phylogenetic trees, which is here realized as the tropical linear space of all
\ ultrametrics. The ${\rm CAT}(0)$-metric of Billera-Holmes-Vogtman arises from
the theory of orthant spaces. While its geodesics can be computed by the
Owen-Provan algorithm, geodesic triangles are complicated. We show that the
dimension of such a triangle can be arbitrarily high. Tropical convexity and
the tropical metric behave better. They exhibit properties desirable for
geometric statistics, such as geodesics of small depth.


  Competition among gametes for fertilization imposes strong selection. For
external fertilizers, this selective pressure extends to eggs for which
spawning conditions can range from sperm limitation (competition among eggs) to
sexual conflict (overabundance of competing sperm toxic to eggs). Yet existing
fertilization models ignore dynamics that can alter the functional nature of
gamete interactions. These factors include attraction of sperm to eggs, egg
crowding effects or other nonlinearities in per capita rates of sperm-egg
interaction. Such processes potentially allow egg concentrations to drastically
affect viable fertilization probabilities. I experimentally tested whether such
egg effects occur using the urchin $\textit{Strongylocentrotus purpuratus}$ and
parameterized a newly derived model of fertilization dynamics and existing
models modified to include such interactions. The experiments revealed that at
low sperm concentrations, eggs compete for sperm while at high sperm
concentrations eggs cooperatively reduce abnormal fertilization (a proxy for
polyspermy). I show that these observations are consistent with declines in the
per capita rate at which sperm and eggs interact as eggs increase in density.
The results suggest a fitness trade-off of egg release during spawning: as
sperm range from scarce to superabundant, interactions among eggs transition
from highly competitive to cooperative in terms of viable fertilization
probabilities.


  Dendritic spines are small membranous structures that protrude from the
neuronal dendrite. Each spine contains a synaptic contact site that may connect
its parent dendrite to the axons of neighboring neurons. Dendritic spines are
markedly distinct in shape and size, and certain types of stimulation prompt
spines to evolve, in fairly predictable fashion, from thin nascent morphologies
to the mushroom-like shapes associated with mature spines. This striking
progression is coincident with the (re)configuration of the neuronal network
during early development, learning and memory formation, and has been
conjectured to be part of the machinery that encodes these processes at the
scale of individual neuronal connections. It is well established that the
structural plasticity of spines is strongly dependent upon the actin
cytoskeleton inside the spine. A general framework that details the precise
role of actin in directing the transitions between the various spine shapes is
lacking. We address this issue, and present a quantitative, model-based
scenario for spine plasticity validated using realistic and physiologically
relevant parameters. Our model points to a crucial role for the actin
cytoskeleton. In the early stages of spine formation, the interplay between the
elastic properties of the spine membrane and the protrusive forces generated in
the actin cytoskeleton propels the incipient spine. In the maturation stage,
actin remodeling in the form of the combined dynamics of branched and bundled
actin is required to form mature, mushroom-like spines. Our model identifies
additional factors that plausibly aid the stabilization and maintenance of
spine morphology. Taken together, our model provides unique insights into the
fundamental role of actin remodeling and polymerization forces during spine
formation and maturation.


  Rewarding cooperation is in many ways expected behaviour from social players.
However, strategies that promote antisocial behaviour are also surprisingly
common, not just in human societies, but also among eusocial insects and
bacteria. Examples include sanctioning of individuals who behave prosocially,
or rewarding of freeriders who do not contribute to collective enterprises. We
therefore study the public goods game with antisocial and prosocial pool
rewarding in order to determine the potential negative consequences on the
effectiveness of positive incentives to promote cooperation. Contrary to a
naive expectation, we show that the ability of defectors to distribute rewards
to their like does not deter public cooperation as long as cooperators are able
to do the same. Even in the presence of antisocial rewarding the spatial
selection for cooperation in evolutionary social dilemmas is enhanced. Since
the administration of rewards to either strategy requires a considerable degree
of aggregation, cooperators can enjoy the benefits of their prosocial
contributions as well as the corresponding rewards. Defectors when aggregated,
on the other hand, can enjoy antisocial rewards, but due to their lack of
contributions to the public good they ultimately succumb to their inherent
inability to secure a sustainable future. Strategies that facilitate the
aggregation of akin players, even if they seek to promote antisocial behaviour,
thus always enhance the long-term benefits of cooperation.


  Calmodulin (CaM) is a ubiquitous calcium binding protein consisting of two
structurally similar domains with distinct stabilities, binding affinities, and
flexibilities. We present coarse grained simulations that suggest the mechanism
for the domain's allosteric transitions between the open and closed
conformations depend on subtle differences in the folded state topology of the
two domains. Throughout a wide temperature range, the simulated transition
mechanism of the N-terminal domain (nCaM) follows a two-state transition
mechanism while domain opening in the C-terminal domain (cCaM) involves
unfolding and refolding of the tertiary structure. The appearance of the
unfolded intermediate occurs at a higher temperature in nCaM than it does in
cCaM consistent with nCaM's higher thermal stability. Under approximate
physiological conditions, the simulated unfolded state population of cCaM
accounts for 10% of the population with nearly all of the sampled transitions
(approximately 95%) unfolding and refolding during the conformational change.
Transient unfolding significantly slows the domain opening and closing rates of
cCaM. This potentially influences the mechanism of calcium binding to each
domain.


  It is challenging to apply the tenets of individuality to filamentous fungi:
a fungal mycelium can contain millions of genetically diverse but totipotent
nuclei, each capable of founding new mycelia. Moreover a single mycelium can
potentially stretch over kilometers and it is unlikely that its distant parts
share resources or have the same fitness. Here we directly measure how a single
mycelium of the model ascomycete {\it Neurospora crassa} is patterned into
Reproductive Units (RUs); subpopulations of nuclei that propagate together as
spores, and function as reproductive individuals. The density of RUs is
sensitive to the geometry of growth; we detect 50-fold smaller RUs when mycelia
had expanding frontiers than when they were constrained to grow in one
direction only. RUs fragmented further when the mycelial network was perturbed.
In mycelia with expanding frontiers RU composition was strongly influenced by
the distribution of genotypes early in development. Our results provide a
concept of fungal individuality that is directly connected to reproductive
potential, and therefore to theories of how fungal individuals adapt and evolve
over time. Our data show that the size of reproductive individuals is a dynamic
and environment-dependent property, even within apparently totally connected
fungal mycelia.


  The recent discovery of zero-determinant strategies for the iterated
Prisoner's Dilemma sparked a surge of interest in the surprising fact that a
player can exert unilateral control over iterated interactions. These
remarkable strategies, however, are known to exist only in games in which
players choose between two alternative actions such as "cooperate" and
"defect." Here we introduce a broader class of autocratic strategies by
extending zero-determinant strategies to iterated games with more general
action spaces. We use the continuous Donation Game as an example, which
represents an instance of the Prisoner's Dilemma that intuitively extends to a
continuous range of cooperation levels. Surprisingly, despite the fact that the
opponent has infinitely many donation levels from which to choose, a player can
devise an autocratic strategy to enforce a linear relationship between his or
her payoff and that of the opponent even when restricting his or her actions to
merely two discrete levels of cooperation. In particular, a player can use such
a strategy to extort an unfair share of the payoffs from the opponent.
Therefore, although the action space of the continuous Donation Game dwarfs
that of the classical Prisoner's Dilemma, players can still devise relatively
simple autocratic and, in particular, extortionate strategies.


  Recent years have witnessed a surge of biological interest in the minimum
spanning tree (MST) problem for its relevance to automatic model construction
using the distances between data points. Despite the increasing use of MST
algorithms for this purpose, the goodness-of-fit of an MST to the data is often
elusive because no quantitative criteria have been developed to measure it.
Motivated by this, we provide a necessary and sufficient condition to ensure
that a metric space on n points can be represented by a fully labeled tree on n
vertices, and thereby determine when an MST preserves all pairwise distances
between points in a finite metric space.


  The dynamics of contact networks and epidemics of infectious diseases often
occur on comparable time scales. Ignoring one of these time scales may provide
an incomplete understanding of the population dynamics of the infection
process. We develop an individual-based approximation for the
susceptible-infected-recovered epidemic model applicable to arbitrary dynamic
networks. Our framework provides, at the individual-level, the probability flow
over time associated with the infection dynamics. This computationally
efficient framework discards the correlation between the states of different
nodes, yet provides accurate results in approximating direct numerical
simulations. It naturally captures the temporal heterogeneities and
correlations of contact sequences, fundamental ingredients regulating the
timing and size of an epidemic outbreak. Using real-life data, we show that the
static network model overestimates the reproduction number but underestimates
the infection potential of super-spreading individuals. The high accuracy of
our approximation further allows us to detect the index individual of an
epidemic outbreak.


  A new class of formal latent-variable stochastic processes called hidden
quantum models (HQM's) is defined in order to clarify the theoretical
foundations of ion channel signal processing. HQM's are based on quantum
stochastic processes which formalize time-dependent observation. They allow the
calculation of autocovariance functions which are essential for
frequency-domain signal processing. HQM's based on a particular type of
observation protocol called independent activated measurements are shown to to
be distributionally equivalent to hidden Markov models yet without an
underlying physical Markov process. Since the formal Markov processes are
non-physical, the theory of activated measurement allows merging energy-based
Eyring rate theories of ion channel behavior with the more common
phenomenological Markov kinetic schemes to form energy-modulated quantum
channels. Using the simplest quantum channel model consistent with neuronal
membrane voltage-clamp experiments, activation eigenenergies are calculated for
the Hodgkin-Huxley K+ and Na+ ion channels. It is also shown that maximizing
entropy under constrained activation energy yields noise spectral densities
approximating $S(f) \sim 1/f^\alpha$, thus offering a biophysical explanation
for the ubiquitous $1/f$-type in neurological signals.


  In this work we explore the temporal dynamics of spatial heterogeneity during
the process of tumorigenesis from healthy tissue. We utilize a spatial
stochastic process model of mutation accumulation and clonal expansion in a
structured tissue to describe this process. Under a two-step tumorigenesis
model, we first derive estimates of a non-spatial measure of diversity:
Simpson's Index, which is the probability that two individuals sampled at
random from the population are identical, in the premalignant population. We
next analyze two new measures of spatial population heterogeneity. In
particular we study the typical length scale of genetic heterogeneity during
the carcinogenesis process and estimate the extent of a surrounding
premalignant clone given a clinical observation of a premalignant point biopsy.
This evolutionary framework contributes to a growing literature focused on
developing a better understanding of the spatial population dynamics of cancer
initiation and progression. Although initially motivated by understanding
questions in cancer, these results can be applied more generally to help
understand the dynamics of heterogeneity and diversity in a variety of
spatially structured, evolving populations.


  The diversity of the immune repertoire is initially generated by random
rearrangements of the receptor gene during early T and B cell development.
Rearrangement scenarios are composed of random events -- choices of gene
templates, base pair deletions and insertions -- described by probability
distributions. Not all scenarios are equally likely, and the same receptor
sequence may be obtained in several different ways. Quantifying the
distribution of these rearrangements is an essential baseline for studying the
immune system diversity. Inferring the properties of the distributions from
receptor sequences is a computationally hard problem, requiring enumerating
every possible scenario for every sampled receptor sequence. We present a
Hidden Markov model, which accounts for all plausible scenarios that can
generate the receptor sequences. We developed and implemented a method based on
the Baum-Welch algorithm that can efficiently infer the parameters for the
different events of the rearrangement process. We tested our software tool on
sequence data for both the alpha and beta chains of the T cell receptor. To
test the validity of our algorithm, we also generated synthetic sequences
produced by a known model, and confirmed that its parameters could be
accurately inferred back from the sequences. The inferred model can be used to
generate synthetic sequences, to calculate the probability of generation of any
receptor sequence, as well as the theoretical diversity of the repertoire. We
estimate this diversity to be $\approx 10^{23}$ for human T cells. The model
gives a baseline to investigate the selection and dynamics of immune
repertoires.


  Motility is characteristic of life, but a common basis for movement has
remained to be identified. Diverse systems in motion shift between two states
depending on interactions that turnover at the rate of an applied cycle of
force. Although one phase of the force cycle terminates the decay of the most
recent state, continuation of the cycle of force regenerates the original decay
process in a recursive cycle. By completing a cycle, kinetic energy is
transformed into probability of sustaining the most recent state and the system
gains a frame of reference for discrete transitions having static rather than
time-dependent probability. The probability of completing a recursive cycle is
computed with a Markov chain comprised of two equilibrium states and a kinetic
intermediate. Given rate constants for the reactions, a random walk reproduces
bias and recurrence times of walking motor molecules and bacterial flagellar
switching with unrivaled fidelity.


  Base-pockets (non-complementary base-pairs) in a double-stranded DNA play a
crucial role in biological processes. Because of thermal fluctuations, it can
lower the stability of DNA, whereas, in case of DNA aptamer, small molecules
e.g. adenosinemonophosphate(AMP), adenosinetriphosphate(ATP) etc, form
additional hydrogen bonds with base-pockets termed as "binding-pockets", which
enhance the stability. Using the Langevin Dynamics simulations of coarse
grained model of DNA followed by atomistic simulations, we investigated the
influence of base-pocket and binding-pocket on the stability of DNA aptamer.
Striking differences have been reported here for the separation induced by
temperature and force, which require further investigation by single molecule
experiments.


  As protein folding is a NP-complete problem, artificial intelligence tools
like neural networks and genetic algorithms are used to attempt to predict the
3D shape of an amino acids sequence. Underlying these attempts, it is supposed
that this folding process is predictable. However, to the best of our
knowledge, this important assumption has been neither proven, nor studied. In
this paper the topological dynamic of protein folding is evaluated. It is
mathematically established that protein folding in 2D hydrophobic-hydrophilic
(HP) square lattice model is chaotic as defined by Devaney. Consequences for
both structure prediction and biology are then outlined.


  How cells reliably infer information about their environment is a
fundamentally important question. While sensing and signaling generally start
with cell-surface receptors, the degree of accuracy with which a cell can
measure external ligand concentration with even the simplest device - a single
receptor - is surprisingly hard to pin down. Recent studies provide conflicting
results for the fundamental physical limits. Comparison is made difficult as
different studies either suggest different readout mechanisms of the
ligand-receptor occupancy, or differ on how ligand diffusion is implemented.
Here we critically analyse these studies and present a unifying perspective on
the limits of sensing, with wide-ranging biological implications.


  Our work is motivated by an interest in constructing a protein-protein
interaction network that captures key features associated with Parkinson's
disease. While there is an abundance of subnetwork construction methods
available, it is often far from obvious which subnetwork is the most suitable
starting point for further investigation. We provide a method to assess whether
a subnetwork constructed from a seed list (a list of nodes known to be
important in the area of interest) differs significantly from a randomly
generated subnetwork. The proposed method uses a Monte Carlo approach. As
different seed lists can give rise to the same subnetwork, we control for
redundancy by constructing a minimal seed list as the starting point for the
significance test. The null model is based on random seed lists of the same
length as a minimum seed list that generates the subnetwork, in this random
seed list the nodes have (approximately) the same degree distribution as the
nodes in the minimum seed list. We use this null model to select subnetworks
which deviate significantly from random on an appropriate set of statistics and
might capture useful information for a real world protein-protein interaction
network.


  Darwinian evolution can be modeled in general terms as a flow in the space of
fitness (i.e. reproductive rate) distributions. In the diffusion approximation,
Tsimring et al. have showed that this flow admits "fitness wave" solutions:
Gaussian-shape fitness distributions moving towards higher fitness values at
constant speed. Here we show more generally that evolving fitness distributions
are attracted to a one-parameter family of distributions with a fixed parabolic
relationship between skewness and kurtosis. Unlike fitness waves, this
statistical pattern encompasses both positive and negative (a.k.a. purifying)
selection and is not restricted to rapidly adapting populations. Moreover we
find that the mean fitness of a population under the selection of pre-existing
variation is a power-law function of time, as observed in microbiological
evolution experiments but at variance with fitness wave theory. At the
conceptual level, our results can be viewed as the resolution of the "dynamic
insufficiency" of Fisher's fundamental theorem of natural selection. Our
predictions are in good agreement with numerical simulations.


  Interventions are made in networks to change the network or its values in a
desired way. The intervention strategies evaluated in the study described here
use network sampling designs to find units to which interventions are applied.
An intervention applied to a network node or link can change a value associated
with that unit. Over time the effect of the intervention can have an effect on
the population that goes beyond the sample units to which it is directly
applied. This paper describes the methods used for this study. These include a
variety of link-tracing sampling designs in networks, a number of types of
interventions, and a temporal spatial network model in which the intervention
strategies are evaluated. An intervention strategy is associated with an agent
and different intervention strategies interact and adapt to each other over
time. Some preliminary results are summarized regarding potential intervention
strategies to help alleviate the HIV epidemic.


  As the understanding of the importance of social contact networks in the
spread of infectious diseases has increased, so has the interest in
understanding the feedback process of the disease altering the social network.
While many studies have explored the influence of individual epidemiological
parameters and/or underlying network topologies on the resulting disease
dynamics, we here provide a systematic overview of the interactions between
these two influences on population-level disease outcomes. We show that the
sensitivity of the population-level disease outcomes to the combination of
epidemiological parameters that describe the disease are critically dependent
on the topological structure of the population's contact network. We introduce
a new metric for assessing disease-driven structural damage to a network as a
population-level outcome. Lastly, we discuss how the expected individual-level
disease burden is influenced by the complete suite of epidemiological
characteristics for the circulating disease and the ongoing process of network
compromise. Our results have broad implications for prediction and mitigation
of outbreaks in both natural and human populations.


  Many membrane-bound molecules in cells form small clusters. It has been
hypothesized that these clusters convert an analog extracellular signal into a
digital intracellular signal and that this conversion increases signaling
fidelity. However, the mechanism by which clusters digitize a signal and the
subsequent effects on fidelity remain poorly understood. Here we demonstrate
using a stochastic model of cooperative cluster formation that sufficient
cooperation leads to digital signaling. We show that despite reducing the
number of output states, which decreases fidelity, digitization also reduces
noise in the system, which increases fidelity. The tradeoff between these
effects leads to an optimal cluster size that agrees with experimental
measurements.


  Assume that $N_m(x)$ denotes the density of the population at a point $x$ at
the beginning of the reproductive season in the $m$th year. We study the
following impulsive reaction-diffusion model for any $m\in \mathbb Z^+$
\begin{eqnarray*}\label{}
  \ \ \ \ \ \left\{ \begin{array}{lcl}
  u^{(m)}_t = div(A\nabla u^{(m)}-a u^{(m)}) + f(u^{(m)}) \quad \text{for} \ \
(x,t)\in\Omega\times (0,1]
  u^{(m)}(x,0)=g(N_m(x)) \quad \text{for} \ \ x\in \Omega
  N_{m+1}(x):=u^{(m)}(x,1) \quad \text{for} \ \ x\in \Omega \end{array}\right.
  \end{eqnarray*} for functions $f,g$, a drift $a$ and a diffusion matrix $A$
and $\Omega\subset \mathbb R^n$. Study of this model requires a simultaneous
analysis of the differential equation and the recurrence relation. When
boundary conditions are hostile we provide critical domain results showing how
extinction versus persistence of the species arises, depending on the size and
geometry of the domain. We show that there exists an {\it extreme volume size}
such that if $|\Omega|$ falls below this size the species is driven extinct,
regardless of the geometry of the domain. To construct such extreme volume
sizes and critical domain sizes, we apply Schwarz symmetrization rearrangement
arguments, the classical Rayleigh-Faber-Krahn inequality and the spectrum of
uniformly elliptic operators. The critical domain results provide qualitative
insight regarding long-term dynamics for the model. Lastly, we provide
applications of our main results to certain biological reaction-diffusion
models regarding marine reserve, terrestrial reserve, insect pest outbreak and
population subject to climate change.


  To characterize the thermodynamical equilibrium of DNA chains interacting
with a solution of non-specific binding proteins, a Flory-Huggins free energy
model was implemented. We explored the dependence on DNA and protein
concentrations of the DNA collapse. For physiologically relevant values of the
DNA-protein affinity, this collapse gives rise to a biphasic regime with a
dense and a dilute phase; the corresponding phase diagram was computed. Using
an approach based on Hamiltonian paths, we show that the dense phase has either
a molten globule or a crystalline structure, depending on the DNA bending
rigidity, which is influenced by the ionic strength. These results are valid at
the thermodynamical equilibrium and should therefore be consistent with many
biological processes, whose characteristic timescales range typically from 1 ms
to 10 s. Our model may thus be applied to biological phenomena that involve
DNA-binding proteins, such as DNA condensation with crystalline order, which
occurs in some bacteria to protect their chromosome from detrimental factors;
or transcription initiation, which occurs in clusters called transcription
factories that are reminiscent of the dense phase characterized in this study.


  Living organisms are inherently out-of-equilibrium systems. We employ new
developments in stochastic energetics and rely on a minimal microscopic model
to predict the amount of mechanical energy dissipated by such dynamics. Our
model includes complex rheological effects and nonequilibrium stochastic
forces. By performing active microrheology and tracking micron-sized vesicles
in the cytoplasm of living oocytes, we provide unprecedented measurements of
the spectrum of dissipated energy. We show that our model is fully consistent
with the experimental data, and we use it to offer predictions for the
injection and dissipation energy scales involved in active fluctuations.


  It is unknown if the current size and shape of protected areas (PA) in Costa
Rica favors retention of viable populations of wildlife. The western border of
Carara National Park (CNP) and its surroundings were characterized at a
landscape level in order to consider their implications for biodiversity
conservation. 54 400 hectares of Costa Rican Central Pacific tree cover were
analyzed (1997, 2000, 2005) to calculate the size of the PAs, mean shape index,
mean patch fractal dimension and edge density. Aside from that, we performed
eight habitat evaluations on the western border of the CNP (costanera-sur
highway) to assess tree composition, regeneration, litter, horizontal
obstruction, and canopy cover. It was determined that the PAs were between 36
and 5 242 hectares long. We observed that the tree cover increased from 21
231,8 hectares in 1997 to 29 006,9 hectares in 2000, and decreased to 26 933,4
hectares in 2005. We found out that most of the tree cover (2005) belongs to
CNP and ZP Cerros de Turrubares, but both of them present high values of edge
density and mean shape index, suggesting that they are susceptible to
alteration and fragmentation. The four evaluated points have a similar
successional stage. In order to maintain the potential of CNP as an area for
biodiversity conservation, we recommend reducing the maximum speed limit in the
region adjacent to the park. It is also important to establish frequent
monitoring of the tree cover and promote reforestation programs to create
corridors that stop the area's fragmentation and facilitate biodiversity
conservation.


  Onychophorans, or velvet worms, are poorly known and rare animals. Here we
report the discovery of a new species that is also the largest onychophoran
found so far, a 22cm long female from the Caribbean coastal forest of Costa
Rica. Specimens were examined with Scanning Electron Microscopy; Peripatus
solorzanoi sp. nov., is diagnosed as follows: primary papillae convex and
conical with rounded bases, with more than 18 scale ranks. Apical section
large, spherical, with a basal diameter of at least 20 ranks. Apical piece with
6-7 scale ranks. Outer blade 1 principal tooth, 1 accessory tooth, 1 vestigial
accessory tooth (formula: 1/1/1); inner blade 1 principal tooth, 1 accessory
tooth, 1 rudimentary accessory tooth, 9 to 10 denticles (formula: 1/1/1/9-10).
Accessory tooth blunt in both blades. Four pads in the fourth and fifth
oncopods; 4th. pad arched. The previously unknown mechanism by which
onychophorans weave their adhesive is simple: muscular action produces a
swinging movement of the adhesive-spelling organs; as a result, the streams
cross in mid air, weaving the net. Like all onychophorans, P. solorzanoi is a
rare species: active protection of the habitat of the largest onychophoran ever
described, is considered urgent.


  A quantitative analysis of several onychophoran characteristics shows that in
habitats with lower rain levels females reproduce at an older age, are more
fecund and tend to have reproductive diapause where rain does not exceed a mean
of 200 cm/year. These habitat characteristics are associated with the southern
family Peripatopsidae. Sex ratio and parental investment per young are not
correlated with general environmental conditions. A comparison of 72 species
showed that larger species are often more variable in morphometry, but species
with the longest females do not always have the longest males. Larger Peripatus
acacioi females (Peripatidae: Brazil) produce more and heavier off spring.
Intrapopulation morphology was studied in 12 peripatid species for which
samples of between II and 798 individuals were available. In general, within
populations the females are more variable than males in length and weight, but
similarly variable in the number of legs. The number of legs has a low
variability (1.73-2.45%). length is intermediate (22.4-25.3%) and weight is
very variable (49.41-75.17%). When sexes are compared within a population,
females can have 14-8.9 % more leg pairs, and be 47-63 % heavier and 26 %
longer than males.


  The populations of pitvipers from south western Costa Rica, have
traditionally been identified as Bothriechis schlegelii (Berthold). However, in
1954 E. H. Taylor described one specimen from the area as a new subspecies, B.
schlegelii supraciliaris. Werman returned supraciliaris to synonymy with
schlegelii four decades later. However, morphometry and color pattern in a SW
Costa Rica population (25 specimens) differ from those of specimens (N=57) from
other parts of Costa Rica and from descriptions of South American specimens.
Here the epithet Bothriechis schlegelii supraciliaris Taylor 1954. is
reestablished as a valid taxon and elevated to specific rank as B.
supraciliaris stat.nov. It is closely related to B. schlegelii from which it
differs by its color patterns based on a uniform ground color with polymorphic
dorsal designs, and its lower counts of ventral and caudal scales.


  Discrete and Continuous Dynamics is the first in a series of articles on
Network Models for Epidemiology. This project began in the Fall quarter of 2014
in my continuous modeling course. Since then, it has taken off and turned into
a series of articles, which I hope to compile into a single report. The purpose
of the report is to explore mathematical epidemiology. In this article, we
discuss the historical approach to disease modeling with compartmental models.
We discuss the issues and benefits of using network models. We build a discrete
dynamical system to describe infection and recovery of individuals in the
population. Lastly, we detail the computational scheme for iterating this
model.


  Consider a system evolving according to an absorbing discrete-time Markov
chain with known transition matrix. The state of the system is observed at two
points in time, separated by an unknown number of generations. We are
interested in calculating the expected elapsed time and its variance. We
provide a novel, exact solution, which is computable from the fundamental
matrix of a related absorbing Markov chain. This solution may be useful in
population genetics for computing the expected age of a segregating allele
without requiring diffusion approximation.


  I hypothesize that re-occurring prior experience of complex systems mobilizes
a fast response, whose attractor is encoded by their strongly connected network
core. In contrast, responses to novel stimuli are often slow and require the
weakly connected network periphery. Upon repeated stimulus, peripheral network
nodes remodel the network core that encodes the attractor of the new response.
This "core-periphery learning" theory reviews and generalizes the heretofore
fragmented knowledge on attractor formation by neural networks,
periphery-driven innovation and a number of recent reports on the adaptation of
protein, neuronal and social networks. The coreperiphery learning theory may
increase our understanding of signaling, memory formation, information encoding
and decision-making processes. Moreover, the power of network periphery-related
'wisdom of crowds' inventing creative, novel responses indicates that
deliberative democracy is a slow yet efficient learning strategy developed as
the success of a billion-year evolution.


  Successful adaptation helped the emergence of complexity. Alternating
plastic- and rigid-like states were recurrently considered to play a role in
adaptive processes. However, this extensive knowledge remained fragmented. In
this paper I describe plasticity-rigidity cycles as a general adaptation
mechanism operating in molecular assemblies, assisted protein folding, cellular
differentiation, learning, memory formation, creative thinking, as well as the
organization of social groups and ecosystems. Plasticity-rigidity cycles enable
a novel understanding of aging, exploration/exploitation trade-off and
evolvability, as well as help the design of efficient interventions in medicine
and in crisis management of financial and biological ecosystems.


  Simulation of biomolecular networks is now indispensable for studying
biological systems, from small reaction networks to large ensembles of cells.
Here we present a novel approach for stochastic simulation of networks embedded
in the dynamic environment of the cell and its surroundings. We thus sample
trajectories of the stochastic process described by the chemical master
equation with time-varying propensities. A comparative analysis shows that
existing approaches can either fail dramatically, or else can impose
impractical computational burdens due to numerical integration of reaction
propensities, especially when cell ensembles are studied. Here we introduce the
Extrande method which, given a simulated time course of dynamic network inputs,
provides a conditionally exact and several orders-of-magnitude faster
simulation solution. The new approach makes it feasible to demonstrate, using
decision-making by a large population of quorum sensing bacteria, that
robustness to fluctuations from upstream signaling places strong constraints on
the design of networks determining cell fate. Our approach has the potential to
significantly advance both understanding of molecular systems biology and
design of synthetic circuits.


  We determine the impact of resource renewal on the lifetime of a forager that
depletes its environment and starves if it wanders too long without eating. In
the framework of the minimal starving random walk model with resource renewal,
there are three universal classes of behavior as a function of the renewal
time. For sufficiently rapid renewal, foragers are immortal, while foragers
have a finite lifetime otherwise. In one dimension, there is a third regime,
for sufficiently slow renewal, in which the lifetime of the forager is
independent of the renewal time. We outline an enumeration method to determine
the mean lifetime of the forager in the mortal regime.


  A disease that caused by dengue virus (DENV) has become the major health
problem of the world. Nowadays, no effective treatment is available to overcome
the disease due to the level of dengue virus pathogeneses. A novel treatment
method such as antiviral drug is highly necessary for coping with the dengue
disease. Envelope protein is one of the non-structural proteins of DENV, which
engaged in the viral fusion process. The fusion process is mediated by the
conformational change in the protein structure from dimer to trimer state. The
previous research showed the existing cavity on the dimer structure of the
envelope protein. The existing ligand could get into cavity of the envelope
protein, stabilize the dimer structure or hamper the transition of dimer
protein into trimer. In this fashion, the fusion process can be prevented. The
aim of this research is designing the cyclic peptide with prolin-prolin bond as
fusion inhibitor of DENV envelope protein through molecular docking and
molecular dynamics simulation. The screening of 3,883 cyclic peptides, each of
them connected by prolin-prolin bond, through molecular docking resulted in
five best ligands. The pharmacological and toxicity character of these five
ligands were analised in silico. The result showed that PYRRP was the best
ligand. PAWRP was also chosen as the best ligand because it showed good
affinity with protein cavity. Stability of ligand-protein complex was analyzed
by molecular dynamics simulation. The result showed that PYRRP ligand was able
to support the stability of DENV envelope protein dimer structure at 310 K and
312 K. While PAWRP ligand actively formed complex with the DENV envelope
protein at 310 K compared to 312 K. Thus the PYRRP ligand has a potential to be
developed as DENV fusion inhibitor.


  A new approach to understanding evolution [Val09], namely viewing it through
the lens of computation, has already started yielding new insights, e.g.,
natural selection under sexual reproduction can be interpreted as the
Multiplicative Weight Update (MWU) Algorithm in coordination games played among
genes [CLPV14]. Using this machinery, we study the role of mutation in changing
environments in the presence of sexual reproduction. Following [WVA05], we
model changing environments via a Markov chain, with the states representing
environments, each with its own fitness matrix. In this setting, we show that
in the absence of mutation, the population goes extinct, but in the presence of
mutation, the population survives with positive probability. On the way to
proving the above theorem, we need to establish some facts about dynamics in
games. We provide the first, to our knowledge, polynomial convergence bound for
noisy MWU in a coordination game. Finally, we also show that in static
environments, sexual evolution with mutation converges, for any level of
mutation.


  The three dimensional structure of DNA in the nucleus (chromatin) plays an
important role in many cellular processes. Recent experimental advances have
led to high-throughput methods of capturing information about chromatin
conformation on genome-wide scales. New models are needed to quantitatively
interpret this data at a global scale. Here we introduce the use of tools from
topological data analysis to study chromatin conformation. We use persistent
homology to identify and characterize conserved loops and voids in contact map
data and identify scales of interaction. We demonstrate the utility of the
approach on simulated data and then look data from both a bacterial genome and
a human cell line. We identify substantial multiscale topology in these
datasets.


  Ordinary differential equations (ODEs) are widely used to model biological,
(bio-)chemical and technical processes. The parameters of these ODEs are often
estimated from experimental data using ODE-constrained optimisation. This
article proposes a simple simulation-based approach for solving optimisation
problems with steady state constraints relying on an ODE. This simulation-based
optimisation method is tailored to the problem structure and exploits the local
geometry of the steady state manifold and its stability properties. A
parameterisation of the steady state manifold is not required. We prove local
convergence of the method for locally strictly convex objective functions.
Effciency and reliability of the proposed method are demonstrated in two
examples. The proposed method demonstrated better convergence properties than
existing general purpose methods and a significantly higher number of converged
starts per time.


  The understanding of kinase structure is mostly based on protein
crystallography, which is limited by the requirement to trap molecules within a
crystal lattice. Therefore, characterisations of the conformational dynamics of
the activation loop in solution are important to enhance the understanding of
molecular processes related to diseases and to support the discovery of small
molecule kinase inhibitors. In this work, we demonstrated that long molecular
dynamics simulations exhaustively sampled all the conformational space of the
activation loop of the Aurora-A kinase and of the methane-thiosulfonate spin
label, introduced into the activation loop for the electron paramagnetic
measurements. MD was used to determine structural fluctuations, order
parameters and rotational correlation times of the motion of the activation
loop and of the MTSL. Theoretical data obtained were used as input for the
calculation of the room temperature 9 GHz continuous wave EPR of the Aurora-A
kinase in solution and the comparison between simulated and experimental date
revealed that the motion of the protein and spin label occurred on comparable
timescales. This work is a starting point for deeper experimental and
theoretical studies of the rotational and translational diffusion properties of
the Aurora-A kinase protein related to its biological activity.


  Deprivation of essential nutrients can have stark consequences for many
processes in a cell. We consider amino acid starvation, which can result in
bottlenecks in mRNA translation when ribosomes stall due to lack of resources,
i.e. tRNAs charged with the missing amino acid. Recent experiments also show
less obvious effects such as increased charging of other (non-starved) tRNA
species and selective charging of isoaccepting tRNAs. We present a mechanism
which accounts for these observations, and shows that production of some
proteins can actually increase under starvation. One might assume that such
responses could only be a result of sophisticated control pathways, but here we
show that these effects can occur naturally due to changes in the supply and
demand for different resources, and that control can be accomplished through
selective use of rare codons. We develop a model for translation which includes
the dynamics of the charging and use of aa-tRNAs, explicitly taking into
account the effect of specific codon sequences. This constitutes a new control
mechanism in gene regulation which emerges at the community level, i.e., via
resources used by all ribosomes.


  Biophysicists are modeling conformations of interphase chromosomes, often
basing the strengths of interactions between segments distant on the genetic
map on contact frequencies determined experimentally. Here, instead, we develop
a fitting-free, minimal model: bivalent red and green "transcription factors"
bind to cognate sites in runs of beads ("chromatin") to form molecular bridges
stabilizing loops. In the absence of additional explicit forces, molecular
dynamic simulations reveal that bound "factors' spontaneously cluster -- red
with red, green with green, but rarely red with green -- to give structures
reminiscent of transcription factories. Binding of just two transcription
factors (or proteins) to active and inactive regions of human chromosomes
yields rosettes, topological domains, and contact maps much like those seen
experimentally. This emergent "bridging-induced attraction" proves to be a
robust, simple, and generic force able to organize interphase chromosomes at
all scales.


  If rock beats scissors and scissors beat paper, one might assume that rock
beats paper too. But this is not the case for intransitive relationships that
make up the famous rock-paper-scissors game. However, the sole presence of
paper might prevent rock from beating scissors, simply because paper beats
rock. This is the blueprint for the rock-paper-scissors game with protection
spillovers, which has recently been introduced as a new paradigm for
biodiversity in well-mixed microbial populations. Here we study the game in
structured populations, demonstrating that protection spillovers give rise to
spatial patterns that are impossible to observe in the classical
rock-paper-scissors game. We show that the spatiotemporal dynamics of the
system is determined by the density of stable vortices, which may ultimately
transform to frozen states, to propagating waves, or to target waves with
reversed propagation direction, depending further on the degree and type of
randomness in the interactions among the species. If vortices are rare, the
fixation to waves and complex oscillatory solutions is likelier. Moreover,
annealed randomness in interactions favors the emergence of target waves, while
quenched randomness favors collective synchronization. Our results demonstrate
that protection spillovers may fundamentally change the dynamics of cyclic
dominance in structured populations, and they outline the possibility of
programming pattern formation in microbial populations.


  Frequencies of $k$-mers in sequences are sometimes used as a basis for
inferring phylogenetic trees without first obtaining a multiple sequence
alignment. We show that a standard approach of using the squared-Euclidean
distance between $k$-mer vectors to approximate a tree metric can be
statistically inconsistent. To remedy this, we derive model-based distance
corrections for orthologous sequences without gaps, which lead to consistent
tree inference. The identifiability of model parameters from $k$-mer
frequencies is also studied. Finally, we report simulations showing the
corrected distance out-performs many other $k$-mer methods, even when sequences
are generated with an insertion and deletion process. These results have
implications for multiple sequence alignment as well, since $k$-mer methods are
usually the first step in constructing a guide tree for such algorithms.


  The correlated motion of flocks is an instance of global order emerging from
local interactions. An essential difference with analogous ferromagnetic
systems is that flocks are active: animals move relative to each other,
dynamically rearranging their interaction network. The effect of this
off-equilibrium element is well studied theoretically, but its impact on actual
biological groups deserves more experimental attention. Here, we introduce a
novel dynamical inference technique, based on the principle of maximum entropy,
which accodomates network rearrangements and overcomes the problem of slow
experimental sampling rates. We use this method to infer the strength and range
of alignment forces from data of starling flocks. We find that local bird
alignment happens on a much faster timescale than neighbour rearrangement.
Accordingly, equilibrium inference, which assumes a fixed interaction network,
gives results consistent with dynamical inference. We conclude that bird
orientations are in a state of local quasi-equilibrium over the interaction
length scale, providing firm ground for the applicability of statistical
physics in certain active systems.


  With recent advances in high-throughput cell biology the amount of cellular
biological data has grown drastically. Such data is often modeled as graphs
(also called networks) and studying them can lead to new insights into
molecule-level organization. A possible way to understand their structure is by
analysing the smaller components that constitute them, namely network motifs
and graphlets. Graphlets are particularly well suited to compare networks and
to assess their level of similarity but are almost always used as small
undirected graphs of up to five nodes, thus limiting their applicability in
directed networks. However, a large set of interesting biological networks such
as metabolic, cell signaling or transcriptional regulatory networks are
intrinsically directional, and using metrics that ignore edge direction may
gravely hinder information extraction. The applicability of graphlets is
extended to directed networks by considering the edge direction of the
graphlets. We tested our approach on a set of directed biological networks and
verified that they were correctly grouped by type using directed graphlets.
However, enumerating all graphlets in a large network is a computationally
demanding task. Our implementation addresses this concern by using a
state-of-the-art data structure, the g-trie, which is able to greatly reduce
the necessary computation. We compared our tool, gtrieScanner, to other
state-of-the art methods and verified that it is the fastest general tool for
graphlet counting.


  Simple elastic network models of DNA were developed to reveal the
structure-dynamics relationships for several nucleotide sequences. First, we
propose a simple all-atom elastic network model of DNA that can explain the
profiles of temperature factors for several crystal structures of DNA. Second,
we propose a coarse-grained elastic network model of DNA, where each nucleotide
is described only by one node. This model could effectively reproduce the
detailed dynamics obtained with the all-atom elastic network model according to
the sequence-dependent geometry. Through normal-mode analysis for the
coarse-grained elastic network model, we exhaustively analyzed the dynamic
features of a large number of long DNA sequences, approximately $\sim 150$ bp
in length. These analyses revealed positive correlations between the
nucleosome-forming abilities and the inter-strand fluctuation strength of
double-stranded DNA for several DNA sequences.


  The rise of multicellularity in the early evolution of life represents a
major challenge for evolutionary biology. Guidance for finding answers has
emerged from disparate fields, from phylogenetics to modelling and synthetic
biology, but little is known about the potential origins of multicellular
aggregates before genetic programs took full control of developmental
processes. Such aggregates should involve spatial organisation of
differentiated cells and the modification of flows and concentrations of
metabolites within well defined boundaries. Here we show that, in an
environment where limited nutrients and toxic metabolites are introduced, a
population of cells capable of stochastic differentiation and differential
adhesion can develop into multicellular aggregates with a complex internal
structure. The morphospace of possible patterns is shown to be very rich,
including proto-organisms that display a high degree of organisational
complexity, far beyond simple heterogeneous populations of cells. Our findings
reveal that there is a potentially enormous richness of organismal complexity
between simple mixed cooperators and embodied living organisms.


  Motivated by the famous Waddington's epigenetic landscape metaphor in
developmental biology, biophysicists and applied mathematicians made different
proposals to realize this metaphor in a rationalized way. We adopt
comprehensive perspectives to systematically investigate three different but
closely related realizations in recent literature: namely the potential
landscape theory from the steady state distribution of stochastic differential
equations (SDEs), the quasi-potential from the large deviation theory, and the
construction through SDE decomposition and A-type integral.The connections
among these theories are established in this paper. We demonstrate that the
quasi-potential is the zero noise limit of the potential landscape. We also
show that the potential function in the third proposal coincides with the
quasi-potential. The most probable transition path by minimizing the
Onsager-Machlup or Freidlin-Wentzell action functional is discussed as well.
Furthermore, we compare the difference between local and global quasi-potential
through the exchange of limit order for time and noise amplitude. As a
consequence of such explorations, we arrive at the existence result for the SDE
decomposition while deny its uniqueness in general cases. It is also clarified
that the A-type integral is more appropriate to be applied to the decomposed
SDEs rather than the original one. Our results contribute to a better
understanding of existing landscape theories for biological systems.


  Notable recent works have focused on the multi-layer properties of coevolving
diseases. We point out that very similar systems play an important role in
population ecology. Specifically we study a meta food-web model that was
recently proposed by Pillai et al. This model describes a network of species
connected by feeding interactions, which spread over a network of spatial
patches. Focusing on the essential case, where the network of feeding
interactions is a chain, we develop an analytical approach for the computation
of the degree distributions of colonized spatial patches for the different
species in the chain. This framework allows us to address ecologically relevant
questions. Considering configuration model ensembles of spatial networks, we
find that there is an upper bound for the fraction of patches that a given
species can occupy, which depends only on the networks mean degree. For a given
mean degree there is then an optimal degree distribution that comes closest to
the upper bound. Notably scale-free degree distributions perform worse than
more homogeneous degree distributions if the mean degree is sufficiently high.
Because species experience the underlying network differently the optimal
degree distribution for one particular species is generally not the optimal
distribution for the other species in the same food web. These results are of
interest for conservation ecology, where, for instance, the task of selecting
areas of old-growth forest to preserve in an agricultural landscape, amounts to
the design of a patch network.


  Interactions between a protein and a ligand are often accompanied by a
redistribution of the population of thermally accessible conformations. This
dynamic response of the protein's functional energy landscape enables a protein
to modulate binding affinities and control binding sensitivity to ligand
concentration. In this paper, we investigate the structural origins of binding
affinity and allosteric cooperativity of binding two calcium ions to each
domain of calmodulin (CaM) through simulations of a simple coarse-grained
model. In this model, the protein's conformational transitions between open and
closed conformational ensembles are simulated explicitly and ligand binding and
unbinding is treated implicitly within the Grand Canonical Ensemble. Ligand
binding is cooperative because the binding sites are coupled through a shift in
the dominant conformational ensemble upon binding. The classic
Monod-Wyman-Changeux model of allostery with appropriate binding free energy to
the open and closed ensembles accurately describes the simulated binding
thermodynamics. The simulations predict that the two domains of CaM have
distinct binding affinity and cooperativity. In particular, C-terminal domain
binds calcium with higher affinity and greater cooperativity than the
N-terminal domain. From a structural point of view, the affinity of an
individual binding loop depends sensitively on the loop's structural
compatibility with the ligand in the bound ensemble, as well as the
conformational flexibility of the binding site in the unbound ensemble.


  We present a method for estimating epidemic parameters in network-based
stochastic epidemic models when the total number of infections is assumed to be
small. We illustrate the method by reanalyzing the data from the 2014
Democratic Republic of the Congo (DRC) Ebola outbreak described in Maganga et
al. (2014).


  The proper positioning of the mitotic spindle is crucial for asymmetric cell
division and generating cell diversity during development. Proper position in
the single-cell embryo of Caenorhabditis elegans is achieved initially by the
migration and rotation of the pronuclear complex (PNC) and its two associated
centrosomal arrays of microtubules (MTs). We present here the first systematic
theoretical study of how these $O(1000)$ centrosomal microtubules (MTs)
interact through the immersing cytoplasm, the cell periphery and PNC, and with
each other, to achieve proper position. This study is made possible through our
development of a highly efficient and parallelized computational framework that
accounts explicitly for long-ranged hydrodynamic interactions (HIs) between the
MTs, while also capturing their flexibility, dynamic instability, and
interactions with molecular motors and boundaries. First, we show through
direct simulation that previous estimates of the PNC drag coefficient, based on
either ignoring or partially including HIs, lead to misprediction of the active
forces and time-scales of migration. We then directly study the dynamics of PNC
migration under various force-transduction models, including the pushing or
pulling of MTs at the cortex, and the pulling of MTs by cytoplasmically-bound
force generators. While achieving proper position and orientation on
physiologically reasonable time-scales does not uniquely choose a model, we
find that each model produces a different signature in its induced cytoplasmic
flow and MT conformations. We suggest then that cytoplasmic flows and MT
conformations can be used to differentiate between mechanisms and to determine
their contribution to the migration process.


  Understanding the relationship between genomic variation and variation in
phenotypes for quantitative traits such as physiology, yield, fitness or
behavior, will provide important insights for both predicting adaptive
evolution and for breeding schemes. A particular question is whether the
genetic variation that influences quantitative phenotypes is typically the
result of one or two mutations of large effect, or multiple mutations of small
effect. In this paper we explore this issue using the wild model legume
Medicago truncatula. We show that phenotypes, such as quantitative disease
resistance, can be well-predicted using genome-wide patterns of admixture, from
which it follows that there must be many mutations of small effect. Our
findings prove the potential of our novel 'whole-genome modeling' -WhoGEM-
method and experimentally validate, for the first time, the infinitesimal model
as a mechanism for adaptation of quantitative phenotypes in plants. This
insight can accelerate breeding and biomedicine research programs.


  Fractional-order SIR models have become increasingly popular in the
literature in recent years, however unlike the standard SIR model, they often
lack a derivation from an underlying stochastic process. Here we derive a
fractional-order infectivity SIR model from a stochastic process that
incorporates a time-since-infection dependence on the infectivity of
individuals. The fractional derivative appears in the generalised master
equations of a continuous time random walk through SIR compartments, with a
power-law function in the infectivity. We show that this model can also be
formulated as an infection-age structured Kermack-McKendrick
integro-differential SIR model. Under the appropriate limit the fractional
infectivity model reduces to the standard ordinary differential equation SIR
model.


  Living organisms survive and multiply even though they have uncertain and
incomplete information about their environment and imperfect models to predict
the consequences of their actions. Bayesian models have been proposed to face
this challenge. Indeed, Bayesian inference is a way to do optimal reasoning
when only uncertain and incomplete information is available. Various
perceptive, sensory-motor, and cognitive functions have been successfully
modeled this way. However, the biological mechanisms allowing animals and
humans to represent and to compute probability distributions are not known. It
has been proposed that neurons and assemblies of neurons could be the
appropriate scale to search for clues to probabilistic reasoning. In contrast,
in this paper, we propose that interacting populations of macromolecules and
diffusible messengers can perform probabilistic computation. This suggests that
probabilistic reasoning, based on cellular signaling pathways, is a fundamental
skill of living organisms available to the simplest unicellular organisms as
well as the most complex brains.


  The Moran process on graphs is a popular model to study the dynamics of
evolution in a spatially structured population. Exact analytical solutions for
the fixation probability and time of a new mutant have been found for only a
few classes of graphs so far. Simulations are time-expensive and many
realizations are necessary, as the variance of the fixation times is high. We
present an algorithm that numerically computes these quantities for arbitrary
small graphs by an approach based on the transition matrix. The advantage over
simulations is that the calculation has to be executed only once. Building the
transition matrix is automated by our algorithm. This enables a fast and
interactive study of different graph structures and their effect on fixation
probability and time. We provide a fast implementation in C with this note
https://github.com/hindersin/efficientFixation. Our code is very flexible, as
it can handle two different update mechanisms (Birth-death or death-Birth), as
well as arbitrary directed or undirected graphs.


  Most species are structured and influenced by processes that either increased
or reduced gene flow between populations. However, most population genetic
inference methods ignore population structure and reconstruct a history
characterized by population size changes under the assumption that species
behave as panmictic units. This is potentially problematic since population
structure can generate spurious signals of population size change. Moreover,
when the model assumed for demographic inference is misspecified, genomic data
will likely increase the precision of misleading if not meaningless parameters.
In a context of model uncertainty (panmixia \textit{versus} structure) genomic
data may thus not necessarily lead to improved statistical inference.
  We consider two haploid genomes and develop a theory which explains why any
demographic model (with or without population size changes) will necessarily be
interpreted as a series of changes in population size by inference methods
ignoring structure. We introduce a new parameter, the IICR (inverse
instantaneous coalescence rate), and show that it is equivalent to a population
size only in panmictic models, and mostly misleading for structured models. We
argue that this general issue affects all population genetics methods ignoring
population structure. We take the PSMC method as an example and show that it
infers population size changes that never took place. We apply our approach to
human genomic data and find a reduction in gene flow at the start of the
Pleistocene, a major increase throughout the Middle-Pleistocene, and an abrupt
disconnection preceding the emergence of modern humans.


  Somatic mutations in cancer cells affect various genomic elements disrupting
important cell functions. In particular, mutations in DNA binding sites
recognized by transcription factors can alter regulator binding affinities and
expression of target genes. A number of promoter mutations have been linked
with an increased risk of cancer, mutations in binding sites of selected
transcription factors have been found under positive selection. However,
negative selection of mutations in coding regions is elusive and significance
of negative selection in non-coding regions remains controversial.
  Here we present analysis of transcription factors with binding sites
co-localized with non-coding variants. To avoid statistical bias we account for
mutation signatures of different cancer types. For many transcription factors,
including multiple members of FOX, HOX, and NR families, we show that human
cancers accumulate fewer mutations than expected by chance that increase or
decrease affinity of binding motifs. Such conservation of motifs is even more
exhibited in DNase accessible regions.
  Our data demonstrate negative selection against binding sites alterations and
suggest that this selection pressure protects cancer cells from rewiring of
regulatory circuits. Further analysis of transcription factors and the
respective conserved binding motifs can reveal cell regulatory pathways crucial
for the survivability of various human cancers.


  Three-dimensional protein structures usually contain regions of local order,
called secondary structure, such as $\alpha$-helices and $\beta$-sheets.
Secondary structure is characterized by the local rotational state of the
protein backbone, quantified by two dihedral angles called $\phi$ and $\psi$.
Particular types of secondary structure can generally be described by a single
(diffuse) location on a two-dimensional plot drawn in the space of the angles
$\phi$ and $\psi$, called a Ramachandran plot. By contrast, a
recently-discovered nanomaterial made from peptoids, structural isomers of
peptides, displays a secondary-structure motif corresponding to two regions on
the Ramachandran plot [Mannige et al., Nature 526, 415 (2015)]. In order to
describe such `higher-order' secondary structure in a compact way we introduce
here a means of describing regions on the Ramachandran plot in terms of a
single Ramachandran number, ${\mathcal{R}}$, which is a structurally meaningful
combination of $\phi$ and $\psi$. We show that the potential applications of
${\mathcal{R}}$ are numerous: it can be used to describe the geometric content
of protein structures, and can be used to draw diagrams that reveal, at a
glance, the frequency of occurrence of regular secondary structures and
disordered regions in large protein datasets. We propose that ${\mathcal{R}}$
might be used as an order parameter for protein geometry for a wide range of
applications.


  Motivation: DNA data is transcribed into single-stranded RNA, which folds
into specific molecular structures. In this paper we pose the question to what
extent sequence- and structure-information correlate. We view this correlation
as structural semantics of sequence data that allows for a different
interpretation than conventional sequence alignment. Structural semantics could
enable us to identify more general embedded "patterns" in DNA and RNA
sequences. Results: We compute the partition function of sequences with respect
to a fixed structure and connect this computation to the mutual information of
a sequence-structure pair for RNA secondary structures. We present a Boltzmann
sampler and obtain the a priori probability of specific sequence patterns. We
present a detailed analysis for the three PDB-structures, 2JXV (hairpin), 2N3R
(3-branch multi-loop) and 1EHZ (tRNA). We localize specific sequence patterns,
contrast the energy spectrum of the Boltzmann sampled sequences versus those
sequences that refold into the same structure and derive a criterion to
identify native structures. We illustrate that there are multiple sequences in
the partition function of a fixed structure, each having nearly the same mutual
information, that are nevertheless poorly aligned. This indicates the
possibility of the existence of relevant patterns embedded in the sequences
that are not discoverable using alignments.


  We study the ultrasensitivity of multisite binding processes where ligand
molecules can bind to several binding sites, considering more particularly
recent models involving complex chemical reactions in phosphorylation systems
such as allosteric phosphorylation processes, or substrate-catalyst chain
reactions and nucleosome mediated cooperativity. New statistics based formulas
for the Hill coefficient and the effective Hill coefficient are provided and
necessary conditions for a system to be ultrasensitive are exhibited. We then
assume that the binding process is described by a density dependent birth and
death process. We provide precise large deviation results for the steady state
distribution of the process, and show that switch-like ultrasensitive responses
are strongly related to the multi-stability of the associated dynamical system.
Ultrasensitivity occurs if and only if the entropy of the dynamical system has
more than one global minimum for some critical ligand concentration. In this
case, the Hill coefficient is proportional to the number of binding sites, and
the systems is highly ultrasensitive. We also discuss the interpretation of an
extension $I_q$ of the effective Hill coefficient $I_{0.9}$ for which we
recommend the computation of a broad range of values of $q$ instead of just the
standard one corresponding to the 10% to 90% variation in the dose-response. It
is shown that this single choice can sometimes mislead the conclusion by not
detecting ultrasensitivity. This new approach allows a better understanding of
multisite ultrasensitive systems and provides new tools for the design of such
systems.


  Dysregulation in signal transduction pathways can lead to a variety of
complex disorders, including cancer. Computational approaches such as network
analysis are important tools to understand system dynamics as well as to
identify critical components that could be further explored as therapeutic
targets. Here, we performed perturbation analysis of a large-scale signal
transduction model in extracellular environments that stimulate cell death,
growth, motility, and quiescence. Each of the model's components was perturbed
under both loss-of-function and gain-of-function mutations. We identified the
most and least influential components based on the magnitude of their influence
on the rest of the system. Based on the premise that the most influential
components might serve as better drug targets, we characterized them for
biological functions, housekeeping genes, essential genes, and druggable
proteins. Moreover, known cancer drug targets were also classified in
influential components based on the affected components in the network.
Additionally, the systemic perturbation analysis of the model revealed a
network motif of most influential components which affect each other.
Furthermore, our analysis predicted novel combinations of cancer drug targets
with various effects on other most influential components. We found that the
combinatorial perturbation consisting of PI3K inactivation and overactivation
of IP3R1 can lead to increased activity levels of apoptosis-related components
and tumor suppressor genes, suggesting that this combinatorial perturbation may
lead to a better target for decreasing cell proliferation and inducing
apoptosis. Lastly, our results suggest that systematic perturbation analyses of
large-scale computational models may serve as an approach to prioritize and
assess signal transduction components in order to identify novel drug targets
in complex disorders.


  Mixed dispersal syndromes have historically been regarded as bet-hedging
mechanisms that enhance survival in unpredictable environments, ensuring that
some propagules stay in the maternal environment while others can potentially
colonize new sites. However, this entails paying the costs of both dispersal
and non-dispersal. Propagules that disperse are likely to encounter unfavorable
conditions for establishment, while non-dispersing propagules might form
populations of close relatives burdened with inbreeding. Here, we investigate
the conditions under which mixed dispersal syndromes emerge and are
evolutionarily stable, taking into account the risks of both environmental
unpredictability and inbreeding. Using mathematical and computational modeling
we show that high dispersal propensity is favored whenever temporal
environmental unpredictability is low and inbreeding depression high, whereas
mixed dispersal syndromes are adaptive under conditions of high environmental
unpredictability, but more particularly if also inbreeding depression is small.
Although pure dispersers can be selected for under some circumstances, mixed
dispersal provides the optimal strategy under most parameterizations of our
models, indicating that this strategy is likely to be favored under a wide
variety of conditions. Furthermore, populations exhibiting any single phenotype
go inevitably extinct when environmental and genetic costs are high, whilst
mixed strategies can maintain viable populations even under such conditions.
Our models support the hypothesis that the interplay between inbreeding
depression and environmental unpredictability shapes dispersal syndromes, often
resulting in mixed strategies. Moreover, mixed dispersal seems to facilitate
persistence whenever conditions are critical or nearly critical for survival.


  Same-sex sexual behavior is ubiquitous in the animal kingdom, but its
adaptive origins remain a prominent puzzle. Here I suggest the possibility that
same-sex sexual behavior arises as a consequence of the competition between an
evolutionary drive for a wide diversity in traits, which improves the
adaptability of a species, and a drive for sexual dichotomization of traits,
which promotes opposite-sex attraction and increases the rate of reproduction.
A simple analytical "toy model" is proposed for describing this tradeoff. The
model exhibits a number of interesting features, and suggests a simple
mathematical form for describing the sexual orientation continuum.


  Cells generally convert nutrient resources to useful products via energy
transduction. Accordingly, the thermodynamic efficiency of this conversion
process is one of the most essential characteristics of living organisms.
However, although these processes occur under conditions of dynamic metabolism,
most studies of cellular thermodynamic efficiency have been restricted to
examining steady states; thus, the relevance of dynamics to this efficiency has
not yet been elucidated. Here, we develop a simple model of metabolic reactions
with anabolism-catabolism coupling catalysed by enzymes. Through application of
external oscillation in the enzyme abundances, the thermodynamic efficiency of
metabolism was found to be improved. This result is in strong contrast with
that observed in the oscillatory input, in which the efficiency always
decreased with oscillation. This improvement was effectively achieved by
separating the anabolic and catabolic reactions, which tend to disequilibrate
each other, and taking advantage of the temporal oscillations so that each of
the antagonistic reactions could progress near equilibrium. In this case,
anti-phase oscillation between the reaction flux and chemical affinity through
oscillation of enzyme abundances is essential. This improvement was also
confirmed in a model capable of generating autonomous oscillations in enzyme
abundances. Finally, the possible relevance of the improvement in thermodynamic
efficiency is discussed with respect to the potential for manipulation of
metabolic oscillations in microorganisms.


  Gene expression levels exhibit stochastic variations among genetically
identical organisms under the same environmental conditions (called gene
expression "noise" or phenotype "fluctuation"). In yeast and Escherichia coli,
positive correlations have been found between such gene expression noise and
"plasticity" with environmental variations. To determine the universality of
such correlations in both unicellular and multicellular organisms, we focused
on the relationships between gene expression "noise" and "plasticity" in
Arabidopsis thaliana, a multicellular model organism. In recent studies on
yeast and E. coli, only some gene groups with specific properties of promoter
architecture, average expression levels, and functions exhibited strong
noise-plasticity correlations. However, we found strong noise-plasticity
correlations for most gene groups in Arabidopsis; additionally, promoter
architecture, functional essentiality of genes, and circadian rhythm appeared
to have only a weak influence on the correlation strength. The differences in
the characteristics of noise-plasticity correlations may result from
three-dimensional chromosomal structures and/or circadian rhythm.


  We examine immunostaining experimental data for the formation of the strip 2
of $even-skipped$ ($eve$) transcripts on $D.$ $melanogaster$ embryos. An
estimate of the factor converting immunofluorescence intensity units into
molecular numbers is given. The analysis of the $eve$ mRNA's dynamics at the
region of the stripe 2 suggests that the promoter site of the gene has two
distinct regimes: an earlier phase when it is predominantly activated until a
critical time when it becomes mainly repressed. That suggests proposing a
stochastic binary model for gene transcription on $D.$ $melanogaster$ embryos.
Our model has two random variables: the transcripts number and the state of the
source of mRNAs given as active or repressed. We are able to reproduce
available experimental data for the average number of transcripts. An analysis
of the random fluctuations on the number of $eve$ mRNA's and their consequences
on the spatial precision of the stripe 2 is presented. We show that the
position of the anterior/posterior borders fluctuate around their average
position by $\sim 1 \%$ of the embryo length which is similar to what is found
experimentally. The fitting of data by such a simple model suggests that it can
be useful to understand the functions of randomness during developmental
processes.


  We consider the general problem of sensitive and specific discrimination
between biochemical species. An important instance is immune discrimination
between self and not-self, where it is also observed experimentally that
ligands just below discrimination threshold negatively impact response, a
phenomenon called antagonism. We characterize mathematically the generic
properties of such discrimination, first relating it to biochemical adaptation.
Then, based on basic biochemical rules, we establish that, surprisingly,
antagonism is a generic consequence of any strictly specific discrimination
made independently from ligand concentration. Thus antagonism constitutes a
"phenotypic spandrel": a phenotype existing as a necessary by-product of
another phenotype. We exhibit a simple analytic model of discrimination
displaying antagonism, where antagonism strength is linear in distance from
detection threshold. This contrasts with traditional proofreading based models
where antagonism vanishes far from threshold and thus displays an inverted
hierarchy of antagonism compared to simpler models. The phenotypic spandrel
studied here is expected to structure many decision pathways such as immune
detection mediated by TCRs and FC$\epsilon$RIs, as well as endocrine
signalling/disruption.


  Many algorithms for sequence analysis rely on word matching or word
statistics. Often, these approaches can be improved if binary patterns
representing match and don't-care positions are used as a filter, such that
only those positions of words are considered that correspond to the match
positions of the patterns. The performance of these approaches, however,
depends on the underlying patterns. Herein, we show that the overlap complexity
of a pattern set that was introduced by Ilie and Ilie is closely related to the
variance of the number of matches between two evolutionarily related sequences
with respect to this pattern set. We propose a modified hill-climbing algorithm
to optimize pattern sets for database searching, read mapping and
alignment-free sequence comparison of nucleic-acid sequences; our
implementation of this algorithm is called rasbhari. Depending on the
application at hand, rasbhari can either minimize the overlap complexity of
pattern sets, maximize their sensitivity in database searching or minimize the
variance of the number of pattern-based matches in alignment-free sequence
comparison. We show that, for database searching, rasbhari generates pattern
sets with slightly higher sensitivity than existing approaches. In our Spaced
Words approach to alignment-free sequence comparison, pattern sets calculated
with rasbhari led to more accurate estimates of phylogenetic distances than the
randomly generated pattern sets that we previously used. Finally, we used
rasbhari to generate patterns for short read classification with CLARK-S. Here
too, the sensitivity of the results could be improved, compared to the default
patterns of the program. We integrated rasbhari into Spaced Words; the source
code of rasbhari is freely available at http://rasbhari.gobics.de/


  Simple growth mechanisms have been proposed to explain the emergence of
seemingly universal network structures. The widely-studied model of
preferential attachment assumes that new nodes are more likely to connect to
highly connected nodes. Preferential attachment explains the emergence of
scale-free degree distributions within complex networks. Yet, it is
incompatible with many network systems, particularly bipartite systems in which
two distinct types of agents interact. For example, the addition of new links
in a host-parasite system corresponds to the infection of hosts by parasites.
Increasing connectivity is beneficial to a parasite and detrimental to a host.
Therefore, the overall network connectivity is subject to conflicting
pressures. Here, we propose a stochastic network growth model of conflicting
attachment, inspired by a particular kind of parasite-host interactions: that
of viruses interacting with microbial hosts. The mechanism of network growth
includes conflicting preferences to network density as well as costs involved
in modifying the network connectivity according to these preferences. We find
that the resulting networks exhibit realistic patterns commonly observed in
empirical data, including the emergence of nestedness, modularity, and
nested-modular structures that exhibit both properties. We study the role of
conflicting interests in shaping network structure and assess opportunities to
incorporate greater realism in linking growth process to pattern in systems
governed by antagonistic and mutualistic interactions.


  The most recent Ebola outbreak in West Africa highlighted critical weaknesses
in the medical infrastructure of the affected countries, including effective
diagnostics tools, sufficient isolation wards, and enough medical personnel.
Here, we develop and analyze a mathematical model to assess the impact of early
diagnosis of pre-symptomatic individuals on the transmission dynamics of Ebola
virus disease in West Africa in scenarios where Ebola may remain at low levels
in the population. Our findings highlight the importance of implementing
integrated control measures of early diagnosis and isolation. The mathematical
analysis shows a threshold where early diagnosis of pre-symptomatic
individuals, combined with a sufficient level of effective isolation, can lead
to an epidemic control of Ebola virus disease. That is, the local erradication
of the disease or the effective management of the disease at low levels of
endemicity.