  Stochasticity in gene expression can give rise to fluctuations in protein
levels and lead to phenotypic variation across a population of genetically
identical cells. Recent experiments indicate that bursting and feedback
mechanisms play important roles in controlling noise in gene expression and
phenotypic variation. A quantitative understanding of the impact of these
factors requires analysis of the corresponding stochastic models. However, for
stochastic models of gene expression with feedback and bursting, exact
analytical results for protein distributions have not been obtained so far.
Here, we analyze a model of gene expression with bursting and feedback
regulation and obtain exact results for the corresponding protein steady-state
distribution. The results obtained provide new insights into the role of
bursting and feedback in noise regulation and optimization. Furthermore, for a
specific choice of parameters, the system studied maps on to a two-state
biochemical switch driven by a bursty input noise source. The analytical
results derived thus provide quantitative insights into diverse cellular
processes involving noise in gene expression and biochemical switching.


  We consider a continuous mathematical description of a population of ants and
simulate numerically their foraging behavior using a system of partial
differential equations of chemotaxis type. We show that this system accurately
reproduces observed foraging behavior, especially spontaneous trail formation
and efficient removal of food sources. We show through numerical experiments
that trail formation is correlated with efficient food removal. Our results
illustrate the emergence of trail formation from simple modeling principles.


  Emerging integrative analysis of genomic and anatomical imaging data which
has not been well developed, provides invaluable information for the holistic
discovery of the genomic structure of disease and has the potential to open a
new avenue for discovering novel disease susceptibility genes which cannot be
identified if they are analyzed separately. A key issue to the success of
imaging and genomic data analysis is how to reduce their dimensions. Most
previous methods for imaging information extraction and RNA-seq data reduction
do not explore imaging spatial information and often ignore gene expression
variation at genomic positional level. To overcome these limitations, we extend
functional principle component analysis from one dimension to two dimension
(2DFPCA) for representing imaging data and develop a multiple functional linear
model (MFLM) in which functional principal scores of images are taken as
multiple quantitative traits and RNA-seq profile across a gene is taken as a
function predictor for assessing the association of gene expression with
images. The developed method has been applied to image and RNA-seq data of
ovarian cancer and KIRC studies. We identified 24 and 84 genes whose
expressions were associated with imaging variations in ovarian cancer and KIRC
studies, respectively. Our results showed that many significantly associated
genes with images were not differentially expressed, but revealed their
morphological and metabolic functions. The results also demonstrated that the
peaks of the estimated regression coefficient function in the MFLM often
allowed the discovery of splicing sites and multiple isoform of gene
expressions.


  Dengue is a vector-borne disease transmitted from an infected human to an
Aedes mosquito, during a blood-meal. Dengue is still a major public health
problem. A model for the disease transmission is presented, composed by human
and mosquitoes compartments. The aim is to simulate the effects of seasonality,
on the vectorial capacity and, consequently, on the disease development. Using
entomological information about the mosquito behavior under different
temperatures and rainfall, simulations are carried out and the repercussions
analyzed. The basic reproduction number of the model is given, as well as a
sensitivity analysis of model's parameters. Finally, an optimal control problem
is proposed and solved, illustrating the difficulty of making a trade-off
between reduction of infected individuals and costs with insecticide.


  Given a group-based Markov model on a tree, one can compute the vertex
representation of a polytope which describes the associated toric variety. The
half-space representation, however, is not easily computable. In the case of
$\mathbb{Z}_2$ or $\mathbb{Z}_2\times\mathbb{Z}_2$, these polytopes have
applications in the field of phylogenetics. We provide a half-space
representation for the $m$-claw tree where $G=\mathbb{Z}_2\times\mathbb{Z}_2$,
which corresponds to the Kimura-3 model of evolution.


  A classic measure of ecological stability describes the tendency of a
community to return to equilibrium after small perturbation. While many
advances show how the network structure of these communities severely
constrains such tendencies, few if any of these advances address one of the
most fundamental properties of network structure: heterogeneity among nodes
with different numbers of links. Here we systematically explore this property
of "degree heterogeneity" and find that its effects on stability systematically
vary with different types of interspecific interactions. Degree heterogeneity
is always destabilizing in ecological networks with both competitive and
mutualistic interactions while its effects on networks of predator-prey
interactions such as food webs depend on prey contiguity, i.e., the extent to
which the species consume an unbroken sequence of prey in community niche
space. Increasing degree heterogeneity stabilizes food webs except those with
the most contiguity. These findings help explain previously unexplained
observations that food webs are highly but not completely contiguous and, more
broadly, deepens our understanding of the stability of complex ecological
networks with important implications for other types of dynamical systems.


  Bacterial resistance to antibiotic treatment is a huge concern: introduction
of any new antibiotic is shortly followed by the emergence of resistant
bacterial isolates in the clinic. This issue is compounded by a severe lack of
new antibiotics reaching the market. The significant rise in clinical
resistance to antibiotics is especially problematic in nosocomial infections,
where already vulnerable patients may fail to respond to treatment, causing
even greater health concern. A recent focus has been on the development of
anti-virulence drugs as a second line of defence in the treatment of
antibiotic-resistant infections. This treatment, which weakens bacteria by
reducing their virulence rather than killing them, should allow infections to
be cleared through the body's natural defence mechanisms. In this way there
should be little to no selective pressure exerted on the organism and, as such,
a predominantly resistant population would be unlikely to emerge. However, much
controversy surrounds this approach with many believing it would not be
powerful enough to clear existing infections, restricting its potential
application to prophylaxis. We have developed a mathematical model that
provides a theoretical framework to reveal the circumstances under which
anti-virulence drugs may or may not be successful. We demonstrate that by
harnessing and combining the advantages of antibiotics with those provided by
anti-virulence drugs, given infection-specific parameters, it is possible to
identify treatment strategies that would efficiently clear bacterial
infections, while preventing the emergence of resistant subpopulations. Our
findings strongly support the continuation of research into anti-virulence
drugs and demonstrate that their applicability may reach beyond infection
prevention.


  The classical method of determining the atomic structure of complex molecules
by analyzing diffraction patterns is currently undergoing drastic developments.
Modern techniques for producing extremely bright and coherent X-ray lasers
allow a beam of streaming particles to be intercepted and hit by an ultrashort
high energy X-ray beam. Through machine learning methods the data thus
collected can be transformed into a three-dimensional volumetric intensity map
of the particle itself. The computational complexity associated with this
problem is very high such that clusters of data parallel accelerators are
required.
  We have implemented a distributed and highly efficient algorithm for
inversion of large collections of diffraction patterns targeting clusters of
hundreds of GPUs. With the expected enormous amount of diffraction data to be
produced in the foreseeable future, this is the required scale to approach real
time processing of data at the beam site. Using both real and synthetic data we
look at the scaling properties of the application and discuss the overall
computational viability of this exciting and novel imaging technique.


  A microtubule (MT) is a hollow tube of approximately 25 nm diameter. The two
