  Background. Biological networks have growing importance for the
interpretation of high-throughput "omics" data. Statistical and combinatorial
methods allow to obtain mechanistic insights through the extraction of smaller
subnetwork modules. Further enrichment analyses provide set-based annotations
of these modules.
  Results. We present eXamine, a set-oriented visual analysis approach for
annotated modules that displays set membership as contours on top of a
node-link layout. Our approach extends upon Self Organizing Maps to
simultaneously lay out nodes, links, and set contours.
  Conclusions. We implemented eXamine as a freely available Cytoscape app.
Using eXamine we study a module that is activated by the virally-encoded
G-protein coupled receptor US28 and formulate a novel hypothesis about its
functioning.


  Collective behavior of proteins on biomembranes is usually studied within the
spontaneous curvature model. Here we consider an alternative phenomenological
approach, which accounts consistently for partial ordering of proteins as well
as the anchoring forces exerted on a membrane by layer of proteins. We show
analytically that such anisotropic interactions can drive membrane bending,
resulting in non-trivial equilibrium morphologies. The predicted instabilities
can advance our conceptual understanding of physical mechanisms behind
collective phenomena in biological systems, in particular those with inherent
anisotropy.


  Highly symmetric nano-shells are found in many biological systems, such as
clathrin cages and viral shells. Several studies have shown that symmetric
shells appear in nature as a result of the free energy minimization of a
generic interaction between their constituent subunits. We examine the physical
basis for the formation of symmetric shells, and using a minimal model we
demonstrate that these structures can readily grow from identical subunits
under non equilibrium conditions. Our model of nano-shell assembly shows that
the spontaneous curvature regulates the size of the shell while the mechanical
properties of the subunit determines the symmetry of the assembled structure.
Understanding the minimum requirements for the formation of closed nano-shells
is a necessary step towards engineering of nano-containers, which will have far
reaching impact in both material science and medicine.


  Pseudomonas aeruginosa is a metabolically flexible member of the
Gammaproteobacteria. Under anaerobic conditions and the presence of nitrate, P.
aeruginosa can perform (complete) denitrification, a respiratory process of
dissimilatory nitrate reduction to nitrogen gas via nitrite ($NO_{2}$), nitric
oxide ($NO$) and nitrous oxide ($N_{2}O$). This study focuses on understanding
the influence of environmental conditions on bacterial denitrification
performance, using a mathematical model of a metabolic network in P.
aeruginosa. To our knowledge, this is the first mathematical model of
denitrification for this bacterium. Analysis of the long-term behavior of the
network under changing concentration levels of oxygen ($O_{2}$), nitrate
($NO_{3}$), and phosphate ($PO_{4}$) suggests that $PO_{4}$ concentration
strongly affects denitrification performance. The model provides five
predictions on denitrification activity of P. aeruginosa under various
environmental conditions, and these predictions are either experimentally
validated or supported by pertinent biological literature. One motivation for
this study is to capture the effect of $PO_{4}$ on a denitrification metabolic
network of P. aeruginosa in order to shed light on mechanisms for greenhouse
gas $N_{2}O$ accumulation during seasonal oxygen depletion in aquatic
environments such as Lake Erie (Laurentian Great Lakes, USA).


  We introduce a model of dyadic social interactions and establish its
correspondence with relational models theory (RMT), a theory of human social
relationships. RMT posits four elementary models of relationships governing
human interactions, singly or in combination: Communal Sharing, Authority
Ranking, Equality Matching, and Market Pricing. To these are added the limiting
cases of asocial and null interactions, whereby people do not coordinate with
reference to any shared principle. Our model is rooted in the observation that
each individual in a dyadic interaction can do either the same thing as the
other individual, a different thing or nothing at all. To represent these three
possibilities, we consider two individuals that can each act in one out of
three ways toward the other: perform a social action X or Y, or alternatively
do nothing. We demonstrate that the relationships generated by this model
aggregate into six exhaustive and disjoint categories. We propose that four of
these categories match the four relational models, while the remaining two
correspond to the asocial and null interactions defined in RMT. We generalize
our results to the presence of N social actions. We infer that the four
relational models form an exhaustive set of all possible dyadic relationships
based on social coordination. Hence, we contribute to RMT by offering an answer
to the question of why there could exist just four relational models. In
addition, we discuss how to use our representation to analyze data sets of
dyadic social interactions, and how social actions may be valued and matched by
the agents.


  Every day tens of thousands of sequence searches and sequence alignment
queries are submitted to webservers. The capitalized word "BLAST" become a
verb, describing the act of performing sequence search and alignment. However,
if one needs to search for sequences that contain, for example, two hydrophobic
and three polar residues at five given positions, the query formation on the
most frequently used webservers will be difficult. Some servers support the
formation of queries with regular expressions, but most of the users are
unfamiliar with their syntax. Here we present an intuitive, easily applicable
webserver, the Protein Sequence Analysis server, that allows the formation of
multiple choice queries by simply drawing the residues to their positions; if
more than one residue are drawn to the same position, then they will be nicely
stacked on the user interface, indicating the multiple choice at he given
position. This computer-game like interface is natural and intuitive, and the
coloring of the residues makes possible to form queries requiring not just
certain amino acids in the given positions, but also small nonpolar, negatively
charged, hydrophobic, positively charged, or polar ones. The webserver is
available at http://psa.pitgroup.org.


  Background: Simulated nucleotide or amino acid sequences are frequently used
to assess the performance of phylogenetic reconstruction methods. BEAST, a
Bayesian statistical framework that focuses on reconstructing time-calibrated
molecular evolutionary processes, supports a wide array of evolutionary models,
but lacked matching machinery for simulation of character evolution along
phylogenies.
  Results: We present a flexible Monte Carlo simulation tool, called piBUSS,
that employs the BEAGLE high performance library for phylogenetic computations
within BEAST to rapidly generate large sequence alignments under complex
evolutionary models. piBUSS sports a user-friendly graphical user interface
(GUI) that allows combining a rich array of models across an arbitrary number
of partitions. A command-line interface mirrors the options available through
the GUI and facilitates scripting in large-scale simulation studies. Analogous
to BEAST model and analysis setup, more advanced simulation options are
supported through an extensible markup language (XML) specification, which in
addition to generating sequence output, also allows users to combine simulation
and analysis in a single BEAST run.
  Conclusions: piBUSS offers a unique combination of flexibility and
ease-of-use for sequence simulation under realistic evolutionary scenarios.
Through different interfaces, piBUSS supports simulation studies ranging from
modest endeavors for illustrative purposes to complex and large-scale
assessments of evolutionary inference procedures. The software aims at
implementing new models and data types that are continuously being developed as
part of BEAST/BEAGLE.


  In a recent paper, Wilmes et al. demonstrated a qualitative integration of
omics data streams to gain a mechanistic understanding of cyclosporine A
toxicity. One of their major conclusions was that cyclosporine A strongly
activates the nuclear factor (erythroid-derived 2)-like 2 pathway (Nrf2) in
renal proximal tubular epithelial cells exposed in vitro. We pursue here the
analysis of those data with a quantitative integration of omics data with a
differential equation model of the Nrf2 pathway. That was done in two steps:
(i) Modeling the in vitro pharmacokinetics of cyclosporine A (exchange between
cells, culture medium and vial walls) with a minimal distribution model. (ii)
Modeling the time course of omics markers in response to cyclosporine A
exposure at the cell level with a coupled PK-systems biology model. Posterior
statistical distributions of the parameter values were obtained by Markov chain
Monte Carlo sampling. Data were well simulated, and the known in vitro toxic
effect EC50 was well matched by model predictions. The integration of in vitro
pharmacokinetics and systems biology modeling gives us a quantitative insight
into mechanisms of cyclosporine A oxidative-stress induction, and a way to
predict such a stress for a variety of exposure conditions.


  The Shiga toxins comprise a family of related protein toxins secreted by
certain types of bacteria. Shigella dysenteriae, some strain of Escherichia
coli and other bacterias can express toxins which caused serious complication
during the infection. Shiga toxin and the closely related Shiga-like toxins
represent a group of very similar cytotoxins that may play an important role in
diarrheal disease and hemolytic-uremic syndrome. The outbreaks caused by this
toxin raised serious public health crisis and caused economic losses. These
toxins have the same biologic activities and according to recent studies also
share the same binding receptor, globotriosyl ceramide (Gb3). Rapid detection
of food contamination is therefore relevant for the containment of food-borne
pathogens. The conventional methods to detect pathogens, such as
microbiological and biochemical identification are time-consuming and
laborious. The immunological or nucleic acid-based techniques require extensive
sample preparation and are not amenable to miniaturization for on-site
detection. In the present are necessary of techniques of rapid identification,
simple and sensitive which can be employed in the countryside with
minimally-sophisticated instrumentation. Biosensors have shown tremendous
promise to overcome these limitations and are being aggressively studied to
provide rapid, reliable and sensitive detection platforms for such
applications.


  Multisite protein phosphorylation plays a prominent role in intracellular
processes like signal transduction, cell-cycle control and nuclear signal
integration. Many proteins are phosphorylated in a sequential and distributive
way at more than one phosphorylation site. Mathematical models of $n$-site
sequential distributive phosphorylation are therefore studied frequently. In
particular, in {\em Wang and Sontag, 2008,} it is shown that models of $n$-site
sequential distributive phosphorylation admit at most $2n-1$ steady states.
Wang and Sontag furthermore conjecture that for odd $n$, there are at most $n$
and that, for even $n$, there are at most $n+1$ steady states. This, however,
is not true: building on earlier work in {\em Holstein et.al., 2013}, we
present a scalar determining equation for multistationarity which will lead to
parameter values where a $3$-site system has $5$ steady states and parameter
values where a $4$-site system has $7$ steady states. Our results therefore are
counterexamples to the conjecture of Wang and Sontag. We furthermore study the
inherent geometric properties of multistationarity in $n$-site sequential
distributive phosphorylation: the complete vector of steady state ratios is
determined by the steady state ratios of free enzymes and unphosphorylated
protein and there exists a linear relationship between steady state ratios of
phosphorylated protein.


  Yeast cells produce daughter cells through a DNA replication and mitosis
