  While the characteristics of the driven translocation for asymptotically long
polymers are well understood, this is not the case for finite-sized polymers,
which are relevant for real-world experiments and simulation studies. Most
notably, the behavior of the exponent $\alpha$, which describes the scaling of
the translocation time with polymer length, when the driving force $f_p$ in the
pore is changed, is under debate. By Langevin dynamics simulations of regular
and modified translocation models using the freely-jointed-chain polymer model
we find that a previously reported incomplete model, where the {\it trans} side
and fluctuations were excluded, gives rise to characteristics that are in stark
contradiction with those of the complete model, for which $\alpha$ increases
with $f_p$. Our results suggest that contribution due to fluctuations is
important. We construct a minimal model where dynamics is completely excluded
to show that close alignment with a full translocation model can be achieved.
Our findings set very stringent requirements for a minimal model that is
supposed to describe the driven polymer translocation correctly.


  Protein-protein interactions are fundamental to many biological processes.
Experimental screens have identified tens of thousands of interactions and
structural biology has provided detailed functional insight for select 3D
protein complexes. An alternative rich source of information about protein
interactions is the evolutionary sequence record. Building on earlier work, we
show that analysis of correlated evolutionary sequence changes across proteins
identifies residues that are close in space with sufficient accuracy to
determine the three-dimensional structure of the protein complexes. We evaluate
prediction performance in blinded tests on 76 complexes of known 3D structure,
predict protein-protein contacts in 32 complexes of unknown structure, and
demonstrate how evolutionary couplings can be used to distinguish between
interacting and non-interacting protein pairs in a large complex. With the
current growth of sequence databases, we expect that the method can be
generalized to genome-wide elucidation of protein-protein interaction networks
and used for interaction predictions at residue resolution.


  The different cell types in a living organism acquire their identity through
the process of cell differentiation in which the multipotent progenitor cells
differentiate into distinct cell types. Experimental evidence and analysis of
large-scale microarray data establish the key role played by a two-gene motif
in cell differentiation in a number of cell systems. The two genes express
transcription factors which repress each other's expression and autoactivate
their own production. A number of theoretical models have recently been
proposed based on the two-gene motif to provide a physical understanding of how
cell differentiation occurs. In this paper, we study a simple model of cell
differentiation which assumes no cooperativity in the regulation of gene
expression by the transcription factors. The latter repress each other's
activity directly through DNA binding and indirectly through the formation of
heterodimers. We specifically investigate how deterministic processes combined
with stochasticity contribute in bringing about cell differentiation. The
deterministic dynamics of our model give rise to a supercritical pitchfork
bifurcation from an undifferentiated stable steady state to two differentiated
stable steady states. The stochastic dynamics of our model are studied using
the approaches based on the Langevin equations and the linear noise
approximation. The simulation results provide a new physical understanding of
recent experimental observations. We further propose experimental measurements
of quantities like the variance and the lag-1 autocorrelation function in
protein fluctuations as the early signatures of an approaching bifurcation
point in the cell differentiation process.


  The fitness of somatic cells of metazoan, the ability of proliferation and
survival, depends on microenvironment. In somatic evolution, a mutated cell in
a tissue clonally expands abnormally because of its high fitness as normal
cells in a corresponding microenvironment. In this study, we propose the cancer
cell hypothesis that cancer cells are the mutated cells with two
characteristics: clonal expansion and damaging the microenvironment through the
behaviours such as producing more poison in metabolism than normal cells. This
model provides an explanation for the nature of invasion and metastasis, which
are still controversial. In addition, we theoretically reasoned out that normal
cells have almost the highest fitness in healthy microenvironments as a result
of long-term organic evolution. This inspires a new kind of therapy of cancer,
which improving microenvironment to make cancer cells lower in fitness than
normal cells and then halt the growth of tumours. This general therapy relies
on a mechanism differing from chemotherapy and targeted therapy.


  The conduction and selectivity of calcium/sodium ion channels are described
in terms of ionic Coulomb blockade, a phenomenon based on charge discreteness
and an electrostatic model of an ion channel. This novel approach provides a
unified explanation of numerous observed and modelled conductance and
selectivity phenomena, including the anomalous mole fraction effect and
discrete conduction bands. Ionic Coulomb blockade and resonant conduction are
similar to electronic Coulomb blockade and resonant tunnelling in quantum dots.
The model is equally applicable to other nanopores.


  We present a machine learning framework for modeling protein dynamics. Our
approach uses L1-regularized, reversible hidden Markov models to understand
large protein datasets generated via molecular dynamics simulations. Our model
is motivated by three design principles: (1) the requirement of massive
scalability; (2) the need to adhere to relevant physical law; and (3) the
necessity of providing accessible interpretations, critical for both cellular
biology and rational drug design. We present an EM algorithm for learning and
introduce a model selection criteria based on the physical notion of
convergence in relaxation timescales. We contrast our model with standard
methods in biophysics and demonstrate improved robustness. We implement our
algorithm on GPUs and apply the method to two large protein simulation datasets
generated respectively on the NCSA Bluewaters supercomputer and the
Folding@Home distributed computing network. Our analysis identifies the
conformational dynamics of the ubiquitin protein critical to cellular
signaling, and elucidates the stepwise activation mechanism of the c-Src kinase
protein.


  Protein-protein interaction networks associated with diseases have gained
prominence as an area of research. We investigate algebraic and topological
indices for protein-protein interaction networks of 11 human cancers derived
from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We find a
strong correlation between relative automorphism group sizes and topological
network complexities on the one hand and five year survival probabilities on
the other hand. Moreover, we identify several protein families (e.g. PIK, ITG,
AKT families) that are repeated motifs in many of the cancer pathways.
Interestingly, these sources of symmetry are often central rather than
peripheral. Our results can aide in identification of promising targets for
anti-cancer drugs. Beyond that, we provide a unifying framework to study
protein-protein interaction networks of families of related diseases (e.g.
neurodegenerative diseases, viral diseases, substance abuse disorders).


  In the present paper we study a tri-trophic Lotka-Volterra food web model
with omnivory and predator switching. We observe that if the intensity of
predator switching increases the chaotic behavior of the omnivory model will be
reduced. We conclude that the predator switching mechanism enhances the
stability of an otherwise chaotic system.


  The world in which we are living is a huge network of networks and should be
described by interdependent networks. The interdependence between networks
significantly affects the evolutionary dynamics of cooperation on them.
Meanwhile, due to the diversity and complexity of social and biological
systems, players on different networks may not interact with each other by the
same way, which should be described by multiple models in evolutionary game
theory, such as the Prisoner's Dilemma and Snowdrift Game. We therefore study
the evolutionary dynamics of cooperation on two interdependent networks playing
different games respectively. We clearly evidence that, with the increment of
network interdependence, the evolution of cooperation is dramatically promoted
on the network playing Prisoner's Dilemma. The cooperation level of the network
playing Snowdrift Game reduces correspondingly, although it is almost
invisible. In particular, there exists an optimal intermediate region of
network interdependence maximizing the growth rate of the evolution of
cooperation on the network playing Prisoner's Dilemma. Remarkably, players
contacting with other network have advantage in the evolution of cooperation
than the others on the same network.


  Trait differences between species may be attributable to natural selection.
However, quantifying the strength of evidence for selection acting on a
particular trait is a difficult task. Here we develop a population-genetic test
for selection acting on a quantitative trait which is based on multiple-line
crosses. We show that using multiple lines increases both the power and the
scope of selection inference. First, a test based on three or more lines
detects selection with strongly increased statistical significance, and we show
explicitly how the sensitivity of the test depends on the number of lines.
Second, a multiple-line test allows to distinguish different lineage-specific
selection scenarios. Our analytical results are complemented by extensive
numerical simulations. We then apply the multiple-line test to QTL data on
floral character traits in plant species of the Mimulus genus and on
photoperiodic traits in different maize strains, where we find a signatures of
lineage-specific selection not seen in a two-line test.


  The distributions of many proteins in rod-shaped bacteria are far from
homogenous. Often they accumulate at the cell poles or in the cell center. At
the same time, the copy number of proteins in a single cell is relatively small
making the patterns noisy. To explore limits to protein patterns due to
molecular noise, we studied a generic mechanism for spontaneous polar protein
assemblies in rod-shaped bacteria, which is based on cooperative binding of
proteins to the cytoplasmic membrane. For mono-polar assemblies, we find that
the switching time between the two poles increases exponentially with the cell
length and with the protein number.


  Cooperation is fundamental to the evolution of human society. We regularly
observe cooperative behaviour in everyday life and in controlled experiments
with anonymous people, even though standard economic models predict that they
should deviate from the collective interest and act so as to maximise their own
individual payoff. However, there is typically heterogeneity across subjects:
some may cooperate, while others may not. Since individual factors promoting
cooperation could be used by institutions to indirectly prime cooperation, this
heterogeneity raises the important question of who these cooperators are. We
have conducted a series of experiments to study whether benevolence, defined as
a unilateral act of paying a cost to increase the welfare of someone else
beyond one's own, is related to cooperation in a subsequent one-shot anonymous
Prisoner's dilemma. Contrary to the predictions of the widely used inequity
aversion models, we find that benevolence does exist and a large majority of
people behave this way. We also find benevolence to be correlated with
cooperative behaviour. Finally, we show a causal link between benevolence and
cooperation: priming people to think positively about benevolent behaviour
makes them significantly more cooperative than priming them to think
malevolently. Thus benevolent people exist and cooperate more.


  We consider systems with two competing species whose actions are completely
