  Protected areas are an important conservation measure. However, there are
controversial findings regarding whether closed areas are beneficial for
species and habitat conservation as well as landings. Species dispersal is
acknowledged as a key factor for the design and impacts of closed areas. A
series of agent based models using random diffusion to model fish dispersal
were run before and after habitat protection. All results were normalised
without the protected habitat in each scenario to detect the relative
difference after closing an area, all else being equal. Results show that
landings of species with short dispersal ranges will take longer to reach the
levels of pre Marine Protected Areas (MPAs) establishment than landings of
species with long dispersal ranges. Further the establishment of an MPA
generates a higher relative population source within the MPA for species with
low dispersal abilities than for species with high dispersal abilities. Results
derived here show that there exists a win-win feasible scenario that maximises
both fish biomass as well as fish catches.


  Motivation: To assess the quality of a protein model, i.e. to estimate how
close it is to its native structure, using no other information than the
structure of the model has been shown to be useful for structure prediction.
The state of the art method, ProQ2, is based on a machine learning approach
that uses a number of features calculated from a protein model. Here, we
examine if these features can be exchanged with energy terms calculated from
Rosetta and if a combination of these terms can improve the quality assessment.
Results: When using the full atom energy function from Rosetta in ProQRosFA the
QA is on par with our previous state-of-the-art method, ProQ2. The method based
on the low-resolution centroid scoring function, ProQRosCen, performs almost as
well and the combination of all the three methods, ProQ2, ProQRosFA and
ProQCenFA into ProQ3 show superior performance over ProQ2. Availability: ProQ3
is freely available on BitBucket at https://bitbucket.org/ElofssonLab/proq3


  Motivation: De Bruijn graphs have been proposed as a data structure to
facilitate the analysis of related whole genome sequences, in both a population
and comparative genomic settings. However, current approaches do not scale well
to many genomes of large size (such as mammalian genomes). Results: In this
paper, we present TwoPaCo, a simple and scalable low memory algorithm for the
direct construction of the compacted de Bruijn graph from a set of complete
genomes. We demonstrate that it can construct the graph for 100 simulated human
genomes in less then a day and eight real primates in less than two hours, on a
typical shared-memory machine. We believe that this progress will enable novel
biological analyses of hundreds of mammalian-sized genomes. Availability: Our
code and data is available for download from github.com/medvedevgroup/TwoPaCo
Contact: ium125@psu.edu


  Many ion channels spontaneously switch between different levels of activity.
Although this behaviour known as modal gating has been observed for a long time
it is currently not well understood. Despite the fact that appropriately
representing activity changes is essential for accurately capturing time course
data from ion channels, systematic approaches for modelling modal gating are
currently not available. In this paper, we develop a modular approach for
building such a model in an iterative process. First, stochastic switching
between modes and stochastic opening and closing within modes are represented
in separate aggregated Markov models. Second, the continuous-time hierarchical
Markov model, a new modelling framework proposed here, then enables us to
combine these components so that in the integrated model both mode switching as
well as the kinetics within modes are appropriately represented. A mathematical
analysis reveals that the behaviour of the hierarchical Markov model naturally
depends on the properties of its components. We also demonstrate how a
hierarchical Markov model can be parameterised using experimental data and show
that it provides a better representation than a previous model of the same data
set. Because evidence is increasing that modal gating reflects underlying
molecular properties of the channel protein, it is likely that biophysical
processes are better captured by our new approach than in earlier models.


  We investigate the interplay between post-translational folding and escape of
two small single-domain proteins at the ribosomal exit tunnel by using Langevin
dynamics with coarse-grained models. It is shown that at temperatures lower or
near the temperature of the fastest folding, folding proceeds concomitantly
with the escape process, resulting in vectorial folding and enhancement of
foldability of nascent proteins. The concomitance between the two processes,
however, deteriorates as temperature increases. Our folding simulations as well
as free energy calculation by using umbrella sampling show that, at low
temperatures, folding at the tunnel follows one or two specific pathways
without kinetic traps. It is shown that the escape time can be mapped to a
one-dimensional diffusion model with two different regimes for temperatures
above and below the folding transition temperature. Attractive interactions
between amino acids and attractive sites on the tunnel wall lead to a free
energy barrier along the escape route of protein. It is suggested that this
barrier slows down the escape process and consequently promotes correct folding
of the released nascent protein.


  Despite the common misconception of nearly static organisms, plants do
interact continuously with the environment and with each other. It is fair to
assume that during their evolution they developed particular features to
overcome problems and to exploit possibilities from environment. In this paper
we introduce various quantitative measures based on recent advancements in
complex network theory that allow to measure the effective similarities of
various species. By using this approach on the similarity in fruit-typology
ecological traits we obtain a clear plant classification in a way similar to
traditional taxonomic classification. This result is not trivial, since a
similar analysis done on the basis of diaspore morphological properties do not
provide any clear parameter to classify plants species. Complex network theory
can then be used in order to determine which feature amongst many can be used
to distinguish scope and possibly evolution of plants. Future uses of this
approach range from functional classification to quantitative determination of
plant communities in nature.


  We study the conformational freedom of a protein made by two rigid domains
connected by a flexible linker. The conformational freedom is represented as an
unknown probability distribution on the space of allowed states. A new
algorithm for the calculation of the Maximum Allowable Probability is proposed,
which can be extended to any type of measurements. In this paper we use Pseudo
Contact Shifts and Residual Dipolar Coupling. We reconstruct a single central
tendency in the distribution and discuss in depth the results.


  3 transmembrane and 4 transmembrane helices models are suggested for the
human vitamin K epoxide reductase (VKOR). In this study, we investigate the
stability of the human 3 transmembrane/4 transmembrane VKOR models employing a
coarse-grained normal mode analysis and molecular dynamics simulation. Based on
the analysis of the mobility of each transmembrane domain, we suggest that the
3 transmembrane human VKOR model is more stable than the 4 transmembrane human
VKOR model.


  The mechanism of threshold elongation (overstretching) of DNA macromolecules
under the action of external force is studied within the framework of
phenomenological approach. When considering the task it is taken into account
that double-stranded DNA is a polymorphic macromolecule with a set of
metastable states. Accordingly to proposed mechanism, DNA threshold elongation
can take place as a cooperative structural transition of macromolecule to
metastable form, which is stabilized by the external force. For the description
of DNA overstretching, the model included external (stretching) and internal
(conformation) displacement components is constructed. As assumed, both
components are coupled on the pathway of double helix transformations. It is
shown, that under force action DNA deformation proceeds in two stages. First
one is the restructuring of the double helix, which leads to formation of
conformational bistability in DNA macromolecule. On the second stage, the
conformational transition and the deformation induced by it consistently cover
the macromolecule as a threshold process. The comparison of calculated
characteristics of DNA overstretching process with experimental data shows good
agreement. Analysis of the results obtained, together with the available
literature data allow to conclude, that overstretching transition in the double
helix is dynamical process and can spread in DNA chain. In the same time, in
DNA with A$\cdot$T-rich content due to large dissipation the overstretching
process leads to force induced melting transition, and should have nearly
static character.


  Isoform quantification is an important goal of RNA-seq experiments, yet it
remains prob- lematic for genes with low expression or several isoforms. These
difficulties may in principle be ameliorated by exploiting correlated
experimental designs, such as time series or dosage response experiments. Time
series RNA-seq experiments, in particular, are becoming in- creasingly popular,
yet there are no methods that explicitly leverage the experimental design to
improve isoform quantification. Here we present DICEseq, the first isoform
quantification method tailored to correlated RNA-seq experiments. DICEseq
explicitly models the corre- lations between different RNA-seq experiments to
aid the quantification of isoforms across experiments. Numerical experiments on
simulated data sets show that DICEseq yields more accurate results than
state-of-the-art methods, an advantage that can become considerable at low
coverage levels. On real data sets, our results show that DICEseq provides
substan- tially more reproducible and robust quantifications, increasing the
correlation of estimates from replicate data sets by up to 10% on genes with
low or moderate expression levels (bot- tom third of all genes). Furthermore,
DICEseq permits to quantify the trade-off between temporal sampling of RNA and
depth of sequencing, frequently an important choice when planning experiments.
Our results have strong implications for the design of RNA-seq ex- periments,
and offer a novel tool for improved analysis of such data sets. Python code is
freely available at http://diceseq.sf.net.


  Epidemic control is of great importance for human society. Adjusting
interacting partners is an effective individualized control strategy.
Intuitively, it is done either by shortening the interaction time between
susceptible and infected individuals or by increasing the opportunities for
contact between susceptible individuals. Here, we provide a comparative study
on these two control strategies by establishing an epidemic model with
non-uniform stochastic interactions. It seems that the two strategies should be
similar, since shortening the interaction time between susceptible and infected
individuals somehow increases the chances for contact between susceptible
individuals. However, analytical results indicate that the effectiveness of the
former strategy sensitively depends on the infectious intensity and the
combinations of different interaction rates, whereas the latter one is quite
robust and efficient. Simulations are shown in comparison with our analytical
predictions. Our work may shed light on the strategic choice of disease
control.


  Overtreatment of early-stage low-risk prostate cancer (PC) patients