cycle associated with checkpoints and governed by the cell cycle regulatory
network. To ensure genome stability and genetic information inheritance, this
regulatory network must be dynamically robust against various fluctuations.
Here we construct a simplified cell cycle model for a budding yeast to
investigate the underlying mechanism that ensures robustness in this process
containing sequential tasks (DNA replication and mitosis). We first establish a
three-variable model and select a parameter set that qualitatively describes
the yeast cell cycle process. Then, through nonlinear dynamic analysis, we
demonstrate that the yeast cell cycle process is an excitable system driven by
a sequence of saddle-node bifurcations with ghost effects. We further show that
the yeast cell cycle trajectory is globally attractive with modularity in both
state and parameter space, while the convergent manifold provides a suitable
control state for cell cycle checkpoints. These results not only highlight a
regulatory mechanism for executing successive cell cycle processes, but also
provide a possible strategy for the synthetic network design of sequential-task
processes.


  The uniform sampling of convex polytopes is an interesting computational
problem with many applications in inference from linear constraints, but the
performances of sampling algorithms can be affected by ill-conditioning. This
is the case of inferring the feasible steady states in models of metabolic
networks, since they can show heterogeneous time scales . In this work we focus
on rounding procedures based on building an ellipsoid that closely matches the
sampling space, that can be used to define an efficient hit-and-run (HR) Markov
Chain Monte Carlo. In this way the uniformity of the sampling of the convex
space of interest is rigorously guaranteed, at odds with non markovian methods.
We analyze and compare three rounding methods in order to sample the feasible
steady states of metabolic networks of three models of growing size up to
genomic scale. The first is based on principal component analysis (PCA), the
second on linear programming (LP) and finally we employ the lovasz ellipsoid
method (LEM). Our results show that a rounding procedure is mandatory for the
application of the HR in these inference problem and suggest that a combination
of LEM or LP with a subsequent PCA perform the best. We finally compare the
distributions of the HR with that of two heuristics based on the Artificially
Centered hit-and-run (ACHR), gpSampler and optGpSampler. They show a good
agreement with the results of the HR for the small network, while on genome
scale models present inconsistencies.


  The activation of the dodecameric Ca2+/calmodulin dependent kinase II
(CaMKII) holoenzyme is critical for the formation of memory. We now report that
CaMKII has a remarkable property, which is that activation of the holoenzyme
triggers the exchange of subunits with other holoenzymes, including unactivated
ones, enabling the calcium-independent phosphorylation and activation of new
subunits. We show, using a single-molecule TIRF microscopy technique, that the
exchange process is triggered by the activation of CaMKII, with the exchange
rate being modulated by phosphorylation of two residues in the
calmodulin-binding segment, Thr 305 and Thr 306. Based on these results, and on
the analysis of molecular dynamics simulations, we suggest that the
phosphorylated regulatory segment of CaMKII interacts with the central hub of
the holoenzyme and weakens its integrity, thereby promoting exchange. Our
results have implications for an earlier idea that subunit exchange in CaMKII
may have relevance for the storage of information resulting from brief
coincident stimuli during neuronal signaling.


  Growth in brain volume is one of the most spectacular changes in the hominid
lineage. The anthropological community agrees on that point. No consensus,
however, has been reached on selection pressures contributing to that growth.
In that respect Martin (1984) can be invoked. In his review of size
relationships among primates he stated that despite the relationship between
brain size, body size and feeding behavior no single interpretation could be
provided that revealed the causality of such relationship.
  This paper deals with one specific aspect of hominid brain growth; the fact
that for most of the hominid period, growth in brain volume was exponential in
character. To the best of our knowledge, no attempt has been made to identify a
selection mechanism that can facilitate just the exponential features of that
growth (as distinct from any of its other characteristics). It is broadly
accepted that the dynamics of this growth were peculiar. Growth was very fast,
or even rapid in the evolutionary scale of time. The most profound evidence of
that opinion was expressed by Haldane that this dramatic increase in brain size
was the most rapid evolutionary change known to him.


  In a molecular communication network, transmitters may encode information in
concentration or frequency of signalling molecules. When the signalling
molecules reach the receivers, they react, via a set of chemical reactions or a
molecular circuit, to produce output molecules. The counts of output molecules
over time is the output signal of the receiver. The aim of this paper is to
investigate the impact of different reaction types on the information
transmission capacity of molecular communication networks. We realise this aim
by using a general molecular circuit model. We derive general expressions of
mean receiver output, and signal and noise spectra. We use these expressions to
investigate the information transmission capacities of a number of molecular
circuits.


  A dendritic spine is a very small structure (~0.1 {\mu}m$^3$) of a neuron
that processes input timing information. Why are spines so small? Here, we
provide functional reasons; the size of spines is optimal for information
coding. Spines code input timing information by the probability of Ca$^{2+}$
increases, which makes robust and sensitive information coding possible. We
created a stochastic simulation model of input timing-dependent Ca$^{2+}$
increases in a cerebellar Purkinje cell's spine. Spines used probability coding
of Ca$^{2+}$ increases rather than amplitude coding for input timing detection
via stochastic facilitation by utilizing the small number of molecules in a
spine volume, where information per volume appeared optimal. Probability coding
of Ca$^{2+}$ increases in a spine volume was more robust against input
fluctuation and more sensitive to input numbers than amplitude coding of
Ca$^{2+}$ increases in a cell volume. Thus, stochasticity is a strategy by
which neurons robustly and sensitively code information.


  Marine viruses shape the structure of the microbial community. They are,
thus, a key determinant of the most important biogeochemical cycles in the
planet. Therefore, a correct description of the ecological and evolutionary
behavior of these viruses is essential to make reliable predictions about their
role in marine ecosystems. The infection cycle, for example, is indistinctly
modeled in two very different ways. In one representation, the process is
described including explicitly a fixed delay between infection and offspring
release. In the other, the offspring are released at exponentially distributed
times according to a fixed release rate. By considering obvious quantitative
differences pointed out in the past, the latter description is widely used as a
simplification of the former. However, it is still unclear how the dichotomy
"delay versus rate description" affects long-term predictions of host-virus
interaction models. Here, we study the ecological and evolutionary implications
of using one or the other approaches, applied to marine microbes. To this end,
we use mathematical and eco-evolutionary computational analysis. We show that
the rate model exhibits improved competitive abilities from both ecological and
evolutionary perspectives in steady environments. However, rate-based
descriptions can fail to describe properly long-term microbe-virus
interactions. Moreover, additional information about trade-offs between
life-history traits is needed in order to choose the most reliable
representation for oceanic bacteriophage dynamics. This result affects deeply
most of the marine ecosystem models that include viruses, especially when used
to answer evolutionary questions.


  The observation that, through a titration mechanism, microRNAs (miRNAs) can
act as mediators of effective interactions among their common targets
(competing endogenous RNAs or ceRNAs) has brought forward the idea ('ceRNA
hypothesis') that RNAs can regulate each other in extended 'cross-talk'
networks. Such an ability might play a major role in post-transcriptional
regulation (PTR) in shaping a cell's protein repertoire. Recent work focusing
on the emergent properties of the cross-talk networks has emphasized the high
flexibility and selectivity that may be achieved at stationarity. On the other
hand, dynamical aspects, possibly crucial on the relevant time scales, are far
less clear. We have carried out a dynamical study of the ceRNA hypothesis on a
model of PTR. Sensitivity analysis shows that ceRNA cross-talk is dynamically
extended, i.e. it may take place on time scales shorter than those required to
achieve stationairity even in cases where no cross-talk occurs in the steady
state, and is possibly amplified. Besides, in case of large, transfection-like
perturbations the system may develop strongly non-linear, threshold response.
Finally, we show that the ceRNA effect provides a very efficient way for a cell
to achieve fast positive shifts in the level of a ceRNA when necessary. These
results indicate that competition for miRNAs may indeed provide an elementary
mechanism to achieve system-level regulatory effects on the transcriptome over
physiologically relevant time scales.


  A proof is presented that gene regulatory networks (GRNs) based solely on
transcription factors cannot control the development of complex multicellular
life. GRNs alone cannot explain the evolution of multicellular life in the
Cambrian Explosion. Networks are based on addressing systems which are used to
construct network links. The more complex the network the greater the number of
links and the larger the required address space. It has been assumed that
combinations of transcription factors generate a large enough address space to
form GRNs that are complex enough to control the development of complex
multicellular life. However, it is shown in this article that transcription
factors do not have sufficient combinatorial power to serve as the basis of an
addressing system for regulatory control of genomes in the development of
complex organisms. It is proven that given $n$ transcription factor genes in a
genome and address combinations of length $k$ then there are at most $n/k$
k-length transcription factor addresses in the address space. The complexity of
embryonic development requires a corresponding complexity of control
information in the cell and its genome. Therefore, a different addressing
system must exist to form the complex control networks required for complex
control systems. It is postulated that a new type of network evolved based on
an RNA-DNA addressing system that utilized and subsumed the extant GRNs. These
new developmental control networks are called CENES (for Control genes). The
evolution of these new higher networks would explain how the Cambrian Explosion
was possible. The architecture of these higher level networks may in fact be
universal (modulo syntax) in the genomes of all multicellular life.


  The quasi-equilibrium approximation is acceptable when molecular interactions
are fast enough compared to circuit dynamics, but is no longer allowed when
cellular activities are governed by rare events. A typical example is the
lactose operon (lac), one of the most famous paradigms of transcription
regulation, for which several theories still coexist to describe its behaviors.
The lac system is generally analyzed by using equilibrium constants,
contradicting single-event hypotheses long suggested by Novick and Weiner
(1957). Enzyme induction as an all-or-none phenomenon. Proc. Natl. Acad. Sci.
USA 43, 553-566) and recently refined in the study of (Choi et al., 2008. A
stochastic single-molecule event triggers phenotype switching of a bacterial
cell. Science 322, 442-446). In the present report, a lac repressor
(LacI)-mediated DNA immunoprecipitation experiment reveals that the natural
LacI-lac DNA complex built in vivo is extremely tight and long-lived compared
to the time scale of lac expression dynamics, which could functionally
disconnect the abortive expression bursts and forbid using the standard modes
of lac bistability. As alternatives, purely kinetic mechanisms are examined for
their capacity to restrict induction through: (i) widely scattered derepression
related to the arrival time variance of a predominantly backward asymmetric
random walk and (ii) an induction threshold arising in a single window of
derepression without recourse to nonlinear multimeric binding and Hill
functions. Considering the complete disengagement of the lac repressor from the
lac promoter as the probabilistic consequence of a transient stepwise
mechanism, is sufficient to explain the sigmoidal lac responses as functions of
time and of inducer concentration. This sigmoidal shape can be misleadingly
interpreted as a phenomenon of equilibrium cooperativity classically used to
explain bistability, but which has been reported to be weak in this system.


  Cells sense and predict their environment via energy-dissipating pathways.