symmetric, with same mobility, reproduction and competition rates. Numerical
implementations of the model in two and three-dimensional space show that
regions of single species are formed by spontaneous symmetry breaking. We
propose a theoretical formalism for describing the static profile of the
interfaces of empty spaces separating domains with different species. We
compute the topological properties of the interfaces and show that these
theoretical functions are useful to the understanding of the dynamics of the
network. Finally, we compare the theoretical functions with results from the
numerical implementation of the mean field equations and verify that our model
fits well the properties of interfaces.


  In this paper we propose a method and discuss its computational
implementation as an integrated tool for the analysis of viral genetic
diversity on data generated by high-throughput sequencing. Most methods for
viral diversity estimation proposed so far are intended to take benefit of the
longer reads produced by some NGS platforms in order to estimate a population
of haplotypes. Our goal here is to take advantage of distinct virtues of a
certain kind of NGS platform - the platform SOLiD (Life Technologies) is an
example - that has not received much attention due to the short length of its
reads, which renders haplotype estimation very difficult. However, this kind of
platform has a very low error rate and extremely deep coverage per site and our
method is designed to take advantage of these characteristics. We propose to
measure the populational genetic diversity through a family of multinomial
probability distributions indexed by the sites of the virus genome, each one
representing the populational distribution of the diversity per site. The
implementation of the method focuses on two main optimization strategies: a
read mapping/alignment procedure that aims at the recovery of the maximum
possible number of short-reads; the estimation of the multinomial parameters
through a Bayesian approach, which, unlike simple frequency counting, allows
one to take into account the prior information of the control population within
the inference of a posterior experimental condition and provides a natural way
to separate signal from noise, since it automatically furnishes Bayesian
confidence intervals. The methods described in this paper have been implemented
as an integrated tool called Tanden (Tool for Analysis of Diversity in Viral
Populations).


  In order to gain an understanding of the effectiveness of phylogenetic Markov
chain Monte Carlo (MCMC), it is important to understand how quickly the
empirical distribution of the MCMC converges to the posterior distribution. In
this paper we investigate this problem on phylogenetic tree topologies with a
metric that is especially well suited to the task: the subtree
prune-and-regraft (SPR) metric. This metric directly corresponds to the minimum
number of MCMC rearrangements required to move between trees in common
phylogenetic MCMC implementations. We develop a novel graph-based approach to
analyze tree posteriors and find that the SPR metric is much more informative
than simpler metrics that are unrelated to MCMC moves. In doing so we show
conclusively that topological peaks do occur in Bayesian phylogenetic
posteriors from real data sets as sampled with standard MCMC approaches,
investigate the efficiency of Metropolis-coupled MCMC (MCMCMC) in traversing
the valleys between peaks, and show that conditional clade distribution (CCD)
can have systematic problems when there are multiple peaks.


  As a model of soft confinement for polymers, we investigated equilibrium
shapes of a flexible vesicle that contains a phase-separating polymer solution.
To simulate such a system, we combined the phase field theory (PFT) for the
vesicle and the self-consistent field theory (SCFT) for the polymer solution.
We observed a transition from a symmetric prolate shape of the vesicle to an
asymmetric pear shape induced by the domain structure of the enclosed polymer
solution. Moreover, when a non-zero spontaneous curvature of the vesicle is
introduced, a re-entrant transition between the prolate and the dumbbell shapes
of the vesicle is observed. This re-entrant transition is explained by
considering the competition between the loss of conformational entropy and that
of translational entropy of polymer chains due to the confinement by the
deformable vesicle. This finding is in accordance with the recent experimental
result reported by Terasawa, et al.


  Mosquito-borne diseases cause significant public health burden, mostly in
tropical and sub-tropical regions, and are widely emerging or re-emerging in
areas where previously absent. Understanding, predicting, and mitigating the
spread of mosquito-borne disease in diverse populations and geographies are
ongoing modeling challenges. We propose a hybrid network-patch model for the
spread of mosquito-borne pathogens that accounts for the movement of
individuals through mosquito habitats and responds to environmental factors
such as rainfall and temperature. Our approach extends the capabilities of
existing agent-based models for individual movement developed to predict the
spread of directly transmitted pathogens in populations. To extend to
mosquito-borne disease, agent-based models are coupled with differential
equations representing `clouds' of mosquitoes in geographic patches that
account for mosquito ecology, including heterogeneity in mosquito density,
emergence rates, and extrinsic incubation period. We illustrate the method by
adapting an agent-based model for human movement across a network to
mosquito-borne disease. We investigated the importance of heterogeneity in
mosquito population dynamics and host movement on pathogen transmission,
motivating the utility of detailed models of individual behavior. We observed
that the total number of infected people is greater in heterogeneous patch
models with one high risk patch and high or medium human movement than it would
be in a random mixing homogeneous model. Our hybrid agent-based/differential
equation model is able to quantify the importance of the heterogeneity in
predicting the spread and invasion of mosquito-borne pathogens. Mitigation
strategies can be more effective when guided by realistic models created by
extending the capabilities of existing agent-based models to include
vector-borne diseases.


  Any organism is embedded in an environment that changes over time. The
timescale for and statistics of environmental change, the precision with which
the organism can detect its environment, and the costs and benefits of
particular protein expression levels all will affect the suitability of
different strategies-such as constitutive expression or graded response-for
regulating protein levels in response to environmental inputs. We propose a
general framework-here specifically applied to the enzymatic regulation of
metabolism in response to changing concentrations of a basic nutrient-to
predict the optimal regulatory strategy given the statistics of fluctuations in
the environment and measurement apparatus, respectively, and the costs
associated with enzyme production. We use this framework to address three
fundamental questions: (i) when a cell should prefer thresholding to a graded
response; (ii) when there is a fitness advantage to implementing a Bayesian
decision rule; and (iii) when retaining memory of the past provides a selective
advantage. We specifically find that: (i) relative convexity of enzyme
expression cost and benefit influences the fitness of thresholding or graded
responses; (ii) intermediate levels of measurement uncertainty call for a
sophisticated Bayesian decision rule; and (iii) in dynamic contexts,
intermediate levels of uncertainty call for retaining memory of the past.
Statistical properties of the environment, such as variability and correlation
times, set optimal biochemical parameters, such as thresholds and decay rates
in signaling pathways. Our framework provides a theoretical basis for
interpreting molecular signal processing algorithms and a classification scheme
that organizes known regulatory strategies and may help conceptualize
heretofore unknown ones.


  We develop a theory of generalist predation showing how alternative prey
species are affected by changes in both mean abundance and variability
(coefficient of variation) of their predator's primary prey. The theory is
motivated by the indirect effects of cyclic rodent populations on
ground-breeding birds, and developed through progressive analytic
simplifications of an empirically-based model. It applies nonetheless to many
other systems where primary prey have fast life-histories and can become
locally superabundant, which facilitates impact on alternative prey species. In
contrast to classic apparent competition theory based on symmetric
interactions, our results suggest that predator effects on alternative prey
should generally decrease with mean primary prey abundance, and increase with
primary prey variability (low to high CV) - unless predators have strong
aggregative responses, in which case these results can be reversed.
Approximations of models including predator dynamics (general numerical
response with possible delays) confirm these results but further suggest that
negative temporal correlation between predator and primary prey is harmful to
alternative prey. We find in general that predator numerical responses are
crucial to predict the response of ecosystems to changes in key prey species
exhibiting outbreaks, and extend the apparent competition/mutualism theory to
asymmetric interactions.


  The ~4-Mbp basic genome shared by 32 independent isolates of E. coli
representing considerable population diversity has been approximated by
whole-genome multiple-alignment and computational filtering designed to remove
mobile elements and highly variable regions. Single nucleotide polymorphisms
(SNPs) in the 496 basic-genome pairs are identified and clonally inherited
stretches are distinguished from those acquired by horizontal transfer (HT) by
sharp discontinuities in SNP density. The six least diverged genome-pairs each
have only one or two HT stretches, each occupying 42-115-kbp of basic genome
and containing at least one gene cluster known to confer selective advantage.
At higher divergences, the typical mosaic pattern of interspersed clonal and HT
stretches across the entire basic genome are observed, including likely
fragmented integrations across a restriction barrier. A simple model suggests
that individual HT events are of the order of 10-kbp and are the chief
contributor to genome divergence, bringing in almost 12 times more SNPs than
point mutations. As a result of continuing horizontal transfer of such large
segments, 400 out of the 496 strain-pairs beyond genomic divergence of share
virtually no genomic material with their common ancestor. We conclude that the
active and continuing horizontal transfer of moderately large genomic fragments
is likely to be mediated primarily by a co evolving population of phages that
distribute random genome fragments throughout the population by generalized
transduction, allowing efficient adaptation to environmental changes.


  Transcription factors (TFs) interact with a multitude of binding sites on DNA
and partner proteins inside cells. We investigate how nonspecific
binding/unbinding to such decoy binding sites affects the magnitude and
time-scale of random fluctuations in TF copy numbers arising from stochastic
gene expression. A stochastic model of TF gene expression, together with decoy
site interactions is formulated. Distributions for the total (bound and
unbound) and free (unbound) TF levels are derived by analytically solving the
chemical master equation under physiologically relevant assumptions. Our
results show that increasing the number of decoy binding sides considerably
reduces stochasticity in free TF copy numbers. The TF autocorrelation function
reveals that decoy sites can either enhance or shorten the time-scale of TF
fluctuations depending on model parameters. To understand how noise in TF
abundances propagates downstream, a TF target gene is included in the model.
Intriguingly, we find that noise in the expression of the target gene decreases
with increasing decoy sites for linear TF-target protein dose-responses, even
in regimes where decoy sites enhance TF autocorrelation times. Moreover,
counterintuitive noise transmissions arise for nonlinear dose-responses. In
summary, our study highlights the critical role of molecular sequestration by
decoy binding sites in regulating the stochastic dynamics of TFs and target
proteins at the single-cell level.


  In this study, through various tests one theoretical model is presented to
describe the relationship between promoter strength and its nucleotide
sequence. Our analysis shows that, promoter strength is greatly influenced by
nucleotide groups with three adjacent nucleotides in its sequence. Meanwhile,
nucleotides in different regions of promoter sequence have different effects on
promoter strength. Based on experimental data for {\it E. coli} promoters, our
calculations indicate, nucleotides in -10 region, -35 region, and the
discriminator region of promoter sequence are more important than those in
spacing region for determining promoter strength. With model parameter values
obtained by fitting to experimental data, four promoter libraries are
theoretically built for the corresponding experimental environments under which
data for promoter strength in gene expression has been measured previously.


  We consider a class of nonlocal reaction-diffusion problems, referred to as
replicator-mutator equations in evolutionary genetics. By using explicit
changes of unknown function, we show that they are equivalent to the heat
equation and, therefore, compute their solution explicitly. Based on this, we
then prove that, in the case of beneficial mutations in asexual populations,
solutions dramatically depend on the tails of the initial data: they can be
global, become extinct in finite time or, even, be defined for no positive
time. In the former case, we prove that solutions are accelerating, and in many
cases converge for large time to some universal Gaussian profile. This sheds
light on the biological relevance of such models.


  We consider the model of spatially distributed population consisting of two
species with "\textsl{predator\,--\,prey}" interaction; each of the species
occupies two stations. Transfer of individuals between the stations (migration)
is not random and yields the maximization of a net reproduction of each
species. Besides, each species implements reflexive behavior strategy to
determine the optimal migration flow.


  We present the exact solution of a microscopic statistical mechanical model
for the transformation of a long polypeptide between an unstructured coil
conformation and an $\alpha$-helix conformation. The polypeptide is assumed to
be adsorbed to the interface between a polar and a non-polar environment such
as realized by water and the lipid bilayer of a membrane. The interfacial
coil-helix transformation is the first stage in the folding process of helical
membrane proteins. Depending on the values of model parameters, the
conformation changes as a crossover, a discontinuous transition, or a
continuous transition with helicity in the role of order parameter. Our model
is constructed as a system of statistically interacting quasiparticles that are
activated from the helix pseudo-vacuum. The particles represent links between
adjacent residues in coil conformation that form a self-avoiding random walk in
two dimensions. Explicit results are presented for helicity, entropy, heat
capacity, and the average numbers and sizes of both coil and helix segments.


  Self-replicating systems based on information-coding polymers are of crucial
importance in biology. They also recently emerged as a paradigm in material
design on nano- and micro-scales. We present a general theoretical and
numerical analysis of the problem of spontaneous emergence of autocatalysis for
heteropolymers capable of template-assisted ligation driven by cyclic changes
in the environment. Our central result is the existence of the first order
transition between the regime dominated by free monomers and that with a
self-sustaining population of sufficiently long chains. We provide a simple,
mathematically tractable model supported by numerical simulations, which
predicts the distribution of chain lengths and the onset of autocatalysis in
terms of the overall monomer concentration and two fundamental rate constants.
Another key result of our study is the emergence of the kinetically-limited
optimal overlap length between a template and each of its two substrates. The
template-assisted ligation allows for heritable transmission of the information
encoded in chain sequences thus opening up the possibility of long-term memory
and evolvability in such systems.


  Hybrid evolution and horizontal gene transfer (HGT) are processes where
evolutionary relationships may more accurately be described by a reticulated
network than by a tree. In such a network, there will often be several paths
between any two extant species, reflecting the possible pathways that genetic
material may have been passed down from a common ancestor to these species.
These paths will typically have different lengths but an `average distance' can
still be calculated between any two taxa. In this article, we ask whether this
average distance is able to distinguish reticulate evolution from pure
tree-like evolution. We consider two types of reticulation networks:
hybridization networks and HGT networks. For the former, we establish a general
result which shows that average distances between extant taxa can appear
tree-like, but only under a single hybridization event near the root; in all
other cases, the two forms of evolution can be distinguished by average
distances. For HGT networks, we demonstrate some analogous but more intricate
results.


  Antibodies are affinity proteins with a wide spectrum of applications in