represents a significant problem in disease management and has socio-economic
implications. Development of genetic and molecular markers of clinically
significant disease in patients diagnosed with low grade localized PC would
have a major impact in disease management. A gene expression signature (GES) is
reported for lethal PC in biopsy specimens obtained at the time of diagnosis
from patients with Gleason 6 and Gleason 7 tumors in a Swedish watchful waiting
cohort with up to 30 years follow-up. A 98-genes GES identified 89 and 100
percent of all death events 4 years after diagnosis in G7 and G6 patients,
respectively; at 6 years follow-up, 83 and 100 percent of all deaths events
were captured. Remarkably, the 98-genes GES appears to perform successfully in
patients stratification with as little as 2% of cancer cells in a specimen,
strongly indicating that it captures a malignant field effect in prostates
harboring cancer cells of different degrees of aggressiveness. In G6 and G7
tumors from PC patients of age 65 or younger, GES identified 86 percent of all
death events during the entire follow-up period. In G6 and G7 tumors from PC
patients of age 70 or younger, GES identified 90 percent of all death events 6
years after diagnosis. Classification performance of the reported in this study
98-genes GES of lethal PC appeared suitable to meet design and feasibility
requirements of a prospective 4 to 6 years clinical trial, which is essential
for regulatory approval of diagnostic and prognostic tests in clinical setting.
Prospectively validated GES of lethal PC in biopsy specimens of G6 and G7
tumors will help physicians to identify, at the time of diagnosis, patients who
should be considered for exclusion from active surveillance programs and who
would most likely benefit from immediate curative interventions.


  Main goal: The aim of this note is to propose a modeling approach for
assessing the effectiveness of pest biocontrol by natural enemies in
diversified agricultural landscapes including several pesticide-based
management strategies. Our approach combines a stochastic landscape model with
a spatially-explicit model of population dynamics. It enables us to analyze the
effect of the landscape composition (proportion of semi-natural habitat,
non-treated crops, slightly treated crops and conventionally treated crops) on
the effectiveness of pest biocontrol. Effectiveness is measured through
environmental and agronomical descriptors, measuring respectively the impact of
the pesticides on the environment and the average agronomic productivity of the
whole landscape taking into account losses caused by pests.
  Conclusions: The effectiveness of the pesticide, the intensity of the
treatment and the pest intrinsic growth rate are found to be the main drivers
of landscape productivity. The loss in productivity due to a reduced use of
pesticide can be partly compensated by biocontrol. However, the model suggests
that it is not possible to maintain a constant level of productivity while
reducing the use of pesticides, even with highly efficient natural enemies.
Fragmentation of the semi-natural habitats and increased crop rotation tend to
slightly enhance the effectiveness of biocontrol but have a marginal effect
compared to the predation rate by natural enemies.
  This note was written in the framework of the ANR project PEERLESS
"Predictive Ecological Engineering for Landscape Ecosystem Services and
Sustainability"(ANR-12-AGRO-0006).


  Transcription factors (TFs), the key players in transcriptional regulation,
have attracted great experimental attention, yet the functions of most human
TFs remain poorly understood. Recent capabilities in genome-wide protein
binding profiling have stimulated systematic studies of the hierarchical
organization of human gene regulatory network and DNA-binding specificity of
TFs, shedding light on combinatorial gene regulation. We show here that these
data also enable a systematic annotation of the biological functions and
functional diversity of TFs. We compiled a human gene regulatory network for
384 TFs covering the 146,096 TF-target gene relationships, extracted from over
850 ChIP-seq experiments as well as the literature. By integrating this network
of TF-TF and TF-target gene relationships with 3,715 functional concepts from
six sources of gene function annotations, we obtained over 9,000 confident
functional annotations for 279 TFs. We observe extensive connectivity between
transcription factors and Mendelian diseases, GWAS phenotypes, and
pharmacogenetic pathways. Further, we show that transcription factors link
apparently unrelated functions, even when the two functions do not share common
genes. Finally, we analyze the pleiotropic functions of TFs and suggest that
increased number of upstream regulators contributes to the functional
pleiotropy of TFs. Our computational approach is complementary to focused
experimental studies on TF functions, and the resulting knowledge can guide
experimental design for discovering the unknown roles of TFs in human disease
and drug response.


  The problem of image restoration in cryo-EM entails correcting for the
effects of the Contrast Transfer Function (CTF) and noise. Popular methods for
image restoration include `phase flipping', which corrects only for the Fourier
phases but not amplitudes, and Wiener filtering, which requires the spectral
signal to noise ratio. We propose a new image restoration method which we call
`Covariance Wiener Filtering' (CWF). In CWF, the covariance matrix of the
projection images is used within the classical Wiener filtering framework for
solving the image restoration deconvolution problem. Our estimation procedure
for the covariance matrix is new and successfully corrects for the CTF. We
demonstrate the efficacy of CWF by applying it to restore both simulated and
experimental cryo-EM images. Results with experimental datasets demonstrate
that CWF provides a good way to evaluate the particle images and to see what
the dataset contains even without 2D classification and averaging.


  Collective sensing is an emergent phenomenon which enables individuals to
estimate a hidden property of the environment through the observation of social
interactions. Previous work on collective sensing shows that gregarious
individuals obtain an evolutionary advantage by exploiting collective sensing
when competing against solitary individuals. This work addresses the question
of whether collective sensing allows for the emergence of groups from a
population of individuals without predetermined behaviors. It is assumed that
group membership does not lessen competition on the limited resources in the
environment, e.g. groups do not improve foraging efficiency. Experiments are
run in an agent-based evolutionary model of a foraging task, where the fitness
of the agents depends on their foraging strategy. The foraging strategy of
agents is determined by a neural network, which does not require explicit
modeling of the environment and of the interactions between agents. Experiments
demonstrate that gregarious behavior is not the evolutionary-fittest strategy
if resources are abundant, thus invalidating previous findings in a specific
region of the parameter space. In other words, resource scarcity makes
gregarious behavior so valuable as to make up for the increased competition
over the few available resources. Furthermore, it is shown that a population of
solitary agents can evolve gregarious behavior in response to a sudden scarcity
of resources, thus individuating a possible mechanism that leads to gregarious
behavior in nature. The evolutionary process operates on the whole parameter
space of the neural networks, hence these behaviors are selected among an
unconstrained set of behavioral models.


  Network ecology is a rising field of quantitative biology representing
ecosystems as complex networks. A suitable example is parasite spreading:
several parasites may be transmitted among their hosts through different
mechanisms, each one giving rise to a network of interactions. Modelling these
networked, ecological interactions at the same time is still an open challenge.
We present a novel spatially-embedded multiplex network framework for modelling
multi-host infection spreading through multiple routes of transmission. Our
model is inspired by T. cruzi, a parasite transmitted by trophic and vectorial
mechanisms. Our ecological network model is represented by a multiplex in which
nodes represent species populations interacting through a food web and a
parasite contaminative layer at the same time. We modelled an SI dynamics in
two different scenarios: a simple theoretical food web and an empirical one.
Our simulations in both scenarios show that the infection is more widespread
when both the trophic and the contaminative interactions are considered with
equal rates. This indicates that trophic and contaminative transmission may
have additive effects in real ecosystems. We also find that the ratio of
vectors-to-host in the community (i) crucially influences the infection spread,
(ii) regulates a percolating phase transition in the rate of parasite
transmission and (iii) increases the infection rate in hosts. By immunising the
same fractions of predator and prey populations, we show that the multiplex
topology is fundamental in outlining the role that each host species plays in
parasite transmission in a given ecosystem. We also show that the multiplex
models provide a richer phenomenology in terms of parasite spreading dynamics
compared to mono-layer models. Our work opens new challenges and provides new
quantitative tools for modelling multi-channel spreading in networked systems.


  The widely used model for evolutionary relationships is a bifurcating tree
with all taxa/observations placed at the leaves. This is not appropriate if the
taxa have been densely sampled across evolutionary time and may be in a direct
ancestral relationship, or if there is not enough information to fully resolve
all the branching points in the evolutionary tree. In this paper, we present a
fast distance-based agglomeration method called family-joining (FJ) for
constructing so-called generally labeled trees in which taxa may be placed at
internal vertices and the tree may contain polytomies. FJ constructs such trees
on the basis of pairwise distances and a distance threshold. We tested three
methods for threshold selection, FJ-AIC, FJ-BIC and FJ-CV, which minimize
Akaike information criterion, Bayesian information criterion, and
cross-validation error, respectively. When compared with related methods on
simulated data, FJ-BIC was among the best at reconstructing the correct tree
across a wide range of simulation scenarios. FJ-BIC was applied to HIV
sequences sampled from individuals involved in a known transmission chain. The
FJ-BIC tree was found to be compatible with almost all transmission events. On
average, internal branches in the FJ-BIC tree have higher bootstrap support
than branches in the leaf-labeled bifurcating tree constructed using RAxML.
$36\%$ and $25\%$ of the internal branches in the FJ-BIC tree and RAxML tree,
respectively, have bootstrap support greater than $70\%$. To the best of our
knowledge the method presented here is the first attempt at modeling
evolutionary relationships using generally labeled trees.


  The aim of this study is to compare the growth speed of different cell
populations measured by their Malthus parameter. We focus on both the
age-structured and size-structured equations. A first population (of reference)
is composed of cells all aging or growing at the same rate $\bar v$. A second
population (with variability) is composed of cells each aging or growing at a
rate $v$ drawn according to a non-degenerated distribution $\rho$ with mean
$\bar v$. In a first part, analytical answers -- based on the study of an
eigenproblem -- are provided for the age-structured model. In a second part,
numerical answers -- based on stochastic simulations -- are derived for the
size-structured model. It appears numerically that the population with
variability proliferates more slowly than the population of reference (for
experimentally plausible division rates). The decrease in the Malthus parameter
we measure, around 2% for distributions $\rho$ with realistic coefficients of
variations around 15-20\%, is determinant since it controls the {\it
exponential} growth of the whole population.


  Structured population models are a class of general evolution equations which