However, it is unclear whether dissipation helps or harms prediction. Here we
study dissipation and prediction for a minimal sensory module of receptors that
reversibly bind ligand. We find that the module performs short-term prediction
optimally when operating in an adiabatic regime where dissipation vanishes. In
contrast, beyond a critical forecast interval, prediction becomes most precise
in a regime of maximal dissipation, suggesting that dissipative sensing in
biological systems can serve to enhance prediction performance.


  A simple way to get insights about the possible functional motions of a
protein is to perform a normal mode analysis (NMA). Indeed, it has been shown
that low-frequency modes thus obtained are often closely related to domain
motions involved in protein function. Moreover, because protein low-frequency
modes are known to be robust, NMA can be performed using coarse-grained models.
As a consequence, it can be done for large ensembles of conformations as well
as for large systems, like the ribosome, whole virus capsids, etc.
Unexpectedly, on the high-frequency side, modes obtained with cutoff-based
coarse-grained models also seem able to provide useful insights on protein
dynamical properties.


  Policies ensuring that research data are available on public archives are
increasingly being implemented at the government [1], funding agency [2-4], and
journal [5,6] level. These policies are predicated on the idea that authors are
poor stewards of their data, particularly over the long term [7], and indeed
many studies have found that authors are often unable or unwilling to share
their data [8-11]. However, there are no systematic estimates of how the
availability of research data changes with time since publication. We therefore
requested datasets from a relatively homogenous set of 516 articles published
between 2 and 22 years ago, and found that availability of the data was
strongly affected by article age. For papers where the authors gave the status
of their data, the odds of a dataset being extant fell by 17% per year. In
addition, the odds that we could find a working email address for the first,
last or corresponding author fell by 7% per year. Our results reinforce the
notion that, in the long term, research data cannot be reliably preserved by
individual researchers, and further demonstrate the urgent need for policies
mandating data sharing via public archives.


  We are interested in a variation of the SIR (Susceptible/Infected/Recovered)
dynamics on the complete graph, in which infected individuals may only spread
to neighboring susceptible individuals at fixed rate $\lambda>0$ while
recovered individuals may only spread to neighboring infected individuals at
fixed rate $1$. This is also a variant of the so-called chase-escape process
introduced by Kordzakhia and then Bordenave & Ganesan. Our work is the first
study of this dynamics on complete graphs. Starting with one infected and one
recovered individuals on the complete graph with $N+2$ vertices, and stopping
the process when one type of individuals disappears, we study the asymptotic
behavior of the probability that the infection spreads to the whole graph as
$N\rightarrow\infty$ and show that for $\lambda\in (0,1)$ (resp. for
$\lambda>1$), the infection dies out (resp. does not die out) with probability
tending to one as $N\rightarrow\infty$, and that the probability that the
infection dies out tends to $1/2$ for $\lambda=1$. We also establish limit
theorems concerning the asymptotic state of the system in all regimes and show
that two additional phase transitions occur in the subcritical phase
$\lambda\in (0,1)$: at $\lambda=1/2$ the behavior of the expected number of
remaining infected individuals changes, while at $\lambda=(\sqrt {5}-1)/2$ the
behavior of the expected number of remaining recovered individuals changes. We
also study the outbreak sizes of the infection, and show that the outbreak
sizes are small if $\lambda \in(0,1/2]$, exhibit a power-law behavior for
$1/2<\lambda<1$, and are pandemic for $\lambda\geq 1$. Our method relies on
different couplings: we first couple the dynamics with two independent Yule
processes by using an Athreya-Karlin embedding, and then we couple the Yule
processes with Poisson processes thanks to Kendall's representation of Yule
processes.


  We study an atomic signaling game under stochastic evolutionary dynamics.
There is a finite number of players who repeatedly update from a finite number
of available languages/signaling strategies. Players imitate the most fit
agents with high probability or mutate with low probability. We analyze the
long-run distribution of states and show that, for sufficiently small mutation
probability, its support is limited to efficient communication systems. We find
that this behavior is insensitive to the particular choice of evolutionary
dynamic, a property that is due to the game having a potential structure with a
potential function corresponding to average fitness. Consequently, the model
supports conclusions similar to those found in the literature on language
competition. That is, we show that efficient languages eventually predominate
the society while reproducing the empirical phenomenon of linguistic drift. The
emergence of efficiency in the atomic case can be contrasted with results for
non-atomic signaling games that establish the non-negligible possibility of
convergence, under replicator dynamics, to states of unbounded efficiency loss.


  As the chief informational molecule of life, DNA is subject to extensive
physical manipulations. The energy required to deform double-helical DNA
depends on sequence, and this mechanical code of DNA influences gene
regulation, such as through nucleosome positioning. Here we examine the
sequence-dependent flexibility of DNA in bacterial transcription
factor-mediated looping, a context for which the role of sequence remains
poorly understood. Using a suite of synthetic constructs repressed by the Lac
repressor and two well-known sequences that show large flexibility differences
in vitro, we make precise statistical mechanical predictions as to how DNA
sequence influences loop formation and test these predictions using in vivo
transcription and in vitro single-molecule assays. Surprisingly,
sequence-dependent flexibility does not affect in vivo gene regulation. By
theoretically and experimentally quantifying the relative contributions of
sequence and the DNA-bending protein HU to DNA mechanical properties, we reveal
that bending by HU dominates DNA mechanics and masks intrinsic
sequence-dependent flexibility. Such a quantitative understanding of how
mechanical regulatory information is encoded in the genome will be a key step
towards a predictive understanding of gene regulation at single-base pair
resolution.


  Multistable gene regulatory systems sustain different levels of gene
expression under identical external conditions. Such multistability is used to
encode phenotypic states in processes including nutrient uptake and persistence
in bacteria, fate selection in viral infection, cell cycle control, and
development. Stochastic switching between different phenotypes can occur as the
result of random fluctuations in molecular copy numbers of mRNA and proteins
arising in transcription, translation, transport, and binding. However, which
component of a pathway triggers such a transition is generally not known. By
linking single-cell experiments on the lactose-uptake pathway in E. coli to
molecular simulations, we devise a general method to pinpoint the particular
fluctuation driving phenotype switching and apply this method to the transition
between the uninduced and induced states of the lac genes. We find that the
transition to the induced state is not caused only by the single event of
lac-repressor unbinding, but depends crucially on the time period over which
the repressor remains unbound from the lac-operon. We confirm this notion in
strains with a high expression level of the repressor (leading to shorter
periods over which the lac-operon remains unbound), which show a reduced
switching rate. Our techniques apply to multi-stable gene regulatory systems in
general and allow to identify the molecular mechanisms behind stochastic
transitions in gene regulatory circuits.


  In addition to the biologically active monomer of the protein Insulin
circulating in human blood, the molecule also exists in dimeric and hexameric
forms that are used as storage. The Insulin monomer contains two distinct
surfaces, namely the dimer forming surface (DFS) and the hexamer forming
surface (HFS) that are specifically designed to facilitate the formation of the
dimer and the hexamer, respectively. In order to characterize the structural
and dynamical behaviour of interfacial water molecules near these two surfaces
(DFS and HFS), we performed atomistic molecular dynamics simulations of Insulin
with explicit water. Dynamical characterization reveals that the structural
relaxation of the hydrogen bonds formed between the residues of DFS and the
interfacial water molecules is faster than those formed between water and that
of the HFS. Furthermore, the residence times of water molecules in the protein
hydration layer for both the DFS and HFS are found to be significantly higher
than those for some of the other proteins studied so far, such as HP-36 and
lysozyme. The surface topography and the arrangement of amino acid residues
work together to organize the water molecules in the hydration layer in order
to provide them with a preferred orientation. HFS having a large polar solvent
accessible surface area and a convex extensive nonpolar region, drives the
surrounding water molecules to acquire predominantly a clathrate-like
structure. In contrast, near the DFS, the surrounding water molecules acquire
an inverted orientation owing to the flat curvature of hydrophobic surface and
interrupted hydrophilic residual alignment. We have followed escape trajectory
of several such quasi-bound water molecules from both the surfaces and
constructed free energy surfaces of these water molecules.These free energy
surfaces reveal the differences between the two hydration layers.


  Virophages are viruses that rely on the replication machinery of other