analytical and therapeutic biology. Proteins showing specific recognition for a
chosen molecular target can be isolated and their encoding sequence identified
in vitro from a large and diverse library by phage display selection. In this
work, we show that this standard biochemical technique rapidly yields a
collection of antibody protein binders for an inorganic target of major
technological importance: crystalline metallic gold surfaces. 21 distinct
anti-gold antibody proteins emerged from a large random library of antibodies
and were sequenced. The systematic statistical analysis of all the protein
sequences reveals a strong occurrence of arginine in anti-gold antibodies,
which corroborates recent molecular dynamics predictions on the crucial role of
arginine in protein/gold interactions. Once tethered to small gold
nanoparticles using histidine tag chemistry, the selected antibodies could
drive the self-assembly of the colloids onto the surface of single crystalline
gold platelets as a first step towards programmable protein-driven construction
of complex plasmonic architectures. Electrodynamic simulations based on the
Green Dyadic Method suggest that the antibody-driven assembly demonstrated here
could be exploited to significantly modify the plasmonic modal properties of
the gold platelets. Our work shows that molecular biology tools can be used to
design the interaction between fully folded proteins and inorganic surfaces
with potential applications in the bottom-up construction of plasmonic hybrid
nanomaterials.


  The inherent probabilistic nature of the biochemical reactions, and low copy
number of species can lead to stochasticity in gene expression across identical
cells. As a result, after induction of gene expression, the time at which a
specific protein count is reached is stochastic as well. Therefore events
taking place at a critical protein level will see stochasticity in their
timing. First-passage time (FPT), the time at which a stochastic process hits a
critical threshold, provides a framework to model such events. Here, we
investigate stochasticity in FPT. Particularly, we consider events for which
controlling stochasticity is advantageous. As a possible regulatory mechanism,
we also investigate effect of auto-regulation, where the transcription rate of
gene depends on protein count, on stochasticity of FPT. Specifically, we
investigate for an optimal auto-regulation which minimizes stochasticity in
FPT, given fixed mean FPT and threshold.
  For this purpose, we model the gene expression at a single cell level. We
find analytic formulas for statistical moments of the FPT in terms of model
parameters. Moreover, we examine the gene expression model with
auto-regulation. Interestingly, our results show that the stochasticity in FPT,
for a fixed mean, is minimized when the transcription rate is independent of
protein count. Further, we discuss the results in context of lysis time of an
\textit{E. coli} cell infected by a $\lambda$ phage virus. An optimal lysis
time provides evolutionary advantage to the $\lambda$ phage, suggesting a
possible regulation to minimize its stochasticity. Our results indicate that
there is no auto-regulation of the protein responsible for lysis. Moreover,
congruent to experimental evidences, our analysis predicts that the expression
of the lysis protein should have a small burst size.


  Motivation: Gene set testing is typically performed in a supervised context
to quantify the association between groups of genes and a clinical phenotype.
In many cases, however, a gene set-based interpretation of genomic data is
desired in the absence of a phenotype variable. Although methods exist for
unsupervised gene set testing, they predominantly compute enrichment relative
to clusters of the genomic variables with performance strongly dependent on the
clustering algorithm and number of clusters. Results: We propose a novel
method, spectral gene set enrichment (SGSE), for unsupervised competitive
testing of the association between gene sets and empirical data sources. SGSE
first computes the statistical association between gene sets and principal
components (PCs) using our principal component gene set enrichment (PCGSE)
method. The overall statistical association between each gene set and the
spectral structure of the data is then computed by combining the PC-level
p-values using the weighted Z-method with weights set to the PC variance scaled
by Tracey-Widom test p-values. Using simulated data, we show that the SGSE
algorithm can accurately recover spectral features from noisy data. To
illustrate the utility of our method on real data, we demonstrate the superior
performance of the SGSE method relative to standard cluster-based techniques
for testing the association between MSigDB gene sets and the variance structure
of microarray gene expression data. Availability:
http://cran.r-project.org/web/packages/PCGSE/index.html Contact:
rob.frost@dartmouth.edu or jason.h.moore@dartmouth.edu


  Understanding the pathways by which viral capsid proteins assemble around
their genomes could identify key intermediates as potential drug targets. In
this work we use computer simulations to characterize assembly over a wide
range of capsid protein-protein interaction strengths and solution ionic
strengths. We find that assembly pathways can be categorized into two classes,
in which intermediates are either predominantly ordered or disordered. Our
results suggest that estimating the protein-protein and the protein-genome
binding affinities may be sufficient to predict which pathway occurs.
Furthermore, the calculated phase diagrams suggest that knowledge of the
dominant assembly pathway and its relationship to control parameters could
identify optimal strategies to thwart or redirect assembly to block infection.
Finally, analysis of simulation trajectories suggests that the two classes of
assembly pathways can be distinguished in single molecule fluorescence
correlation spectroscopy or bulk time resolved small angle x-ray scattering
experiments.


  A prototype for a web application was designed and implemented as a guide to
be used by clinicians when designing the best drug therapy for a specific
cancer patient, given biological data derived from the patients tumor tissue
biopsy. A representation of the patients metabolic pathways is displayed as a
graph in the application, with nodes as substrates and products and edges as
enzymes. The top metabolically active sub- paths in the pathway, ranked using
an algorithm based on both the patients biological data and the graph topology,
are also displayed and can be individually highlighted to examine potential
enzymatic sites to be disrupted by a drug.


  Shapes of interacting RNA complexes are studied using a filtration via their