are widely used in the study of biological systems. Many theoretical methods
are available for establishing existence and stability of steady states of
general evolution equations. However, except for very special cases, finding an
analytical form of stationary solutions for evolution equations is a
challenging task. In the present paper, we develop a numerical framework for
computing approximations to stationary solutions of general evolution
equations, which can also be used to produce existence and stability regions
for steady states. In particular, we use the Trotter-Kato Theorem to
approximate the infinitesimal generator of an evolution equation on a finite
dimensional space, which in turn reduces the evolution equation into a system
of ordinary differential equations. Consequently, we approximate and study the
asymptotic behavior of stationary solutions. We illustrate the convergence of
our numerical framework by applying it to a linear Sinko-Streifer structured
population model for which the exact form of the steady state is known. To
further illustrate the utility of our approach, we apply our framework to
nonlinear population balance equation, which is an extension of well-known
Smoluchowksi coagulation-fragmentation model to biological populations. We also
demonstrate that our numerical framework can be used to gain insight about the
theoretical stability of the stationary solutions of the evolution equations.
Furthermore, the open source Python program that we have developed for our
numerical simulations is freely available from our Github repository
(github.com/MathBioCU).


  Cellular decision making is based on regulatory circuits that associate
signal thresholds to specific physiological actions. This transmission of
information is subjected to molecular noise what can decrease its fidelity.
Here, we show instead how such intrinsic noise enhances information transfer in
the presence of multiple circuit copies. The result is due to the contribution
of noise to the generation of autonomous responses by each copy, which are
altogether associated with a common decision. Moreover, factors that correlate
the responses of the redundant units (extrinsic noise or regulatory cross-talk)
contribute to reduce fidelity, while those that further uncouple them
(heterogeneity within the copies) can lead to stronger information gain.
Overall, our study emphasizes how the interplay of signal thresholding,
redundancy, and noise influences the accuracy of cellular decision making.
Understanding this interplay provides a basis to explain collective cell
signaling mechanisms, and to engineer robust decisions with noisy genetic
circuits.


  Contemporary techniques in biology produce readouts for large numbers of
genes simultaneously, the typical example being differential gene expression
measurements. Moreover, those genes are often richly annotated using GO terms
that describe gene function and that can be used to summarize the results of
the genome-scale experiments. However, making sense of such GO enrichment
analyses may be challenging. For instance, overrepresented GO functions in a
set of differentially expressed genes are typically output as a flat list, a
format not adequate to capture the complexities of the hierarchical structure
of the GO annotation labels.
  In this chapter, we survey the various methods to visualize large,
difficult-to-interpret lists of GO terms. We catalogue their availability
(web-based or standalone), the main principles they employ in summarizing large
lists of GO terms, and the visualization styles they support. These brief
commentaries on each software are intended as a helpful inventory, rather than
comprehensive descriptions of the underlying algorithms. Instead, we show
examples of their use and suggest that the choice of an appropriate
visualization tool may be crucial to the utility of GO in biological discovery.


  The paper presents the importance of research which characterizes the natural
forest structure for the forest management. The lessons learned in these
particular forest ecosystems can be integrated by the forest management
objectives, in order to increase the sustainability of this type of resources.
The project NATFORMAN was focused on the structure of the natural forest, thus
research methodologies and modern technology (such as Field-Map) investigation
and determination were used in order to record information on forest structural
parameters. The results obtained refer to these structural parameters and to
the possibility of transferring such information in practice, in order to
achieve forest sustainable management.


  More than 300,000 new cases worldwide are being diagnosed with oral cancer
annually. Complexity of oral cancer renders designing drug targets very
difficult. We analyse protein-protein interaction network for the normal and
oral cancer tissue and detect crucial changes in the structural properties of
the networks in terms of the interactions of the hub proteins and the
degree-degree correlations. Further analysis of the spectra of both the
networks, while exhibiting universal statistical behavior, manifest distinction
in terms of the zero degeneracy, providing insight to the complexity of the
underlying system.


  A re-refinement of 4xan, hen egg white lysozyme (HEWL) with carboplatin
crystallised in NaBr solution, has been made (Tanley et al 2016). This follows
our Response article (Tanley et al 2015) to the Critique article of Shabalin et
al 2015, suggesting the need for corrections to some solute molecule
interpretations of electron density in 4xan and removal of an organic moiety as
a ligand to the platinum ion coordinated to His15. We note Shabalin et al
(2015) model of a chlorine in that density and a close by bromine at partial
occupancy to explain the shape. However, as the bromide concentration is in
huge excess over chloride (by 20 fold), we think that the 4yem Shabalin et al
2015 interpretation highly unlikely, but nevertheless we still cannot offer an
explanation for that shape, confirming our earlier analysis described in Tanley
et al (2014). Following Shabalin et al (2015) reprocessing of the raw
diffraction data for 4g4a, we also redid the diffraction data processing for
4xan to a higher resolution using EVAL (Schreurs et al 2010) concluding in
favour of 1.3 Angstrom as the resolution limit and which is the basis for our
revised PDB file for 4xan (5HMJ). It is very interesting that there is extra
X-ray diffraction data from 1.47 to 1.30 Angstrom resolution e.g. with
<I/sigma(I)> =0.39 and CC1/2 = 0.181 in the final shell (1.30 to 1.322
Angstrom). In this arXiv article we document in detail our different solvent
and split occupancy side chain electron density interpretations as evidence for
our statement of approach in our Response article (Tanley et al 2015). Our
critical re-examination includes comparisons based on the 4xan diffraction data
images reprocessing with three different software packages so as to evaluate
the possibility of variations in electron density interpretations due to that.
Overall our finalised model (PDB code 5HMJ) is now improved over 4xan.


  Changes in the frequency, duration, and severity of climate extremes are
forecast to occur under global climate change. The impacts of climate extremes
on forest productivity and health are complicated by potential interactions
with disturbance events and stand dynamics. The effects of stand dynamics on
forest responses to climate and disturbance are particularly important given
forest characteristics driven by stand dynamics can be modified through forest
management with the goal of increasing forest resistance and resilience to
climate change. We develop a hierarchical Bayesian state-space model allowing
climate effects on tree growth to vary over time and in relation to climate
extremes, disturbance events, and stand dynamics. We apply the model to a
dendrochronology dataset comprising measurements from forest stands of varying
composition, structure, and development stage in northeastern Minnesota.
Results indicate average forest growth was most sensitive to variables
describing climatic water deficit. Forest growth responses to water deficit
were partitioned into responses driven by climatic threshold exceedances and
interactions with forest tent caterpillar defoliation. Forest growth was both
resistant and resilient to climate extremes with the majority of forest growth
responses occurring after multiple climatic threshold exceedances or insect
defoliation events. Forest growth was most sensitive to water deficit during
periods of high stem density following major regeneration events when average
inter-tree competition was high. Results suggest that forest growth resistance
and resilience to interactions between climate extremes and insect defoliation
can be increased through management steps such as thinning to reduce
competition during early stages of stand development and small-group selection
harvests to maintain forest structures characteristic of older, mature stands.


  Adaptation of bacteria occurs predominantly via horizontal gene transfer
(HGT). While it is widely recognized that horizontal acquisitions frequently
encompass multiple genes, it is unclear what the size distribution of
successfully transferred DNA segments looks like and what evolutionary forces
shape this distribution. Here, we identified 1790 gene family pairs that were
consistently co-gained on the same branches across a phylogeny of 53 E. coli
strains. We estimated a lower limit of their genomic distances at the time they
were transferred to their host genomes; this distribution shows a sharp upper
bound at 30 kb. The same gene-pairs can have larger distances (up to 70 kb) in
other genomes. These more distant pairs likely represent recent acquisitions
via transduction that involve the co-transfer of excised prophage genes, as
they are almost always associated with intervening phage-associated genes. The
observed distribution of genomic distances of co-transferred genes is much
broader than expected from a model based on the co-transfer of genes within
operons; instead, this distribution is highly consistent with the size
distribution of supra-operonic clusters (SOCs), groups of co-occurring and
co-functioning genes that extend beyond operons. Thus, we propose that SOCs
form a basic unit of horizontal gene transfer.


  The hominid-specific non-LTR retrotransposon termed SINE VNTR Alu (SVA) is
the youngest of the transposable elements in the human genome. The propagation
of the most ancient SVA type A took place about thirteen millions years ago
ago, and the youngest SVA types appeared in the human genome after the
chimpanzee divergence. Functional enrichment analysis of genes associated with
SVA insertions demonstrated their strong link to multiple ontological
categories attributed to brain function and the disorders. SVA types that
expanded their presence in the human genome at different stages of hominoid
life history were also associated with progressively evolving behavioural
features that indicated a potential impact of SVA propagation on a cognitive
ability of a modern human. The SVA-associated genes were highly cross-linked in
functional networks suggesting an accumulative impact of functional alterations
potentially caused by SVA insertions. Our analysis uncovered a potential role
of SVAs in evolution of human CNS and especially emergence of functional trends
relevant to social and parental behaviour. It also supports models which
explain in part how brain function can be modulated by both the immune and
reproductive systems based on the gene expression patterns and gene pathways
potentially altered by SVA insertions.


  Planar curves with periodically varying curvature arise in the natural
sciences as the result of a wide variety of periodic processes. The total
curvature of a periodic arc in such curves constrains their symmetry. It is
shown how the total curvature can be computed without reparameterizing the
curve to unit speed. The use of the total curvature of the periodic arcs is
demonstrated through a series of four examples from various branches of
science. Insights gained from these examples are applied to improve the
modeling of spiral wave meander.


  We present a method of dimensional reduction for the general Markov model of