viruses to reproduce within eukaryotic hosts. Two different modes of
coinfection have been posited based on experimental observation. In one mode,
the virophage and virus enter the host independently. In the other mode, the
virophage adheres to the virus so both virophage and virus enter the host
together. Here we ask: what are the ecological effects of these different modes
of coinfection? In particular, what ecological effects are common to both
infection modes, and what are the differences particular to each mode? We
develop a pair of biophysically motivated ODE models of viral-host population
dynamics, corresponding to dynamics arising from each mode of infection. We
find both modes of coinfection allow for the coexistence of the virophage,
virus, and host either at a stable fixed point or through cyclical dynamics. In
both models, virophage tend to be the most abundant population and their
presence always reduces the viral abundance and increases the host abundance.
However, we do find qualitative differences between models. For example, via
extensive sampling of biologically relevant parameter space, we only observe
bistability when the virophage and virus enter the host together. We discuss
how such differences may be leveraged to help identify modes of infection in
natural environments from population level data.


  Population structures can be crucial determinants of evolutionary processes.
For the Moran process on graphs certain structures suppress selective pressure,
while others amplify it (Lieberman et al. 2005 Nature 433 312-316).
Evolutionary amplifiers suppress random drift and enhance selection. Recently,
some results for the most powerful known evolutionary amplifier, the superstar,
have been invalidated by a counter example (D\'iaz et al. 2013 Proc. R. Soc. A
469 20130193). Here we correct the original proof and derive improved upper and
lower bounds, which indicate that the fixation probability remains close to
$1-1/(r^4 H)$ for population size $N\to\infty$ and structural parameter
$H\gg1$. This correction resolves the differences between the two
aforementioned papers. We also confirm that in the limit $N,H\to\infty$
superstars remain capable of eliminating random drift and hence of providing
arbitrarily strong selective advantages to any beneficial mutation. In
addition, we investigate the robustness of amplification in superstars and find
that it appears to be a fragile phenomenon with respect to changes in the
selection or mutation processes.


  Cancer is increasingly perceived as a systems-level, network phenomenon. The
major trend of malignant transformation can be described as a two-phase
process, where an initial increase of network plasticity is followed by a
decrease of plasticity at late stages of tumor development. The fluctuating
intensity of stress factors, like hypoxia, inflammation and the either
cooperative or hostile interactions of tumor inter-cellular networks, all
increase the adaptation potential of cancer cells. This may lead to the bypass
of cellular senescence, and to the development of cancer stem cells. We propose
that the central tenet of cancer stem cell definition lies exactly in the
indefinability of cancer stem cells. Actual properties of cancer stem cells
depend on the individual "stress-history" of the given tumor. Cancer stem cells
are characterized by an extremely large evolvability (i.e. a capacity to
generate heritable phenotypic variation), which corresponds well with the
defining hallmarks of cancer stem cells: the possession of the capacity to
self-renew and to repeatedly re-build the heterogeneous lineages of cancer
cells that comprise a tumor in new environments. Cancer stem cells represent a
cell population, which is adapted to adapt. We argue that the high evolvability
of cancer stem cells is helped by their repeated transitions between plastic
(proliferative, symmetrically dividing) and rigid (quiescent, asymmetrically
dividing, often more invasive) phenotypes having plastic and rigid networks.
Thus, cancer stem cells reverse and replay cancer development multiple times.
We describe network models potentially explaining cancer stem cell-like
behavior. Finally, we propose novel strategies including combination therapies
and multi-target drugs to overcome the Nietzschean dilemma of cancer stem cell
targeting: "what does not kill me makes me stronger".


  Genomic datasets generated with massively parallel sequencing methods have
the potential to propel systematics in new and exciting directions, but
selecting appropriate markers and methods is not straightforward. We applied
two approaches with particular promise for systematics, restriction site
associated DNA sequencing (RAD-Seq) and sequence capture (Seq-cap) of
ultraconserved elements (UCEs), to the same set of samples from a non-model,
Neotropical bird. We found that both RAD-Seq and Seq-cap produced genomic
datasets containing thousands of loci and SNPs and that the inferred population
assignments and species trees were concordant between datasets. However,
model-based estimates of demographic parameters differed between datasets,
particularly when we estimated the parameters using a method based on allele
frequency spectra. The differences we observed may result from differences in
assembly, alignment, and filtering of sequence data between methods, and our
findings suggest that caution is warranted when using allele frequencies to
estimate parameters from low-coverage sequencing data. We further explored the
differences between methods using simulated Seq-cap- and RAD-Seq-like datasets.
Analyses of simulated data suggest that increasing the number of loci from 500
to 5000 increased phylogenetic concordance factors and the accuracy and
precision of demographic parameter estimates, but increasing the number of loci
past 5000 resulted in minimal gains. Increasing locus length from 64 bp to 500
bp improved phylogenetic concordance factors and minimal gains were observed
with loci longer than 500 bp, but locus length did not influence the accuracy
and precision of demographic parameter estimates. We discuss our results
relative to the diversity of data collection methods available, and we provide
advice for harnessing next-generation sequencing for systematics research.


  Simulation testing is an important approach to evaluating fishery stock
assessment methods. In the last decade, the fisheries stock assessment modeling
framework Stock Synthesis (SS3) has become widely used around the world.
However, there lacks a generalized and scriptable framework for SS3 simulation
testing. Here, we introduce ss3sim, an R package that facilitates reproducible,
flexible, and rapid end-to-end simulation testing with SS3. ss3sim requires an
existing SS3 model configuration along with plain-text control files describing
alternative population dynamics, fishery properties, sampling scenarios, and
assessment approaches. ss3sim then generates an underlying 'truth' from a
specified operating model, samples from that truth, modifies and runs an
estimation model, and synthesizes the results. The simulations can be run in
parallel, reducing runtime, and the source code is free to be modified under an
open-source MIT license. ss3sim is designed to explore structural differences
between the underlying truth and assumptions of an estimation model, or between
multiple estimation model configurations. For example, ss3sim can be used to
answer questions about model misspecification, retrospective patterns, and the
relative importance of different types of fisheries data. We demonstrate the
software with an example, discuss how ss3sim complements other simulation
software, and outline specific research questions that ss3sim could address.


  The diversity of the arthropod fauna from the pacific coast in southern
Mexico was investigated in order to study and compare the biodiversity of two
different regions. One is located within the Biosphere reserve "La encrucijada"
and the second one on the edges of this reserve at a non-protected region. Both
regions are also ecologically different in terms of humidity and temperature.
To collect the material seven different trap methods were used, arthropods were
classified at order level, Diptera also at family level and of the Diptera
family Psychodidae, genera were also identified. Twenty-eight orders have been
identified, showing some regional differences in composition and abundance. As
a general conclusion we observe that no big statistical differences can be
observed at order or family level (only Diptera) between the protected area and
the non-protected area. This difference shows an unexpected greater diversity
in the non-protected area. At genus level (only Psychodidae) however, the
diversity of moth flies was slightly greater in the preserved area. We conclude
that climatic variations are probably playing an important role at low
systematic level (order and family) and that human intervention (e.g.
agriculture) becomes more important at genus level.


  Pyrrolizidine Alkaloids (PAs) are a group of naturally occurring alkaloids
that are produced by plants as a defense mechanism against insect herbivores.
The analytical methodologies employed for their detection have come a long way
since the first analytical experiment and in the last 30 years had an enormous
development, both technological and experimental. It is notorious that before
the generalization of certain technologies, especially in a post-war
atmosphere, most scientific researches relied on what it is today thin-layer
chromatography. Nevertheless this technique was not sufficient for accurately
measure quantities and unambiguously identify compounds, therefore
spectroscopic techniques arose as well as chromatographic techniques. While the
first never really coped with PAs analysis requirements the latter, either as
gas or liquid chromatography allowed the analysis of complex sample matrices.
Simultaneously, nuclear magnetic resonance also suffered significant
developments while mass spectrometry has become an attractive technique due to
increasingly higher maximum resolutions. The observed tendency in recent years,
in pyrrolizidine detection and quantification, as well as in many other areas,
is that hyphenated techniques are the chosen methods. A large number of papers
report multihyphenated methodologies, and the overwhelming majority relies on
gas or liquid chromatography.


  We sequenced genomes from a $\sim$7,000 year old early farmer from Stuttgart
in Germany, an $\sim$8,000 year old hunter-gatherer from Luxembourg, and seven
$\sim$8,000 year old hunter-gatherers from southern Sweden. We analyzed these
data together with other ancient genomes and 2,345 contemporary humans to show
that the great majority of present-day Europeans derive from at least three
highly differentiated populations: West European Hunter-Gatherers (WHG), who
contributed ancestry to all Europeans but not to Near Easterners; Ancient North
Eurasians (ANE), who were most closely related to Upper Paleolithic Siberians
and contributed to both Europeans and Near Easterners; and Early European
Farmers (EEF), who were mainly of Near Eastern origin but also harbored
WHG-related ancestry. We model these populations' deep relationships and show
that EEF had $\sim$44% ancestry from a "Basal Eurasian" lineage that split
prior to the diversification of all other non-African lineages.


  The fundamental processes of biological development are governed by multiple
signaling molecules that create non-uniform concentration profiles known as
morphogen gradients. It is widely believed that the establishment of morphogen
gradients is a result of complex processes that involve diffusion and
degradation of locally produced signaling molecules. We developed a
multi-dimensional discrete-state stochastic approach for investigating the
corresponding reaction-diffusion models. It provided a full analytical
description for stationary profiles and for important dynamic properties such
as local accumulation times, variances and mean first-passage times. The role
of discreteness in developing of morphogen gradients is analyzed by comparing
with available continuum descriptions. It is found that the continuum models
prediction about multiple time scales near the source region in two-dimensional
and three-dimensional systems is not supported in our analysis. Using ideas
that view the degradation process as an effective potential, the effect of
dimensionality on establishment of morphogen gradients is also discussed. In
addition, we investigated how these reaction-diffusion processes are modified
with changing the size of the source region.


  Multiple phenotypic states often arise in a single cell with different