topological genus. A shape of an RNA complex is obtained by (iteratively)
collapsing stacks and eliminating hairpin loops. This shape-projection
preserves the topological core of the RNA complex and for fixed topological
genus there are only finitely many such shapes.Our main result is a new
bijection that relates the shapes of RNA complexes with shapes of RNA
structures.This allows to compute the shape polynomial of RNA complexes via the
shape polynomial of RNA structures. We furthermore present a linear time
uniform sampling algorithm for shapes of RNA complexes of fixed topological
genus.


  The fitness landscape - the mapping between genotypes and fitness -
determines properties of the process of adaptation. Several small genetic
fitness landscapes have recently been built by selecting a handful of
beneficial mutations and measuring fitness of all combinations of these
mutations. Here we generate several testable predictions for the properties of
these landscapes under Fisher's geometric model of adaptation (FGMA). When far
from the fitness optimum, we analytically compute the fitness effect of
beneficial mutations and their epistatic interactions. We show that epistasis
may be negative or positive on average depending on the distance of the
ancestral genotype to the optimum and whether mutations were independently
selected or co-selected in an adaptive walk. Using simulations, we show that
genetic landscapes built from FGMA are very close to an additive landscape when
the ancestral strain is far from the optimum. However, when close to the
optimum, a large diversity of landscape with substantial ruggedness and sign
epistasis emerged. Strikingly, landscapes built from different realizations of
stochastic adaptive walks in the same exact conditions were highly variable,
suggesting that several realizations of small genetic landscapes are needed to
gain information about the underlying architecture of the global adaptive
landscape.


  Bayesian model averaging (BMA) is the state of the art approach for
overcoming model uncertainty. Yet, especially on small data sets, the results
yielded by BMA might be sensitive to the prior over the models. Credal Model
Averaging (CMA) addresses this problem by substituting the single prior over
the models by a set of priors (credal set). Such approach solves the problem of
how to choose the prior over the models and automates sensitivity analysis. We
discuss various CMA algorithms for building an ensemble of logistic regressors
characterized by different sets of covariates. We show how CMA can be
appropriately tuned to the case in which one is prior-ignorant and to the case
in which instead domain knowledge is available. CMA detects prior-dependent
instances, namely instances in which a different class is more probable
depending on the prior over the models. On such instances CMA suspends the
judgment, returning multiple classes. We thoroughly compare different BMA and
CMA variants on a real case study, predicting presence of Alpine marmot burrows
in an Alpine valley. We find that BMA is almost a random guesser on the
instances recognized as prior-dependent by CMA.


  The KaiABC circadian clock from cyanobacteria is the only known three-protein
oscillatory system which can be reconstituted outside the cell and which
displays sustained periodic dynamics in various molecular state variables.
Despite many recent experimental and theoretical studies there are several open
questions regarding the central mechanism(s) responsible for creating this ~24
hour clock in terms of molecular assembly/disassembly of the proteins and
site-dependent phosphorylation and dephosphorylation of KaiC monomers.
Simulations of protein-protein interactions and phosphorylation reactions
constrained by analytical fits to partial reaction experimental data support
the central mechanism of oscillation as KaiB-induced KaiA sequestration in
KaiABC complexes associated with the extent of Ser431 phosphorylation in KaiC
hexamers. A simple two-state deterministic model in terms of the degree of
phosphorylation of Ser431 and Thr432 sites alone can reproduce the previously
observed circadian oscillation in the four population monomer phosphoforms in
terms of waveform, amplitude and phase. This suggests that a cyclic
phosphorylation scheme (involving cooperativity between adjacent Ser431 and
Thr432 sites) is not necessary for creating oscillations. Direct simulations of
the clock predict the minimum number of serine-only monomer subunits associated
with KaiA sequestration and release, highlight the role of monomer exchange in
rapid synchronization, and predict the average number of KaiA dimers
sequestered per KaiC hexamer.


  In this PhD thesis, a novel method to determine protein structures using
chemical shifts is presented.


  Boolean networks are discrete dynamical systems for modeling regulation and
signaling in living cells. We investigate a particular class of Boolean
functions with inhibiting inputs exerting a veto (forced zero) on the output.
We give analytical expressions for the sensitivity of these functions and
provide evidence for their role in natural systems. In an intracellular signal
transduction network [Helikar et al., PNAS (2008)], the functions with veto are
over-represented by a factor exceeding the over-representation of threshold
functions and canalyzing functions in the same system. In Boolean networks for
control of the yeast cell cycle [Fangting Li et al., PNAS (2004), Davidich et
al., PLoS One (2009)], none or minimal changes to the wiring diagrams are
necessary to formulate their dynamics in terms of the veto functions introduced
here.


  The convergent interests of different scientific disciplines, from
biochemistry to electronics, toward the investigation of protein electrical
properties, has promoted the development of a novel bailiwick, the so called
proteotronics. The main aim of proteotronics is to propose and achieve
innovative electronic devices, based on the selective action of specific
proteins. This paper gives a sketch of the fields of applications of
proteotronics, by using as significant example the detection of a specific
odorant molecule carried out by an olfactory receptor. The experiment is
briefly reviewed and its theoretical interpretation given. Further experiments
are envisioned and expected results discussed in the perspective of an
experimental validation.


  The transcriptional regulation of gene expression is orchestrated by complex
networks of interacting genes. Increasing evidence indicates that these
transcriptional regulatory networks (TRNs) in bacteria have an inherently
hierarchical architecture, although the design principles and the specific
advantages offered by this type of organization have not yet been fully
elucidated. In this study, we focussed on the hierarchical structure of the TRN
of the gram-positive bacterium Bacillus subtilis and performed a comparative
analysis with the TRN of the gram-negative bacterium Escherichia coli. Using a
graph-theoretic approach, we organized the transcription factors (TFs) and
sigma-factors in the TRNs of B. subtilis and E. coli into three hierarchical
levels (Top, Middle and Bottom) and studied several structural and functional
properties across them. In addition to many similarities, we found also
specific differences, explaining the majority of them with variations in the
distribution of sigma-factors across the hierarchical levels in the two
organisms. We then investigated the control of target metabolic genes by
transcriptional regulators to characterize the differential regulation of three
distinct metabolic subsystems (catabolism, anabolism and central energy
metabolism). These results suggest that the hierarchical architecture that we
observed in B. subtilis represents an effective organization of its TRN to
achieve flexibility in the response to diverse stimuli.


  The issue of the nucleation and slow closure mechanisms of non superhelical
stress-induced denaturation bubbles in DNA is tackled using coarse-grained
MetaDynamics and Brownian simulations. A minimal mesoscopic model is used where
the double helix is made of two interacting bead-spring rotating strands with a
prescribed torsional modulus in the duplex state. We demonstrate that
timescales for the nucleation (resp. closure) of an approximately 10 base-pair
bubble, in agreement with experiments, are associated with the crossing of a
free-energy barrier of $22~k_{\rm B}T$ (resp. $13~k_{\rm B}T$) at room
temperature $T$. MetaDynamics allows us to reconstruct accurately the
free-energy landscape, to show that the free-energy barriers come from the
difference in torsional energy between the bubble and duplex states, and thus
to highlight the limiting step, a collective twisting, that controls the
nucleation/closure mechanism, and to access opening time scales on the
millisecond range. Contrary to small breathing bubbles, these more than
4~base-pair bubbles are of biological relevance, for example when a preexisting
state of denaturation is required by specific DNA-binding proteins.


  We develop a detailed theoretical framework for various types of
transcription factor gene oscillators. We further demonstrate that one can
build genetic-oscillators which are tunable and robust against perturbations in
the critical control parameters by coupling two or more independent
Goodwin-Griffith oscillators through either -OR- or -AND- type logic. Most of
the coupled oscillators constructed in the literature so far seem to be of -OR-
type. When there are transient perturbations in one of the -OR- type
coupled-oscillators, then the overall period of the system remains constant
(period-buffering) whereas in case of -AND- type coupling the overall period of
the system moves towards the perturbed oscillator. Though there is a
period-buffering, the amplitudes of oscillators coupled through -OR- type logic
are more sensitive to perturbations in the parameters associated with the
promoter state dynamics than -AND- type. Further analysis shows that the period
of -AND- type coupled dual-feedback oscillators can be tuned without conceding
on the amplitudes. Using these results we derive the basic design principles
governing the robust and tunable synthetic gene oscillators without
compromising on their amplitudes.


  Most new mutations are deleterious and are eventually eliminated by natural
selection. But in an adapting population, the rapid amplification of beneficial
mutations can hinder the removal of deleterious variants in nearby regions of
the genome, altering the patterns of sequence evolution. Here, we analyze the
interactions between beneficial "driver" mutations and linked deleterious
"passengers" during the course of adaptation. We derive analytical expressions
for the substitution rate of a deleterious mutation as a function of its
fitness cost, as well as the reduction in the beneficial substitution rate due
to the genetic load of the passengers. We find that the fate of each
deleterious mutation varies dramatically with the rate and spectrum of
beneficial mutations, with a non-monotonic dependence on both the population
size and the rate of adaptation. By quantifying this dependence, our results
allow us to estimate which deleterious mutations will be likely to fix, and how
many of these mutations must arise before the progress of adaptation is
significantly reduced.


  Auto-regulation, a process wherein a protein negatively regulates its own