sequence evolution on a phylogenetic tree. We show that taking certain linear
combinations of the associated random variables (site pattern counts) reduces
the dimensionality of the model from exponential in the number of extant taxa,
to quadratic in the number of taxa, while retaining the ability to
statistically identify phylogenetic divergence events. A key feature is the
identification of an invariant subspace which depends only bilinearly on the
model parameters, in contrast to the usual multi-linear dependence in the full
space. We discuss potential applications including the computation of split
(edge) weights on phylogenetic trees from observed sequence data.


  While the individuals chosen for a genome-wide association study (GWAS) may
not be closely related to each other, there can be distant (cryptic)
relationships that confound the evidence of disease association. These cryptic
relationships violate the GWAS assumption regarding the independence of the
subjects' genomes, and failure to account for these relationships results in
both false positives and false negatives.
  This paper presents a method to correct for these cryptic relationships. We
accurately detect distant relationships using an expectation maximization (EM)
algorithm for finding the identity coefficients from genotype data with know
prior knowledge of relationships. From the identity coefficients we compute the
kinship coefficients and employ a kinship-corrected association test.
  We assess the accuracy of our EM kinship estimation algorithm. We show that
on genomes simulated from a Wright-Fisher pedigree, our method converges
quickly and requires only a relatively small number of sites to be accurate. We
also demonstrate that our kinship coefficient estimates outperform
state-of-the-art covariance-based approaches and PLINK's kinship estimate.
  To assess the kinship-corrected association test, we simulated individuals
from deep pedigrees and drew one site to recessively determine the disease
status. We applied our EM algorithm to estimate the kinship coefficients and
ran a kinship-adjusted association test. Our approach compares favorably with
the state-of-the-art and far out-performs a na\"ive association test.
  We advocate use of our method to detect cryptic relationships and for
correcting association tests. Not only is our model easy to interpret due to
the use of identity states as latent variables, but also inference provides
state-of-the-art accuracy.


  Proteins change their charge state through protonation and redox reactions as
well as through binding charged ligands. The free energy of these reactions are
dominated by solvation and electrostatic energies and modulated by protein
conformational relaxation in response to the ionization state changes. Although
computational methods for calculating these interactions can provide very
powerful tools for predicting protein charge states, they include several
critical approximations of which users should be aware. This chapter discusses
the strengths, weaknesses, and approximations of popular computational methods
for predicting charge states and understanding their underlying electrostatic
interactions. The goal of this chapter is to inform users about applications
and potential caveats of these methods as well as outline directions for future
theoretical and computational research.


  Structural relationships among proteins are important in the study of their
evolution as well as in drug design and development. The protein 3D structure
has been shown to be effective with respect to classifying proteins. Prior work
has shown that the Double Centroid Reduced Representation (DCRR) model is a
useful geometric representation for protein structure with respect to visual
models, reducing the quantity of modeled information for each amino acid, yet
retaining the most important geometrical and chemical features of each: the
centroids of the backbone and of the side-chain. DCRR has not yet been applied
in the calculation of geometric structural similarity. Meanwhile,
multi-dimensional indexing (MDI) of protein structure combines protein
structural analysis with distance metrics to facilitate structural similarity
queries and is also used for clustering protein structures into related groups.
In this respect, the combination of geometric models with MDI has been shown to
be effective. Prior work, notably Distance and Density based Protein Indexing
(DDPIn), applies MDI to protein models based on the geometry of the CA
backbone. DDPIn distance metrics are based on radial and density functions that
incorporate spherical-based metrics, and the indices are built from metric tree
(M-tree) structures. This work combines DCRR with DDPIn for the development of
new DCRR centroid-based metrics: spherical binning distance and inter-centroid
spherical distance. The use of DCRR models will provide additional significant
structural information via the inclusion of side-chain centroids. Additionally,
the newly developed distance metric functions combined with DCRR and M-tree
indexing attempt to improve upon the performance of prior work (DDPIn), given
the same data set, with respect to both individual k-nearest neighbor (kNN)
search queries as well as clustering all proteins in the index.


  Cancer Biotherapy (CB), as opposed to cancer chemotherapy, is the use of
macromolecular, biological agents instead of organic chemicals or drugs to
treat cancer. Biological agents usually have higher selectivity and have less
toxic side effects than chemical agents. The I.S.B.T.C., being the only major
information database for CB, seems lacking in some crucial information on
various cancer biotherapy regimens. It is thus necessary to have a
comprehensive curated CB database. The database accessible to cancer patients
and also should be a sounding board for scientific ideas by cancer researchers.
The database/web server has information about main families of cancer
biotherapy regimens to date, namely, Protein Kinase Inhibitors, Ras Pathway
Inhibitors, Cell-Cycle Active Agents, MAbs (monoclonal antibodies), ADEPT
(Antibody-Directed Enzyme Pro-Drug Therapy), Cytokines, Anti-Angiogenesis
Agents, Cancer Vaccines, Cell-based Immunotherapeutics, Gene Therapy,
Hematopoietic Growth Factors, Retinoids, and CAAT. For each biotherapy regimen,
we will extract the following attributes in populating the database: Cancer
type, Gene/s and gene product/s involved, Gene sequence, Organs affected,
Reference papers, Clinical phase/stage, Survival rate, Clinical test center
locations, Cost, Patient blogs, Researcher blogs, and Future work. The database
will be accessible to public through a website and had FAQs for making it
understandable to the laymen and discussion page for researchers to express
their views and ideas. In addition to information about the biotherapy
regimens, the website will link to other biologically significant databases
like structural proteomics, metabolomics, glycomics, and lipidomics databases,
as well as to news around the world regarding cancer therapy results. The
database attributes would be regularly updated for novel attributes as
discoveries are made.


  Based on overall 3D structure, proteins may be grouped into two broad
categories, namely, globular proteins (spheroproteins), and elongated or
rod-shaped proteins (RSP). The former comprises a significant majority of
proteins. This work concerns the second category. Unlike a spheroprotein, an
RSP possesses a conspicuous axis along its longest dimension. To take advantage
of this symmetry element, we decided to represent RSPs using cylindrical
coordinates, (rho, theta, z), with the z-axis as the main axis, and one tip of
the protein at the origin. A "tip" is one of two extreme points in the protein
lying along the protein axis along its longest dimension. We first identify the
two tips, T1 and T2, of the RSP using a protein graphics software, then
determine their (Cartesian) coordinates, (h, k, l) and (m, n, o), respectively.
Arbitrarily selecting T1 as the tip at the origin, we translate the protein by
subtracting (h, k, l) from all structural coordinates. We then find the angle
alpha between vector T1-T2 and the positive z-axis by computing the scalar
product of vectors T1- T2 and OP where P is an arbitrary point along the
positive z-axis. We typically use (0, 0, p) where p is a suitable positive
number. Then we compute the cross product of the two vectors to determine the
axis about which we should rotate vector T1-T2 to make it coincide with the
positive z-axis. We use a matrix form of Rodrigues' formula to do the rotation.
We then apply the Cartesian to cylindrical coordinate transformation equations
to the system. We have applied the above transformation to 15 RSPs: 1QCE, 2JJ7,
2KPE, 3K2A, 3LHP, 2LOE, 2L3H, 2L1P, 1KSG, 1KSJ, 1KSH, 2KOL, 2KZG, 2KPF and
3MQC. We have also created a web server that can take the PDB coordinate file
of an RSP and output its cylindrical coordinates. The URL of our web server
will be announced publicly in due course.


  We have developed a web server that implements two complementary methods to
quantify the depth of ligand binding site (LBS) in protein-ligand complexes:
the "secant plane" (SP) and "tangent sphere" (TS) methods. The protein
molecular centroid (global centroid, GC), and the LBS centroid (local centroid,
LC) are first determined. The SP is the plane passing through the LC and normal
to the line passing through the LC and the GC. The "exterior side" of the SP is
the side opposite GC. The TS is the sphere with center at GC and tangent to the
SP at LC. The percentage of protein atoms inside the TS (TS index) and on the
exterior side of the SP (SP index), are complementary measures of LBS depth.
The SPi is directly proportional to LBS depth while the TSi is inversely
proportional. We tested the two methods using a test set of 67
well-characterized protein-ligand structures (Laskowski, et al. 1996), as well
as that of an artificial protein in the form of a grid of points in the overall
shape of a sphere and in which LBS of any depth can be specified. Results from
both the SP and TS methods agree well with reported data (ibid.), and results
from the artificial case confirm that both methods are suitable measures of LBS
depth. The web server may be used in two modes. In the "ligand mode", user
inputs the protein PDB coordinates as well as those of the ligand. The "LBS
mode" is the same as the former, except that the ligand coordinates are assumed
to be unavailable; hence the user inputs what s/he believes to be the
coordinates of the LBS amino acid residues. In both cases, the web server
outputs the SP and TS indices. LBS depth is usually directly related to the
amount of conformational change a protein undergoes upon ligand binding -
ability to quantify it could allows meaningful comparison of protein
flexibility and dynamics. The URL of our web server will be announced publicly
in due course.


  Protein structure representation is an important tool in structural biology.