gene-expression states that undergo transcription regulation with positive
feedback. Recent experiments have shown that at least in E. coli, the gene
state switching can be neither extremely slow nor exceedingly rapid as many
previous theoretical treatments assumed. Rather it is in the intermediate
region which is difficult to handle mathematically.Under this condition, from a
full chemical-master-equation description we derive a model in which the
protein copy-number, for a given gene state, follow a deterministic mean-field
description while the protein synthesis rates fluctuate due to stochastic
gene-state switching. The simplified kinetics yields a nonequilibrium landscape
function, which, similar to the energy function for equilibrium fluctuation,
provides the leading orders of fluctuations around each phenotypic state, as
well as the transition rates between the two phenotypic states. This rate
formula is analogous to Kramers theory for chemical reactions. The resulting
behaviors are significantly different from the two limiting cases studied
previously.


  The transfer of electrons and holes along DNA dimers, trimers and polymers is
described at the base-pair level, using the relevant on-site energies of the
base-pairs and the hopping parameters between successive base-pairs. The
temporal and spatial evolution of carriers along a $N$ base-pair DNA segment is
determined, solving a system of $N$ coupled differential equations. Useful
physical quantities are calculated including the pure mean carrier transfer
rate $k$, the inverse decay length $\beta$ used for exponential fit ($k = k_0
\textrm{exp}(-\beta d)$) of the transfer rate as a function of the charge
transfer distance $d = N \times$ 3.4 {\AA} and the exponent $\eta$ used for a
power law fit ($k = k_0' N^{-\eta}$) of the transfer rate as function of the
number of monomers $N$. Among others, the electron and hole transfer along the
polymers poly(dG)-poly(dC), poly(dA)-poly(dT), GCGCGC..., ATATAT... is studied.
$\beta$ ($\eta$) falls in the range $\approx$ 0.2 - 2 {\AA}$^{-1}$ (1.7 - 17),
$k_0$ ($k_0'$) is usually $\approx 10^{-2}$-10$^{-1}$ ($ 10^{-2}$-10$^{-1}$)
PHz although, generally, it falls in the wider range $\approx 10^{-4}$-10
($10^{-4}$-10$^3$) PHz. The results are compared with past predictions and
experiments. Our approach illustrates to which extent a specific DNA segment
can serve as an efficient medium for charge transfer.


  As many emerging and re-emerging infectious diseases are associated with food
animals, the relationship between available healthy food sources and population
health and social stability has become evident. A recent example of the
importance of this relationship was observed during the current flu pandemic.
This recent pandemic brought attention to novel target groups of susceptible
people at the interface of the animal and human populations. Veterinarians,
producers and processors are uniquely exposed to emerging zoonoses. Therefore
these individuals may serve as key sentinels and allow efficient evaluation of
the effectiveness of zoonoses prophylaxis and control, including evaluation of
the cost-effectiveness in the broader view. We also suggest some valuable
approaches for rapid diagnosis of emerging and re-emerging infectious diseases
and supportive systemic research which may address related ethical questions.
We also highly recommend more research investigations characterizing this
human/animal zoonosis interface, a potentially productive target for emerging
disease diagnosis and control.


  Predicting protein 3D structure from amino acid sequence remains as a
challenge in the field of computational biology. If protein structure
homologues are not found, one has to construct 3D structural conformations from
the very beginning by the so-called ab initio approach, using some empirical
energy functions. A successful algorithm in this category, Rosetta, creates an
ensemble of decoy conformations by assembling selected best short fragments of
known protein structures and then recognizes the native state as the highly
populated one with a very low energy. Typically, an energy function is a
combination of a variety of terms characterizing different structural features,
say hydrophobic interactions, van der Waals force, hydrogen bonding, etc. It is
critical for an energy function to be capable to distinguish native-like
conformations from non-native ones and to drive most initial conformations
assembled from fragments to a native-like one in a conformation search process.
In this paper we propose a linear programming algorithm to optimize weighting
of a total of 14 energy terms used in Rosetta. We reverse the Monte Carlo
process of Rosetta to approach native-like conformations to a process
generating from the native state an ensemble of initial conformations most
relevant to the native state. Intuitively, an ideal weighting scheme would
result in a large "basin of attraction" of the native structure, which leads to
an objective function for the linear programming. We have examined the proposal
on several benchmark proteins, and the experimental results suggest that the
optimized weights enlarge the attraction basin of the native state and improve
the quality of the predicted native states as well. In addition, a comparison
of optimal weighting schema for proteins of different classes indicates that in
different protein classes energy terms may have different effects.


  Typically, contagion strength is modeled by a transmission rate $\lambda$,
whereby all nodes in a network are treated uniformly in a mean-field
approximation. However, local agents react differently to the same contagion
based on their local characteristics. Following our recent work [EPL
\textbf{99}, 58002 (2012)], we investigate contagion spreading models with
local dynamics on complex networks. We therefore quantify contagions by their
quality, $0 \leq \alpha \leq 1$, and follow their spreading as their
transmission condition (fitness) is evaluated by local agents. We choose
various deterministic local rules. Initial spreading with exponential
quality-dependent time scales is followed by a stationary state with a
prevalence depending on the quality of the contagion. We also observe various
interesting phenomena, for example, high prevalence without the participation
of the hubs. This is in sharp contrast with the general belief that hubs play a
central role in a typical spreading process. We further study the role of
network topology in various models and find that as long as small-world effect
exists, the underlying topology does not contribute to the final stationary
state but only affects the initial spreading velocity.


  X-ray crystallography is the predominant method for obtaining atomic-scale
information about biological macromolecules. Despite the success of the
technique, obtaining well diffracting crystals still critically limits going
from protein to structure. In practice, the crystallization process proceeds
through knowledge-informed empiricism. Better physico-chemical understanding
remains elusive because of the large number of variables involved, hence little
guidance is available to systematically identify solution conditions that
promote crystallization. To help determine relationships between macromolecular
properties and their crystallization propensity, we have trained statistical
models on samples for 182 proteins supplied by the Northeast Structural
Genomics consortium. Gaussian processes, which capture trends beyond the reach
of linear statistical models, distinguish between two main physico-chemical
mechanisms driving crystallization. One is characterized by low levels of side
chain entropy and has been extensively reported in the literature. The other
identifies specific electrostatic interactions not previously described in the
crystallization context. Because evidence for two distinct mechanisms can be
gleaned both from crystal contacts and from solution conditions leading to
successful crystallization, the model offers future avenues for optimizing
crystallization screens based on partial structural information. The
availability of crystallization data coupled with structural outcomes analyzed
through state-of-the-art statistical models may thus guide macromolecular
crystallization toward a more rational basis.


  Determining the full complement of protein-coding genes is a key goal of
genome annotation. The most powerful approach for confirming protein coding
potential is the detection of cellular protein expression through peptide mass
spectrometry experiments. Here we map the peptides detected in 7 large-scale
proteomics studies to almost 60% of the protein coding genes in the GENCODE
annotation the human genome. We find that conservation across vertebrate
species and the age of the gene family are key indicators of whether a peptide
will be detected in proteomics experiments. We find peptides for most highly
conserved genes and for practically all genes that evolved before bilateria. At
the same time there is almost no evidence of protein expression for genes that
have appeared since primates, or for genes that do not have any protein-like
features or cross-species conservation. We identify 19 non-protein-like
features such as weak conservation, no protein features or ambiguous
annotations in major databases that are indicators of low peptide detection
rates. We use these features to describe a set of 2,001 genes that are
potentially non-coding, and show that many of these genes behave more like
non-coding genes than protein-coding genes. We detect peptides for just 3% of
these genes. We suggest that many of these 2,001 genes do not code for proteins
under normal circumstances and that they should not be included in the human
protein coding gene catalogue. These potential non-coding genes will be revised
as part of the ongoing human genome annotation effort.


  Autophagy is a conserved biological stress response in mammalian cells that
is responsible for clearing damaged proteins and organelles from the cytoplasm
and recycling their contents via the lysosomal pathway. In cases of mild
stress, autophagy acts as a survival mechanism, while in cases of severe stress
cells may switch to programmed cell death. Understanding the decision process
that moves a cell from autophagy to apoptosis is important since abnormal
regulation of autophagy occurs in many diseases, including cancer. To integrate
existing knowledge about this decision process into a rigorous, analytical
framework, we built a mathematical model of cell fate decisions mediated by
autophagy. Our dynamical model is consistent with existing quantitative
measurements of autophagy and apoptosis in rat kidney proximal tubular cells
responding to cisplatin-induced stress.


  The problem of equivalence in terms of multistability properties between gene
regulatory network models of different dimensionality has been recently
addressed by Schittler et al. (2013). The authors in that work proposed
construction rules for a high-dimensional dynamical system, when given a
low-dimensional dynamical system and the high-dimensional network structure.
However, the proof therein was restricted to the class of systems for which all
internal feedback loops can be broken by a loopbreaking approach yielding a
single-input single-output (SISO) system. In this report, we present the
generalization of the proof to systems with any number of internal feedback
loops, which will be broken by a generalized loopbreaking approach resulting in
a multi-input multi-output (MIMO) system. This generalization of the method
renders the construction method applicable to a broad class of gene regulatory
network models, thus promoting the transfer of results from core motif models
to more realistic, high-dimensional models of gene regulation. We demonstrate
the potential and value of our method by applying it to an example of a gene
regulatory network in mesenchymal stem cell differentiation.


  Site-specific transcription factors (TFs) bind to their target sites on the
DNA, where they regulate the rate at which genes are transcribed. Bacterial TFs
undergo facilitated diffusion (a combination of 3D diffusion around and 1D
random walk on the DNA) when searching for their target sites. Using computer
simulations of this search process, we show that the organisation of the
binding sites, in conjunction with TF copy number and binding site affinity,
plays an important role in determining not only the steady state of promoter
occupancy, but also the order at which TFs bind. These effects can be captured
by facilitated diffusion-based models, but not by standard thermodynamics. We
show that the spacing of binding sites encodes complex logic, which can be
derived from combinations of three basic building blocks: switches, barriers
and clusters, whose response alone and in higher orders of organisation we
characterise in detail. Effective promoter organizations are commonly found in
the E. coli genome and are highly conserved between strains. This will allow
studies of gene regulation at a previously unprecedented level of detail, where
our framework can create testable hypothesis of promoter logic.


  This review is a tutorial for scientists interested in the problem of protein
structure prediction, particularly those interested in using coarse-grained
molecular dynamics models that are optimized using lessons learned from the
energy landscape theory of protein folding. We also present a review of the
results of the AMH/AMC/AMW/AWSEM family of coarse-grained molecular dynamics
protein folding models to illustrate the points covered in the first part of
the article. Accurate coarse-grained structure prediction models can be used to
investigate a wide range of conceptual and mechanistic issues outside of
protein structure prediction; specifically, the paper concludes by reviewing
how AWSEM has in recent years been able to elucidate questions related to the
unusual kinetic behavior of artificially designed proteins, multidomain protein
misfolding, and the initial stages of protein aggregation.


  In this survey paper, we will present a number of core algorithmic questions
concerning several transitive reduction problems on network that have
applications in network synthesis and analysis involving cellular processes.
Our starting point will be the so-called minimum equivalent digraph problem, a
classic computational problem in combinatorial algorithms. We will subsequently
consider a few non-trivial extensions or generalizations of this problem
motivated by applications in systems biology. We will then discuss the
applications of these algorithmic methodologies in the context of three major
biological research questions: synthesizing and simplifying signal transduction
networks, analyzing disease networks, and measuring redundancy of biological
networks.


  Associative learning allows animals to rapidly adapt to changes in the
environment. Whether and what aspects of such acquired traits may be
transmittable across generations remains unclear. Using prolonged olfactory
training and subsequent two-forced choice testing in Drosophila melanogaster,
it is observed that certain aspects of learned behavior were transmitted from
parents to offspring. Offspring of parents exposed to distinct odors during
both aversive and appetitive conditioning displayed a heightened sensitivity to
those same odors. The conditioned responses associated with those odors,
however, were not transmitted to the offspring as they displayed a constitutive
preference to the parent-exposed stimuli irrespective of whether they were
associated with aversive or appetitive training. Moreover, the degree to which
the offspring preferred the conditioned stimuli markedly varied from
odor-to-odor. These findings suggest that heightened sensitivities to certain
salient stimuli in the environment, but not their associated conditioned
behaviors, may be transmittable from parents to offspring. Such
trans-generational adaptations may influence animal traits over short
evolutionary time-scales.


  Metastasis is the process by which cells from a primary tumor disperse and
form new tumors at distant anatomical locations. The treatment and prevention
of metastatic cancer remains an extremely challenging problem. This work
introduces a novel biologically motivated objective function to the radiation
optimization community that takes into account metastatic risk instead of the
status of the primary tumor. In this work, we consider the problem of
developing fractionated irradiation schedules that minimize production of
metastatic cancer cells while keeping normal tissue damage below an acceptable
level. A dynamic programming framework is utilized to determine the optimal
fractionation scheme. We evaluated our approach on a breast cancer case using
the heart and the lung as organs-at-risk (OAR). For small tumor $\alpha/\beta$
values, hypo-fractionated schedules were optimal, which is consistent with
standard models. However, for relatively larger $\alpha/\beta$ values, we found
the type of schedule depended on various parameters such as the time when
metastatic risk was evaluated, the $\alpha/\beta$ values of the OARs, and the
normal tissue sparing factors. Interestingly, in contrast to standard models,
hypo-fractionated and semi-hypo-fractionated schedules (large initial doses
with doses tapering off with time) were suggested even with large tumor
$\alpha$/$\beta$ values. Numerical results indicate potential for significant
reduction in metastatic risk.


  Interactions between proteins are hard to decipher. Protein-protein
interactions are difficult problem to address because they are not based on
differences in charge type like protein-DNA or protein-lipid interactions. In
this manuscript we present a few methods aimed at recognizing contact and their
neighboring first layer residues. The methods are based on the simplified
analysis of normal modes via the Gaussian Network Model (GNM). The methods
adjust the number of modes used in analysis on the basis of the comparison of
predicted and expected number of targets (contact and first layer residues).
The number of modes used in weighted sum calculation is increased or decreased
if the number of predictions falls out of range for a given protein sequence
length. The methods are able to recognize more than 50% of targets on average.
At the same time they generate slightly more than 40% of false positives. They
are also fairly successful in recognizing near native decoys of the Vakser
decoy sets. The results depicted here show that kinetically active residues are
important in protein-protein interactions. The GNM analysis depicted also shows
that the two interacting chains usually exhibit opposite behavior. While the
binding surface of one of the partners is rigid and stable, the interfacial
area of the other partner is more flexible.


  Cadherin adhesion receptors are fundamental determinants of tissue
organization in health and disease. Increasingly, we have come to appreciate
that classical cadherins exert their biological actions through active
cooperation with the contractile actin cytoskeleton. Rather than being passive
resistors of detachment forces, cadherins can regulate the assembly and
mechanics of the contractile apparatus itself. Moreover, coordinate spatial
patterning of adhesion and contractility is emerging as a determinant of
morphogenesis. Here we review recent developments in cadherins and actin
cytoskeleton cooperativity, by focusing on E-cadherin adhesive patterning in
the epithelia. Next, we discuss the underlying principles of cellular
rearrangement during Drosophila germband extension and epithelial cell
extrusion, as models of how planar and apical-lateral patterns of contractility
organizes tissue architecture.


  We discuss the approach to investigation of molecular machines using systems
of integro--differential ultrametric (p-adic) reaction--diffusion equations
with drift. This approach combines the features of continuous and discrete
dynamic models. We apply this model to investigation of actomyosin molecular
motor.
  The introduced system of equations is solved analytically using p-adic
wavelet theory. We find explicit stationary solutions and behavior in the
relaxation regime.


  Bacteria are known to exchange genetic information by horizontal gene
transfer. Since the frequency of homologous recombination depends on the
similarity of recombining segments, several studies examined whether this could
lead to the emergence of subspecies. Most of them used simulations of fixed
size Wright-Fisher populations, in which the genetic drift should be taken into
account. Here, we use non-linear Markov chains to describe a bacterial
population evolving under mutation and recombination processes. We consider a
population as a space of genomes and a probability measure on it. Thus, the
genetic drift is not assumed. We prove that under these conditions the
emergence of subspecies is impossible, so the genetic drift is a necessary
driving force of bacterial speciation.


  Motivation: Analysis of relationships of drug structure to biological
response is key to understanding off-target and unexpected drug effects, and
for developing hypotheses on how to tailor drug thera-pies. New methods are
required for integrated analyses of a large number of chemical features of
drugs against the corresponding genome-wide responses of multiple cell models.
Results: In this paper, we present the first comprehensive multi-set analysis
on how the chemical structure of drugs impacts on ge-nome-wide gene expression
across several cancer cell lines (CMap database). The task is formulated as
searching for drug response components across multiple cancers to reveal shared
effects of drugs and the chemical features that may be responsible. The
com-ponents can be computed with an extension of a very recent ap-proach called
Group Factor Analysis (GFA). We identify 11 compo-nents that link the
structural descriptors of drugs with specific gene expression responses
observed in the three cell lines, and identify structural groups that may be
responsible for the responses. Our method quantitatively outperforms the
limited earlier studies on CMap and identifies both the previously reported
associations and several interesting novel findings, by taking into account
multiple cell lines and advanced 3D structural descriptors. The novel
observations include: previously unknown similarities in the effects induced by
15-delta prostaglandin J2 and HSP90 inhibitors, which are linked to the 3D
descriptors of the drugs; and the induction by simvastatin of leukemia-specific
anti-inflammatory response, resem-bling the effects of corticosteroids.


  A range of a priori hypotheses about the evolution of modern and archaic
genomes are further evaluated and tested. In addition to the well-known
splits/introgressions involving Neanderthal genes into out-of- Africa people,
or Denisovan genes into Oceanians, a further series of archaic splits and
hypotheses proposed in Waddell et al. (2011) are considered in detail. These
include signals of Denisovans with something markedly more archaic and possibly
something more archaic into Papuans as well. These are compared and contrasted
with some well-advertised introgressions such as Denisovan genes across East
Asia, archaic genes into San or non-tree mixing between Oceanians, East Asians
and Europeans. The general result is that these less appreciated and surprising
archaic splits have just as much or more support in genome sequence data.
Further, evaluation confirms the hypothesis that archaic genes are much rarer
on modern X chromosomes, and may even be near totally absent, suggesting strong
selection against their introgression. Modeling of relative split weights
allows an inference of the proportion of the genome the Denisovan seems to have
gotten from an older archaic, and the best estimate is around 2%. Using a mix
of quantitative and qualitative morphological data and novel phylogenetic
methods, robust support is found for multiple distinct middle Pleistocene
lineages. Of these, fossil hominids such as SH5, Petralona, and Dali, in
particular, look like prime candidates for contributing pre-Neanderthal/Modern
archaic genes to Denisovans, while the Jinniu-Shan fossil looks like the best
candidate for a close relative of the Denisovan. That the Papuans might have
received some truly archaic genes appears a good possibility and they might
even be from Homo erectus.


  This article focuses on the development of the method for the genetic
classification of agamospermous reproduction types in plants using sugar beet
as an example. The classification feasibility is ensured by the use of isozymes
as genetic markers allowing the identification of all three phenotypic classes
in the progeny of individual heterozygous diploid plant and is based on
different phenotypic class ratios in the progenies obtained by meiotic and
mitotic agamospermy. The data indicate that for sugar beet meiotic agamospermy
is the more typical since 13 of 15 explored progenies were classified as those
produced by meiotic agamospermy and only 2 as produced by mitotic agamospermy.


  Power-law scalings are ubiquitous to physical phenomena undergoing a
continuous phase transition. The classic Susceptible-Infectious-Recovered (SIR)
model of epidemics is one such example where the scaling behavior near a
critical point has been studied extensively. In this system the distribution of
outbreak sizes scales as $P(n) \sim n^{-3/2}$ at the critical point as the
system size $N$ becomes infinite. The finite-size scaling laws for the outbreak
size and duration are also well understood and characterized. In this work, we
report scaling laws for a model with SIR structure coupled with a constant
force of infection per susceptible, akin to a `reservoir forcing'. We find that
the statistics of outbreaks in this system are fundamentally different than
those in a simple SIR model. Instead of fixed exponents, all scaling laws
exhibit tunable exponents parameterized by the dimensionless rate of external
forcing. As the external driving rate approaches a critical value, the scale of
the average outbreak size converges to that of the maximal size, and above the
critical point, the scaling laws bifurcate into two regimes. Whereas a simple
SIR process can only exhibit outbreaks of size $\mathcal{O}(N^{1/3})$ and
$\mathcal{O}(N)$ depending on whether the system is at or above the epidemic
threshold, a driven SIR process can exhibit a richer spectrum of outbreak sizes
that scale as $O(N^{\xi})$ where $\xi \in (0,1] \backslash \{2/3\}$ and
$\mathcal{O}((N/\log N)^{2/3})$ at the multi-critical point.


  The process of desertification is usually modeled as a first order
transition, where a change of an external parameter (e.g. precipitation) leads
to a catastrophic bifurcation followed by an ecological regime shift. However,
vegetation elements like shrubs and trees undergo a stochastic birth-death
process with an absorbing state; such a process supports a second order
continuous transition with no hysteresis. We present a numerical study of a
minimal model that supports bistability and catastrophic shift on spatial
domain with demographic noise and an absorbing state. When the external
parameter varies adiabatically the transition is continuous and the front
velocity renormalizes to zero at the extinction transition. Below the
transition one may identify three modes of desertification: accumulation of
local catastrophes, desert invasion and global collapse. A catastrophic regime
shift occurs as a dynamical hysteresis, when the pace of environmental
variations is too fast. We present some empirical evidence, suggesting that the
mid-holocene desertification of the Sahara was, indeed, continuous.


  We analyze the problem of the helix-coil transition in explicit solvents
analytically by using spin-based models incorporating two different mechanisms
of solvent action: explicit solvent action through the formation of
solvent-polymer hydrogen bonds that can compete with the intrinsic
intra-polymer hydrogen bonded configurations (competing interactions) and
implicit solvent action, where the solvent-polymer interactions tune biopolymer
configurations by changing the activity of the solvent (non-competing
interactions). The overall spin Hamiltonian is comprised of three terms: the
background \emph{in vacuo} Hamiltonian of the "Generalized Model of Polypeptide
Chain" type and two additive terms that account for the two above mechanisms of
solvent action. We show that on this level the solvent degrees of freedom can
be {\sl explicitly} and {\sl exactly} traced over, the ensuing effective
partition function combining all the solvent effects in a unified framework. In
this way we are able to address helix-coil transitions for polypeptides,
proteins, and DNA, with different buffers and different external constraints.
Our spin-based effective Hamiltonian is applicable for treatment of such
diverse phenomena as cold denaturation, effects of osmotic pressure on the cold
and warm denaturation, complicated temperature dependence of the hydrophobic
effect as well as providing a conceptual base for understanding the behavior of
Intrinsically Disordered Proteins and their analogues.


  We review recent results obtained from simple individual-based models of
biological competition in which birth and death rates of an organism depend on
the presence of other competing organisms close to it. In addition the
individuals perform random walks of different types (Gaussian diffusion and
L\'{e}vy flights). We focus on how competition and random motions affect each
other, from which spatial instabilities and extinctions arise. Under suitable
conditions, competitive interactions lead to clustering of individuals and
periodic pattern formation. Random motion has a homogenizing effect and then
delays this clustering instability. When individuals from species differing in
their random walk characteristics are allowed to compete together, the ones
with a tendency to form narrower clusters get a competitive advantage over the
others. Mean-field deterministic equations are analyzed and compared with the
outcome of the individual-based simulations.


  We consider the self-assembly of composite structures from a group of
nanocomponents, each consisting of particles within an $N$-atom system.
Self-assembly pathways and rates for nanocomposites are derived via a
multiscale analysis of the classical Liouville equation. From a reduced
statistical framework, rigorous stochastic equations for population levels of
beginning, intermediate, and final aggregates are also derived. It is shown
that the definition of an assembly type is a self-consistency criterion that
must strike a balance between precision and the need for population levels to
be slowly varying relative to the time scale of atomic motion. The deductive
multiscale approach is complemented by a qualitative notion of multicomponent
association and the ensemble of exact atomic-level configurations consistent
with them. In processes such as viral self-assembly from proteins and RNA or
DNA, there are many possible intermediates, so that it is usually difficult to
predict the most efficient assembly pathway. However, in the current study,
rates of assembly of each possible intermediate can be predicted. This avoids
the need, as in a phenomenological approach, for recalibration with each new
application. The method accounts for the feedback across scales in space and
time that is fundamental to nanosystem self-assembly. The theory has
applications to bionanostructures, geomaterials, engineered composites, and
nanocapsule therapeutic delivery systems.


  The effects of ionic ordering in DNA water solutions are studied by
conductivity experiments. The conductivity measurements are performed for the
solutions of DNA with KCl salt in the temperature range from 28 to 70 C. Salt
concentration vary from 0 to 2 M. The conductivity of solutions without DNA but
with the same concentration of KCl salt are also performed. The results show
that in case of salt free solution of DNA the melting process of the double
helix is observed, while in case of DNA solution with added salt the
macromolecule denaturation is not featured. For salt concentrations lower than
some critical one (0.4 M) the conductivity of DNA solution is higher than the
conductivity of KCl water solution without DNA. Starting from the critical
concentration the conductivity of KCl solution is higher than the conductivity
of DNA solution with added salt. For description of the experimental data
phenomenological model is elaborated basing on electrolyte theory. In framework
of the developed model a mechanism of counterion ordering is introduced.
According to this mechanism under the low salt concentrations electrical
conductivity of the system is caused by counterions of DNA ion-hydrate shell.
Increasing the amount of salt to the critical concentration counterions
condense on DNA polyanion. Further increase of salt concentration leads to the
formation of DNA-salt complexes that decreases the conductivity of the system.


  Summary: Longer sequencing reads, with at least 200 bases per template are
now common. While traditional aligners have adopted new strategies to improve
the mapping of longer reads, aligners specific to bisulfite-sequencing were
optimized when much shorter reads were the norm. We sought to perform the first
comparison using longer reads to determine which aligners were most accurate
and efficient and to evaluate a novel software tool, bwa-meth, built on a
traditional mapper that supports insertions, deletions and clipped alignments.
We gauge accuracy by comparing the number of on and off-target reads from a
targeted sequencing project and by simulations. Availability and
Implementation: The benchmarking scripts and the bwa-meth software are
available at https://github/com/brentp/bwa-meth/ under the MIT License.


  We present a novel system of equations to describe the evolution of
self-organized structured societies (biological or human) composed of several
trait groups. The suggested approach is based on the combination of ideas
employed in the theory of biological populations, system theory, and utility
theory. The evolution equations are defined as utility rate equations, whose
parameters are characterized by the utility of each group with respect to the
society as a whole and by the mutual utilities of groups with respect to each
other. We analyze in detail the cases of two groups (cooperators and defectors)
and of three groups (cooperators, defectors, and regulators) and find that, in
a self-organized society, neither defectors nor regulators can overpass the
maximal fractions of about 10% each. This is in agreement with the data for bee
and ant colonies. The classification of societies by their distance from
equilibrium is proposed. We apply the formalism to rank the countries according
to the introduced metric quantifying their relative stability, which depends on
the cost of defectors and regulators as well as their respective population
fractions. We find a remarkable concordance with more standard economic ranking
based, for instance, on GDP per capita.


  Epigenetic histone modifications play an important role in the maintenance of
different cell phenotypes. The exact molecular mechanism for inheritance of the
modification patterns over cell generations remains elusive. We construct a
Potts-type model based on experimentally observed nearest-neighbor enzyme
lateral interactions and nucleosome covalent modification state biased enzyme
recruitment. The model can lead to effective nonlocal interactions among
nucleosomes suggested in previous theoretical studies, and epigenetic memory is
robustly inheritable against stochastic cellular processes.


  Statistical methods for analyzing large-scale biomolecular data are
commonplace in computational biology. A notable example is phenotype prediction
from gene expression data, for instance, detecting human cancers,
differentiating subtypes and predicting clinical outcomes. Still, clinical
applications remain scarce. One reason is that the complexity of the decision
rules that emerge from standard statistical learning impedes biological
understanding, in particular, any mechanistic interpretation. Here we explore
decision rules for binary classification utilizing only the ordering of
expression among several genes; the basic building blocks are then two-gene
expression comparisons. The simplest example, just one comparison, is the TSP
classifier, which has appeared in a variety of cancer-related discovery
studies. Decision rules based on multiple comparisons can better accommodate
class heterogeneity, and thereby increase accuracy, and might provide a link
with biological mechanism. We consider a general framework ("rank-in-context")
for designing discriminant functions, including a data-driven selection of the
number and identity of the genes in the support ("context"). We then specialize
to two examples: voting among several pairs and comparing the median expression
in two groups of genes. Comprehensive experiments assess accuracy relative to
other, more complex, methods, and reinforce earlier observations that simple
classifiers are competitive.


  The bioactivity of nanoparticles crucially depends on their ability to cross
biomembranes. Recent simulations indicate the cooperative wrapping and
internalization of spherical nanoparticles in tubular membrane structures. In
this article, we systematically investigate the energy gain of this cooperative
wrapping by minimizing the energies of the rotationally symmetric shapes of the
membrane tubes and of membrane segments wrapping single particles. We find that
the energy gain for the cooperative wrapping of nanoparticles in membrane tubes
relative to their individual wrapping as single particles strongly depends on
the ratio of the particle radius and the range of the particle-membrane
adhesion potential. For a potential range of the order of one nanometer, the
cooperative wrapping in tubes is highly favorable for particles with a radius
of tens of nanometers and intermediate adhesion energies, but not for particles
that are significantly larger.


  Whole genome sequencing of pathogens from multiple hosts in an epidemic
offers the potential to investigate who infected whom with unparalleled
resolution, potentially yielding important insights into disease dynamics and
the impact of control measures. We considered disease outbreaks in a setting
with dense genomic sampling, and formulated stochastic epidemic models to
investigate person-to-person transmission, based on observed genomic and
epidemiological data. We constructed models in which the genetic distance
between sampled genotypes depends on the epidemiological relationship between
the hosts. A data augmented Markov chain Monte Carlo algorithm was used to
sample over the transmission trees, providing a posterior probability for any
given transmission route. We investigated the predictive performance of our
methodology using simulated data, demonstrating high sensitivity and
specificity, particularly for rapidly mutating pathogens with low
transmissibility. We then analyzed data collected during an outbreak of
methicillin-resistant Staphylococcus aureus in a hospital, identifying probable
transmission routes and estimating epidemiological parameters. Our approach
overcomes limitations of previous methods, providing a framework with the
flexibility to allow for unobserved infection times, multiple independent
introductions of the pathogen, and within-host genetic diversity, as well as
allowing forward simulation.


  Standard models of population dynamics focus on the the interaction,
survival, and extinction of the competing species individually. Real ecological
systems, however, are characterized by an abundance of species (or strategies,
in the terminology of evolutionary-game theory) that form intricate, complex
interaction networks. The description of the ensuing dynamics may be aided by
studying associations of certain strategies rather than individual ones. Here
we show how such a higher-level description can bear fruitful insight.
Motivated from different strains of colicinogenic Escherichia coli bacteria, we
investigate a four-strategy system which contains a three-strategy cycle and a
neutral alliance of two strategies. We find that the stochastic, spatial model
exhibits a mobility-dependent selection of either the three-strategy cycle or
of the neutral pair. We analyze this intriguing phenomenon numerically and
analytically.


  To illustrate the power of the biophysical approach in solving important
problems in life science, I present here one of our current research projects
as an example. We have developed special biophotonic techniques to study the
dynamic properties of signaling proteins in a single living cell. Such a study
allowed us to gain new insight into the signaling mechanism that regulates
programmed cell death.


  Summary: Genotype networks are a method used in systems biology to study the
innovability of a given phenotype, determining whether the phenotype is robust
to mutations, and how do the genotypes associated to it are distributed in the
genotype space. Here we developed VCF2Networks, a tool to apply this method to
population genetics data, and in particular to single Nucleotide Variants data
encoded in the Variant Call file Format (VCF). A complete summary of the
properties of the genotype network that can be calculated by VCF2Networks is
given in the Supplementary Materials 1.
  Availability and Implementation: The home page of the project is
https://bitbucket.org/dalloliogm/vcf2networks . VCF2Networks is also available
directly from the Python Package Index (PyPI), under the name vcf2networks.


  Observation of phenotypic diversity in a population of genetically identical
cells is often linked to the stochastic nature of chemical reactions involved
in gene regulatory networks. We investigate the distribution of population
averaged gene expression levels as a function of population, or sample, size
for several stochastic gene expression models to find out to what extent
population averaged quantities reflect the underlying mechanism of gene
expression. We consider three basic gene regulation networks corresponding to
transcription with and without gene state switching and translation. Using
analytical expressions for the probability generating function of observables
and Large Deviation Theory, we calculate the distribution and first two moments
of the population averaged mRNA and protein levels as a function of model
parameters, population size and number of measurements contained in a data set.
We validate our results using stochastic simulations also report exact results
on the asymptotic properties of population averages which show qualitative
differences among different models.


  Genome-wide analysis of distributions of densities of long-range interactions
of human chromosomes with each other, nucleoli, nuclear lamina, and binding
sites of chromatin state regulatory proteins, CTCF and STAT1, identifies
non-random highly correlated patterns of density distributions along the
chromosome length for all these features. Marked co-enrichments and clustering
of all these interactions are detected at discrete genomic regions on selected
chromosomes, which are located within pericentromeric heterochromatin and
designated Centromeric Regions of Interphase Chromatin Homing (CENTRICH).
CENTRICH manifest 199-716-fold higher density of inter-chromosomal binding
sites compared to genome-wide or chromosomal averages (p =
2.10E-101-1.08E-292). Sequence alignment analysis shows that CENTRICH represent
unique DNA sequences of 3.9 to 22.4 Kb in size which are: 1) associated with
nucleolus; 2) exhibit highly diverse set of DNA-bound chromatin state
regulators, including marked enrichment of CTCF and STAT1 binding sites; 3)
bind multiple intergenic disease-associated genomic loci (IDAGL) with
documented long-range enhancer activities and established links to increased
risk of developing epithelial malignancies and other common human disorders.
Using distances of SNP loci homing sites within genomic coordinates of CENTRICH
as a proxy of likelihood of disease-linked SNP loci binding to CENTRICH, we
demonstrate statistically significant correlations between the probability of
SNP loci binding to CENTRICH and GWAS-defined odds ratios of increased risk of
a disease for cancer, coronary artery disease, and type 2 diabetes. Our
analysis suggests that centromeric sequences and pericentromeric
heterochromatin may play an important role in human cells beyond the critical
functions in chromosome segregation.


  Wolbachia are common endosymbionts of terrestrial arthropods, and are also
found in nematodes: the animal-parasitic filaria, and the plant-parasite
Radopholus similis. Lateral transfer of Wolbachia DNA to the host genome is
common. We generated a draft genome sequence for the strongyloidean nematode
parasite Dictyocaulus viviparus, the cattle lungworm. In the assembly, we
identified nearly 1 Mb of sequence with similarity to Wolbachia. The fragments
were unlikely to derive from a live Wolbachia infection: most were short, and
the genes were disabled through inactivating mutations. Many fragments were
co-assembled with definitively nematode-derived sequence. We found limited
evidence of expression of the Wolbachia-derived genes. The D. viviparus
Wolbachia genes were most similar to filarial strains, and strains from the
host-promiscuous clade F. We conclude that D. viviparus was infected by
Wolbachia in the past. Genome sequence based surveys are a powerful tool for
revealing the genome archaeology of infection and symbiosis.


  Weight gain is a well-established side effect of both conventional and newer
anti-psychotic drugs, but the cause is not well understood. Recent studies
correlate obesity with the presence or absence of particular genetic sequences
in the gut microbiota. We identified strong associations between protein
targets of antipsychotics and microbiota sequences directly related to weight
regulation in human body, leading to a potential metagenomic mechanism of
action. Further experimental study is recommended.


  English version of abstract
  A growing number of scientific researches have been demonstrating that olive
oil operates a crucial role on the prevention of cardiovascular and tumoral
diseases, being related with low mortality and morbidity in populations that
tend to follow a Mediterranean diet. Amongst its minor components, polyphenols
have been subject to several clinical trials that have established health
benefits associated to their antioxidant, antitumoral and anti-atherosclerotic
activity. However the biological activity of polyphenols is dependent not only
of their absorption but also of their metabolization. Bioavailability studies
have demonstrated that after olive oil intake, not only hydroxytyrosol but also
its metabolites can be found in plasma, such as
3,4-dihydroxyphenylacetaldehyde, 3,4-dihydroxyphenylacetic acid, homovanillyl
alcohol and its glucoronides.
  -----
  Portuguese version of abstract
  Um numero crescente de pesquisas cientificas demonstram que o azeite opera um
papel crucial na preven\c{c}ao de doen\c{c}as cardiovasculares e doen\c{c}as
tumorais, estando relacionado com a baixa mortalidade e morbilidade em
popula\c{c}oes que tradicionalmente seguem uma dieta Mediterranica. De entre os
seus componentes minorit\'arios, os polifenois t\^em vindo a ser alvo de
estudos clinicos que demonstraram o seu benef\'icio para a saude pela atividade
antioxidante, anti-tumoral e anti-aterosclerotica que possuem. No entanto a
atividade biol\'ogica destes polifenois \'e dependente nao s\'o da sua
absor\c{c}ao como poder\'a tambem depender da sua metaboliza\c{c}ao. Estudos de
biodisponibilidade demonstraram que ap\'os o consumo de azeite encontram-se no
plasma nao s\'o hidroxitirosol, mas igualmente os seus metabolitos, tais como o
\'acido 3,4-dihidroxifenilacetico, o 3,4-dihidroxifenilacetaldeido, o \'alcool
homovanilico e os respectivos glucoronideos.


  Delay is an important and ubiquitous aspect of many biochemical processes.
For example, delay plays a central role in the dynamics of genetic regulatory
networks as it stems from the sequential assembly of first mRNA and then
protein. Genetic regulatory networks are therefore frequently modeled as
stochastic birth-death processes with delay. Here we examine the relationship
between delay birth-death processes and their appropriate approximating delay
chemical Langevin equations. We prove that the distance between these two
descriptions, as measured by expectations of functionals of the processes,
converges to zero with increasing system size. Further, we prove that the delay
birth-death process converges to the thermodynamic limit as system size tends
to infinity. Our results hold for both fixed delay and distributed delay.
Simulations demonstrate that the delay chemical Langevin approximation is
accurate even at moderate system sizes. It captures dynamical features such as
the spatial and temporal distributions of transition pathways in metastable
systems, oscillatory behavior in negative feedback circuits, and
cross-correlations between nodes in a network. Overall, these results provide a
foundation for using delay stochastic differential equations to approximate the
dynamics of birth-death processes with delay.