production, is a common motif in gene expression networks. Negative feedback in
gene expression plays a critical role in buffering intracellular fluctuations
in protein concentrations around optimal value. Due to the nonlinearities
present in these feedbacks, moment dynamics are typically not closed, in the
sense that the time derivative of the lower-order statistical moments of the
protein copy number depends on high-order moments. Moment equations are closed
by expressing higher-order moments as nonlinear functions of lower-order
moments, a technique commonly referred to as moment closure. Here, we compare
the performance of different moment closure techniques. Our results show that
the commonly used closure method, which assumes a priori that the protein
population counts are normally distributed, performs poorly. In contrast,
conditional derivative matching, a novel closure scheme proposed here provides
a good approximation to the exact moments across different parameter regimes.
In summary our study provides a new moment closure method for studying
stochastic dynamics of genetic negative feedback circuits, and can be extended
to probe noise in more complex gene networks.


  The deep connection between thermodynamics, computation, and information is
now well established both theoretically and experimentally. Here, we extend
these ideas to show that thermodynamics also places fundamental constraints on
statistical estimation and learning. To do so, we investigate the constraints
placed by (nonequilibrium) thermodynamics on the ability of biochemical
signaling networks within cells to estimate the concentration of an external
signal. We show that accuracy is limited by energy consumption, suggesting that
there are fundamental thermodynamic constraints on statistical inference.


  Moment closure on general discrete structures often requires one of the
following: (i) an absence of short closed loops (zero clustering); (ii)
existence of a spatial scale; (iii) ad hoc assumptions. Algebraic methods are
presented to avoid the use of such assumptions for populations based on clumps,
and are applied to both SIR and macroparasite disease dynamics. One approach
involves a series of approximations that can be derived systematically, and
another is exact and based on Lie algebraic methods.


  The driving force behind the origin and evolution of life has been the
thermodynamic imperative of increasing the entropy production of the biosphere
through increasing the global solar photon dissipation rate. In the upper
atmosphere of today, oxygen and ozone derived from life processes are
performing the short wavelength UVC and UVB dissipation. On Earth's surface,
water and organic pigments in water facilitate the near UV and visible photon
dissipation. The first organic pigments probably formed, absorbed, and
dissipated at those photochemically active wavelengths in the UVC that could
have reached Earth's surface during the Archean. Proliferation of these
pigments can be understood as an autocatalytic photochemical process obeying
non-equilibrium thermodynamic directives related to increasing solar photon
dissipation rate. Under these directives, organic pigments would have evolved
over time to increase the global photon dissipation rate by; 1) increasing the
ratio of their effective photon cross sections to their physical size, 2)
decreasing their electronic excited state life times, 3) quenching
non-radiative de-excitation channels (e.g. fluorescence), 4) covering ever more
completely the solar spectrum, and 5) dispersing into an ever greater surface
area of Earth. From knowledge of the evolution of the spectrum of G-type stars,
and considering the most probable history of the transparency of Earths
atmosphere, we construct the most probable surface solar spectrum as a function
of time and compare this with the history of molecular absorption maxima
obtained from the available data in the literature. This comparison supports
the thermodynamic dissipation theory for the origin of life, constrains models
for Earth's early atmosphere, and sheds some new light on the origin of
photosynthesis.


  Cellular automata have been useful artificial models for exploring how
relatively simple rules combined with spatial memory can give rise to complex
emergent patterns. Moreover, studying the dynamics of how rules emerge under
artificial selection for function has recently become a powerful tool for
understanding how evolution can innovate within its genetic rule space.
However, conventional cellular automata lack the kind of state feedback that is
surely present in natural evolving systems. Each new generation of a population
leaves an indelible mark on its environment and thus affects the selective
pressures that shape future generations of that population. To model this
phenomenon, we have augmented traditional cellular automata with
state-dependent feedback. Rather than generating automata executions from an
initial condition and a static rule, we introduce mappings which generate
iteration rules from the cellular automaton itself. We show that these new
automata contain disconnected regions which locally act like conventional
automata, thus encapsulating multiple functions into one structure.
Consequently, we have provided a new model for processes like cell
differentiation. Finally, by studying the size of these regions, we provide
additional evidence that the dynamics of self-reference may be critical to
understanding the evolution of natural language. In particular, the rules of
elementary cellular automata appear to be distributed in the same way as words
in the corpus of a natural language.


  We study stochastic evolution of optional games on simple graphs. There are
two strategies, A and B, whose interaction is described by a general payoff
matrix. In addition there are one or several possibilities to opt out from the
game by adopting loner strategies. Optional games lead to relaxed social
dilemmas. Here we explore the interaction between spatial structure and
optional games. We find that increasing the number of loner strategies (or
equivalently increasing mutational bias toward loner strategies) facilitates
evolution of cooperation both in well-mixed and in structured populations. We
derive various limits for weak selection and large population size. For some
cases we derive analytic results for strong selection. We also analyze strategy
selection numerically for finite selection intensity and discuss combined
effects of optionality and spatial structure.


  We consider a predator-prey population model with prey gathering together for
defense purposes. A transmissible unrecoverable disease affects the prey. We
characterize the system behavior, establishing that ultimately either only the
susceptible prey survive, or the disease becomes endemic, but the predators are
wiped out. Another alternative is that the disease is eradicated, with sound
prey and predators thriving at an equilibrium or through persistent population
oscillations. Finally, the populations can thrive together, with the endemic
disease. The only impossible alternative in these circumstances is predators
thriving just with infected prey. But this follows from the model assumptions,
in that infected prey are too weak to sustain themselves. A mathematical
peculiarity of the model is the singularity-free reformulation, which leads to
three entirely new dependent variables to describe the system. The model is
then extended to encompass the situation in which ingestion of diseased prey is
fatal for the predators.


  We present two new models for interacting populations subject to a
transmissible disease. The novelty lies in the assumption that herd behavior
influences the disease incidence, rather than the demographic description of
the interactions, as in previous related similar models. As it is already known
from other ecoepidemiological situations, the epidemics may affect the system
demographic outcomes.


  Deep transcriptome sequencing has revealed the existence of many transcripts
that lack long or conserved open reading frames and which have been termed long
non-coding RNAs (lncRNAs). Despite the existence of several well-characterized
lncRNAs that play roles in the regulation of gene expression, the vast majority
of them do not yet have a known function. Motivated by the existence of
ribosome profiling data for several species, we have tested the hypothesis that
they may act as a repository for the synthesis of new peptides using data from
human, mouse, zebrafish, fruit fly, Arabidopsis and yeast. The ribosome
protection patterns are consistent with the presence of translated open reading
frames (ORFs) in a very large number of lncRNAs. Most of the ribosome-protected
ORFs are shorter than 100 amino acids and usually cover less than half the
transcript. Ribosome density in these ORFs is high and contrasts sharply with
the 3UTR region, in which very often there is no detectable ribosome binding,
similar to bona fide protein-coding genes. The coding potential of
ribosome-protected ORFs, measured using hexamer frequencies, is significantly
higher than that of randomly selected intronic ORFs and similar to that of
evolutionary young coding sequences. Selective constraints in
ribosome-protected ORFs from lncRNAs are lower than in typical protein-coding
genes but again similar to young proteins. These results strongly suggest that
lncRNAs play an important role in de novo protein evolution.


  We define a new strategy for population games based on techniques from
machine learning and statistical inference that is essentially uninvadable and
can successfully invade (significantly more likely than a neutral mutant)
essentially all known memory-one strategies for the prisoner's dilemma and
other population games, including ALLC (always cooperate), ALLD (always
defect), tit-for-tat (TFT), win-stay-lose-shift (WSLS), and zero determinant
(ZD) strategies, including extortionate and generous strategies. We will refer
to a player using this strategy as an "information player" and the specific
implementation as $IP_0$. Such players use the history of play to identify
opponent's strategies and respond accordingly, and naturally learn to cooperate
with each other.


  In a paper of August 2013, I discussed the so-called SuperSpreader (SS)
epidemic model and emphasized that it has dynamics differing greatly from the
more-familiar uniform (or Poisson) textbook model. In that paper, SARS in 2003
was the representative instance and it was suggested that MERS may be another.
In April 2014, MERS incident cases showed a spectacular spike (going from a
handful in the previous April to more than 260 in that month of 2014)
reminiscent of a figure I published nine months earlier. Here I refit the
two-level and several variant SS models to incident data from January 1,
2013--April 30, 2014 and conclude that MERS will go pandemic (all other factors
remaining the same). In addition, I discuss a number of model-realism and
fitting methodology issues relevant to analysing SS epidemics.


  The fate of every eukaryotic cell subtly relies on the exceptional mechanical
properties of microtubules. Despite significant efforts, understanding their
unusual mechanics remains elusive. One persistent, unresolved mystery is the
formation of long-lived arcs and rings, e.g. in kinesin-driven gliding assays.
To elucidate their physical origin we develop a model of the inner workings of
the microtubule's lattice, based on recent experimental evidence for a
conformational switch of the tubulin dimer. We show that the microtubule
lattice itself coexists in discrete polymorphic states. Curved states can be
induced via a mechanical hysteresis involving torques and forces typical of few
molecular motors acting in unison. This lattice switch renders microtubules not
only virtually unbreakable under typical cellular forces, but moreover provides
them with a tunable response integrating mechanical and chemical stimuli.


  Toxic blooms of Lyngbya majuscula occur in coastal areas worldwide and have
major ecological, health and economic consequences. The exact causes and
combinations of factors which lead to these blooms are not clearly understood.
Lyngbya experts and stakeholders are a particularly diverse group, including
ecologists, scientists, state and local government representatives, community
organisations, catchment industry groups and local fishermen. An integrated
Bayesian network approach was developed to better understand and model this
complex environmental problem, identify knowledge gaps, prioritise future
research and evaluate management options.


  Inputs to signaling pathways can have complex statistics that depend on the