There exists different methods of representing the protein 3D structures and
different biologists favor different methods based on the information they
require. Currently there is no available method of protein 3D structure
representation which captures enough chemical information from the protein
sequence and clearly shows the intra-molecular and the inter-molecular H-bonds
and VDW interactions at the same time. This project aims to reduce the 3D
structure of a protein and display the reduced representation along with
intermolecular and the intra- molecular H-bonds and van der Waals interactions.
A reduced protein representation has a significantly lower atomicity (i.e.,
number of atom coordinates) than one which is in all- atom representation. In
this work, we transform the protein structure from all-atom representation"
(AAR) to double-centroid reduced representation (DCRR), which contains amino
acid backbone (N, CA, C', O) and side chain (CB and beyond) centroid
coordinates instead of atomic coordinates. Another aim of this project is to
develop a visualization interface for the reduced representation. This
interface is implemented in MATLAB and displays the protein in DCRR along with
its inter-molecular, as well as intra-molecular, interaction. Visually, DCRR is
easier to comprehend than AAR. We also developed a Web Server called the
Protein DCRR Web Server wherein users can enter the PDB id or upload a modeled
protein and get the DCRR of that protein. The back end to the Web Server is a
database which has the reduced representation for all the x-ray
crystallographic structure in the PDB.


  A pharmacophore consists of the parts of the structure of the ligand that are
sufficient to express the biological and pharmacological effects of the ligand.
It is usually a substructure of the entire structure of the ligand. Small
organic molecules called ligands or metabolites in the cell form complexes with
biomolecules (usually proteins) to serve different purposes. The sites at which
the ligands bind are known as ligand binding sites, which are essentially
"pockets" which have complementary shapes and patterns of charge distribution
with the ligands. Sometimes a pocket is induced by the ligand itself. If we
study different bound conformations of ligands it is found that they share a
specific three- dimensional pattern that is more or less common and is
responsible for its binding and which is complementary in three-dimensional
geometry and charge distribution pattern with its cognate binding site in the
protein. This work studies the three dimensional structure of the consensus
ligand binding site for the ligand FMN. A training set for the ligand binding
sites was made and a 3D consensus binding site motif was determined for FMN.
The FMN system was studied and its binding sites in its respective regulator
proteins. The ability to identify ligand binding site by scanning the 3D
binding site consensus motif in protein 3D structures is an important step in
drug target discovery. Once a pharmacophore template is found it can also be
used to design other potential molecules that can bind to it and thus serve as
novel drugs.


  Some methods aim to correct or test for relationships or to reconstruct the
pedigree, or family tree. We show that these methods cannot resolve ties for
correct relationships due to identifiability of the pedigree likelihood which
is the probability of inheriting the data under the pedigree model. This means
that no likelihood-based method can produce a correct pedigree inference with
high probability. This lack of reliability is critical both for health and
forensics applications.
  In this paper we present the first discussion of multiple typed individuals
in non-isomorphic pedigrees, $\mathcal{P}$ and $\mathcal{Q}$, where the
likelihoods are non-identifiable, $Pr[G~|~\mathcal{P},\theta] =
Pr[G~|~\mathcal{Q},\theta]$, for all input data $G$ and all recombination rate
parameters $\theta$. While there were previously known non-identifiable pairs,
we give an example having data for multiple individuals.
  Additionally, deeper understanding of the general discrete structures driving
these non-identifiability examples has been provided, as well as results to
guide algorithms that wish to examine only identifiable pedigrees. This paper
introduces a general criteria for establishing whether a pair of pedigrees is
non-identifiable and two easy-to-compute criteria guaranteeing identifiability.
Finally, we suggest a method for dealing with non-identifiable likelihoods: use
Bayes rule to obtain the posterior from the likelihood and prior. We propose a
prior guaranteeing that the posterior distinguishes all pairs of pedigrees.
  Shortened version published as: B. Kirkpatrick. Non-identifiable pedigrees
and a Bayesian solution. Int. Symp. on Bioinformatics Res. and Appl. (ISBRA),
7292:139-152 2012.


  We call a \emph{comb} a map $f:I\to [0,\infty)$, where $I$ is a compact
interval, such that $\{f\ge \varepsilon\}$ is finite for any $\varepsilon$. A
comb induces a (pseudo)-distance $\dtf$ on $\{f=0\}$ defined by $\dtf(s,t) =
\max_{(s\wedge t, s\vee t)} f$. We describe the completion $\bar I$ of
$\{f=0\}$ for this metric, which is a compact ultrametric space called
\emph{comb metric space}.
  Conversely, we prove that any compact, ultrametric space $(U,d)$ without
isolated points is isometric to a comb metric space. We show various examples
of the comb representation of well-known ultrametric spaces: the Kingman
coalescent, infinite sequences of a finite alphabet, the $p$-adic field and
spheres of locally compact real trees. In particular, for a rooted, locally
compact real tree defined from its contour process $h$, the comb isometric to
the sphere of radius $T$ centered at the root can be extracted from $h$ as the
depths of its excursions away from $T$.


  Phylogenetic reconstruction aims at finding plausible hypotheses of the
evolutionary history of genes or species based on genomic sequence information.
The distinction of orthologous genes (genes that having a common ancestry and
diverged after a speciation) is crucial and lies at the heart of many genomic
studies. However, existing methods that rely only on 1:1 orthologs to infer
species trees are strongly restricted to a small set of allowed genes that
provide information about the species tree. The use of larger gene sets that
consist in addition of non-orthologous genes (e.g. so-called paralogous or
xenologous genes) considerably increases the information about the evolutionary
history of the respective species. In this work, we introduce a novel method to
compute species phylogenies based on sequence data including orthologs,
paralogs or even xenologs.


  Molecular networks act as the backbone of cellular activities, providing an
{excellent} opportunity to understand the developmental changes in an organism.
While network data usually constitute only stationary network graphs,
constructing multilayer PPI network may provide clues to the particular
developmental role at each {stage of life} and may unravel the importance of
these developmental changes. The developmental biology model of {Caenorhabditis
elegans} {analyzed} here provides a ripe platform to understand the patterns of
evolution during life stages of an organism. In the present study, the widely
studied network properties exhibit overall similar statistics for all the PPI
layers. Further, the analysis of the degree-degree correlation and spectral
properties not only reveals crucial differences in each PPI layer but also
indicates the presence of the varying complexity among them. The PPI layer of
Nematode life stage exhibits various network properties different to rest of
the PPI layers, indicating the specific role of cellular diversity and
developmental transitions at this stage. The framework presented here provides
a direction to explore and understand developmental changes occurring in
different life stages of an organism.


  Deterministic population growth models with power-law rates can exhibit a
large variety of growth behaviors, ranging from algebraic, exponential to
hyperexponential (finite time explosion). In this setup, selfsimilarity
considerations play a key role, together with two time substitutions. Two
stochastic versions of such models are investigated, showing a much richer
variety of behaviors. One is the Lamperti construction of selfsimilar positive
stochastic processes based on the exponentiation of spectrally positive
processes , followed by an appropriate time change. The other one is based on
stable continuous-state branching processes, given by another Lamperti time
substitution applied to stable spectrally positive processes.


  The well-known Deficiency One Theorem gives structural conditions on a
chemical reaction network under which, for any set of parameter values, the
steady states of the corresponding mass action system may be easily
characterized. It is also known, however, that mass action systems are not
uniquely associated with reaction networks and that some representations may
satisfy the Deficiency One Theorem while others may not. In this paper we
present a mixed-integer linear programming framework capable of determining
whether a given mass action system has a dynamically equivalent or linearly
conjugate representation which has an underlying network satisfying the
Deficiency One Theorem. This extends recent computational work determining
linearly conjugate systems which are weakly reversible and have a deficiency of
zero.


  In multiple scientific and technological applications we face the problem of
having low dimensional data to be justified by a linear model defined in a high
dimensional parameter space. The difference in dimensionality makes the problem
ill-defined: the model is consistent with the data for many values of its
parameters. The objective is to find the probability distribution of parameter
values consistent with the data, a problem that can be cast as the exploration
of a high dimensional convex polytope. In this work we introduce a novel
algorithm to solve this problem efficiently. It provides results that are
statistically indistinguishable from currently used numerical techniques while
its running time scales linearly with the system size. We show that the
algorithm performs robustly in many abstract and practical applications. As
working examples we simulate the effects of restricting reaction fluxes on the
space of feasible phenotypes of a {\em genome} scale E. Coli metabolic network
and infer the traffic flow between origin and destination nodes in a real
communication network.


  In this paper, we study the dynamics of epidemic processes taking place in
adaptive networks of arbitrary topology. We focus our study on the adaptive
susceptible-infected-susceptible (ASIS) model, where healthy individuals are
allowed to temporarily cut edges connecting them to infected nodes in order to
prevent the spread of the infection. In this paper, we derive a closed-form
expression for a lower bound on the epidemic threshold of the ASIS model in
arbitrary networks with heterogeneous node and edge dynamics. For networks with
homogeneous node and edge dynamics, we show that the resulting \blue{lower
bound} is proportional to the epidemic threshold of the standard SIS model over
static networks, with a proportionality constant that depends on the adaptation
rates. Furthermore, based on our results, we propose an efficient algorithm to
optimally tune the adaptation rates in order to eradicate epidemic outbreaks in
arbitrary networks. We confirm the tightness of the proposed lower bounds with
several numerical simulations and compare our optimal adaptation rates with
popular centrality measures.


  Living systems, particularly multicellular systems, often seem hopelessly
complex. But recent studies have suggested that beneath this complexity, there
may be unifying quantitative principles that we are only now starting to
unravel. All cells interact with their environments and with other cells.
Communication among cells is a primary means for cells to interact with each
other. The complexity of these multicellular systems, due to the large numbers
of cells and the diversity of intracellular and intercellular interactions,
makes understanding multicellular systems a daunting task. To overcome this
challenge, we will likely need judicious simplifications and conceptual
frameworks that can reveal design principles that are shared among diverse
multicellular systems. Here we review some recent progress towards developing
such frameworks.


  Cells measure concentrations of external ligands by capturing ligand
molecules with cell surface receptors. The numbers of molecules captured by
different receptors co-vary because they depend on the same extrinsic ligand
fluctuations. However, these numbers also counter-vary due to the intrinsic
stochasticity of chemical processes because a single molecule randomly captured
by a receptor cannot be captured by another. Such structure of receptor
correlations is generally believed to lead to an increase in information about
the external signal compared to the case of independent receptors. We analyze a
solvable model of two molecular receptors and show that, contrary to this
widespread expectation, the correlations have a small and negative effect on
the information about the ligand concentration. Further, we show that
measurements that average over multiple receptors are almost as informative as
those that track the states of every individual one.


  Inside individual cells, expression of genes is stochastic across organisms
ranging from bacterial to human cells. A ubiquitous feature of stochastic
expression is burst-like synthesis of gene products, which drives considerable
intercellular variability in protein levels across an isogenic cell population.
One common mechanism by which cells control such stochasticity is negative
feedback regulation, where a protein inhibits its own synthesis. For a single
gene that is expressed in bursts, negative feedback can affect the burst
frequency or the burst size. In order to compare these feedback types, we study
a piecewise deterministic model for gene expression of a self-regulating gene.
Mathematically tractable steady-state protein distributions are derived and
used to compare the noise suppression abilities of the two feedbacks. Results
show that in the low noise regime, both feedbacks are similar in term of their
noise buffering abilities. Intriguingly, feedback in burst size outperforms the
feedback in burst frequency in the high noise regime. Finally, we discuss
various regulatory strategies by which cells implement feedback to control
burst sizes of expressed proteins at the level of single cells.


  We discuss a mathematical model of migration channel based on the truncated
Waring distribution. The truncated Waring distribution is obtained for a more
general model of motion of substance through a channel containing finite number
of boxes. The model is applied then for case of migrants moving through a
channel consisting of finite number of countries or cities. The number of
migrants in the channel strongly depends on the number of migrants that enter
the channel through the country of entrance. It is shown that if the final
destination country is very popular then large percentage of migrants may
concentrate there.


  The shortest common superstring (SCS) problem has been studied at great
length because of its connections to the de novo assembly problem in
computational genomics. The base problem is APX-complete, but several
generalizations of the problem have also been studied. In particular, previous
results include that SCS with Negative strings (SCSN) is in Log-APX (though
there is no known hardness result) and SCS with Wildcards (SCSW) is
Poly-APX-hard. Here, we prove two new hardness results: (1) SCSN is
Log-APX-hard (and therefore Log-APX-complete) by a reduction from Minimum Set
Cover and (2) SCS with Negative strings and Wildcards (SCSNW) is NPOPB-hard by
a reduction from Minimum Ones 3SAT.


  One of the most common analysis tasks in genomic research is to identify
genes that are differentially expressed (DE) between experimental conditions.
Empirical Bayes (EB) statistical tests using moderated genewise variances have
been very effective for this purpose, especially when the number of biological
replicate samples is small. The EB procedures can however be heavily influenced
by a small number of genes with very large or very small variances. This
article improves the differential expression tests by robustifying the
hyperparameter estimation procedure. The robust procedure has the effect of
decreasing the informativeness of the prior distribution for outlier genes
while increasing its informativeness for other genes. This effect has the
double benefit of reducing the chance that hypervariable genes will be
spuriously identified as DE while increasing statistical power for the main
body of genes. The robust EB algorithm is fast and numerically stable. The
procedure allows exact small-sample null distributions for the test statistics
and reduces exactly to the original EB procedure when no outlier genes are
present. Simulations show that the robustified tests have similar performance
to the original tests in the absence of outlier genes but have greater power
and robustness when outliers are present. The article includes case studies for
which the robust method correctly identifies and downweights genes associated
with hidden covariates and detects more genes likely to be scientifically
relevant to the experimental conditions. The new procedure is implemented in
the limma software package freely available from the Bioconductor repository.


  Classical models of biological invasions assess species spread in homogeneous
landscapes by assuming constant growth rates and random local movement.
Mounting evidence suggests, however, that demographic stochasticity,
environmental heterogeneity and non-random movement of individuals affect
considerably the spread dynamics. Here, we show that the dynamics of biological
invasions are controlled by the spatial heterogeneity of the resource
distribution. We show theoretically that increasing the landscape resource
autocorrelation length causes a reduction in the average speed of species
spread. Demographic stochasticity plays a key role in the slowdown, which is
streghtened when individuals can actively move towards resources. The reduction
in the front propagation speed is verified in laboratory microcosm experiments
with the flagellated protist Euglena gracilis by comparing spread in habitats
characterized by different resource heterogeneity. Our theoretical and
experimental findings highlight the need to account for the intrinsic
stochasticity of population dynamics to describe spread in spatially extended
landscapes, which are inevitably characterized by heterogeneous spatial
distributions of resources controlling vital rates. Our work identifies the
resource autocorrelation length as a key modulator and a simple measure of
landscape susceptibility to biological invasions, with implications for
predicting the characters of biological invasions within naturally
heterogeneous environmental corridors.


  Coregulation of the expression of groups of genes has been extensively
demonstrated empirically in bacterial and eukaryotic systems. Such coregulation
can arise through the use of shared regulatory motifs, which allow the
coordinated expression of modules (and module groups) of functionally related
genes across the genome. Coregulation can also arise through the physical
association of multi-gene complexes through chromosomal looping, which are then
transcribed together. We present a general formalism for modeling coregulation
rules in the framework of Random Boolean Networks (RBN), and develop specific
models for transcription factor networks with modular structure (including
module groups, and multi-input modules (MIM) with autoregulation) and
multi-gene complexes (including hierarchical differentiation between multi-gene
complex members). We develop a mean-field approach to analyse the stability of
large networks incorporating coregulation, and show that autoregulated MIM and
hierarchical gene-complex models can achieve greater stability than networks
without coregulation whose rules have matching activation frequency. We provide
further analysis of the stability of small networks of both kinds through
simulations. We also characterize several general properties of the transients
and attractors in the hierarchical coregulation model, and show using
simulations that the steady-state distribution factorizes hierarchically as a
Bayesian network in a Markov Jump Process analogue of the RBN model.


  Reaction coordinates are widely used throughout chemical physics to model and
understand complex chemical transformations. We introduce a definition of the
natural reaction coordinate, suitable for condensed phase and biomolecular
systems, as a maximally predictive one-dimensional projection. We then show
this criterion is uniquely satisfied by a dominant eigenfunction of an integral
operator associated with the ensemble dynamics. We present a new sparse
estimator for these eigenfunctions which can search through a large candidate
pool of structural order parameters and build simple, interpretable
approximations that employ only a small number of these order parameters.
Example applications with a small molecule's rotational dynamics and
simulations of protein conformational change and folding show that this
approach can filter through statistical noise to identify simple reaction
coordinates from complex dynamics.


  A common idiom in biology education states, "Eyes in the front, the animal
hunts. Eyes on the side, the animal hides." In this paper, we explore one
possible explanation for why predators tend to have forward-facing, high-acuity
visual systems. We do so using an agent-based computational model of evolution,
where predators and prey interact and adapt their behavior and morphology to
one another over successive generations of evolution. In this model, we observe
a coevolutionary cycle between prey swarming behavior and the predator's visual
system, where the predator and prey continually adapt their visual system and
behavior, respectively, over evolutionary time in reaction to one another due
to the well-known "predator confusion effect." Furthermore, we provide evidence
that the predator visual system is what drives this coevolutionary cycle, and
suggest that the cycle could be closed if the predator evolves a hybrid visual
system capable of narrow, high-acuity vision for tracking prey as well as
broad, coarse vision for prey discovery. Thus, the conflicting demands imposed
on a predator's visual system by the predator confusion effect could have led
to the evolution of complex eyes in many predators.


  A fundamental question in protein folding is whether the coil to globule
collapse transition occurs during the initial stages of folding (burst-phase)
or simultaneously with the protein folding transition. Single molecule
fluorescence resonance energy transfer (FRET) and small angle X-ray scattering
(SAXS) experiments disagree on whether Protein L collapse transition occurs
during the burst-phase of folding. We study Protein L folding using a
coarse-grained model and molecular dynamics simulations. The collapse
transition in Protein L is found to be concomitant with the folding transition.
In the burst-phase of folding, we find that FRET experiments overestimate
radius of gyration, $R_g$, of the protein due to the application of Gaussian
polymer chain end-to-end distribution to extract $R_g$ from the FRET
efficiency. FRET experiments estimate $\approx$ 6\AA \ decrease in $R_g$ when
the actual decrease is $\approx$ 3\AA \ on Guanidinium Chloride denaturant
dilution from 7.5M to 1M, and thereby suggesting pronounced compaction in the
protein dimensions in the burst-phase. The $\approx$ 3\AA \ decrease is close
to the statistical uncertainties of the $R_g$ data measured from SAXS
experiments, which suggest no compaction, leading to a disagreement with the
FRET experiments. The transition state ensemble (TSE) structures in Protein L
folding are globular and extensive in agreement with the $\Psi$-analysis
experiments. The results support the hypothesis that the TSE of single domain
proteins depend on protein topology, and are not stabilised by local
interactions alone.


  A crucial step in the early development of multicellular organisms involves
the establishment of spatial patterns of gene expression which later direct
proliferating cells to take on different cell fates. These patterns enable the
cells to infer their global position within a tissue or an organism by reading
out local gene expression levels. The patterning system is thus said to encode
positional information, a concept that was formalized recently in the framework
of information theory. Here we introduce a toy model of patterning in one
spatial dimension, which can be seen as an extension of Wolpert's paradigmatic
"French Flag" model, to patterning by several interacting, spatially coupled
genes subject to intrinsic and extrinsic noise. Our model, a variant of an
Ising spin system, allows us to systematically explore expression patterns that
optimally encode positional information. We find that optimal patterning
systems use positional cues, as in the French Flag model, together with
gene-gene interactions to generate combinatorial codes for position which we
call "Counter" patterns. Counter patterns can also be stabilized against noise
and variations in system size or morphogen dosage by longer-range spatial
interactions of the type invoked in the Turing model. The simple setup proposed
here qualitatively captures many of the experimentally observed properties of
biological patterning systems and allows them to be studied in a single,
theoretically consistent framework.


  On the basis of a mathematical model, we continue the study of the metabolic
Krebs cycle (or the tricarboxilic acid cycle). For the first time, we consider
its consistency and stability, which depend on the dissipation of a
transmembrane potential formed by the respiratory chain in the plasmatic
membrane of a cell. The phase-parametric characteristic of the dynamics of the
ATP level depending on a given parameter is constructed. The scenario of
formation of multiple autoperiodic and chaotic modes is presented. Poincar\'{e}
sections and mappings are constructed. The stability of modes and the
fractality of the obtained bifurcations are studied. The full spectra of
Lyapunov indices, divergences, KS-entropies, horizons of predictability, and
Lyapunov dimensionalities of strange attractors are calculated. Some
conclusions about the structural-functional connections determining the
dependence of the cell respiration cyclicity on the synchronization of the
functioning of the tricarboxilic acid cycle and the electron transport chain
are presented.


  Background: Animal trade plays an important role for the spread of infectious
diseases in livestock populations. As a case study, we consider pig trade in
Germany, where trade actors (agricultural premises) form a complex network. The
central question is how infectious diseases can potentially spread within the
system of trade contacts. We address this question by analyzing the underlying
network of animal movements.
  Methodology/Findings: The considered pig trade dataset spans several years
and is analyzed with respect to its potential to spread infectious diseases.
Focusing on measurements of network-topological properties, we avoid the usage
of external parameters, since these properties are independent of specific
pathogens. They are on the contrary of great importance for understanding any
general spreading process on this particular network. We analyze the system
using different network models, which include varying amounts of information:
(i) static network, (ii) network as a time series of uncorrelated snapshots,
(iii) temporal network, where causality is explicitly taken into account.
  Findings: Our approach provides a general framework for a
topological-temporal characterization of livestock trade networks. We find that
a static network view captures many relevant aspects of the trade system, and
premises can be classified into two clearly defined risk classes. Moreover, our
results allow for an efficient allocation strategy for intervention measures
using centrality measures. Data on trade volume does barely alter the results
and is therefore of secondary importance. Although a static network description
yields useful results, the temporal resolution of data plays an outstanding
role for an in-depth understanding of spreading processes. This applies in
particular for an accurate calculation of the maximum outbreak size.


  This study applies dynamical and statistical modeling techniques to quantify
the proliferation and popularity of trending hashtags on Twitter. Using
time-series data reflecting actual tweets in New York City and San Francisco,
we present estimates for the dynamics (i.e., rates of infection and recovery)
of several hundred trending hashtags using an epidemic modeling framework
coupled with Bayesian Markov Chain Monte Carlo (MCMC) methods. This
methodological strategy is an extension of techniques traditionally used to
model the spread of infectious disease. We demonstrate that in some models,
hashtags can be grouped by infectiousness, possibly providing a method for
quantifying the trendiness of a topic.


  We adapt Forman's discretization of Ricci curvature to the case of undirected
networks, both weighted and unweighted, and investigate the measure in a
variety of model and real-world networks. We find that most nodes and edges in
model and real networks have a negative curvature. Furthermore, the
distribution of Forman curvature of nodes and edges is narrow in random and
small-world networks, while the distribution is broad in scale-free and
real-world networks. In most networks, Forman curvature is found to display
significant negative correlation with degree and centrality measures. However,
Forman curvature is uncorrelated with clustering coefficient in most networks.
Importantly, we find that both model and real networks are vulnerable to
targeted deletion of nodes with highly negative Forman curvature. Our results
suggest that Forman curvature can be employed to gain novel insights on the
organization of complex networks.


  In November 2015, El Salvador reported their first case of Zika virus (Zv)
leading to an explosive outbreak that in just two months had over 6000
suspected cases. Many communities along with national agencies initiated the
process to implement control measures that ranged from vector control and the
use of repellents to the suggestion of avoiding pregnancies for two years, the
latter one, in response to the growing number of microcephaly cases in Brazil.
In our study, we explore the impact of short term mobility between two
idealized interconnected communities where disparities and violence contribute
to the Zv epidemic. Using a Lagrangian modeling approach in a two-patch
setting, it is shown via simulations that short term mobility may be beneficial
in the control of a Zv outbreak when risk is relative low and patch disparities
are not too extreme. However, when the reproductive number is too high, there
seems to be no benefits. This paper is dedicated to the inauguration of the
Centro de Modelamiento Matem\'{a}tico Carlos Castillo-Ch\'{a}vez at Universidad
Francisco Gavidia in San Salvador, El Salvador.


  Staphylococcus aureus is a coccal bacterium, which divides by binary fission.
After division the cells remain attached giving rise to small clusters, with a
characteristic 'bunch of grapes' morphology. S. aureus is an important human
pathogen and this, combined with the increasing prevalence of
antibiotic-resistant strains, such as Methicillin Resistant S. aureus (MRSA),
make it an excellent subject for studies of new methods of antimicrobial
action. Many antibiotics, such as penicillin, prevent S. aureus cell division
and so an understanding of this fundamental process may pave the way to the
identification of novel drugs. We present here a novel image analysis framework
for automated detection and segmentation of cells in S. aureus clusters, and
identification of their cell division planes. We demonstrate the technique on
GFP labelled EzrA, a protein that localises to a mid-cell plane during division
and is involved in regulation of cell size and division. The algorithms may
have wider applicability in detecting morphologically complex structures of
fluorescently-labelled proteins within cells in other cell clusters.


  Flies that walk in a covered planar arena on straight paths avoid colliding
with each other, but which of the two flies stops is not random.
High-throughput video observations, coupled with dedicated experiments with
controlled robot flies have revealed that flies utilize the type of optic flow
on their retina as a determinant of who should stop, a strategy also used by
ship captains to determine which of two ships on a collision course should
throw engines in reverse. We use digital evolution to test whether this
strategy evolves when collision avoidance is the sole penalty. We find that the
strategy does indeed evolve in a narrow range of cost/benefit ratios, for
experiments in which the "regressive motion" cue is error free. We speculate
that these stringent conditions may not be sufficient to evolve the strategy in
real flies, pointing perhaps to auxiliary costs and benefits not modeled in our
study


  Recent magnetic tweezers experiments have reported systematic deviations of
the twist response of double-stranded DNA from the predictions of the twistable
worm-like chain model. Here we show, by means of analytical results and
computer simulations, that these discrepancies can be resolved if a coupling
between twist and bend is introduced. We obtain an estimate of 40 $\pm$ 10 nm
for the twist-bend coupling constant. Our simulations are in good agreement
with high-resolution, magnetic-tweezers torque data. Although the existence of
twist-bend coupling was predicted long ago (Marko and Siggia, Macromolecules
27, 981 (1994)), its effects on the mechanical properties of DNA have been so
far largely unexplored. We expect that this coupling plays an important role in
several aspects of DNA statics and dynamics.


  We present evidence of a novel form of group hunting. Individual sailfish
(Istiophorus platypterus) alternate attacks with other group members on their
schooling prey (Sardinella aurita). While only 24% of attacks result in prey
capture, multiple prey are injured in 95% of attacks, resulting in an increase
of injured fish in the school with the number of attacks. How quickly prey are
captured is positively correlated with the level of injury of the school,
suggesting that hunters can benefit from other conspecifics' attacks on the
prey. To explore this, we built a mathematical model capturing the dynamics of
the hunt. We show that group hunting provides major efficiency gains (prey
caught per unit time) for individuals in groups of up to 70 members. We also
demonstrate that a free riding strategy, where some individuals wait until the
prey are sufficiently injured before attacking, is only beneficial if the cost
of attacking is high, and only then when waiting times are short. Our findings
provide evidence that cooperative benefits can be realised through the
facilitative effects of individuals' hunting actions without spatial
coordination of attacks. Such "proto-cooperation" may be the pre-cursor to more
complex group-hunting strategies.


  Background: The increasing prevalence of food allergies requires development
of specific and sensitive tests capable of identifying the allergen responsible
for the disease. The development of serologic tests that can detect specific
IgE antibodies to allergenic proteins would therefore be highly received.
Results: Here we present two new quantitative immuno-PCR assays for the
sensitive detection of antibodies specific to the shrimp allergen tropomyosin.
Both assays are based on the self-assembling Tus-Ter-lock protein-DNA
conjugation system. Significantly elevated levels of tropomyosin-specific IgE
were detected in sera from patients allergic to shrimp. Conclusions: This is
the first time an allergenic protein has been fused with Tus to enable specific
IgE antibody detection in human sera by quantitative immuno-PCR.


  We would like to show how biological grouping, especially in the case of
species formation, is emerged through a nature of interactive populations with
a number theory. First, we are able to define a species as a $p$-Sylow subgroup
of a particular community in a single niche, confirmed by topological analysis.
We named this model the patch with zeta dominance (PzDom) model. Next, the
topological nature of the system is carefully examined. We confirm the
induction of hierarchy and time through a one-dimensional probability space
with certain topologies. For further clarification of induced fractals
including the relation to renormalization, a theoretical development is
proposed based on a newly identified fact, namely that scaling parameters for
magnetization analogs exactly correspond to imaginary parts of the Riemann zeta
function's nontrivial zeros. In our PzDom model, calculations only require
knowledge of the density of individuals over time.


  Mechanisms leading to speciation are a major focus in evolutionary biology.
In this paper, we present and study a stochastic model of population where
individuals, with type a or A, are equivalent from ecological, demographical
and spatial points of view, and differ only by their mating preference: two
individuals with the same genotype have a higher probability to mate and
produce a viable offspring. The population is subdivided in several patches and
individuals may migrate between them. We show that mating preferences by
themselves, even if they are very small, are enough to entail reproductive
isolation between patches, and we provide the time needed for this isolation to
occur as a function of the population size. Our results rely on a fine study of
the stochastic process and of its deterministic limit in large population,
which is given by a system of coupled nonlinear differential equations.
Besides, we propose several generalisations of our model, and prove that our
findings are robust for those generalisations.


  Mathematical models are often regarded as recent innovations in the
description and analysis of infectious disease outbreaks and epidemics, but
simple models have been in use for projection of epidemic trajectories for more
than a century. We recently described a single equation model (the incidence
decay with exponential adjustment, or IDEA, model) that can be used for short
term forecasting. In the mid-19th century, Dr. William Farr developed a single
equation approach (Farr's law) for epidemic forecasting. We show here that the
two models are in fact identical, and can be expressed in terms of one another,
and also in terms of a susceptible-infectious-removed (SIR) compartmental model
with improving control. This demonstrates that the concept of the reproduction
number, R0, is implicit to Farr's (pre-microbial era) work, and also suggests
that control of epidemics, whether via behavior change or intervention, is as
integral to the natural history of epidemics as is the dynamics of disease
transmission.


  The time taken for gene expression varies not least because proteins vary in
length considerably. This paper uses an abstract, tuneable Boolean regulatory
network model to explore gene expression time variation. In particular, it is
shown how non-uniform expression times can emerge under certain conditions
through simulated evolution. That is, gene expression time variance appears
beneficial in the shaping of the dynamical behaviour of the regulatory network
without explicit consideration of protein function.


  Early estimates of the transmission potential of emerging and re-emerging
infections are increasingly used to inform public health authorities on the
level of risk posed by outbreaks. Existing methods to estimate the reproduction
number generally assume exponential growth in case incidence in the first few
disease generations, before susceptible depletion sets in. In reality,
outbreaks can display sub-exponential (i.e., polynomial) growth in the first
few disease generations, owing to clustering in contact patterns, spatial
effects, inhomogeneous mixing, reactive behavior changes, or other mechanisms.
Here, we introduce the generalized growth model to characterize the early
growth profile of outbreaks and estimate the effective reproduction number,
with no need for explicit assumptions about the shape of epidemic growth. We
demonstrate this phenomenologic approach using analytical results and
simulations from mechanistic models, and provide validation against a range of
empirical disease datasets. Our results suggest that sub-exponential growth in
the early phase of an epidemic is the rule rather the exception. For empirical
outbreaks, the generalized-growth model consistently outperforms the
exponential model for a variety of directly and indirectly transmitted diseases
datasets with model estimates supporting sub-exponential growth dynamics. The
rapid decline in effective reproduction number predicted by analytical results
and observed in real and synthetic datasets within 3-5 disease generations
contrasts with the expectation of invariant reproduction number in epidemics
obeying exponential growth. Overall, our approach promotes a more reliable and
data-driven characterization of the early epidemic phase, which is important
for accurate estimation of the reproduction number and prediction of disease
impact.


  Preventing infectious disease like flu from spreading to large communities is
one of the most important issues for humans. One effective strategy is
voluntary vaccination, however, there is always the temptation for people
refusing to be vaccinated because once herd immunity is achieved, infection
risk is greatly reduced. In this paper, we study the effect of social impact on
the vaccination behavior resulting in preventing infectious disease in
networks. The evolutionary simulation results show that the social impact has
both positive and negative effects on the vaccination behavior. Especially, in
heterogeneous networks, if the vaccination cost is low, the behavior is more
promoted than the case without social impact. In contrast, if the cost is high,
the behavior is reduced compared to the case without social impact. Moreover,
the vaccination behavior is effective in heterogeneous networks more than in
homogeneous networks. This implies that the social impact puts people at risk
in homogeneous networks. We also evaluate the results from the social cost
related to the vaccination policy.


  The study of cellular signalling pathways and their deregulation in disease
states, such as cancer, is a large and extremely complex task. Indeed, these
systems involve many parts and processes but are studied piecewise and their
literatures and data are consequently fragmented, distributed and sometimes--at
least apparently--inconsistent. This makes it extremely difficult to build
significant explanatory models with the result that effects in these systems
that are brought about by many interacting factors are poorly understood.
  The rule-based approach to modelling has shown some promise for the
representation of the highly combinatorial systems typically found in
signalling where many of the proteins are composed of multiple binding domains,
capable of simultaneous interactions, and/or peptide motifs controlled by
post-translational modifications. However, the rule-based approach requires
highly detailed information about the precise conditions for each and every
interaction which is rarely available from any one single source. Rather, these
conditions must be painstakingly inferred and curated, by hand, from
information contained in many papers--each of which contains only part of the
story.
  In this paper, we introduce a graph-based meta-model, attuned to the
representation of cellular signalling networks, which aims to ease this massive
cognitive burden on the rule-based curation process. This meta-model is a
generalization of that used by Kappa and BNGL which allows for the flexible
representation of knowledge at various levels of granularity. In particular, it
allows us to deal with information which has either too little, or too much,
detail with respect to the strict rule-based meta-model. Our approach provides
a basis for the gradual aggregation of fragmented biological knowledge
extracted from the literature into an instance of the meta-model from which we
can define an automated translation into executable Kappa programs.


  Protein allostery requires dynamical structural correlations. Physical origin
of which, however, remain elusive despite intensive studies during last two
decades. Based on analysis of molecular dynamics (MD) simulation trajectories
for ten proteins with different sizes and folds, we found that nonlinear
backbone torsional pair (BTP) correlations, which are spatially more
long-ranged and are mainly executed by loop residues, exist extensively in most
analyzed proteins. Examination of torsional motion for correlated BTPs
suggested that aharmonic torsional state transitions are essential for such
non-linear correlations, which correspondingly occur on widely different and
relatively longer time scales. In contrast, BTP correlations between backbone
torsions in stable $\alpha$ helices and $\beta$ strands are mainly linear and
spatially more short-ranged, and are more likely to associate with intra-well
torsional dynamics. Further analysis revealed that the direct cause of
non-linear contributions are heterogeneous, and in extreme cases canceling,
linear correlations associated with different torsional states of participating
torsions. Therefore, torsional state transitions of participating torsions for
a correlated BTP are only necessary but not sufficient condition for
significant non-linear contributions. These findings implicate a general search
strategy for novel allosteric modulation of protein activities. Meanwhile, it
was suggested that ensemble averaged correlation calculation and static contact
network analysis, while insightful, are not sufficient to elucidate mechanisms
underlying allosteric signal transmission in general, dynamical and time scale
resolved analysis are essential.


  Free energy calculation is critical in predictive tasks such as protein
folding, docking and design. However, rigorous calculation of free energy
change is prohibitively expensive in these practical applications. The minimum
potential energy is therefore widely utilized to approximate free energy. In
this study, based on analysis of extensive molecular dynamics (MD) simulation
trajectories of a few native globular proteins, we found that change of minimum
and corresponding maximum potential energy terms exhibit similar level of
correlation with change of free energy. More importantly, we demonstrated that
change of span (maximum - minimum) of potential energy terms, which engender
negligible additional computational cost, exhibit considerably stronger
correlations with change of free energy than the corresponding change of
minimum and maximum potential energy terms. Therefore, potential energy span
may serve as an alternative efficient approximate free energy proxy.


  Certain biological reactions, such as receptor-ligand binding at cell-cell
interfaces and macromolecules binding to biopolymers, require many smaller
molecules crowding a reaction site to be cleared. Examples include the T cell
interface, a key player in immunological information processing. Diffusion sets
a limit for such cavitation to occur spontaneously, thereby defining a
timescale below which active mechanisms must take over. We consider $N$
independent diffusing particles in a closed domain, containing a sub-region
with $N_{0}$ particles, on average. We investigate the time until the
sub-region is empty, allowing a subsequent reaction to proceed. The first
passage time is computed using an efficient exact simulation algorithm and an
asymptotic approximation in the limit that cavitation is rare. In this limit,
we find that the mean first passage time is sub-exponential, $T \propto
e^{N_{0}}/N_{0}^2$. For the case of T cell receptors, we find that stochastic
cavitation is exceedingly slow, $10^9$ seconds at physiological densities,
however can be accelerated to occur within 5 second with only a four-fold
dilution.


  We investigate the formation of opinion against authority in an authoritarian
society composed of agents with different levels of authority. We explore a
"dissenting" opinion, held by lower-ranking, obedient, or less authoritative
people, spreading in an environment of an "affirmative" opinion held by
authoritative leaders. A real-world example would be a corrupt society where
people revolt against such leaders, but it can be applied to more general
situations. In our model, agents can change their opinion depending on their
authority relative to their neighbors and their own confidence level. In
addition, with a certain probability, agents can override the affirmative
opinion to take the dissenting opinion of a neighbor. Based on analytic
derivation and numerical simulations, we observe that both the network
structure and heterogeneity in authority, and their correlation, significantly
affect the possibility of the dissenting opinion to spread through the
population. In particular, the dissenting opinion is suppressed when the
authority distribution is very heterogeneous and there exists a positive
correlation between the authority and the number of neighbors of people
(degree). Except for such an extreme case, though, spreading of the dissenting
opinion takes place when people have the tendency to override the authority to
hold the dissenting opinion, but the dissenting opinion can take a long time to
spread to the entire society, depending on the model parameters. We argue that
the internal social structure of agents sets the scale of the time to reach
consensus, based on the analysis of the underlying structural properties of
opinion spreading.