ends of the tube are dissimilar and are designated as `plus' and `minus' ends.
Motivated by the collective push and pull exerted by a bundle of MTs during
chromosome segregation in a living cell, we have developed here a much
simplified theoretical model of a bundle of parallel dynamic MTs. The plus-end
of all the MTs in the bundle are permanently attached to a movable `wall' by a
device whose detailed structure is not treated explicitly in our model. The
only requirement is that the device allows polymerization and depolymerization
of each MT at the plus-end. In spite of the absence of external force and
direct lateral interactions between the MTs, the group of polymerizing MTs
attached to the wall create a load force against the group of depolymerizing
MTs and vice-versa; the load against a group is shared equally by the members
of that group. Such indirect interactions among the MTs gives rise to the rich
variety of possible states of collective dynamics that we have identified by
computer simulations of the model in different parameter regimes. The
bi-directional motion of the cargo, caused by the load-dependence of the
polymerization kinetics, is a "proof-of-principle" that the bi-directional
motion of chromosomes before cell division does not necessarily need active
participation of motor proteins.


  Biochemical reactions can happen on different time scales and also the
abundance of species in these reactions can be very different from each other.
Classical approaches, such as deterministic or stochastic approach, fail to
account for or to exploit this multi-scale nature, respectively. In this paper,
we propose a jump-diffusion approximation for multi-scale Markov jump processes
that couples the two modeling approaches. An error bound of the proposed
approximation is derived and used to partition the reactions into fast and slow
sets, where the fast set is simulated by a stochastic differential equation and
the slow set is modeled by a discrete chain. The error bound leads to a very
efficient dynamic partitioning algorithm which has been implemented for several
multi-scale reaction systems. The gain in computational efficiency is
illustrated by a realistically sized model of a signal transduction cascade
coupled to a gene expression dynamics.


  The paper concerns two interacting consumer-resource pairs based on
chemostat-like equations under the assumption that the dynamics of the resource
is considerably slower than that of the consumer. The presence of two different
time scales enables to carry out a fairly complete analysis of the problem.
This is done by treating consumers and resources in the coupled system as
fast-scale and slow-scale variables respectively and subsequently considering
developments in phase planes of these variables, fast and slow, as if they are
independent. When uncoupled, each pair has unique asymptotically stable steady
state and no self-sustained oscillatory behavior (although damped oscillations
about the equilibrium are admitted). When the consumer-resource pairs are
weakly coupled through direct reciprocal inhibition of consumers, the whole
system exhibits self-sustained relaxation oscillations with a period that can
be significantly longer than intrinsic relaxation time of either pair. It is
shown that the model equations adequately describe locally linked
consumer-resource systems of quite different nature: living populations under
interspecific interference competition and lasers coupled via their cavity
losses.


  Experimental and empirical observations on cell metabolism cannot be
understood as a whole without their integration into a consistent systematic
framework. However, the characterization of metabolic flux phenotypes is
typically reduced to the study of a single optimal state, like maximum biomass
yield that is by far the most common assumption. Here we confront optimal
growth solutions to the whole set of feasible flux phenotypes (FFP), which
provides a benchmark to assess the likelihood of optimal and high-growth states
and their agreement with experimental results. In addition, FFP maps are able
to uncover metabolic behaviors, such as aerobic fermentation accompanying
exponential growth on sugars at nutrient excess conditions, that are
unreachable using standard models based on optimality principles. The
information content of the full FFP space provides us with a map to explore and
evaluate metabolic behavior and capabilities, and so it opens new avenues for
biotechnological and biomedical applications.


  An Ebola outbreak of unparalleled size is currently affecting several
countries in West Africa, and international efforts to control the outbreak are
underway. However, the efficacy of these interventions, and their likely impact
on an Ebola epidemic of this size, is unknown. Forecasting and simulation of
these interventions may inform public health efforts. We use existing data from
Liberia and Sierra Leone to parameterize a mathematical model of Ebola and use
this model to forecast the progression of the epidemic, as well as the efficacy
of several interventions, including increased contact tracing, improved
infection control practices, the use of a hypothetical pharmaceutical
intervention to improve survival in hospitalized patients. Model forecasts
until Dec. 31, 2014 show an increasingly severe epidemic with no sign of having
reached a peak. Modeling results suggest that increased contact tracing,
improved infection control, or a combination of the two can have a substantial
impact on the number of Ebola cases, but these interventions are not sufficient
to halt the progress of the epidemic. The hypothetical pharmaceutical
intervention, while impacting mortality, had a smaller effect on the forecasted
trajectory of the epidemic. Near-term, practical interventions to address the
ongoing Ebola epidemic may have a beneficial impact on public health, but they
will not result in the immediate halting, or even obvious slowing of the
epidemic. A long-term commitment of resources and support will be necessary to
address the outbreak.


  Sol argues that innovation propensity is not a specialized adaptation
resulting from targeted selection but an instance of exaptation because
selection cannot act on situations that are only encountered once. In
exaptation, a trait that originally evolved to solve one problem is co-opted to
solve a new problem; thus the trait or traits in question must be necessary and
sufficient to solve the new problem. Sol claims that traits such as persistence
and neophilia, are necessary and sufficient for animal innovation, which is a
matter of trial and error. We suggest that this explanation does not extend to
human innovation, which involves strategy, logic, intuition, and insight, and
requires traits that evolved, not as a byproduct of some other function, but
for the purpose of coming up with adaptive responses to environmental
variability itself. We point to an agent based model that indicates the
feasibility of two such proposed traits: (1) chaining, the ability to construct
complex thoughts from simple ones, and (2) contextual focus, the ability to
shift between convergent and divergent modes of thought. We agree that there is
a sense in which innovation is exaptation--it occurs when an existing object or
behaviour is adapted to new needs or tastes--and refer to a mathematical model
of biological and cultural exaltation. We conclude that much is gained by
comparing and contrasting animal and human innovation.


  A wide range of stochastic processes that model the growth and decline of
populations exhibit a curious dichotomy: with certainty either the population
goes extinct or its size tends to infinity. There is a elegant and classical
theorem that explains why this dichotomy must hold under certain assumptions
concerning the process. In this note, I explore how these assumptions might be
relaxed further in order to obtain the same, or a similar conclusion, and
obtain both positive and negative results.


  The genotype-phenotype map is an essential object in our understanding of
organismal complexity and adaptive properties, determining at once genomic
plasticity and those constraints that may limit the ability of genomes to
attain evolutionary innovations. An exhaustive experimental characterization of
the relationship between genotypes and phenotypes is at present out of reach.
Therefore, several models mimicking that map have been proposed and
investigated, leading to the identification of a number of general features:
genotypes differ in their robustness to mutations, phenotypes are represented
by a broadly varying number of genotypes, and simple point mutations seem to
suffice to navigate the space of genotypes while maintaining a phenotype.
However, most current models address only one level of the map (sequences and
folded structures in RNA or proteins; networks of genes and their dynamical
attractors; sets of chemical reactions and their ability to undergo molecular
catalysis), such that many relevant questions cannot be addressed. Here we
introduce toyLIFE, a multi-level model for the genotype-phenotype map based on
simple genomes and interaction rules from which a complex behavior at upper
levels emerges, remarkably plastic gene regulatory networks and metabolism.
toyLIFE is a tool that permits the investigation of how different levels are
coupled, in particular how and where do mutations affect phenotype or how the
presence of certain metabolites determines the dynamics of toyLIFE gene
regulatory networks. The possibilities of this model are not exhausted by the
results presented in this contribution. It can be easily generalized to
incorporate evolution through mutations that change genome length or through
recombination, to consider gene duplication or deletion, and therefore to
explore further properties of extended genotype-phenotype maps.


  HIF (Hypoxia Inducible Factor) is an oxygen-regulated transcription factor
that mediates the intracellular response to hypoxia in human cells. There is
increasing evidence that cell signaling pathways encode temporal information,
and thus cell fate may be determined by the dynamics of protein levels. We have
developed a mathematical model to describe the transient dynamics of the
HIF-1{\alpha} protein measured in single cells subjected to hypoxic shock. The
essential characteristics of these data are modeled with a system of
differential equations describing the feedback inhibition between HIF-1{\alpha}
and Prolyl Hydroxylases (PHD) oxygen sensors. Heterogeneity in the single-cell
data is accounted for through parameter variation in the model. We previously
identified the PHD2 isoform as the main PHD responsible for controlling the
HIF-1{\alpha} transient response, and make here testable predictions regarding
HIF-1{\alpha} dynamics subject to repetitive hypoxic pulses. The model is
further developed to describe the dynamics of HIF-1{\alpha} in cells cultured
as 3D spheroids, with oxygen dynamics parameterized using experimental
measurements of oxygen within spheroids. We show that the dynamics of
HIF-1{\alpha} and transcriptional targets of HIF-1{\alpha} display a
non-monotone response to the oxygen dynamics. Specifically we demonstrate that
the dynamic transient behavior of HIF-1{\alpha} results in differential
dynamics in transcriptional targets.


  C$_4$ photosynthesis has independently evolved from the ancestral C$_3$
pathway in at least 60 plant lineages, but, as with other complex traits, how
it evolved is unclear. Here we show that the polyphyletic appearance of C$_4$
photosynthesis is associated with diverse and flexible evolutionary paths that
group into four major trajectories. We conducted a meta-analysis of 18 lineages
containing species that use C$_3$, C$_4$, or intermediate C$_3$-C$_4$ forms of
photosynthesis to parameterise a 16-dimensional phenotypic landscape. We then
developed and experimentally verified a novel Bayesian approach based on a
hidden Markov model that predicts how the C$_4$ phenotype evolved. The
alternative evolutionary histories underlying the appearance of C$_4$
photosynthesis were determined by ancestral lineage and initial phenotypic
alterations unrelated to photosynthesis. We conclude that the order of C$_4$
trait acquisition is flexible and driven by non-photosynthetic drivers. This
flexibility will have facilitated the convergent evolution of this complex
trait.


  We study the scheduling problem of a self-replicating factory. We show that
by maintaining a sufficiently large inventory of intermediate metabolites and
catalysts required for self-replication, optimal replication times can be
achieved by a family of random scheduling algorithms that are biochemically
feasible, for which catalysts never idle if they can perform de-novo
bio-synthesis. Optimally scheduled self-replication is facilitated by allowing
several production lines to run in parallel. The excess inventory of catalysts
and substrates decouples these lines, while dynamical balancing tunes average
and variance completion, resulting in an overall universal distribution for the
replication times belonging to the generalized extreme value family. We discuss
biological implications and postulate that bacteria that are tuned by evolution
for fast replication employ this natural scheduling strategy to achieve optimal
asymptotic growth rates by stoichiometrically balancing the amount of work in
progress thus globally controlling the number of parallel basic
self-replicating units within them. Analysis of recently measured data of E.
coli growth in rich media shows data-collapse on a single universal curve
consistent with our prediction, suggesting wild type E. coli optimally schedule
its replication.


  Using the dynamic mean-field approximation of the totally asymmetric simple
exclusion process (TASEP), we investigate the effect of small changes in the
initiation, exit, and elongation rates along the mRNA strand on the steady
state protein translation rate. We focus on two special cases where exact
closed-form expressions for the translation rate sensitivity can be derived. We
discuss the ramifications of our results in the context of functional genomics,
molecular evolution, and synthetic biology.


  Over time, a population acquires neutral genetic substitutions as a
consequence of random drift. A famous result in population genetics asserts
that the rate, $K$, at which these substitutions accumulate in the population
coincides with the mutation rate, $u$, at which they arise in individuals:
$K=u$. This identity enables genetic sequence data to be used as a "molecular
clock" to estimate the timing of evolutionary events. While the molecular clock
is known to be perturbed by selection, it is thought that $K=u$ holds very
generally for neutral evolution. Here we show that asymmetric spatial
population structure can alter the molecular clock rate for neutral mutations,
leading to either $K<u$ or $K>u$. Deviations from $K=u$ occur because mutations
arise unequally at different sites and have different probabilities of fixation
depending on where they arise. If birth rates are uniform across sites, then $K
\leq u$. In general, $K$ can take any value between 0 and $Nu$. Our model can
be applied to a variety of population structures. In one example, we
investigate the accumulation of genetic mutations in the small intestine. In
another application, we analyze over 900 Twitter networks to study the effect
of network topology on the fixation of neutral innovations in social evolution.


  Polymeric filament like type IV Pilus (TFP) can transfer forces in excess of
100pN during their retraction before stalling, powering surface
translocation(twitching). Single TFP level experiments have shown remarkable
nonlinearity in the retraction behavior influenced by the external load as well
as levels of PilT molecular motor protein. This includes reversal of motion
near stall forces when the concentration of the PilT protein is lowered
significantly. In order to explain this behavior, we analyze the coupling of
TFP elasticity and interfacial behavior with PilT kinetics. We model retraction
as reaction controlled and elongation as transport controlled process. The
reaction rates vary with TFP deformation which is modeled as a compound elastic
body consisting of multiple helical strands under axial load. Elongation is
controlled by monomer transport which suffer entrapment due to excess PilT in
the cell periplasm. Our analysis shows excellent agreement with a host of
experimental observations and we present a possible biophysical relevance of
model parameters through a mechano-chemical stall force map


  This paper develops a deterministic differential equations model that
captures the H7N9 virus transmission from live poultry to human via
poultry-human contacts in live poultry markets (LPMs). The virus circulation
among live poultry, which happens but is hard to be detected (since
contaminated live poultry appear to be asymptomatic), is also incorporated in
the model. The time-dependent contact rate between human and live poultry based
on LPMs closing information can be estimated. From data of LPMs closing news,
the contact rate function can be easily estimated. This model could serve as a
rational basis for public health authorities to evaluate the effectiveness of
LPM closing, as well as other interventions according to simple modifications.
Without data about daily cases, I also provide suggestions for some of the
basic parameters that would be a useful fitness parameter set for future
simulation.


  Experiments show that macromolecular crowding modestly reduces the size of
intrinsically disordered proteins (IDPs) even at volume fraction ($\phi$)
similar to that in the cytosol whereas DNA undergoes a coil-to-globule
transition at very small $\phi$. We show using a combination of scaling
arguments and simulations that the polymer size $\overline{R}_g(\phi)$ depends
on $x = \overline{R}_g(0)/D$ where $D$ is the $\phi$-dependent distance between
the crowders. If $x\lesssim \mathcal{O}(1)$, there is only a small decrease in
$\overline{R}_g(\phi)$ as $\phi$ increases. When $x\gg \mathcal{O}(1)$, a
cooperative coil-to-globule transition is induced. Our theory quantitatively
explains a number of experiments.


  Microarray technology is still an important way to assess gene expression in
molecular biology, mainly because it measures expression profiles for thousands
of genes simultaneously, what makes this technology a good option for some
studies focused on systems biology. One of its main problem is complexity of
experimental procedure, presenting several sources of variability, hindering
statistical modeling. So far, there is no standard protocol for generation and
evaluation of microarray data. To mitigate the analysis process this paper
presents an R package, named maigesPack, that helps with data organization.
Besides that, it makes data analysis process more robust, reliable and
reproducible. Also, maigesPack aggregates several data analysis procedures
reported in literature, for instance: cluster analysis, differential
expression, supervised classifiers, relevance networks and functional
classification of gene groups or gene networks.


  Alignment-based sequence similarity searches, while accurate for some type of
sequences, can produce incorrect results when used on more divergent but
functionally related sequences that have undergone the sequence rearrangements
observed in many bacterial and viral genomes. Here, we propose a classification
model that exploits the complementary nature of alignment-based and
alignment-free similarity measures with the aim to improve the accuracy with
which DNA and protein sequences are characterized. Our model classifies
sequences using a combined sequence similarity score calculated by adaptively
weighting the contribution of different sequence similarity measures. Weights
are determined independently for each sequence in the test set and reflect the
discriminatory ability of individual similarity measures in the training set.
Since the similarity between some sequences is determined more accurately with
one type of measure rather than another, our classifier allows different sets
of weights to be associated with different sequences. Using five different
similarity measures we show that our model significantly improves the
classification accuracy over the current composition and alignment based
models, when predicting the taxonomic lineage for both short viral sequence
fragments and complete viral sequences. We also show that our model can be used
effectively for the classification of reads from a real metagenome dataset as
well as protein sequences.


  We propose an improved prediction method of the tertiary structures of
$\alpha$-helical membrane proteins based on the replica-exchange method by
taking into account helix deformations. Our method allows wide applications
because transmembrane helices of native membrane proteins are often distorted.
In order to test the effectiveness of the present method, we applied it to the
structure predictions of glycophorin A and phospholamban. The results were in
accord with experiments.


  Constraints can affect dramatically the behavior of diffusion processes.
Recently, we analyzed a natural and a technological system and reported that
they perform diffusion-like discrete steps displaying a peculiar constraint,
whereby the increments of the diffusing variable are subject to
configuration-dependent bounds. This work explores theoretically some of the
revealing landmarks of such phenomenology, termed "soft bound". At long times,
the system reaches a steady state irreversibly (i.e., violating detailed
balance), characterized by a skewed "shoulder" in the density distribution, and
by a net local probability flux, which has entropic origin. The largest point
in the support of the distribution follows a saturating dynamics, expressed by
the Gompertz law, in line with empirical observations. Finally, we propose a
generic allometric scaling for the origin of soft bounds. These findings shed
light on the impact on a system of such "scaling" constraint and on its
possible generating mechanisms.


  Prion diseases are invariably fatal neurodegenerative diseases that affect
humans and animals. Unlike most other amyloid forming neurodegenerative
diseases, these can be highly infectious. Prion diseases occur in a variety of
species. They include the fatal human neurodegenerative diseases
Creutzfeldt-Jakob Disease (CJD), Fatal Familial Insomnia (FFI),
Gerstmann-Straussler-Scheinker syndrome (GSS), Kuru, the bovine spongiform
encephalopathy (BSE or 'mad-cow' disease) in cattle, the chronic wasting
disease (CWD) in deer and elk, and scrapie in sheep and goats, etc.
Transmission across the species barrier to humans, especially in the case of
BSE in Europe, CWD in North America, and variant CJDs (vCJDs) in young people
of UK, is a major public health concern. Fortunately, scientists reported that
the (central) hydrophobic region of prion proteins (PrP) controls the formation
of diseased prions. This article gives a detailed survey on PrP hydrophobic
region and does molecular dynamics studies of human PrP(110-136) to confirm
some findings from the survey. The structural bioinformatics presented in this
article can be helpful as a reference in three-dimensional images for
laboratory experimental works to study PrP hydrophobic region.


  Many peak detection algorithms have been proposed for ChIP-seq data analysis,
but it is not obvious which method and what parameters are optimal for any
given data set. In contrast, peaks can easily be located by visual inspection
of profile data on a genome browser. We thus propose a supervised machine
learning approach to ChIP-seq data analysis, using annotated regions that
encode an expert's qualitative judgments about which regions contain or do not
contain peaks. The main idea is to manually annotate a small subset of the
genome, and then learn a model that makes consistent predictions on the rest of
the genome. We show how our method can be used to quantitatively calibrate and
benchmark the performance of peak detection algorithms on specific data sets.
We compare several peak detectors on 7 annotated region data sets, consisting
of 2 histone marks, 4 expert annotators, and several different cell types. In
these data the macs algorithm was best for a narrow peak histone profile
(H3K4me3) while the hmcan.broad algorithm was best for a broad histone profile
(H3K36me3). Our benchmark annotated region data sets can be downloaded from a
public website, and there is an R package for computing the annotation error on
GitHub.


  Due to the recent evolution of sequencing techniques, the number of available
genomes is rising steadily, leading to the possibility to make large scale
genomic comparison between sets of close species. An interesting question to
answer is: what is the common functionality genes of a collection of species,
or conversely, to determine what is specific to a given species when compared
to other ones belonging in the same genus, family, etc. Investigating such
problem means to find both core and pan genomes of a collection of species,
\textit{i.e.}, genes in common to all the species vs. the set of all genes in
all species under consideration. However, obtaining trustworthy core and pan
genomes is not an easy task, leading to a large amount of computation, and
requiring a rigorous methodology. Surprisingly, as far as we know, this
methodology in finding core and pan genomes has not really been deeply
investigated. This research work tries to fill this gap by focusing only on
chloroplastic genomes, whose reasonable sizes allow a deep study. To achieve
this goal, a collection of 99 chloroplasts are considered in this article. Two
methodologies have been investigated, respectively based on sequence
similarities and genes names taken from annotation tools. The obtained results
will finally be evaluated in terms of biological relevance.


  E-values have been the dominant statistic for protein sequence analysis for
the past two decades: from identifying statistically significant local sequence
alignments to evaluating matches to hidden Markov models describing protein
domain families. Here we formally show that for "stratified" multiple
hypothesis testing problems, controlling the local False Discovery Rate (lFDR)
per stratum, or partition, yields the most predictions across the data at any
given threshold on the FDR or E-value over all strata combined. For the
important problem of protein domain prediction, a key step in characterizing
protein structure, function and evolution, we show that stratifying statistical
tests by domain family yields excellent results. We develop the first
FDR-estimating algorithms for domain prediction, and evaluate how well
thresholds based on q-values, E-values and lFDRs perform in domain prediction
using five complementary approaches for estimating empirical FDRs in this
context. We show that stratified q-value thresholds substantially outperform
E-values. Contradicting our theoretical results, q-values also outperform
lFDRs; however, our tests reveal a small but coherent subset of domain
families, biased towards models for specific repetitive patterns, for which
FDRs are greatly underestimated due to weaknesses in random sequence models.
Usage of lFDR thresholds outperform q-values for the remaining families, which
have as-expected noise, suggesting that further improvements in domain
predictions can be achieved with improved modeling of random sequences.
Overall, our theoretical and empirical findings suggest that the use of
stratified q-values and lFDRs could result in improvements in a host of
structured multiple hypothesis testing problems arising in bioinformatics,
including genome-wide association studies, orthology prediction, motif
scanning, and multi-microarray analyses.


  The efficient calculation of rare-event kinetics in complex dynamical
systems, such as the rate and pathways of ligand dissociation from a protein,
is a generally unsolved problem. Markov state models can systematically
integrate ensembles of short simulations and thus effectively parallelize the
computational effort, but the rare events of interest still need to be
spontaneously sampled in the data. Enhanced sampling approaches, such as
parallel tempering or umbrella sampling, can accelerate the computation of
equilibrium expectations massively - but sacrifice the ability to compute
dynamical expectations. In this work we establish a principle to combine
knowledge of the equilibrium distribution with kinetics from fast "downhill"
relaxation trajectories using reversible Markov models. This approach is
general as it does not invoke any specific dynamical model, and can provide
accurate estimates of the rare event kinetics. Large gains in sampling
efficiency can be achieved whenever one direction of the process occurs more
rapid than its reverse, making the approach especially attractive for downhill
processes such as folding and binding in biomolecules.


  A transversal topic of my research has been the development and application
of computational methods for DNA sequence analysis. The methods I have been
developing aim at improving our understanding of the regulation processes
happening in normal and cancer cells. This topic connects together the projects
presented in this thesis. Two chapters of the thesis represent major areas of
my research interests: (1) methods for deciphering transcriptional regulation
and their application to answer specific biological questions, and (2) methods
to study the genome structure and their application in cancer studies. The
first chapter predominantly focuses on transcriptional regulation. Here I
describe my contribution to the development of methodology for the discovery of
transcription factor binding sites and the positioning of histone proteins. I
also explain how sequence analysis, in combination with gene expression data,
can allow the identification of direct target genes of a transcription factor
under study, as well as the physical mechanisms of its action. As two examples,
I provide the results of my study of transcriptional regulation by (i)
oncogenic protein EWS-FLI1 in Ewing sarcoma and (ii) oncogenic transcription
factor Spi-1/PU.1 in erythroleukemia. In the second chapter, I describe the
sequence analysis methods aimed at the identification of the genomic
rearrangements in species with existing reference genome. I explain how the
developed methodology can be applied to detect the structure of cancer genomes.
I provide an example of how such an analysis of tumor genomes can result in a
discovery of a new phenomenon: chromothripsis, when hundreds of rearrangements
occur in a single cellular catastrophe. The thesis is concluded by listing the
major challenges in high-throughput sequencing analysis. I also discuss the
current top questions demanding the integration of sequencing data.


  The exceptional reactivity of animal collectives to predatory attacks is
thought to be due to rapid, but local, transfer of information between group
members. These groups turn together in unison and produce escape waves.
However, it is not clear how escape waves are created from local interactions,
nor is it understood how these patterns are shaped by natural selection. By
startling schools of fish with a simulated attack in an experimental arena, we
demonstrate that changes in the direction and speed by a small percentage of
individuals that detect the danger initiate an escape wave. This escape wave
consists of a densely packed band of individuals that causes other school
members to change direction. In the majority of cases this wave passes through
the entire group. We use a simulation model to demonstrate that this mechanism
can, through local interactions alone, produce arbitrarily large escape waves.
In the model, when we set the group density to that seen in real fish schools,
we find that the risk to the members at the edge of the group is roughly equal
to the risk of those within the group. Our experiments and modelling results
provide a plausible explanation for how escape waves propagate in Nature
without centralised control.


  Direct detector device (DDD) cameras have revolutionized single particle
electron cryomicroscopy (cryo-EM). In addition to an improved camera detective
quantum efficiency, acquisition of DDD movies allows for correction of movement
of the specimen, due both to instabilities in the microscope specimen stage and
electron beam-induced movement. Unlike specimen stage drift, beam-induced
movement is not always homogeneous within an image. Local correlation in the
trajectories of nearby particles suggests that beam-induced motion is due to
deformation of the ice layer. Algorithms have already been described that can
correct movement for large regions of frames and for > 1 MDa protein particles.
Another algorithm allows individual < 1 MDa protein particle trajectories to be
estimated, but requires rolling averages to be calculated from frames and fits
linear trajectories for particles. Here we describe an algorithm that allows
for individual < 1 MDa particle images to be aligned without frame averaging or
linear trajectories. The algorithm maximizes the overall correlation of the
shifted frames with the sum of the shifted frames. The optimum in this single
objective function is found efficiently by making use of analytically
calculated derivatives of the function. To smooth estimates of particle
trajectories, rapid changes in particle positions between frames are penalized
in the objective function and weighted averaging of nearby trajectories ensures
local correlation in trajectories. This individual particle motion correction,
in combination with weighting of Fourier components to account for increasing
radiation damage in later frames, can be used to improve 3-D maps from single
particle cryo-EM.


  Above a critical concentration a wide variety of peptides and proteins
self-assemble into amyloid fibrils which entangle to form percolating networks
called hydrogels. Such hydrogels have important applications as biomaterials
and in nanotechnology, but their applicability often depends on their
mechanical properties for which we currently have no predictive capability.
Here we use a peptide model to simulate the formation of amyloid fibril
networks, and couple these to elastic network theory to determine their
mechanical properties. The simulations reveal that the time-dependence of
morphological quantities characterizing the network length scales can be
collapsed onto master curves by using a time scaling function that depends on
the interaction parameter between the peptides. The same scaling function is
used to unveil a universal, non-monotonic dependence of the shear modulus with
time. The obtained insight into the structure-function relationship between the
peptide building blocks, network morphology and network mechanical properties
can aid experimentalists to design amyloid fibril networks with tailored
mechanical properties.


  The dynamics of ecosystem collapse are fundamental to determining how and why
biological communities change through time, as well as the potential effects of
extinctions on ecosystems. Here we integrate depictions of mammals from
Egyptian antiquity with direct lines of paleontological and archeological
evidence to infer local extinctions and community dynamics over a 6000-year
span. The unprecedented temporal resolution of this data set enables
examination of how the tandem effects of human population growth and climate
change can disrupt mammalian communities. We show that the extinctions of
mammals in Egypt were nonrandom, and that destabilizing changes in community
composition coincided with abrupt aridification events and the attendant
collapses of some complex societies. We also show that the roles of species in
a community can change over time, and that persistence is predicted by measures
of species sensitivity, a function of local dynamic stability. Our study is the
first high-resolution analysis of the ecological impacts of environmental
change on predator-prey networks over millennial timescales, and sheds light on
the historical events that have shaped modern animal communities.


  It is known that a gamma ray burst (GRB) originating near the Earth could be
devastating to life. The mechanism of ozone depletion and subsequent increased
UVB exposure is the primary risk, but models also show increased nitrification
culminating in nitric acid rainout. These effects are also expected after
nearby supernovae and extreme solar proton events. In this work we considered
specifically whether the increased nitric acid rainout from such events is a
threat to modern terrestrial ecosystems. We also considered its potential
benefit to early terrestrial Paleozoic ecosystems. We used established critical
loads for nitrogen deposition in ecoregions of Europe and the US and compared
them with previously predicted values of nitric acid rainout from a typical GRB
within our galaxy. The predicted rainout was found to be too low to harm modern
ecosystems, however, it is large compared with probable nitrate flux onto land
prior to the invasion of plants. We suggest that this flux may have contributed
nutrients to this invasion if, as hypothesized, the end-Ordovician extinction
event were initiated by a GRB or other ionizing radiation event.


  Establishing a distance between genomes is a significant problem in
computational genomics, because its solution can be used to establish
evolutionary relationships including phylogeny.
  The "double cut and join" (DCJ) model of chromosomal rearrangement proposed
by Yancopoulos et al. has received attention as it can model inversions,
translocations, fusion and fission on a multichromosomal genome that may
contain both linear and circular chromosomes. In this paper, we realize the DCJ
operator as a group action on the space of multichromosomal genomes. We study
this group action, deriving some properties of the group and finding
group-theoretic analogues for the key results in the DCJ theory.


  Tumor samples are heterogeneous. They consist of different subclones that are
characterized by differences in DNA nucleotide sequences and copy numbers on
multiple loci. Heterogeneity can be measured through the identification of the
subclonal copy number and sequence at a selected set of loci. Understanding
that the accurate identification of variant allele fractions greatly depends on
a precise determination of copy numbers, we develop a Bayesian feature
allocation model for jointly calling subclonal copy numbers and the
corresponding allele sequences for the same loci. The proposed method utilizes
three random matrices, L, Z and w to represent subclonal copy numbers (L),
numbers of subclonal variant alleles (Z) and cellular fractions of subclones in
samples (w), respectively. The unknown number of subclones implies a random
number of columns for these matrices. We use next-generation sequencing data to
estimate the subclonal structures through inference on these three matrices.
Using simulation studies and a real data analysis, we demonstrate how posterior
inference on the subclonal structure is enhanced with the joint modeling of
both structure and sequencing variants on subclonal genomes. Software is
available at http://compgenome.org/BayClone2.


  MEGAHIT is a NGS de novo assembler for assembling large and complex
metagenomics data in a time- and cost-efficient manner. It finished assembling
a soil metagenomics dataset with 252Gbps in 44.1 hours and 99.6 hours on a
single computing node with and without a GPU, respectively. MEGAHIT assembles
the data as a whole, i.e., it avoids pre-processing like partitioning and
normalization, which might compromise on result integrity. MEGAHIT generates 3
times larger assembly, with longer contig N50 and average contig length than
the previous assembly. 55.8% of the reads were aligned to the assembly, which
is 4 times higher than the previous. The source code of MEGAHIT is freely
available at https://github.com/voutcn/megahit under GPLv3 license.


  Populations of isogenic embryonic stem cells or clonal bacteria often exhibit
extensive phenotypic heterogeneity which arises from stochastic intrinsic
dynamics of cells. The internal state of the cell can be transmitted
epigenetically in cell division, leading to correlations in the phenotypic
states of cells related by descent. Therefore, a phenotypic snapshot of a
collection of cells with known genealogical structure, contains information on
phenotypic dynamics. Here we use a model of phenotypic dynamics on a
genealogical tree to define an inference method which allows to extract an
approximate probabilistic description of phenotypic dynamics based on measured
correlations as a function of the degree of kinship. The approach is tested and
validated on the example of Pyoverdine dynamics in P. aeruginosa colonies.
Interestingly, we find that correlations among pairs and triples of distant
relatives have a simple but non-trivial structure indicating that observed
phenotypic dynamics on the genealogical tree is approximately conformal - a
symmetry characteristic of critical behavior in physical systems. Proposed
inference method is sufficiently general to be applied in any system where
lineage information is available.


  Epidemic spreading phenomena are ubiquitous in nature and society. Examples
include the spreading of diseases, information, and computer viruses. Epidemics
can spread by local spreading, where infected nodes can only infect a limited
set of direct target nodes and global spreading, where an infected node can
infect every other node. In reality, many epidemics spread using a hybrid
mixture of both types of spreading. In this study we develop a theoretical
framework for studying hybrid epidemics, and examine the optimum balance
between spreading mechanisms in terms of achieving the maximum outbreak size.
We show the existence of critically hybrid epidemics where neither spreading
mechanism alone can cause a noticeable spread but a combination of the two
spreading mechanisms would produce an enormous outbreak. Our results provide
new strategies for maximising beneficial epidemics and estimating the worst
outcome of damaging hybrid epidemics.


  To analyze high-dimensional systems, many fields in science and engineering
rely on high-level descriptions, sometimes called "macrostates,"
"coarse-grainings," or "effective theories". Examples of such descriptions
include the thermodynamic properties of a large collection of point particles
undergoing reversible dynamics, the variables in a macroeconomic model
describing the individuals that participate in an economy, and the summary
state of a cell composed of a large set of biochemical networks.
  Often these high-level descriptions are constructed without considering the
ultimate reason for needing them in the first place. Here, we formalize and
quantify one such purpose: the need to predict observables of interest
concerning the high-dimensional system with as high accuracy as possible, while
minimizing the computational cost of doing so. The resulting State Space
Compression (SSC) framework provides a guide for how to solve for the {optimal}
high-level description of a given dynamical system, rather than constructing it
based on human intuition alone.
  In this preliminary report, we introduce SSC, and illustrate it with several
information-theoretic quantifications of "accuracy", all with different
implications for the optimal compression. We also discuss some other possible
applications of SSC beyond the goal of accurate prediction. These include SSC
as a measure of the complexity of a dynamical system, and as a way to quantify
information flow between the scales of a system.


  Membrane deformation inside living cells is crucial for the proper shaping of
various intracellular organelles and is necessary during the fission/fusion
processes that allow membrane recycling and transport (e.g. endocytosis).
Proteins that induce membrane curvature play a key role in such processes,
mostly by adsorbing to the membrane and forming a scaffold that deforms the
membrane according to the curvature of the proteins. In this paper we explore
the possibility of membrane tube destabilisation through a pearling mechanism
enabled by the combined effects of the adsorbed curved proteins and the actin
polymerization they may recruit. The pearling instability can furthermore serve
as the initiation for fission of the tube into vesicles. We find that adsorbed
proteins are more likely to stabilise the tubes, while the actin polymerization
can provide the additional constrictive force needed for the robust
instability. We discuss the relevance of the theoretical results to in-vivo and
in-vitro experiments.


  Reproductive success and survival are influenced by wealth in human
populations. Wealth is transmitted to offsprings and strategies of transmission
vary over time and among populations, the main variation being how equally
wealth is transmitted to children. Here we propose a model where we simulate
both the dynamics of wealth in a population and the evolution of a trait that
determines how wealth is transmitted from parents to offspring, in a darwinian
context.


  Sirenia are marine mammals that colonized the European shores up to 2.7 Ma.
Their biodiversity evolution follows the climate evolution of the Cenozoic.
However, several climate events, as well as the global climate trend of this
Era are still struggling to be understood. When considering only Earth
processes, the climate evolution of the Cenozoic is hard to understand. If the
galactic environment is taken into account, some of these climate events, as
well as the global climate trend, became more easily understood. The Milky Way,
through Supernovae, may bring some answers to why Cenozoic climate had this
evolution. With the assumption that SN can induced changes in Earth climate in
long time scales, Sirenia disappearance from Europe would be a side effect of
this process.


  In early HIV infection, the virus population escapes from multiple CD8+ cell
responses. The later an escape mutation emerges, the slower it outgrows its
competition, i. e. the escape rate is lower. This pattern could indicate that
the strength of the CD8+ cell responses is waning, or that later viral escape
mutants carry a larger fitness cost. In this paper, we investigate whether the
pattern of decreasing escape rate could also be caused by genetic interference
among different escape strains. To this end, we developed a mathematical
multi-epitope model of HIV dynamics, which incorporates stochastic effects,
recombination and mutation. We used cumulative linkage disequilibrium measures
to quantify the amount of interference. We found that nearly synchronous,
similarly strong immune responses in two-locus systems enhance the generation
of genetic interference. This effect, combined with densely spaced sampling
times at the beginning of infection, leads to decreasing successive escape rate
estimates, even when there were no selection differences among alleles. These
predictions are supported by experimental data from one HIV-infected patient.
Thus, interference could explain why later escapes are slower. Considering
escape mutations in isolation, neglecting their genetic linkage, conceals the
underlying haplotype dynamics and can affect the estimation of the selective
pressure exerted by CD8+ cells. In systems in which multiple escape mutations
appear, the occurrence of interference dynamics should be assessed by measuring
the linkage between different escape mutations.


  We discuss a theoretical model for bidirectional cargo transport in
biological cells, which is driven by teams of molecular motors and subject to
thermal fluctuations. The model describes explicitly the directed motion of the
molecular motors on the filament. The motor-cargo coupling is implemented via
linear springs. By means of extensive Monte Carlo simulations we show that the
model describes the experimentally observed regimes of anomalous diffusion,
i.e. subdiffusive behavior at short times followed by superdiffusion at
intermediate times. The model results indicate that subdiffuse regime is
induced by thermal fluctuations while the superdiffusive motion is generated by
correlations of the motors' activity. We also tested the efficiency of
bidirectional cargo transport in crowded areas by measuring its ability to pass
barriers with increased viscosity. Our results show a remarkable gain of
efficiency for high viscosities.


  Epigenetic aberrations have profound effects on phenotypic output. Genome
wide methylation alterations are inheritable to pass down the aberrations
through multiple generations. We developed a statistical method, Genome-wide
Identification of Significant Methylation Alteration, GISAIM, to study the
significant transgenerational methylation changes. GISAIM finds the significant
methylation aberrations that are inherited through multiple generations. In a
concrete biological study, we investigated whether exposing pregnant rats (F0)
to a high fat (HF) diet throughout pregnancy or ethinyl estradiol
(EE2)-supplemented diet during gestation days 14 20 affects carcinogen-induced
mammary cancer risk in daughters (F1), granddaughters (F2) and
great-granddaughters (F3). Mammary tumorigenesis was higher in daughters and
granddaughters of HF rat dams, and in daughters, granddaughters and
great-granddaughters of EE2 rat dams. Outcross experiments showed that
increased mammary cancer risk was transmitted to HF granddaughters equally
through the female or male germlines, but is only transmitted to EE2
granddaughters through the female germline. Transgenerational effect on mammary
cancer risk was associated with increased expression of DNA methyltransferases,
and across all three EE2 generations hypo or hyper methylation of the same 375
gene promoter regions in their mammary glands. Our study shows that maternal
dietary estrogenic exposures during pregnancy can increase breast cancer risk
in multiple generations of offspring, and the increase in risk may be inherited
through non-genetic means, possibly involving DNA methylation.


  We study the ensemble velocity of non-processive motor proteins, described
with multiple chemical states. In particular, we discuss the velocity as a
function of ATP concentration. Even a simple model which neglects the
strain-dependence of transition rates, reverse transition rates and
nonlinearities in the elasticity can show interesting functional dependencies,
which deviate significantly from the frequently assumed Michaelis-Menten form.
We discuss how the oder of events in the duty cycle can be inferred from the
measured dependence. The model also predicts the possibility of velocity
reversal at a certain ATP concentration if the duty cycle contains several
conformational changes of opposite directionalities.


  Adaptation to one environment often results in fitness gains and losses in
other conditions. To characterize how these consequences of adaptation depend
on the physical similarity between environments, we evolved 180 populations of
Saccharomyces cerevisiae at different degrees of stress induced by either salt,
temperature, pH, or glucose depletion. We measure how the fitness of clones
adapted to each environment depends on the intensity of the corresponding type
of stress. We find that clones evolved in a given type and intensity of stress
tend to gain fitness in other similar intensities of that stress, and lose
progressively more fitness in more physically dissimilar environments. These
fitness trade-offs are asymmetric: adaptation to permissive conditions incurred
a smaller trade-off in stressful conditions than vice versa. We also find that
fitnesses of clones are highly correlated across similar intensities of stress,
but these correlations decay towards zero in more dissimilar environments. To
interpret these results, we introduce the concept of a joint distribution of
fitness effects of new mutations in multiple environments (the JDFE), which
describes the probability that a given mutation has particular fitness effects
in some set of conditions. We find that our observations are consistent with
JDFEs that are highly correlated between physically similar environments, and
that become less correlated and asymmetric as the environments become more
dissimilar. The JDFE provides a framework for quantifying evolutionary
similarity between conditions, and forms a useful basis for theoretical work
aimed at predicting the outcomes of evolution in fluctuating environments.


  Dengue is a vector-borne disease and 40% of world population is at risk.
Dengue transcends international borders and can be found in tropical and
subtropical regions around the world, predominantly in urban and semi-urban
areas. A model for dengue disease transmission, composed by mutually-exclusive
compartments representing the human and vector dynamics, is presented in this
study. The data is from Madeira, a Portuguese island, where an unprecedented
outbreak was detected on October 2012. The aim of this work is to simulate the
repercussions of the control measures in the fight of the disease.


  According to the World Health Organization, dengue fever is the most
important mosquito-borne disease of humans, and it is currently estimated that
there may be 50 - 100 million yearly dengue infections worldwide. For the
purpose to provide new techniques to public health policies in course, we
introduce a predictive non-linear population dynamics model to describe the
population size of four stages of the development of \emph{Aedes aegypti},
having the coefficients set to be dependent on the rainfall index data. In
spite of the population dynamics of the \emph{Ae. aegypti} be mainly ruled by
the rainfall regime, most models are dedicated exclusively to effects of
temperature and only few models are dedicated to influence of rainfall. Vector
control actions are also implemented in many periods of the year in order to
compare relative efficiency of public health policies. The analysis of
equilibrium and stability was performed. Field rainfall time series data from
the City of Lavras (Minas Gerais, Brazil) was used for the model evaluation.
The model was validated in a comparison with experimental mosquito abundance
data acquired by field health agents. We evaluated and validated an
entomological conjecture that claims that control actions should be performed
during the dry season, instead of the common procedure adopted by vector
control programs, in which those are mainly applied in the rainy season.


  We propose a network metric, edge proximity, ${\cal P}_e$, which demonstrates
the importance of specific edges in a network, hitherto not captured by
existing network metrics. The effects of removing edges with high ${\cal P}_e$
might initially seem inconspicuous but are eventually shown to be very harmful
for networks. Compared to existing strategies, the removal of edges by ${\cal
P}_e$ leads to a remarkable increase in the diameter and average shortest path
length in undirected real and random networks till the first disconnection and
well beyond. ${\cal P}_e$ can be consistently used to rupture the network into
two nearly equal parts, thus presenting a very potent strategy to greatly harm
a network. Targeting by ${\cal P}_e$ causes notable efficiency loss in U.S. and
European power grid networks. ${\cal P}_e$ identifies proteins with essential
cellular functions in protein-protein interaction networks. It pinpoints
regulatory neural connections and important portions of the neural and brain
networks, respectively. Energy flow interactions identified by ${\cal P}_e$
form the backbone of long food web chains. Finally, we scrutinize the potential
of ${\cal P}_e$ in edge controllability dynamics of directed networks.


  The connectome, or the entire connectivity of a neural system represented by
network, ranges various scales from synaptic connections between individual
neurons to fibre tract connections between brain regions. Although the
modularity they commonly show has been extensively studied, it is unclear
whether connection specificity of such networks can already be fully explained
by the modularity alone. To answer this question, we study two networks, the
neuronal network of C. elegans and the fibre tract network of human brains
yielded through diffusion spectrum imaging (DSI). We compare them to their
respective benchmark networks with varying modularities, which are generated by
link swapping to have desired modularity values but otherwise maximally random.
We find several network properties that are specific to the neural networks and
cannot be fully explained by the modularity alone. First, the clustering
coefficient and the characteristic path length of C. elegans and human
connectomes are both higher than those of the benchmark networks with similar
modularity. High clustering coefficient indicates efficient local information
distribution and high characteristic path length suggests reduced global
integration. Second, the total wiring length is smaller than for the
alternative configurations with similar modularity. This is due to lower
dispersion of connections, which means each neuron in C. elegans connectome or
each region of interest (ROI) in human connectome reaches fewer ganglia or
cortical areas, respectively. Third, both neural networks show lower
algorithmic entropy compared to the alternative arrangements. This implies that
fewer rules are needed to encode for the organisation of neural systems.


  We present a coarse-grained two dimensional mechanical model for the
microtubule-tau bundles in neuronal axons in which we remove taus, as can
happen in various neurodegenerative conditions such as Alzheimer's disease,
tauopathies, and chronic traumatic encephalopathy. Our simplified model
includes (i) taus modeled as entropic springs between microtubules, (ii)
removal of taus from the bundles due to phosphorylation, and (iii) a possible
depletion force between microtubules due to these dissociated phosphorylated
taus. We equilibrate upon tau removal using steepest descent relaxation. In the
absence of the depletion force, the transverse rigidity to radial compression
of the bundle falls to zero at about 60% tau occupancy, in agreement with
standard percolation theory results. However, with the attractive depletion
force, spring removal leads to a first order collapse of the bundles over a
wide range of tau occupancies for physiologically realizable conditions. While
our simplest calculations assume a constant concentration of microtubule
intercalants to mediate the depletion force, including a dependence that is
linear in the detached taus yields the same collapse. Applying percolation
theory to removal of taus at microtubule tips, which are likely to be the
protective sites against dynamic instability, we argue that the microtubule
instability can only obtain at low tau occupancy, from 0.06-0.30 depending upon
the tau coordination at the microtubule tips. Hence, the collapse we discover
is likely to be more robust over a wide range of tau occupancies than the
dynamic instability. We suggest in vitro tests of our predicted collapse.


  With novel developments in sequencing technologies, time-sampled data are
becoming more available and accessible. Naturally, there have been efforts in
parallel to infer population genetic parameters from these datasets. Here, we
compare and analyze four recent approaches based on the Wright-Fisher model for
inferring selection coefficients (s) given effective population size (Ne), with
simulated temporal datasets. Furthermore, we demonstrate the advantage of a
recently proposed ABC-based method that is able to correctly infer genome-wide
average Ne from time-serial data, which is then set as a prior for inferring
per-site selection coefficients accurately and precisely. We implement this ABC
method in a new software and apply it to a classical time-serial dataset of the
medionigra genotype in the moth Panaxia dominula. We show that a recessive
lethal model is the best explanation for the observed variation in allele
frequency by implementing an estimator of the dominance ratio (h).


  Certain sequences of peptoid polymers (synthetic analogs of peptides)
assemble into bilayer nanosheets via a nonequilibrium assembly pathway of
adsorption, compression, and collapse at an air-water interface. As with other
large-scale dynamic processes in biology and materials science, understanding
the details of this supramolecular assembly process requires a modeling
approach that captures behavior on a wide range of length and time scales, from
those on which individual sidechains fluctuate to those on which assemblies of
polymers evolve. Here we demonstrate that a new coarse-grained modeling
approach is accurate and computationally efficient enough to do so. Our
approach uses only a minimal number of coarse-grained sites, but retains
independently fluctuating orientational degrees of freedom for each site. These
orientational degrees of freedom allow us to accurately parameterize both
bonded and nonbonded interactions, and to generate all-atom configurations with
sufficient accuracy to perform atomic scattering calculations and to interface
with all-atom simulations. We have used this approach to reproduce all
available experimental X-ray scattering spectra (for stacked nanosheets, and
for peptoids adsorbed at air-water interfaces and in solution), in order to
resolve the microscopic, real-space structures responsible for these
Fourier-space features. By interfacing with all-atom simulations, we have also
laid the foundations for future multiscale simulations of sequence-specific
polymers that communicate in both directions across scales.


  We introduce a method to bring nearly atomistic resolution to coarse-grained
models, and we apply the method to proteins. Using a small number of
coarse-grained sites (about one per eight atoms) but assigning an independent
three-dimensional orientation to each site, we preferentially integrate out
stiff degrees of freedom (bond lengths and angles, as well as dihedral angles
in rings) that are accurately approximated by their average values, while
retaining soft degrees of freedom (unconstrained dihedral angles) mostly
responsible for conformational variability. We demonstrate that our scheme
retains nearly atomistic resolution by mapping all experimental protein
configurations in the Protein Data Bank onto coarse-grained configurations,
then analytically backmapping those configurations back to all-atom
configurations. This roundtrip mapping throws away all information associated
with the eliminated (stiff) degrees of freedom except for their average values,
which we use to construct optimal backmapping functions. Despite the 4:1
reduction in the number of degrees of freedom, we find that heavy atoms move
only 0.051 angstroms on average during the roundtrip mapping, while hydrogens
move 0.179 angstroms on average, an unprecedented combination of efficiency and
accuracy among coarse-grained protein models. We discuss the advantages of such
a high-resolution model for parameterizing effective interactions and
accurately calculating observables through direct or multiscale simulations.


  Percolation offers acknowledged models of random media when the relevant
medium characteristics can be described as a binary feature. However, when
considering habitat modeling in ecology, a natural constraint comes from
nearest-neighbor correlations between the suitable/unsuitable states of the
spatial units forming the habitat. Such constraints are also relevant in the
physics of aggregation where underlying processes may lead to a form of
correlated percolation. However, in ecology, the processes leading to habitat
correlations are in general not known or very complex. As proposed by Hiebeler
[Ecology {\bf 81}, 1629 (2000)], these correlations can be captured in a
lattice model by an observable aggregation parameter $q$, supplementing the
density $p$ of suitable sites. We investigate this model as an instance of
correlated percolation. We analyze the phase diagram of the percolation
transition and compute the cluster size distribution, the pair-connectedness
function $C(r)$ and the correlation function $g(r)$. We find that while $g(r)$
displays a power-law decrease associated with long-range correlations in a wide
domain of parameter values, critical properties are compatible with the
universality class of uncorrelated percolation. We contrast the correlation
structures obtained respectively for the correlated percolation model and for
the Ising model, and show that the diversity of habitat configurations
generated by the Hiebeler model is richer than the archetypal Ising model. We
also find that emergent structural properties are peculiar to the implemented
algorithm, leading to questioning the notion of a well-defined model of
aggregated habitat. We conclude that the choice of model and algorithm have
strong consequences on what insights ecological studies can get using such
models of species habitat.


  Economic incentives to harvest a species usually diminish as its abundance
declines, because harvest costs increase. This prevents harvesting to
extinction. A known exception can occur if consumer demand causes a declining
species' harvest price to rise faster than costs. This threat may affect rare
and valuable species, such as large land mammals, sturgeons, and bluefin tunas.
We analyze a similar but underappreciated threat, which arises when the
geographic area (range) occupied by a species contracts as its abundance
declines. Range contractions maintain the local densities of declining
populations, which facilitates harvesting to extinction by preventing abundance
declines from causing harvest costs to rise. Factors causing such range
contractions include schooling, herding, or flocking behaviors--which,
ironically, can be predator-avoidance adaptations; patchy environments; habitat
loss; and climate change. We use a simple model to identify combinations of
range contractions and price increases capable of causing extinction from
profitable overharvesting, and we compare these to an empirical review. We find
that some aquatic species that school or forage in patchy environments
experience sufficiently severe range contractions as they decline to allow
profitable harvesting to extinction even with little or no price increase; and
some high-value declining aquatic species experience severe price increases.
For terrestrial species, the data needed to evaluate our theory are scarce, but
available evidence suggests that extinction-enabling range contractions may be
common among declining mammals and birds. Thus, factors causing range
contraction as abundance declines may pose unexpectedly large extinction risks
to harvested species.


  Plant diseases often cause serious yield losses in agriculture. A pathogen's
reproductive fitness can be quantified by the basic reproductive number, R0.
Since pathogen transmission between host plants depends on the spatial
separation between them, R0 is strongly influenced by the spatial scales of
pathogen dispersal and the spatial scales of the host population. The basic
reproductive number was found to increase with the field size at small field
sizes and to saturate to a constant value at large field sizes. It reaches a
maximum in quadratic fields and decreases as the field becomes elongated. This
pattern appears to be quite general: it holds for dispersal kernels that
decrease exponentially or faster as well as for "fat-tailed" dispersal kernels
that decrease slower than exponential (i.e. power-law kernels). We used this
approach to estimate R0 in wheat stripe rust (an important pathogen caused by
Puccinia striiformis), since disease gradients for this pathogen were
thoroughly measured over large distances [Sackett and Mundt, Phytopathology,
95, 983 (2005)]. For the two largest datasets, we estimated R0 in the limit of
large fields to be of the order of 50. These estimates are consistent with
independent field observations [Cowger et al. (2005), Phytopathology, 95,
97282; Farber et al. (2013), Phytopathology, 103, 41]. We present a proof of
principle of a novel approach to estimate the basic reproductive number, R0, of
plant pathogens using wheat stripe rust as a case study. We found that the
spatial extent over which R0 changes strongly is quite fine-scaled (about 30 m
of the linear extension of the field). Our results indicate that in order to
optimize the spatial scale of deployment of fungicides or host resistances, the
adjustments should be made at a fine spatial scale.


  In this paper, we present a new approach for analyzing gene expression data
that builds on topological characteristics of time series. Our goal is to
identify cell cycle regulated genes in micro array dataset. We construct a
point cloud out of time series using delay coordinate embeddings. Persistent
homology is utilized to analyse the topology of the point cloud for detection
of periodicity. This novel technique is accurate and robust to noise, missing
data points and varying sampling intervals. Our experiments using Yeast
Saccharomyces cerevisiae dataset substantiate the capabilities of the proposed
method.


  The critical importance of human milk to infants and even human civilization
has been well established. Although the human milk microbiome has received
increasing attention with the expansion of research on the human microbiome,
our understanding of the milk microbiome has been limited to cataloguing OTUs
and computation of community diversity indexes. To the best of our knowledge,
there has been no report on the bacterial interactions within the human milk
microbiome. To bridge this gap, we reconstructed a milk bacterial community
network with the data from Hunt et al (2011), which is the largest 16S-rRNA
sequence data set of human milk microbiome available to date. Our analysis
revealed that the milk microbiome network consists of two disconnected
sub-networks. One sub-network is a fully connected complete graph consisting of
seven genera as nodes and all of its pair-wise interactions among the bacteria
are facilitative or cooperative. In contrast, the interactions in the other
sub-network of 8 nodes are mixed but dominantly cooperative. Somewhat
surprisingly, the only 'non-cooperative' nodes in the second sub-network are
mutually cooperative Staphylococcus and Corynebacterium, genera that include
some opportunistic pathogens. This potentially 'evil' alliance between
Staphylococcus and Corynebacterium could be inhibited by the remaining nodes
who cooperate with one another in the second sub-network. We postulate that the
'confrontation' between the 'evil' alliance and 'benign' alliance in human milk
microbiome should have important health implications to lactating women and
their infants and shifting the balance between the two alliances may be
responsible for dysbiosis of the milk microbiome that permits mastitis. A
related study focusing on ecological analysis was reported at
(http://www.eurekalert.org/pub_releases/2014-09/scp-ahb090214.php).


  Riboswitches, RNA elements found in the untranslated region, regulate gene
expression by binding to target metaboloites with exquisite specificity.
Binding of metabolites to the conserved aptamer domain allosterically alters
the conformation in the downstream expression platform. The fate of gene
expression is determined by the changes in the downstream RNA sequence. As the
metabolite-dependent cotranscriptional folding and unfolding dynamics of
riboswitches is the key determinant of gene expression, it is important to
investigate both the thermodynamics and kinetics of riboswitches both in the
presence and absence of metabolite. Single molecule force experiments that
decipher the free energy landscape of riboswitches from their mechanical
responses, theoretical and computational studies have recently shed light on
the distinct mechanism of folding dynamics in different classes of
riboswitches. Here we first discuss the dynamics of water around riboswitch,
highlighting that water dynamics can enhance the fluctuation of nucleic acid
structure. To go beyond native state fluctuations we used the Self-Organized
Polymer (SOP) model to predict the dynamics of add adenine riboswitch under
mechanical forces. In addition to quantitatively predicting the folding
landscape of add-riboswitch our simulations also explain the difference in the
dynamics between pbuE adenine- and add adenine-riboswitches. In order to probe
the function {\it in vivo} we use the folding landscape to propose a system
level kinetic network model to quantitatively predict how gene expression is
regulated for riboswitches that are under kinetic control.


  Ultrasensitive detection and real-time monitoring of biological processes can
benefit significantly from the improved brightness and photostability of the
popular organic dyes such as cyanines. Here, using a model cyanine dye, Cy3, we
demonstrate that brightness and photostability of the dye is significantly
altered when trapped in a molecular container, e.g. cucurbit[n]urils (CB[n])
and cyclodextrins (CD).Through computational modeling, we predicted that Cy3
forms a stable inclusion complex with three different hosts, CB[7], beta-CD,
and methyl-beta-CD, which was further confirmed by single-molecule diffusion
measurements using fluorescence correlation spectroscopy. The effect of
supramolecular encapsulation on Cy3 photophysical properties was found to be
highly host-specific. Up to three-fold increase in brightness of Cy3 was
observed when the dye was trapped in methyl-beta-CD, due to an increase in both
dye absorption and quantum yield. Steady-state and time-resolved spectroscopy
of the various complexes revealed that host polarizability and restricted
mobility of the dye in the host both contribute to the observed increase in
molecular brightness. Furthermore, entrapment of the dye molecule in CDs
resulted in a marked increase in dye photostability, whereas the dye degraded
more rapidly in CB[7]. These results suggest that the changes in photophysical
properties of the dye afforded by supramolecular encapsulation are highly
dependent on the host molecule. The reported improvement in brightness and
photostability together with the excellent biocompatibility of cyclodextrins
makes supramolecular encapsulation a viable strategy for routine dye
enhancement.


  The multifractal detrended fluctuation analysis of time series is able to
reveal the presence of long-range correlations and, at the same time, to
characterize the self-similarity of the series. The rich information derivable
from the characteristic exponents and the multifractal spectrum can be further
analyzed to discover important insights about the underlying dynamical process.
In this paper, we employ multifractal analysis techniques in the study of
protein contact networks. To this end, initially a network is mapped to three
different time series, each of which is generated by a stationary unbiased
random walk. To capture the peculiarities of the networks at different levels,
we accordingly consider three observables at each vertex: the degree, the
clustering coefficient, and the closeness centrality. To compare the results
with suitable references, we consider also instances of three well-known
network models and two typical time series with pure monofractal and
multifractal properties. The first result of notable interest is that time
series associated to proteins contact networks exhibit long-range correlations
(strong persistence), which are consistent with signals in-between the typical
monofractal and multifractal behavior. Successively, a suitable embedding of
the multifractal spectra allows to focus on ensemble properties, which in turn
gives us the possibility to make further observations regarding the considered
networks. In particular, we highlight the different role that small and large
fluctuations of the considered observables play in the characterization of the
network topology.


  Animal movements have been related to optimal foraging strategies where
self-similar trajectories are central. Most of the experimental studies done so
far have focused mainly on fitting statistical models to data in order to test
for movement patterns described by power-laws. Here we show by analyzing over
half a million movement displacements that isolated termite workers actually
exhibit a range of very interesting dynamical properties --including L\'evy
flights-- in their exploratory behaviour. Going beyond the current trend of
statistical model fitting alone, our study analyses anomalous diffusion and
structure functions to estimate values of the scaling exponents describing
displacement statistics. We evince the fractal nature of the movement patterns
and show how the scaling exponents describing termite space exploration
intriguingly comply with mathematical relations found in the physics of
transport phenomena. By doing this, we rescue a rich variety of physical and
biological phenomenology that can be potentially important and meaningful for
the study of complex animal behavior and, in particular, for the study of how
patterns of exploratory behaviour of individual social insects may impact not
only their feeding demands but also nestmate encounter patterns and, hence,
their dynamics at the social scale.


  This paper describes types of errors arising in a recently proposed method of
incidence estimation from prevalence data. The errors are illustrated by a
simulation study about a hypothetical irreversible disease. In addition, a way
of obtaining error bounds in practical applications of the method is proposed.


  Cancers follow a clonal Darwinian evolution, with fitter subclones replacing
more quiescent cells, ultimately giving rise to macroscopic disease.
High-throughput genomics provides the opportunity to investigate these
processes and determine specific genetic alterations driving disease
progression. Genomic sampling of a patient's cancer provides a molecular
history, represented by a phylogenetic tree. Cohorts of patients represent a
forest of related phylogenetic structures. To extract clinically relevant
information, one must represent and statistically compare these collections of
trees. We propose a framework based on an application of the work by Billera,
Holmes and Vogtmann on phylogenetic tree spaces to the case of unrooted trees
of intra-individual cancer tissue samples. We observe that these tree spaces
are globally nonpositively curved, allowing for statistical inference on
populations of patient histories. A projective tree space is introduced,
permitting visualizations of aggregate evolutionary behavior. Published data
from three types of human malignancies are explored within our framework.


  The design of biologically-inspired wireless communication systems using
bacteria as the basic element of the system is initially motivated by a
phenomenon called \emph{Quorum Sensing}. Due to high randomness in the
individual behavior of a bacterium, reliable communication between two bacteria
is almost impossible. Therefore, we have recently proposed that a population of
bacteria in a cluster is considered as a bio node in the network capable of
molecular transmission and reception. This proposition enables us to form a
reliable bio node out of many unreliable bacteria.
  In this paper, we study the communication between two nodes in such a network
where information is encoded in the concentration of molecules by the
transmitter. The molecules produced by the bacteria in the transmitter node
propagate through the diffusion channel. Then, the concentration of molecules
is sensed by the bacteria population in the receiver node which would decode
the information and output light or fluorescent as a result. The uncertainty in
the communication is caused by all three components of communication, i.e.,
transmission, propagation and reception. We study the theoretical limits of the
information transfer rate in the presence of such uncertainties. Finally, we
consider M-ary signaling schemes and study their achievable rates and
corresponding error probabilities.


  Molecular communication between biological entities is a new paradigm in
communications. Recently, we studied molecular communication between two nodes
formed from synthetic bacteria. Due to high randomness in behavior of bacteria,
we used a population of them in each node. The reliability of such
communication systems depends on both the maximum concentration of molecules
that a transmitter node is able to produce at the receiver node as well as the
number of bacteria in each nodes. This maximum concentration of molecules falls
with distance which makes the communication to the far nodes nearly impossible.
In order to alleviate this problem, in this paper, we propose to use a
molecular relaying node. The relay node can resend the message either by the
different or the same type of molecules as the original signal from the
transmitter. We study two scenarios of relaying. In the first scenario, the
relay node simply senses the received concentration and forwards it to the
receiver. We show that this sense and forward scenario, depending on the type
of molecules used for relaying, results in either increasing the range of
concentration of molecules at the receiver or increasing the effective number
of bacteria in the receiver node. For both cases of sense and forward relaying,
we obtain the resulting improvement in channel capacity. We conclude that
multi-type molecular relaying outperforms the single-type relaying. In the
second scenario, we study the decode and forward relaying for the M-ary
signaling scheme. We show that this relaying strategy increases the reliability
of M-ary communication significantly.


  We study the dynamics of DNA hairpin formation using oxDNA, a
nucleotide-level coarse-grained model of DNA. In particular, we explore the
effects of the loop stacking interactions and non-native base pairing on the
hairpin closing times. We find a non-monotonic variation of the hairpin closing
time with temperature, in agreement with the experimental work of Wallace et
al. [Proc. Nat. Acad. Sci. USA 2001, 98, 5584-5589]. The hairpin closing
process involves the formation of an initial nucleus of one or two bonds
between the stems followed by a rapid zippering of the stem. At high
temperatures, typically above the hairpin melting temperature, an effective
negative activation enthalpy is observed because the nucleus has a lower
enthalpy than the open state. By contrast, at low temperatures, the activation
enthalpy becomes positive mainly due to the increasing energetic cost of
bending a loop that becomes increasingly highly stacked as the temperature
decreases. We show that stacking must be very strong to induce this
experimentally observed behavior, and that the existence of just a few weak
stacking points along the loop can substantially suppress it. Non-native base
pairs are observed to have only a small effect, slightly accelerating hairpin
formation.


  When estimating a phylogeny from a multiple sequence alignment, researchers
often assume the absence of recombination. However, if recombination is
present, then tree estimation and all downstream analyses will be impacted,
because different segments of the sequence alignment support different
phylogenies. Similarly, convergent selective pressures at the molecular level
can also lead to phylogenetic tree incongruence across the sequence alignment.
Current methods for detection of phylogenetic incongruence are not equipped to
distinguish between these two different mechanisms and assume that the
incongruence is a result of recombination or other horizontal transfer of
genetic information. We propose a new recombination detection method that can
make this distinction, based on synonymous codon substitution distances.
Although some power is lost by discarding the information contained in the
nonsynonymous substitutions, our new method has lower false positive
probabilities than the comparable recombination detection method when the
phylogenetic incongruence signal is due to convergent evolution. We apply our
method to three empirical examples, where we analyze: 1) sequences from a
transmission network of the human immunodeficiency virus, 2) tlpB gene
sequences from a geographically diverse set of 38 Helicobacter pylori strains,
and 3) Hepatitis C virus sequences sampled longitudinally from one patient.


  The intricate network of interactions observed in RNA three-dimensional
structures is often described in terms of a multitude of geometrical
properties, including helical parameters, base pairing/stacking, hydrogen
bonding and backbone conformation. We show that a simple molecular
representation consisting in one oriented bead per nucleotide can account for
the fundamental structural properties of RNA. In this framework, canonical
Watson-Crick, non-Watson-Crick base-pairing and base-stacking interactions can
be unambiguously identified within a well-defined interaction shell. We
validate this representation by performing two independent, complementary
tests. First, we use it to construct a sequence-independent, knowledge-based
scoring function for RNA structural prediction, which compares favorably to
fully atomistic, state-of-the-art techniques. Second, we define a metric to
measure deviation between RNA structures that directly reports on the
differences in the base-base interaction network. The effectiveness of this
metric is tested with respect to the ability to discriminate between
structurally and kinetically distant RNA conformations, performing better
compared to standard techniques. Taken together, our results suggest that this
minimalist, nucleobase-centric representation captures the main interactions
that are relevant for describing RNA structure and dynamics.


  The Kutch desert (Great Rann of Kutch, Gujarat, India) is a unique ecosystem:
in the larger part of the year it is a hot, salty desert that is flooded
regularly in the Indian monsoon season. In the dry season, the crystallized
salt deposits form the "white desert" in large regions. The first metagenomic
analysis of the soil samples of Kutch was published in 2013, and the data was
deposited in the NCBI Sequence Read Archive. The sequences were analyzed at the
same time phylogenetically for prokaryotes, especially for bacterial taxa.
  In the present work, we are searching for the DNA sequences of the recently
discovered giant viruses in the soil samples of the Kutch desert. Since most
giant viruses were discovered in biofilms in industrial cooling towers, ocean
water and freshwater ponds, we were surprised to find their DNA sequences in
the soil samples of a seasonally very hot and arid, salty environment.


  Situations where people have to decide between hurting themselves or another
person are at the core of many individual and global conflicts. Yet little is
known about how people behave when facing these situations in the lab. Here we
report a large experiment in which participants could either take $x$ dollars
from another anonymous participant or give $y$ dollars to the same participant.
Depending on the treatments, participants could also exit the game without
making any decision, but paying a cost. Across different protocols and
parameter specifications, we provide evidence of three regularities: (i) when
exiting is allowed and costless, subjects tend to exit the game; (ii) females
are more likely than males to exit the game, but only when the cost is small;
(iii) when exiting is not allowed, altruistic actions are more common than
predicted by the dominant economic models. In particular, against the
predictions of every dominant economic model, about one sixth of the subjects
show hyper-altruistic tendencies, that is, they prefer giving $y$ rather than
taking $x>y$. In doing so, our findings shed light on human decision-making in
conflictual situations and suggest that economic models should be revised in
order to take into account hyper-altruistic behaviour.


  In this work we performed a numerical study of an epidemic model that mimics
the endemic state of whooping cough in the pre-vaccine era. We considered a
stochastic SIR model on dynamical networks that involve local and global
contacts among individuals and analyzed the influence of the network properties
on the characterization of the quasi-stationary state. We computed probability
density functions (PDF) for infected fraction of individuals and found that
they are well fitted by gamma functions, excepted the tails of the
distributions that are q-exponentials. We also computed the fluctuation power
spectra of infective time series for different networks. We found that network
effects can be partially absorbed by rescaling the rate of infective contacts
of the model. An explicit relation between the effective transmission rate of
the disease and the correlation of susceptible individuals with their infective
nearest neighbours was obtained. This relation quantifies the known screening
of infective individuals observed in these networks. We finally discuss the
goodness and limitations of the SIR model with homogeneous mixing and
parameters taken from epidemiological data to describe the dynamic behaviour
observed in the networks studied.


  Leading Ebola subtypes exhibit a wide mortality range, here explained at the
molecular level by using fractal hydropathic scaling of amino acid sequences
based on protein self-organized criticality. Specific hydrophobic features in
the hydrophilic mucin-like domain suffice to account for the wide mortality
range. Significance statement: Ebola virus is spreading rapidly in Africa. The
connection between protein amino acid sequence and mortality is identified
here.