environment and on the behavioral response to previous stimuli. Such behavioral
feedback is particularly important in navigation. Successful navigation relies
on proper coupling between sensors, which gather information during motion, and
actuators, which control behavior. Because reorientation conditions future
inputs, behavioral feedback can place sensors and actuators in an operational
regime different from the resting state. How then can organisms maintain proper
information transfer through the pathway while navigating diverse environments?
In bacterial chemotaxis, robust performance is often attributed to the zero
integral feedback control of the sensor, which guarantees that activity returns
to resting state when the input remains constant. While this property provides
sensitivity over a wide range of signal intensities, it remains unclear how
other parameters affect chemotactic performance, especially when considering
that the swimming behavior of the cell determines the input signal. Using
analytical models and simulations that incorporate recent experimental
evidences about behavioral feedback and flagellar motor adaptation we identify
an operational regime of the pathway that maximizes drift velocity for various
environments and sensor adaptation rates. This optimal regime is outside the
dynamic range of the motor response, but maximizes the contrast between run
duration up and down gradients. In steep gradients, the feedback from
chemotactic drift can push the system through a bifurcation. This creates a
non-chemotactic state that traps cells unless the motor is allowed to adapt.
Although motor adaptation helps, we find that as the strength of the feedback
increases individual phenotypes cannot maintain the optimal operational regime
in all environments, suggesting that diversity could be beneficial.


  Dispersal of species is a fundamental ecological process in the evolution and
maintenance of biodiversity. Limited control over ecological parameters has
hindered progress in understanding of what enables species to colonise new
area, as well as the importance of inter-species interactions. Such control is
necessary to construct reliable mathematical models of ecosystems. In our work,
we studied dispersal in the context of bacterial range expansions and
identified the major determinants of species coexistence for a bacterial model
system of three Escherichia coli strains (toxin producing, sensitive, and
resistant). Genetic engineering allowed us to tune strain growth rates and to
design different ecological scenarios (cyclic and hierarchical). We found that
coexistence of all strains depended on three strongly interdependent factors:
composition of inoculum, relative strain growth rates, and effective toxin
range. Robust agreement between our experiments and a thoroughly calibrated
computational model enabled us to extrapolate these intricate interdependencies
in terms of phenomenological biodiversity laws. Our mathematical analysis also
suggested that cyclic dominance between strains is not a prerequisite for
coexistence in competitive range expansions. Instead, robust three-strain
coexistence required a balance between growth rates and either a reduced
initial ratio of the toxin-producing strain, or a sufficiently short toxin
range.


  We studied a relations between the triplet frequency composition of
mitochondria genomes, and the phylogeny of their bearers. First, the clusters
in 63dimensional space were developed due to $K$-means. Second, the clade
composition of those clusters has been studied. It was found that genomes are
distributed among the clusters very regularly, with strong correlation to
taxonomy. Strong co-evolution manifests through this correlation: the proximity
in frequency space was determined over the mitochondrion genomes, while the
proximity in taxonomy was determined morphologically.


  Background: Several sources of noise obfuscate the identification of single
nucleotide variation (SNV) in next generation sequencing data. For instance,
errors may be introduced during library construction and sequencing steps. In
addition, the reference genome and the algorithms used for the alignment of the
reads are further critical factors determining the efficacy of variant calling
methods. It is crucial to account for these factors in individual sequencing
experiments.
  Results: We introduce a simple data-adaptive model for variant calling. This
model automatically adjusts to specific factors such as alignment errors. To
achieve this, several characteristics are sampled from sites with low mismatch
rates, and these are used to estimate empirical log-likelihoods. These
likelihoods are then combined to a score that typically gives rise to a mixture
distribution. From these we determine a decision threshold to separate
potentially variant sites from the noisy background.
  Conclusions: In simulations we show that our simple proposed model is
competitive with frequently used much more complex SNV calling algorithms in
terms of sensitivity and specificity. It performs specifically well in cases
with low allele frequencies. The application to next-generation sequencing data
reveals stark differences of the score distributions indicating a strong
influence of data specific sources of noise. The proposed model is specifically
designed to adjust to these differences.


  We describe a conceptual design of a distributed classifier formed by a
population of genetically engineered microbial cells. The central idea is to
create a complex classifier from a population of weak or simple classifiers. We
create a master population of cells with randomized synthetic biosensor
circuits that have a broad range of sensitivities towards chemical signals of
interest that form the input vectors subject to classification. The randomized
sensitivities are achieved by constructing a library of synthetic gene circuits
with randomized control sequences (e.g. ribosome-binding sites) in the front
element. The training procedure consists in re-shaping of the master population
in such a way that it collectively responds to the "positive" patterns of input
signals by producing above-threshold output (e.g. fluorescent signal), and
below-threshold output in case of the "negative" patterns. The population
re-shaping is achieved by presenting sequential examples and pruning the
population using either graded selection/counterselection or by
fluorescence-activated cell sorting (FACS). We demonstrate the feasibility of
experimental implementation of such system computationally using a realistic
model of the synthetic sensing gene circuits.


  In this paper we study the interplay between epidemic spreading and risk
perception on multiplex networks. The basic idea is that the effective
infection probability is affected by the perception of the risk of being
infected, which we assume to be related to the fraction of infected neighbours,
as introduced by Bagnoli et al., PRE 76:061904 (2007). We re-derive previous
results using a self-organized method, that automatically gives the percolation
threshold in just one simulation. We then extend the model to multiplex
networks considering that people get infected by contacts in real life but
often gather information from an information networks, that may be quite
different from the real ones. The similarity between the real and information
networks determine the possibility of stopping the infection for a sufficiently
high precaution level: if the networks are too different there is no mean of
avoiding the epidemics.


  We introduce a simple ecological model describing the spatial organization of
two interacting populations whose individuals are indifferent to conspecifics
and avoid the proximity to heterospecifics. At small population densities
$\Phi$ a non-trivial structure is observed where clusters of individuals
arrange into a rhomboidal bipartite network with an average degree of four. For
$\Phi\rightarrow0$ the length scale, order parameter and susceptibility of the
network exhibit power-law divergences compatible with hyper-scaling, suggesting
the existence of a zero density - non-trivial - critical point. At larger
densities a critical threshold $\Phi_{c}$ is identified above which the
evolution toward a partially ordered configuration is prevented and the system
becomes jammed in a fully mixed state.


  The design and synthesis of complex and large mimicked biochemical networks
de novo is an unsolved problem in synthetic biology. To address this limitation
without resorting to ad hoc computations and experiments, a predictive
mathematical theory is required to reduce these complex chemical networks into
natural physico-chemical expressions. Here we provide a theory that offers a
physico-chemical expression for a large chemical network that is almost
arbitrarily both nonlinear and complex. Unexpectedly, the theory demonstrates
that such networks can be decomposed into reactions based solely on the Hill
equation, a simple chemical logic gate. This theory, analogous to implemented
electrical logic gates or functional algorithms in a computer, is proposed for
implementing regulated sequences of functional chemical reactions, such as
mimicked genes, transcriptional regulation, signal transduction, protein
interaction, and metabolic networks, into an artificial designed chemical
network.


  Motivated by studies on the recurrent properties of animal and human
mobility, we introduce a path-dependent random walk model with long range
memory for which not only the mean square displacement (MSD) can be obtained
exactly in the asymptotic limit, but also the propagator. The model consists of
a random walker on a lattice, which, at a constant rate, stochastically
relocates at a site occupied at some earlier time. This time in the past is
chosen randomly according to a memory kernel, whose temporal decay can be
varied via an exponent parameter. In the weakly non-Markovian regime, memory
reduces the diffusion coefficient from the bare value. When the mean backward
jump in time diverges, the diffusion coefficient vanishes and a transition to
an anomalous subdiffusive regime occurs. Paradoxically, at the transition, the
process is an anti-correlated L\'evy flight. Although in the subdiffusive
regime the model exhibits some features of the continuous time random walk with
infinite mean waiting time, it belongs to another universality class. If memory
is very long-ranged, a second transition takes place to a regime characterized
by a logarithmic growth of the MSD with time. In this case the process is
asymptotically Gaussian and effectively described as a scaled Brownian motion
with a diffusion coefficient decaying as 1/t.


  The cytoskeleton is regulated by a plethora of enzymes that influence the
stability and dynamics of cytoskeletal filaments. Molecular motors of the
kinesin-8 protein family depolymerise microtubules in a length-dependent
manner, and experimental and theoretical evidence suggest a role for kinesin-8
in the dynamic regulation of microtubules. However, so far the detailed
molecular mechanisms how these molecular motors interact with the growing
microtubule tip remain elusive. Here we investigate two interaction scenarios
for kinesin-8 and the microtubule tip. We give a comprehensive analysis of
regimes where length-regulation is possible and characterise how the stationary
length depends on the biochemical rates and the bulk concentrations of the
various proteins. For a neutral scenario, where microtubules grow irrespective
of whether the microtubule tip is occupied by a molecular motor, length
regulation is possible only for a narrow range of biochemical rates and limited
to small polymerisation rates. In contrast, for an inhibition scenario, where
the presence of a motor at the microtubule tip inhibits microtubule growth, the
regime of length regulation is extremely broad and includes high growth rates.
These results also apply to situations where polymerising enzymes, like
XMAP215, and kinesin-8 mutually exclude each other from the microtubule tip. We
also investigate the stochastic dynamics of the two scenarios. While for the
neutral scenario length is tightly controlled, length dynamics is intermittent
for the inhibition scenario and exhibits extended periods of microtubule growth
and shrinkage, reminiscent of microtubule dynamic instability. On a broader
perspective, the set of models established in this work quite generally
suggests that mutual exclusion of molecules at the ends of cytoskeletal
filaments is an important factor for filament dynamics and regulation.


  Escherichia coli and Bacillus subtilis are two of the best-studied
prokaryotic model organisms. Previous analyses of their transcriptional
regulatory networks have shown that they exhibit high plasticity during
evolution and suggested that both converge to scale-free-like structures.
Nevertheless, beyond this suggestion, no analyses have been carried out to
identify the common systems-level components and principles governing these
organisms. Here we show that these two phylogenetically distant organisms
follow a set of common novel biologically consistent systems principles
revealed by the mathematically and biologically founded natural decomposition
approach. The discovered common functional architecture is a diamond-shaped,
matryoshka-like, three-layer (coordination, processing, and integration)
hierarchy exhibiting feedback, which is shaped by four systems-level
components: global transcription factors (global TFs), locally autonomous
modules, basal machinery and intermodular genes. The first mathematical
criterion to identify global TFs, the $\kappa$-value, was reassessed on B.
subtilis and confirmed its high predictive power by identifying all the
previously reported, plus three potential, master regulators and eight sigma
factors. The functional conserved cores of modules, basal cell machinery, and a
set of non-orthologous common physiological global responses were identified
via both orthologous genes and non-orthologous conserved functions. This study
reveals novel common systems principles maintained between two phylogenetically
distant organisms and provides a comparison of their lifestyle adaptations. Our
results shed new light on the systems-level principles and the fundamental
functions required by bacteria to sustain life.


  Background: Previous studies have used different methods in an effort to
extract the modular organization of transcriptional regulatory networks (TRNs).
However, these approaches are not natural, as they try to cluster highly
connected genes into a module or locate known pleiotropic transcription factors
in lower hierarchical layers. Here, we unravel the TRN of Escherichia coli by
separating it into its key elements, thus revealing its natural organization.
We also present a mathematical criterion, based on the topological features of
the TRN, to classify the network elements into one of two possible classes:
hierarchical or modular genes.
  Results: We found that modular genes are clustered into physiologically
correlated groups validated by a statistical analysis of the enrichment of the
functional classes. Hierarchical genes encode transcription factors responsible
for coordinating module responses based on general interest signals.
Hierarchical elements correlate highly with the previously studied global
regulators, suggesting that this could be the first mathematical method to
identify global regulators. We identified a new element in TRNs never described
before: intermodular genes. These are structural genes that integrate, at the
promoter level, signals coming from different modules, and therefore from
different physiological responses. Using the concept of pleiotropy, we have
reconstructed the hierarchy of the network and discuss the role of feedforward
motifs in shaping the hierarchical backbone of the TRN.
  Conclusions: This study sheds new light on the design principles underpinning
the organization of TRNs, showing a novel nonpyramidal architecture composed of
independent modules globally governed by hierarchical transcription factors,
whose responses are integrated by intermodular genes.


  We consider the lifetimes of metastable states in bistable evolutionary games
(coordination games), and examine how they are affected by spatial structure. A
semiclassical approximation based on a path integral method is applied to
stochastic evolutionary game dynamics with and without spatial structure, and
the lifetimes of the metastable states are evaluated. It is shown that the
population dependence of the lifetimes is qualitatively different in these two
models. Our result indicates that spatial structure can accelerate the
transitions between metastable states.


  We introduce a variational approximation to the microscopic dynamics of rare
conformational transitions of macromolecules. Within this framework it is
possible to simulate on a small computer cluster reactions as complex as
protein folding, using state of the art all-atom force fields in explicit
solvent. We test this method against molecular dynamics (MD) simulations of the
folding of an alpha- and a beta-protein performed with the same all-atom force
field on the Anton supercomputer. We find that our approach yields results
consistent with those of MD simulations, at a computational cost orders of
magnitude smaller.


  Dengue is a viral infection that can produce a severe fever and symptoms that
may require hospitalization. It is transmitted between humans by the
urban-adapted, day-biting Aedes mosquitoes and is therefore a particular
problem in towns and cities. To explore this risk, one has assessed the
potential levels of exposure by a climate-driven model for dengue risk
transmission in Brazil and records of its seasonal variation at the key sites.
Like the weather, it is unworkable to forecast the precise situation with
regard to dengue in Brazil in 2014. One can, however, make informed guesses on
the basis of averaged records of dengue in previous years. For the areas around
the World Cup stadiums, these records show that the main dengue season will
have passed before the World Cup is held in June and July. Unfortunately, the
risk remains, even this is low but not negliegible, during these months in the
Brazilian north and northeast. But, in fact, the risk of an outbreak of dengue
fever during the upcoming soccer World Cup in Brazil is not serious enough to
warrant a high alert in the host cities, according to a reliable early warning
system for the disease.


  Clustering is the propensity of nodes that share a common neighbour to be
connected. It is ubiquitous in many networks but poses many modelling
challenges. Clustering typically manifests itself by a higher than expected
frequency of triangles, and this has led to the principle of constructing
networks from such building blocks. This approach has been generalised to
networks being constructed from a set of more exotic subgraphs. As long as
these are fully connected, it is then possible to derive mean-field models that
approximate epidemic dynamics well. However, there are virtually no results for
non-fully connected subgraphs. In this paper, we provide a general and
automated approach to deriving a set of ordinary differential equations, or
mean-field model, that describes, to a high degree of accuracy, the expected
values of system-level quantities, such as the prevalence of infection. Our
approach offers a previously unattainable degree of control over the
arrangement of subgraphs and network characteristics such as classical node
degree, variance and clustering. The combination of these features makes it
possible to generate families of networks with different subgraph compositions
while keeping classical network metrics constant. Using our approach, we show
that higher-order structure realised either through the introduction of loops
of different sizes or by generating clustered networks based on different
subgraphs, leads to significant differences in epidemic dynamics despite
controlling for basic network metrics.


  The reconstruction of a central tendency `species tree' from a large number
of conflicting gene trees is a central problem in systematic biology. Moreover,
it becomes particularly problematic when taxon coverage is patchy, so that not
all taxa are present in every gene tree. Here, we list four apparently
desirable properties that a method for estimating a species tree from gene
trees could have (the strongest property states that building a species tree
from input gene trees and then pruning leaves gives a tree that is the same as,
or more resolved than, the tree obtained by first removing the taxa from the
input trees and then building the species tree). We show that while it is
technically possible to simultaneously satisfy these properties when taxon
coverage is complete, they cannot all be satisfied in the more general
supertree setting. In part two, we discuss a concordance-based consensus method
based on Baum's `plurality clusters', and an extension to concordance
supertrees.


  Persistent Organic Pollutants represent a global ecological concern due to
their ability to accumulate in organisms and to spread species-by-species via
feeding connections. In this work we focus on the estimation and simulation of
the bioaccumulation dynamics of persistent pollutants in the marine ecosystem,
and we apply the approach for reconstructing a model of PCBs bioaccumulation in
the Adriatic sea, estimated after an extensive review of trophic and PCBs
concentration data on Adriatic species. Our estimations evidence the occurrence
of PCBs biomagnification in the Adriatic food web, together with a strong
dependence of bioaccumulation on trophic dynamics and external factors like
fishing activity.


  We study general properties of the leading eigenvalue $\overline{w}(q)$ of
Eigen's evolutionary matrices depending on the probability $q$ of faithful
reproduction. This is a linear algebra problem that has various applications in
theoretical biology, including such diverse fields as the origin of life,
evolution of cancer progression, and virus evolution. We present the exact
expressions for $\overline{w}(q),\overline{w}'(q),\overline{w}''(q)$ for
$q=0,0.5,1$ and prove that the absolute minimum of $\overline{w}(q)$, which
always exists, belongs to the interval $[0,0.5]$. For the specific case of a
single peaked landscape we also find lower and upper bounds on
$\overline{w}(q)$, which are used to estimate the critical mutation rate, after
which the distribution of the types of individuals in the population becomes
almost uniform. This estimate is used as a starting point to conjecture another
estimate, valid for any fitness landscape, and which is checked by numerical
calculations. The last estimate stresses the fact that the inverse dependence
of the critical mutation rate on the sequence length is not a generally valid
fact. Therefore, the discussions of the error threshold applied to biological
systems must take this fact into account.


  One reason why classical phylogenetic reconstruction methods fail to
correctly infer the underlying topology is because they assume oversimplified
models. In this paper we propose a topology reconstruction method consistent
with the most general Markov model of nucleotide substitution, which can also
deal with data coming from mixtures on the same topology. It is based on an
idea of Eriksson on using phylogenetic invariants and provides a system of
weights that can be used as input of quartet-based methods. We study its
performance on real data and on a wide range of simulated 4-taxon data (both
time-homogeneous and nonhomogeneous, with or without among-site rate
heterogeneity, and with different branch length settings). We compare it to the
classical methods of neighbor-joining (with paralinear distance), maximum
likelihood (with different underlying models), and maximum parsimony. Our
results show that this method is accurate and robust, has a similar performance
to ML when data satisfies the assumptions of both methods, and outperforms all
methods when these are based on inappropriate substitution models or when both
long and short branches are present. If alignments are long enough, then it
also outperforms other methods when some of its assumptions are violated.


  Elucidating the demographic and phylogeographic histories of species provides
insight into the processes responsible for generating biological diversity, and
genomic datasets are now permitting the estimation of histories and demographic
parameters with unprecedented accuracy. We used a genomic single nucleotide
polymorphism (SNP) dataset generated using a RAD-Seq method to investigate the
historical demography and phylogeography of a widespread lowland Neotropical
bird (Xenops minutus). As expected, we found that prominent landscape features
that act as dispersal barriers, such as Amazonian rivers and the Andes
Mountains, are associated with the deepest phylogeographic breaks, and also
that isolation by distance is limited in areas between these barriers. In
addition, we inferred positive population growth for most populations and
detected evidence of historical gene flow between populations that are now
physically isolated. Even with genomic estimates of historical demographic
parameters, we found the prominent diversification hypotheses to be untestable.
We conclude that investigations into the multifarious processes shaping species
histories, aided by genomic datasets, will provide greater resolution of
diversification in the Neotropics, but that future efforts should focus on
understanding the processes shaping the histories of lineages rather than
trying to reconcile these histories with landscape and climatic events in Earth
history.


  The scientific enterprise depends critically on the preservation of and open
access to published data. This basic tenet applies acutely to phylogenies
(estimates of evolutionary relationships among species). Increasingly,
phylogenies are estimated from increasingly large, genome-scale datasets using
increasingly complex statistical methods that require increasing levels of
expertise and computational investment. Moreover, the resulting phylogenetic
data provide an explicit historical perspective that critically informs
research in a vast and growing number of scientific disciplines. One such use
is the study of changes in rates of lineage diversification (speciation -
extinction) through time. As part of a meta-analysis in this area, we sought to
collect phylogenetic data (comprising nucleotide sequence alignment and tree
files) from 217 studies published in 46 journals over a 13-year period. We
document our attempts to procure those data (from online archives and by direct
request to corresponding authors), and report results of analyses (using
Bayesian logistic regression) to assess the impact of various factors on the
success of our efforts. Overall, complete phylogenetic data for ~60% of these
studies are effectively lost to science. Our study indicates that phylogenetic
data are more likely to be deposited in online archives and/or shared upon
request when: (1) the publishing journal has a strong data-sharing policy; (2)
the publishing journal has a higher impact factor, and; (3) the data are
requested from faculty rather than students. Although the situation appears
dire, our analyses suggest that it is far from hopeless: recent initiatives by
the scientific community -- including policy changes by journals and funding
agencies -- are improving the state of affairs.


  Experiments indicate that unbinding rates of proteins from DNA can depend on
the concentration of proteins in nearby solution. Here we present a theory of
multi-step replacement of DNA-bound proteins by solution-phase proteins. For
four different kinetic scenarios we calculate the depen- dence of protein
unbinding and replacement rates on solution protein concentration. We find (1)
strong effects of progressive 'rezipping' of the solution-phase protein onto
DNA sites liberated by 'unzipping' of the originally bound protein; (2) that a
model in which solution-phase proteins bind non-specifically to DNA can
describe experiments on exchanges between the non specific DNA- binding
proteins Fis-Fis and Fis-HU; (3) that a binding specific model describes
experiments on the exchange of CueR proteins on specific binding sites.


  Widely used models in genetics include the Wright-Fisher diffusion and its
moment dual, Kingman's coalescent. Each has a multilocus extension but under
neither extension is the sampling distribution available in closed-form, and
their computation is extremely difficult. In this paper we derive two new
multilocus population genetic models, one a diffusion and the other a
coalescent process, which are much simpler than the standard models, but which
capture their key properties for large recombination rates. The diffusion model
is based on a central limit theorem for density dependent population processes,
and we show that the sampling distribution is a linear combination of moments
of Gaussian distributions and hence available in closed-form. The coalescent
process is based on a probabilistic coupling of the ancestral recombination
graph to a simpler genealogical process which exposes the leading dynamics of
the former. We further demonstrate that when we consider the sampling
distribution as an asymptotic expansion in inverse powers of the recombination
parameter, the sampling distributions of the new models agree with the standard
ones up to the first two orders.


  A mathematical model of coastal forest growth is proposed to describe and
test the effects of salinity and inundation in the long-term growth performance
and carbon sequestration of monospecific mangrove ($Rhizophora\; mucronata$)
plantation in the Philippines. We used allometry in expressing the mangrove
growth equation, and stochasticity in scheduling population-level events that
drive the development of the mangrove forest. Analysis of the model unveils an
index, $\xi$, that could be used in assessing a strategy which could promote
optimum carbon-stock accumulation in the long run. If initial plot is
configured such that $\xi > 1$, the $R.\; mucronata$ plantations could achieve
an above-ground biomass per hectare (AGB) of $1000\mbox{ t/ha}$, or about
$500\mbox{ tC/ha}$, in approximately $200$ to $250$ years post planting. In
contrast, the current restoration strategy implemented in the Philippines
corresponds to the case that $\xi <1$. Consequently, the restored mangroves
could not achieve stable growth without the support of costly human assistance
such as frequent replanting. Rather, through that typical strategy and in the
absence of assistance, the AGB decreases with time until all trees die.
Mangrove restoration could therefore be planned strategically to mitigate
costly and wasteful implementation. The proposed index $\xi$ thus serves as an
early indicator for the progress or demise of restored mangroves.


  We study translocation dynamics of a driven compressible semi-flexible chain
consisting of alternate blocks of stiff ($S$) and flexible ($F$) segments of
size $m$ and $n$ respectively for different chain length $N$ in two dimension
(2D). The free parameters in the model are the bending rigidity $\kappa_b$
which controls the three body interaction term, the elastic constant $k_F$ in
the FENE (bond) potential between successive monomers, as well as the segmental
lengths $m$ and $n$ and the repeat unit $p$ ($N=m_pn_p$) and the solvent
viscosity $\gamma$. We demonstrate that due to the change in entropic barrier
and the inhomogeneous viscous drag on the chain backbone a variety of scenarios
are possible amply manifested in the waiting time distribution of the
translocating chain. These information can be deconvoluted to extract the
mechanical properties of the chain at various length scales and thus can be
used to nanopore based methods to probe bio-molecules, such as DNA, RNA and
proteins.


  Statistical coupling analysis (SCA) is a method for analyzing multiple
sequence alignments that was used to identify groups of coevolving residues
termed "sectors". The method applies spectral analysis to a matrix obtained by
combining correlation information with sequence conservation. It has been
asserted that the protein sectors identified by SCA are functionally
significant, with different sectors controlling different biochemical
properties of the protein. Here we reconsider the available experimental data
and note that it involves almost exclusively proteins with a single sector. We
show that in this case sequence conservation is the dominating factor in SCA,
and can alone be used to make statistically equivalent functional predictions.
Therefore, we suggest shifting the experimental focus to proteins for which SCA
identifies several sectors. Correlations in protein alignments, which have been
shown to be informative in a number of independent studies, would then be less
dominated by sequence conservation.


  The nature of the phase transition in a system of self-propelling particles
has been extensively studied during the last few decades. A theoretical model
was proposed by T. Vicsek, {\it et. al.} [Phys. Rev. Lett. {\bf 75}, 1226
(1995)] with a simple rule for updating the direction of motion of each
particle. Based on the Vicsek's model (VM) [1], in this work, we consider a
group of animals as particles moving freely on a two-dimensional space. Due to
the fact that the viewable area of animals depends on the species, we consider
the motion of each individual within an angle $\varphi=\Phi/2$ ($\Phi$ is
called angle of view) of a circle centered at its position, of radius $R$. We
obtained a phase diagram in the space ($\varphi$, $\eta_c$) with $\eta_c$ being
the critical noise. We show that, the phase transition exists only in the case
of a wide view's angle $\varphi \geq 0.5\pi$. The flocking of animals is an
universal behavior of the species of prey, but not the one of the predator. Our
simulation results are in good agreement with experimental observation [2].


  Assignment of protein backbone resonances is most routinely carried out using
triple resonance three dimensional NMR experiments involving amide 1H and 15N
resonances. However for intrinsically unstructured proteins, alpha-helical
proteins or proteins containing several disordered fragments, the assignment
becomes problematic because of high degree of backbone shift degeneracy. In
this backdrop, a novel reduced dimensionality (RD) experiment
-(5,3)D-hNCO-CANH- is presented to facilitate (and/or to validate) the
sequential backbone resonance assignment in such proteins. The proposed 3D NMR
experiment makes use of the modulated amide 15N chemical shifts (resulting from
the joint sampling along both its indirect dimensions) to resolve the ambiguity
involved in connecting the neighboring amide resonances (i.e. HiNi and
Hi-1Ni-1) for overlapping amide NH peaks. The experiment -encoding 5D spectral
information- leads to a conventional 3D spectrum with significantly reduced
spectral crowding and complexity. The improvisation is based on the fact that
the linear combinations of intra-residue and inter-residue backbone chemical
shifts along both the co-evolved indirect dimensions span a wider spectral
range and produce better peak dispersion than the individual shifts themselves.
Taken together, the experiment -in combination with routine triple resonance 3D
NMR experiments involving backbone amide (1H and 15N) and carbon (13C-alpha and
13C') chemical shifts- will serve as a powerful complementary tool to achieve
the nearly complete assignment of protein backbone resonances in a time
efficient manner. The performance of the experiment and application of the
method have been demonstrated here using a 15.4 kDa size folded protein and a
12 kDa size unfolded protein.


  We show that a rate of conditional Shannon entropy reduction, characterizing
the learning of an internal process about an external process, is bounded by
the thermodynamic entropy production. This approach allows for the definition
of an informational efficiency that can be used to study cellular information
processing. We analyze three models of increasing complexity inspired by the E.
coli sensory network, where the external process is an external ligand
concentration jumping between two values. We start with a simple model for
which ATP must be consumed so that a protein inside the cell can learn about
the external concentration. With a second model for a single receptor we show
that the rate at which the receptor learns about the external environment can
be nonzero even without any dissipation inside the cell since chemical work
done by the external process compensates for this learning rate. The third
model is more complete, also containing adaptation. For this model we show
inter alia that a bacterium in an environment that changes at a very slow
time-scale is quite inefficient, dissipating much more than it learns. Using
the concept of a coarse-grained learning rate, we show for the model with
adaptation that while the activity learns about the external signal the option
of changing the methylation level increases the concentration range for which
the learning rate is substantial.