  Many biological characteristics of evolutionary interest are not scalar
variables but continuous functions. Given a dataset of function-valued traits
generated by evolution, we develop a practical statistical approach to infer
ancestral function-valued traits, and estimate the generative evolutionary
process. We do this by combining dimension reduction and phylogenetic Gaussian
process regression, a nonparametric procedure which explicitly accounts for
known phylogenetic relationships. We test the methods' performance on simulated
function-valued data generated from a stochastic evolutionary model. The
methods are applied assuming that only the phylogeny and the function-valued
traits of taxa at its tips are known. Our method is robust and applicable to a
wide range of function-valued data, and also offers a phylogenetically aware
method for estimating the autocorrelation of function-valued traits.


  The conditions that can lead to the exploitative depletion of a shared
resource, i.e, the tragedy of the commons, can be reformulated as a game of
prisoner's dilemma: while preserving the common resource is in the best
interest of the group, over-consumption is in the interest of each particular
individual at any given point in time. One way to try and prevent the tragedy
of the commons is through infliction of punishment for over-consumption and/or
encouraging under-consumption, thus selecting against over-consumers. Here, the
effectiveness of various punishment functions in an evolving consumer-resource
system is evaluated within a framework of a parametrically heterogeneous system
of ordinary differential equations (ODEs). Conditions leading to the
possibility of sustainable coexistence with the common resource for a subset of
cases are identified analytically using adaptive dynamics; the effects of
punishment on heterogeneous populations with different initial composition are
evaluated using the Reduction theorem for replicator equations. Obtained
results suggest that one cannot prevent the tragedy of the commons through
rewarding of under-consumers alone - there must also be an implementation of
some degree of punishment that increases in a non-linear fashion with respect
to over-consumption and which may vary depending on the initial distribution of
clones in the population.


  Gene innovation is a key mechanism on the evolution and phenotypic diversity
of life forms. There is a need for tools able to study gene innovation across
an increasingly large number of genomic sequences to maximally capitalise our
understanding of biological systems. Here we present
Comparative-Phylostratigraphy, an open-source software suite that enables to
time the emergence of new genes across evolutionary time and to correlate
patterns of gene emergence with species traits simultaneously across whole
genomes from multiple species. Such a comparative strategy is a new powerful
tool for starting to dissect the relationship between gene innovation and
phenotypic diversity. We describe and showcase our method by analysing recently
published ant genomes. This new methodology identified significant bouts of new
gene evolution in ant clades, that are associated with shifts in life-history
traits. Our method allows easy integration of new genomic data as it becomes
available, and thus will be a valuable analytical tool for evolutionary
biologists interested in explaining the evolution of diversity of life at the
level of the genes.


  Using a tensorial approach, we show how to construct a one-one correspondence
between pattern probabilities and edge parameters for any group-based model.
This is a generalisation of the "Hadamard conjugation" and is equivalent to
standard results that use Fourier analysis. In our derivation we focus on the
connections to group representation theory and emphasize that the inversion is
possible because, under their usual definition, group-based models are defined
for abelian groups only. We also argue that our approach is elementary in the
sense that it can be understood as simple matrix multiplication where matrices
are rectangular and indexed by ordered-partitions of varying sizes.


  One of the fundamental cellular processes governed by genetic regulatory
networks in cells is the transition among different states under the intrinsic
and extrinsic noise. Based on a two-state genetic switching model with positive
feedback, we develop a framework to understand the metastability in gene
expressions. This framework is comprised of identifying the transition path,
reconstructing the global quasi-potential energy landscape, analyzing the
uphill and downhill transition paths, etc. It is successfully utilized to
investigate the stability of genetic switching models and fluctuation
properties in different regimes of gene expression with positive feedback. The
quasi-potential energy landscape, which is the rationalized version of
Waddington potential, provides a quantitative tool to understand the
metastability in more general biological processes with intrinsic noise.


  A new algorithm, "HiER-leap", is derived which improves on the computational
properties of the ER-leap algorithm for exact accelerated simulation of
stochastic chemical kinetics. Unlike ER-leap, HiER-leap utilizes a hierarchical
or divide-and-conquer organization of reaction channels into tightly coupled
"blocks" and is thereby able to speed up systems with many reaction channels.
Like ER-leap, HiER-leap is based on the use of upper and lower bounds on the
reaction propensities to define a rejection sampling algorithm with inexpensive
early rejection and acceptance steps. But in HiER-leap, large portions of
intra-block sampling may be done in parallel. An accept/reject step is used to
synchronize across blocks. This method scales well when many reaction channels
are present and has desirable asymptotic properties. The algorithm is exact,
parallelizable and achieves a significant speedup over SSA and ER-leap on
certain problems. This algorithm offers a potentially important step towards
efficient in silico modeling of entire organisms.


  The thesis consists of two projects. In the first project, we present a
software that analyses RNA secondary structures and compares them. The goal of
this software is to find the differences between two secondary structures
(experimental or predicted) in order to improve or compare algorithms for
predicting secondary structures. Then, a comparison between secondary
structures predicted by the Vienna package to those found experimentally is
presented and cases in which there exists a difference between the prediction
and the experimental structure are identified. As the differences originate
mainly from faces and hydrogen bonds that are not allowed by the Vienna
package, it is suggested that prediction may be improved by integrating them
into the software. In the second project we calculate the free energy of an
interior loop using Monte-Carlo simulation. We first present a semi-coarse
grained model for interior loops of RNA, and the energy model for the different
interactions. We then introduce the Monte-Carlo simulations and the method of
Parallel Tempering which enables good sampling of configuration space by
simulating a system simultaneously at several temperatures. Next we present
Thermodynamic Integration, which is a method for calculating free energy
differences. Then, We introduce a method that calculates the free energy
significantly faster since we need to use only one parameter, $T$. To implement
this method, we had to reach a regime in which the partition functions of the
two systems are equal, which isn't satisfied for systems that have different
entropies in the high temperature limit, so a solution to this problem had to
be found. Free energy values calculated for various interior loops are shown,
and may, if verified or with a more realistic modelling, be integrated into the
Vienna package and supply an alternative to the experiments done.


  The contribution to an organism's phenotype from one genetic locus may depend
upon the status of other loci. Such epistatic interactions among loci are now
recognized as fundamental to shaping the process of adaptation in evolving
populations. Although little is known about the structure of epistasis in most
organisms, recent experiments with bacterial populations have concluded that
antagonistic interactions abound and tend to de-accelerate the pace of
adaptation over time. Here, we use a broad class of mathematical fitness
landscapes to examine how natural selection biases the mutations that
substitute during evolution based on their epistatic interactions. We find
that, even when beneficial mutations are rare, these biases are strong and
change substantially throughout the course of adaptation. In particular,
epistasis is less prevalent than the neutral expectation early in adaptation
and much more prevalent later, with a concomitant shift from predominantly
antagonistic interactions early in adaptation to synergistic and sign epistasis
later in adaptation. We observe the same patterns when re-analyzing data from a
recent microbial evolution experiment. Since these biases depend on the
population size and other parameters, they must be quantified before we can
hope to use experimental data to infer an organism's underlying fitness
landscape or to understand the role of epistasis in shaping its adaptation. In
particular, we show that when the order of substitutions is not known to an
experimentalist, then standard methods of analysis may suggest that epistasis
retards adaptation when in fact it accelerates it.


  The Genographic Project is an international effort using genetic data to
chart human migratory history. The project is non-profit and non-medical, and
through its Legacy Fund supports locally led efforts to preserve indigenous and
traditional cultures. In its second phase, the project is focusing on markers
from across the entire genome to obtain a more complete understanding of human
genetic variation. Although many commercial arrays exist for genome-wide SNP
genotyping, they were designed for medical genetic studies and contain
medically related markers that are not appropriate for global population
genetic studies. GenoChip, the Genographic Project's new genotyping array, was
designed to resolve these issues and enable higher-resolution research into
outstanding questions in genetic anthropology. We developed novel methods to
identify AIMs and genomic regions that may be enriched with alleles shared with
ancestral hominins. Overall, we collected and ascertained AIMs from over 450
populations. Containing an unprecedented number of Y-chromosomal and mtDNA SNPs
and over 130,000 SNPs from the autosomes and X-chromosome, the chip was
carefully vetted to avoid inclusion of medically relevant markers. The GenoChip
results were successfully validated. To demonstrate its capabilities, we
compared the FST distributions of GenoChip SNPs to those of two commercial
arrays for three continental populations. While all arrays yielded similarly
shaped (inverse J) FST distributions, the GenoChip autosomal and X-chromosomal
distributions had the highest mean FST, attesting to its ability to discern
subpopulations. The GenoChip is a dedicated genotyping platform for genetic
anthropology and promises to be the most powerful tool available for assessing
population structure and migration history.


  High-throughput shotgun sequence data makes it possible in principle to
accurately estimate population genetic parameters without confounding by SNP
ascertainment bias. One such statistic of interest is the proportion of
heterozygous sites within an individual's genome, which is informative about
inbreeding and effective population size. However, in many cases, the available
sequence data of an individual is limited to low coverage, preventing the
confident calling of genotypes necessary to directly count the proportion of
heterozygous sites. Here, we present a method for estimating an individual's
genome-wide rate of heterozygosity from low-coverage sequence data, without an
intermediate step calling genotypes. Our method jointly learns the shared
allele distribution between the individual and a panel of other individuals,
together with the sequencing error distributions and the reference bias. We
show our method works well, first by its performance on simulated sequence
data, and secondly on real sequence data where we obtain estimates using low
coverage data consistent with those from higher coverage. We apply our method
to obtain estimates of the rate of heterozygosity for 11 humans from diverse
world-wide populations, and through this analysis reveal the complex dependency
of local sequencing coverage on the true underlying heterozygosity, which
complicates the estimation of heterozygosity from sequence data. We show
filters can correct for the confounding by sequencing depth. We find in
practice that ratios of heterozygosity are more interpretable than absolute
estimates, and show that we obtain excellent conformity of ratios of
heterozygosity with previous estimates from higher coverage data.


  The growing interest for comparing protein internal dynamics owes much to the
realization that protein function can be accompanied or assisted by structural
fluctuations and conformational changes. Analogously to the case of functional
structural elements, those aspects of protein flexibility and dynamics that are
functionally oriented should be subject to evolutionary conservation.
Accordingly, dynamics-based protein comparisons or alignments could be used to
detect protein relationships that are more elusive to sequence and structural
alignments. Here we provide an account of the progress that has been made in
recent years towards developing and applying general methods for comparing
proteins in terms of their internal dynamics and advance the understanding of
the structure-function relationship.


  Multi-agent systems have been studied in various contexts of both application
and theory. We take Dynamic Epistemic Logic (DEL), one of the formalisms
designed to reason about such systems, as the foundation of the language we
will build.
  BioAmbient calculus is an extension of \pi-calculus, developed largely for
applications to biomolecular systems. It deals with ambients and their ability
to communicate and to execute concurrent processes while moving.
  In this paper we combine the formalism of Dynamic Epistemic Logic together
with the formalism of BioAmbient Calculus in order to reason about knowledge
maintained and gained upon process transitions. The motivation lies in
developing a language that captures locally available information through
assignment of knowledge, with potential application to biological systems as
well as social, virtual, and others.
  We replace the ambients of BioAmbient Calculus with agents, to which we
attribute knowledge, and explore the parallels of this treatment. The resulting
logic describes the information flow governing mobile structured agents,
organized hierarchically, whose architecture (and local information) may change
due to actions such as communication, merging (of two agents), entering (of an
agent into the inner structure of another agent) and exiting (of an agent from
the structure of another). We show how the main axioms of DEL must be altered
to accommodate the informational effects of the agents' dynamic architecture.


  Despite numerous mass extinctions in the Phanerozoic eon, the overall trend
in biodiversity evolution was not blocked and the life has never been wiped
out. Almost all possible catastrophic events (large igneous province, asteroid
impact, climate change, regression and transgression, anoxia, acidification,
sudden release of methane clathrate, multi-cause etc.) have been proposed to
explain the mass extinctions. However, we should, above all, clarify at what
timescale and at what possible levels should we explain the mass extinction?
Even though the mass extinctions occurred at short-timescale and at the species
level, we reveal that their cause should be explained in a broader context at
tectonic timescale and at both the molecular level and the species level. The
main result in this paper is that the Phanerozoic biodiversity evolution has
been explained by reconstructing the Sepkoski curve based on climatic, eustatic
and genomic data. Consequently, we point out that the P-Tr extinction was
caused by the tectonically originated climate instability. We also clarify that
the overall trend of biodiversification originated from the underlying genome
size evolution, and that the fluctuation of biodiversity originated from the
interactions among the earth's spheres. The evolution at molecular level had
played a significant role for the survival of life from environmental
disasters.


  We propose the first, to our knowledge, coarse-grained modeling strategy for
peptides where the effect of changes of the pH can be efficiently described.
The idea is based on modeling the effects of the pH value on the main driving
interactions. We use reference data from atomistic simulations and experimental
databases and transfer its main physical features to the coarse-grained
resolution according the principle of "consistency across the scales". The
coarse-grained model is refined by finding a set of parameters that, when
applied to peptides with different sequences and experimental properties,
reproduces the experimental and atomistic data of reference. We use the such
parameterized model for performing several numerical tests to check its
transferability to other systems and to prove the universality of the related
modeling strategy. We have tried systems with rather different response to pH
variations, showing a highly satisfactory performance of the model.


  Rubella is a completely immunizing and mild infection in children.
Understanding its behavior is of considerable public health importance because
of Congenital Rubella Syndrome, which results from infection with rubella
during early pregnancy and may entail a variety of birth defects. The recurrent
dynamics of rubella are relatively poorly resolved, and appear to show
considerable diversity globally. Here, we investigate the behavior of a
stochastic seasonally forced susceptible-infected-recovered model to
characterize the determinants of these dynamics and illustrate patterns by
comparison with measles. We perform a systematic analysis of spectra of
stochastic fluctuations around stable attractors of the corresponding
deterministic model and compare them with spectra from full stochastic
simulations in large populations. This approach allows us to quantify the
effects of demographic stochasticity and to give a coherent picture of measles
and rubella dynamics, explaining essential differences in the recurrent
patterns exhibited by these diseases. We discuss the implications of our
findings in the context of vaccination and changing birth rates as well as the
persistence of these two childhood infections.


  Self-sustaining autocatalytic chemical networks represent a necessary, though
not sufficient condition for the emergence of early living systems. These
networks have been formalised and investigated within the framework of RAF
theory, which has led to a number of insights and results concerning the
likelihood of such networks forming. In this paper, we extend this analysis by
focussing on how {\em small} autocatalytic networks are likely to be when they
first emerge. First we show that simulations are unlikely to settle this
question, by establishing that the problem of finding a smallest RAF within a
catalytic reaction system is NP-hard. However, irreducible RAFs (irrRAFs) can
be constructed in polynomial time, and we show it is possible to determine in
polynomial time whether a bounded size set of these irrRAFs contain the
smallest RAFs within a system. Moreover, we derive rigorous bounds on the sizes
of small RAFs and use simulations to sample irrRAFs under the binary polymer
model. We then apply mathematical arguments to prove a new result suggested by
those simulations: at the transition catalysis level at which RAFs first form
in this model, small RAFs are unlikely to be present. We also investigate
further the relationship between RAFs and another formal approach to
self-sustaining and closed chemical networks, namely chemical organisation
theory (COT).


  We report the experimental verification of noise-enhanced logic behaviour in
an electronic analog of a synthetic genetic network, composed of two repressors
and two constitutive promoters. We observe good agreement between circuit
measurements and numerical prediction, with the circuit allowing for robust
logic operations in an optimal window of noise. Namely, the input-output
characteristics of a logic gate is reproduced faithfully under moderate noise,
which is a manifestation of the phenomenon known as Logical Stochastic
Resonance. The two dynamical variables in the system yield complementary logic
behaviour simultaneously. The system is easily morphed from AND/NAND to OR/NOR
logic.


  Modeling of ion transport via plasma membrane needs identification and
quantitative understanding of the involved processes. Brief characterization of
main ion transport systems of a yeast cell (Pma1, Ena1, TOK1, Nha1, Trk1, Trk2,
non-selective cation conductance) and determining the exact number of molecules
of each transporter per a typical cell allow us to predict the corresponding
ion flows. In this review a comparison of ion transport in small yeast cell and
several animal cell types is provided. The importance of cell volume to surface
ratio is emphasized. The role of cell wall and lipid rafts is discussed in
respect to required increase in spatial and temporal resolution of
measurements. Conclusions are formulated to describe specific features of ion
transport in a yeast cell. Potential directions of future research are outlined
based on the assumptions.


  We apply a recently-developed coarse-grained model of DNA, designed to
capture the basic physics of nanotechnological DNA systems, to the study of a
`burnt-bridges' DNA motor consisting of a single-stranded cargo that steps
processively along a track of single-stranded stators. We demonstrate that the
model is able to simulate such a system, and investigate the sensitivity of the
stepping process to the spatial separation of stators, finding that an
increased distance can suppress successful steps due to the build up of
unfavourable tension. The mechanism of suppression suggests that varying the
distance between stators could be used as a method for improving
signal-to-noise ratios for motors that are required to make a decision at a
junction of stators.


  In the study of amyloid fiber formation, the fiber length distribution plays
a central role. From this microscopic distribution, all macroscopic statistical
quantities, like the fiber number concentration P(t) and mass concentration
M(t), can be easily achieved. However the inverse problem is usually
non-trivial. In this paper, a practical method is introduced to show how to
reconstruct the fiber length distribution function based on the knowledge of
P(t) and M(t). Compared to directly solving mass-action equations, our method
speeds up the calculation by at least ten thousand times (from days to
seconds). And the accuracy is also quite satisfactory. Furthermore our method
provides an interesting linkage between the fluorescence imaging and more
detailed TEM studies, which may inspire wider applications in this field.


  We propose a method of estimating the Total Loss of Healthy Life Years based
on the first exit time theory for a stochastic process, the resulting Health
State Function and the Deterioration Function estimated as the curvature of the
health state function. We have done many applications in UK and Scotland and
Sweden supporting our theory. Furthermore it was proven that both the WHO and
EU estimates of the healthy life expectancy can result from our method. The WHO
system takes into account the severe and moderate causes in estimating the loss
of healthy life years; instead the EU system calculates the total loss of
healthy life years. For both cases our methodology provides both estimators
from only death and population data. The advantages of our method are
straightforward. We do not need survey data to make the calculations. The
resulting estimates should be used to test and improve the existing survey
based methodologies. While the WHO and EU systems tend to approach each other
differences continue to appear based on the methodology of the related surveys
and the analysis of data. Two main schools are working to these directions one
based on USA the Institute for Health Metrics and Evaluation (IHME) headed by
Christopher J.L. Murray and contributors in all over the world and the European
Health and Life Expectancy Information System (EHLEIS) with Jean-Marie Robine
and a team from the EU member states. Keywords: Deterioration, Loss of healthy
years, HALE, DALE, World Health Organization, WHO, European Union, EU, EHEMU,
IHME, EHLEIS, Healthy life expectancy, Life expectancy.


  We present results of a forecast initiated Week 49 (beginning December 9,
2012) of the 2012-2013 influenza season for municipalities in the United
States. The forecast was made on December 14, 2012. Results from forecasts
initiated the two previous weeks (Weeks 47 and 48) are also presented. Also
results from the forecast generated with the SIRS model without AH forcing (no
AH) are shown


  Recent experimental data indicate that the elastic wormlike rod model of DNA
that works well on long length scales may break down on shorter scales relevant
to biology. According to Noy and Golestanian (Phys. Rev. Lett. 109, 228101,
2012) molecular dynamics (MD) simulations predict DNA rigidity close to
experimental data and confirm one scenario of such breakdown, namely, that for
lengths of a few helical turns, DNA dynamics exhibit long-range bending and
stretching correlations. Earlier studies using similar forcefields concluded
that (i) MD systematically overestimate the DNA rigidity, and (ii) no
deviations from the WLR model are detectable. Here it is argued that the data
analysis in the above mentioned paper was incorrect and that the earlier
conclusions are valid.


  We present a minimal motif model for transmembrane cell signaling. The model
assumes signaling events taking place in spatially distributed nanoclusters
regulated by a birth/death dynamics. The combination of these spatio-temporal
aspects can be modulated to provide a robust and high-fidelity response
behavior without invoking sophisticated modeling of the signaling process as a
sequence of cascade reactions and fine-tuned parameters. Our results show that
the fact that the distributed signaling events take place in nanoclusters with
a finite lifetime regulated by local production is sufficient to obtain a
robust and high-fidelity response.


  Mapping expression Quantitative Trait Loci (eQTLs) represents a powerful and
widely-adopted approach to identifying putative regulatory variants and linking
them to specific genes. Up to now eQTL studies have been conducted in a
relatively narrow range of tissues or cell types. However, understanding the
biology of organismal phenotypes will involve understanding regulation in
multiple tissues, and ongoing studies are collecting eQTL data in dozens of
cell types. Here we present a statistical framework for powerfully detecting
eQTLs in multiple tissues or cell types (or, more generally, multiple
subgroups). The framework explicitly models the potential for each eQTL to be
active in some tissues and inactive in others. By modeling the sharing of
active eQTLs among tissues this framework increases power to detect eQTLs that
are present in more than one tissue compared with "tissue-by-tissue" analyses
that examine each tissue separately. Conversely, by modeling the inactivity of
eQTLs in some tissues, the framework allows the proportion of eQTLs shared
across different tissues to be formally estimated as parameters of a model,
addressing the difficulties of accounting for incomplete power when comparing
overlaps of eQTLs identified by tissue-by-tissue analyses. Applying our
framework to re-analyze data from transformed B cells, T cells and fibroblasts
we find that it substantially increases power compared with tissue-by-tissue
analysis, identifying 63% more genes with eQTLs (at FDR=0.05). Further the
results suggest that, in contrast to previous analyses of the same data, the
majority of eQTLs detectable in these data are shared among all three tissues.


  Motivation: The rapid growth in genome-wide association studies (GWAS) in
plants and animals has brought about the need for a central resource that
facilitates i) performing GWAS, ii) accessing data and results of other GWAS,
and iii) enabling all users regardless of their background to exploit the
latest statistical techniques without having to manage complex software and
computing resources.
  Results: We present easyGWAS, a web platform that provides methods, tools and
dynamic visualizations to perform and analyze GWAS. In addition, easyGWAS makes
it simple to reproduce results of others, validate findings, and access larger
sample sizes through merging of public datasets.
  Availability: Detailed method and data descriptions as well as tutorials are
available in the supplementary materials. easyGWAS is available at
http://easygwas.tuebingen.mpg.de/.
  Contact: dominik.grimm@tuebingen.mpg.de


  Cued and un-cued cell polarization is a fundamental mechanism in cell
biology. As an alternative to the classical Turing bifurcation, it has been
proposed that the cell polarity might onset by means of the well-known
phenomenon of wave-pinning (Gamba et al, PNAS, 2005). A particularly simple and
elegant model of wave-pinning has been proposed by Edelstein-Keshet and
coworkers (Biop. J., 2008). However, biomolecular networks do communicate with
other networks as well as with the external world. As such, their dynamics has
to be considered as perturbed by extrinsic noises. These noises may have both a
spatial and a temporal correlation, but any case they must be bounded to
preserve the biological meaningfulness of the perturbed parameters. Here we
numerically show that the inclusion of external spatio-temporal bounded
perturbations may sometime destroy the polarized state. The polarization loss
depends on both the extent of temporal and spatial correlations, and on the
kind of adopted noise. Namely, independently of the specific model of noise, an
increase of the spatial correlation induces an increase of the probability of
polarization. However, if the noise is spatially homogeneous then the
polarization is lost in the majority of cases. On the contrary, an increase of
the temporal autocorrelation of the noise induces an effect that depends on the
noise model.


  Mitochondrial outer membrane permeabilization (MOMP) is one of the most
important points, in majority of apoptotic signaling cascades. Decision
mechanism controlling whether the MOMP occurs or not, is formed by an interplay
between members of the Bcl-2 family. To understand the role of individual
members of this family within the MOMP regulation, we constructed a boolean
network-based mathematical model of interactions between the Bcl-2 proteins.
Results of computational simulations reveal the existence of the potentially
malign configurations of activities of the Bcl-2 proteins, blocking the
occurrence of MOMP, independently of the incoming stimuli. Our results suggest
role of the antiapoptotic protein Mcl-1 in relation to these configurations. We
demonstrate here, the importance of the Bid and Bim according to activation of
effectors Bax and Bak, and the irreversibility of this activation. The model
further shows the distinct requirements for effectors activation, where the
antiapoptic protein Bcl-w is seemingly a key factor preventing the Bax
activation. We believe that this work may help to describe the functioning of
the Bcl-2 regulation of MOMP better, and hopefully provide some contribution
regarding the anti-cancer drug development research.


  Metacommunity theory is considered a promising approach for explaining
species diversity and food web complexity. Recently Pillai et al. proposed a
simple modeling framework for the dynamics of food webs at the metacommunity
level. Here, we employ this framework to compute general conditions for the
persistence of complex food webs in metacommunities. The persistence conditions
found depend on the connectivity of the resource patches and the structure of
the assembled food web, thus linking the underlying spatial patch-network and
the species interaction network. We find that the persistence of omnivores is
more likely when it is feeding on (a) prey on low trophic levels, and (b) prey
on similar trophic levels.


  Construction of synthetic genetic networks requires the assembly of DNA
fragments encoding functional biological parts in a defined order. Yet this may
become a time-consuming procedure. To address this technical bottleneck, we
have created a series of Gateway shuttle vectors and an integration vector,
which facilitate the assembly of artificial genes and their expression in the
budding yeast Saccharomyces cerevisiae. Our method enables the rapid
construction of an artificial gene from a promoter and an open reading frame
(ORF) cassette by one-step recombination reaction in vitro. Furthermore, the
plasmid thus created can readily be introduced into yeast cells to test the
assembled gene's functionality. As flexible regulatory components of a
synthetic genetic network, we also created new versions of the
tetracycline-regulated transactivators tTA and rtTA by fusing them to the
auxin-inducible degron (AID). Using our gene assembly approach, we made yeast
expression vectors of these engineered transactivators, AIDtTA and AIDrtTA and
then tested their functions in yeast. We showed that these factors can be
regulated by doxycycline and degraded rapidly after addition of auxin to the
medium. Taken together, the method for combinatorial gene assembly described
here is versatile and would be a valuable tool for yeast synthetic biology.


  An important question in molecular evolution is whether an amino acid that
occurs at a given position makes an independent contribution to fitness, or
whether its effect depends on the state of other loci in the organism's genome,
a phenomenon known as epistasis. In a recent letter to Nature, Breen et al.
(2012) argued that epistasis must be "pervasive throughout protein evolution"
because the observed ratio between the per-site rates of non-synonymous and
synonymous substitutions (dN/dS) is much lower than would be expected in the
absence of epistasis. However, when calculating the expected dN/dS ratio in the
absence of epistasis, Breen et al. assumed that all amino acids observed in a
protein alignment at any particular position have equal fitness. Here, we relax
this unrealistic assumption and show that any dN/dS value can in principle be
achieved at a site, without epistasis. Furthermore, for all nuclear and
chloroplast genes in the Breen et al. dataset, we show that the observed dN/dS
values and the observed patterns of amino acid diversity at each site are
jointly consistent with a non-epistatic model of protein evolution.


  In recent work discussing model choice for continuous-time Markov chains, we
have argued that it is important that the Markov matrices that define the model
are closed under matrix multiplication (Sumner 2012a, 2012b). The primary
requirement is then that the associated set of rate matrices form a Lie
algebra. For the generic case, this connection to Lie theory seems to have
first been made by Johnson (1985), with applications for specific models given
in Bashford (2004) and House (2012). Here we take a different perspective:
given a model that forms a Lie algebra, we apply existing Lie theory to gain
additional insight into the geometry of the associated Markov matrices. In this
short note, we present the simplest case possible of 2x2 Markov matrices. The
main result is a novel decomposition of 2x2 Markov matrices that parameterises
the general Markov model as a perturbation away from the binary-symmetric
model. This alternative parameterisation provides a useful tool for visualising
the binary-symmetric model as a submodel of the general Markov model.


  Recently, a large number of studies have been carried out on the early
signatures of sudden regime shifts in systems as diverse as ecosystems,
financial markets, population biology and complex diseases. Signatures of
regime shifts in gene expression dynamics are less systematically investigated.
In this paper, we consider sudden regime shifts in the gene expression dynamics
described by a fold-bifurcation model involving bistability and hysteresis. We
consider two alternative models, Models 1 and 2, of competence development in
the bacterial population B. subtilis and determine some early signatures of the
regime shifts between competence and noncompetence. We use both deterministic
and stochastic formalisms for the purpose of our study. The early signatures
studied include the critical slowing down as a transition point is approached,
rising variance and the lag-1 autocorrelation function, skewness and a ratio of
two mean first passage times. Some of the signatures could provide the
experimental basis for distinguishing between bistability and excitability as
the correct mechanism for the development of competence.


  Diverse complex dynamical systems are known to exhibit abrupt regime shifts
at bifurcation points of the saddle-node type. The dynamics of most of these
systems, however, have a stochastic component resulting in noise driven regime
shifts even if the system is away from the bifurcation points. In this paper,
we propose a new quantitative measure, namely, the propensity transition point
as an indicator of stochastic regime shifts. The concepts and the methodology
are illustrated for the one-variable May model, a well-known model in ecology
and the genetic toggle, a two-variable model of a simple genetic circuit. The
general applicability and usefulness of the method for the analysis of regime
shifts is further demonstrated in the case of the mycobacterial switch to
persistence for which experimental data are available.


  Phytoplankton are tiny floating plants (algae) living in oceans. In the
process of photosynthesis, phytoplankton produces half of the world's oxygen.
Moreover, by primary production, death and sinking, they transport carbon from
the ocean's surface layer to marine sediments. There are many species of
phytoplankton that can be distinguished according to morphology.
  In this paper, we investigate the generalised Truscott-Brindley model of the
dynamics of zoologically defined interacting populations which have spatial
structure. Specifically, we consider conjointly marine phytoplankton and
zooplankton populations, and model them as an excitable medium. The resolution
using the Boubaker polynomials expansion scheme (BPES) along with stability
analysis is carried out.


  Phylogenetic studies show that horizontal gene transfer (HGT) is a
significant contributor to genetic variability of prokaryotes, and was perhaps
even more abundant during early evolution. Hitherto, research of natural HGT
has mainly focused on three mechanisms: conjugation, natural competence, and
viral transduction. This paper discusses the feasibility of a fourth such
mechanism - cell electroporation and/or electrofusion triggered by atmospheric
electrostatic discharges (lightnings). A description of electroporation as a
phenomenon is followed by a review of experimental evidence that
electroporation of prokaryotes in aqueous environments can result in release of
non-denatured DNA, as well as uptake of DNA from the surroundings and
transformation. Similarly, a description of electrofusion is followed by a
review of experiments showing that prokaryotes devoid of cell wall can
electrofuse into hybrids expressing the genes of their both precursors. Under
sufficiently fine-tuned conditions, electroporation and electrofusion are
efficient tools for artificial transformation and hybridization, respectively,
but the quantitative analysis developed here shows that conditions for
electroporation-based DNA release, DNA uptake and transformation, as well as
for electrofusion are also present in many natural aqueous environments exposed
to lightnings. Electroporation is thus a plausible contributor to natural HGT
among prokaryotes, and could have been particularly important during the early
evolution, when the other mechanisms might have been scarcer or nonexistent. In
modern prokaryotes, natural absence of the cell wall is rare, but it is
reasonable to assume that the wall has formed during evolution, and at least
prior to this, electrofusion could also have contributed to HGT. The concluding
section outlines several guidelines for assessment of the feasibility of
lightning-triggered HGT.


  Gynandromorphs are creatures where at least two different body sections are a
different sex. Bilateral gynandromorphs are half male and half female. Here we
develop a theory of gynandromorph ontogeny based on developmental control
networks. The theory explains the embryogenesis of all known variations of
gynandromorphs found in multicellular organisms. The theory also predicts a
large variety of more subtle gynandromorphic morphologies yet to be discovered.
The network theory of gynandromorph development has direct relevance to
understanding sexual dimorphism (differences in morphology between male and
female organisms of the same species) and medical pathologies such as
hemihyperplasia (asymmetric development of normally symmetric body parts in a
unisexual individual). The network theory of gynandromorphs brings up
fundamental open questions about developmental control in ontogeny. This in
turn suggests a new theory of the origin and evolution of species that is based
on cooperative interactions and conflicts between developmental control
networks in the haploid genomes and epigenomes of potential sexual partners for
reproduction. This network-based theory of the origin of species is a
paradigmatic shift in our understanding of evolutionary processes that goes
beyond gene-centered theories.


  A spatial three level food web model with a closed nutrient cycle is
presented and analyzed via Monte Carlo simulations. The time evolution of the
model reveals two asymptotic states: an absorbing one with all species being
extinct, and a coexisting one, in which concentrations of all species are
non-zero. There are two possible ways for the system to reach the absorbing
state. In some cases the densities increase very quickly at the beginning of a
simulation and then decline slowly and almost monotonically. In others, well
pronounced peaks in the $R$, $C$ and $D$ densities appear regularly before the
extinction. Those peaks correspond to density outbursts (waves) traveling
through the system. We investigate the mechanisms leading to the waves. In
particular, we show that the percolation of the detritus (i. e. the
accumulation of nutrients) is necessary for the emergence of the waves.
Moreover, our results corroborate the hypothesis that top-level predators play
an essential role in maintaining the stability of a food web (top-down
control).


  Various biological sensory systems exhibit a response to a relative change of
the stimulus, often referred to as fold-change detection. In the last few years
fold-change detecting mechanisms, based on transcriptional networks, have been
proposed. Here we present fold-change detecting mechanism, based on
protein-protein interactions, consisting of two interacting proteins. This
mechanism, in contrast to previously proposed mechanisms, does not consume
chemical energy and is not subject to transcriptional and translational noise.
We show by analytical and numerical calculations, that the mechanism can have a
fast, precise and efficient response for parameters that are relevant to
eukaryotic cells.


  We present results of a forecast initiated following assimilation of
observations for week Week 50 (i.e. the forecast begins December 16, 2012) of
the 2012-2013 influenza season for municipalities in the United States. The
forecast was made on December 21, 2012. Results from forecasts initiated the
three previous weeks (Weeks 47-49) are also presented. Also results from
forecasts generated with an SIRS model without absolute humidity forcing (no
AH) are shown.


  During the application of mass-action equation models to the study of amyloid
fiber formation, time-consuming numerical calculations constitute a major
bottleneck when no analytical solution is available. To conquer this
difficulty, here an alternative efficient method is introduced for the
fragmentation-only model. It includes two basic steps: (1) simulate
close-formed time-evolution equations for the number concentration P(t) derived
from the moment-closure method; (2) reconstruct the detailed fiber length
distribution based on the knowledge of moments obtained in the first step.
Compared to direct solution, current method speeds up the calculation by at
least ten thousand times. The accuracy is also quite satastifactory if suitable
forms of approximate distribution fucntion is taken. Further application to
PI264-b-PFS48 micelles study performed by Guerin et al. confirms our method is
very promising for the real-time analysis of the experimental data on fiber
fragmentation.


  Accumulation of standardized data collections is opening up novel
opportunities for holistic characterization of genome function. The limited
scalability of current preprocessing techniques has, however, formed a
bottleneck for full utilization of contemporary microarray collections. While
short oligonucleotide arrays constitute a major source of genome-wide profiling
data, scalable probe-level preprocessing algorithms have been available only
for few measurement platforms based on pre-calculated model parameters from
restricted reference training sets. To overcome these key limitations, we
introduce a fully scalable online-learning algorithm that provides tools to
process large microarray atlases including tens of thousands of arrays. Unlike
the alternatives, the proposed algorithm scales up in linear time with respect
to sample size and is readily applicable to all short oligonucleotide
platforms. This is the only available preprocessing algorithm that can learn
probe-level parameters based on sequential hyperparameter updates at small,
consecutive batches of data, thus circumventing the extensive memory
requirements of the standard approaches and opening up novel opportunities to
take full advantage of contemporary microarray data collections. Moreover,
using the most comprehensive data collections to estimate probe-level effects
can assist in pinpointing individual probes affected by various biases and
provide new tools to guide array design and quality control. The implementation
is freely available in R/Bioconductor at
http://www.bioconductor.org/packages/devel/bioc/html/RPA.html


  The effect of a finite torque on semiflexible polymers in a confined
environment is investigated. It is shown how a new lengthscale appears in the
strongly confined limit. The influence of a torque on the extension of
biopolymers in nanochannels is also touched upon and it is argued that the
presence of a torque has a strong influence on the dimensions of nanochannels
needed to prevent hairpins.


  This paper is concerned with conditions for the existence of oscillations in
gene regulatory networks with negative cyclic feedback, where time delays in
transcription, translation and translocation process are explicitly considered.
The primary goal of this paper is to propose systematic analysis tools that are
useful for a broad class of cyclic gene regulatory networks, and to provide
novel biological insights. To this end, we adopt a simplified model that is
suitable for capturing the essence of a large class of gene regulatory
networks. It is first shown that local instability of the unique equilibrium
state results in oscillations based on a Poincare-Bendixson type theorem. Then,
a graphical existence condition, which is equivalent to the local instability
of a unique equilibrium, is derived. Based on the graphical condition, the
existence condition is analytically presented in terms of biochemical
parameters. This allows us to find the dimensionless parameters that primarily
affect the existence of oscillations, and to provide biological insights. The
analytic conditions and biological insights are illustrated with two existing
biochemical networks, Repressilator and the Hes7 gene regulatory networks.


  We study the Susceptible-Infected-Recovered model in complex networks,
considering that not all individuals in the population interact in the same way
between them. This heterogeneity between contacts is modeled by a continuous
disorder. In our model the disorder represents the contact time or the
closeness between individuals. We find that the duration time of an epidemic
has a crossover with the system size, from a power law regime to a logarithmic
regime depending on the transmissibility related to the strength of the
disorder. Using percolation theory, we find that the duration of the epidemic
scales as the average length of the branches of the infection. Our theoretical
findings, supported by simulations, explains the crossover between the two
regimes.


  This document is part of a series of near real-time weekly influenza
forecasts made during the 2012-2013 influenza season. Here we present results
of a forecast initiated following assimilation of observations for Week 51
(i.e. the forecast begins December 23, 2012) for municipalities in the United
States. The forecast was made on December 28, 2012. Results from forecasts
initiated the four previous weeks (Weeks 47-50) are also presented. Predictions
generated with an alternate SIRS model, run without absolute humidity forcing
(no AH), are also presented.


  The purpose of this article is to look at how information criteria, such as
AIC and BIC, relate to the g%SD fit criterion derived in Waddell et al. (2007,
2010a). The g%SD criterion measures the fit of data to model based on a
normalized weighted root mean square percentage deviation between the observed
data and model estimates of the data, with g%SD = 0 being a perfectly fitting
model. However, this criterion may not be adjusting for the number of
parameters in the model comprehensively. Thus, its relationship to more
traditional measures for maximizing useful information in a model, including
AIC and BIC, are examined. This results in an extended set of fit criteria
including g%AIC and g%BIC. Further, a broader range of asymptotically most
powerful fit criteria of the power divergence family, which includes maximum
likelihood (or minimum G^2) and minimum X^2 modeling as special cases, are used
to replace the sum of squares fit criterion within the g%SD criterion. Results
are illustrated with a set of genetic distances looking particularly at a range
of Jewish populations, plus a genomic data set that looks at how Neanderthals
and Denisovans are related to each other and modern humans. Evidence that Homo
erectus may have left a significant fraction of its genome within the Denisovan
is shown to persist with the new modeling criteria.


  Relevant events in a three state illness-death model (IDM) of a chronic
disease are the diagnosis of the disease and death with or without the disease.
In this article a simulation framework for populations moving in the IDM is
presented. The simulation is closely related to the concept of Lexis diagrams
in event history analysis. Details of the implementation and an example of a
hypothetical disease are described.


  This paper shows that differentiating the lifetimes of two phenotypes
independently from their fertility can lead to a qualitative change in the
equilibrium of a population: since survival and reproduction are distinct
functional aspects of an organism, this observation contributes to extend the
population-genetical characterisation of biological function. To support this
statement a mathematical relation is derived to link the lifetime ratio
T_1/T_2, which parametrizes the different survival ability of two phenotypes,
with population variables that quantify the amount of neutral variation
underlying a population's phenotypic distribution.


  We present a framework for the design of optimal assembly algorithms for
shotgun sequencing under the criterion of complete reconstruction. We derive a
lower bound on the read length and the coverage depth required for
reconstruction in terms of the repeat statistics of the genome. Building on
earlier works, we design a de Brujin graph based assembly algorithm which can
achieve very close to the lower bound for repeat statistics of a wide range of
sequenced genomes, including the GAGE datasets. The results are based on a set
of necessary and sufficient conditions on the DNA sequence and the reads for
reconstruction. The conditions can be viewed as the shotgun sequencing analogue
of Ukkonen-Pevzner's necessary and sufficient conditions for Sequencing by
Hybridization.


  Prion diseases caused by the conversion from a soluble normal cellular prion
protein into insoluble abnormally folded infectious prions, are invariably
fatal and highly infectious degenerative diseases that affect a wide variety of
mammalian species. The palindrome and the Glycine-rich conserved segment in the
hydrophobic region 109-136 control the conversion from normal prion protein to
form into diseased prions. This paper gives detailed reviews on the 109-136
region and presents the studies of its 3D structures and structural dynamics.


  This paper is concerned with stochastic reaction-diffusion kinetics governed
by the reaction-diffusion master equation. Specifically, the primary goal of
this paper is to provide a mechanistic basis of Turing pattern formation that
is induced by intrinsic noise. To this end, we first derive an approximate
reaction-diffusion system by using linear noise approximation. We show that the
approximated system has a certain structure that is associated with a coupled
dynamic multi-agent system. This observation then helps us derive an efficient
computation tool to examine the spatial power spectrum of the intrinsic noise.
We numerically demonstrate that the result is quite effective to analyze
noise-induced Turing pattern. Finally, we illustrate the theoretical mechanism
behind the noise-induced pattern formation with a H2 norm interpretation of the
multi-agent system.


  How to determine the community structure of complex networks is an open
question. It is critical to establish the best strategies for community
detection in networks of unknown structure. Here, using standard synthetic
benchmarks, we show that none of the algorithms hitherto developed for
community structure characterization perform optimally. Significantly,
evaluating the results according to their modularity, the most popular measure
of the quality of a partition, systematically provides mistaken solutions.
However, a novel quality function, called Surprise, can be used to elucidate
which is the optimal division into communities. Consequently, we show that the
best strategy to find the community structure of all the networks examined
involves choosing among the solutions provided by multiple algorithms the one
with the highest Surprise value. We conclude that Surprise maximization
precisely reveals the community structure of complex networks.


  We have studied an evolutionary game with spatially arranged players who can
choose one of the two strategies (named cooperation and defection for social
dilemmas) when playing with their neighbors. In addition to the application of
the usual strategies in the present model the players are also characterized by
one of the two extreme personal features representing the egoist or fraternal
behavior. During the evolution each player can modify both her own strategy
and/or personal feature via a myopic update process in order to improve her
utility. The results of numerical simulations and stability analysis are
summarized in phase diagrams representing a wide scale of spatially ordered
distribution of strategies and personal features when varying the payoff
parameters. In most of the cases only two of the four possible options prevail
and may form sublattice ordered spatial structure. The evolutionary advantage
of the fraternal attitude is demonstrated within a large range of payoff
parameters including the region of prisoner's dilemma where egoist defectors
and fraternal cooperators form a role-separating chessboard-like pattern.


  Over the last few years, several computational techniques have been devised
to recover protein complexes from the protein interaction (PPI) networks of
organisms. These techniques model "dense" subnetworks within PPI networks as
complexes. However, our comprehensive evaluations revealed that these
techniques fail to reconstruct many 'gold standard' complexes that are "sparse"
in the networks (only 71 recovered out of 123 known yeast complexes embedded in
a network of 9704 interactions among 1622 proteins). In this work, we propose a
novel index called Component-Edge (CE) score to quantitatively measure the
notion of "complex derivability" from PPI networks. Using this index, we
theoretically categorize complexes as "sparse" or "dense" with respect to a
given network. We then devise an algorithm SPARC that selectively employs
functional interactions to improve the CE scores of predicted complexes, and
thereby elevates many of the "sparse" complexes to "dense". This empowers
existing methods to detect these "sparse" complexes. We demonstrate that our
approach is effective in reconstructing significantly many complexes missed
previously (104 recovered out of the 123 known complexes or ~47% improvement).


  A device is investigated that continuously and directly converts light into
mechanical energy, using polymers and photodissociation. A polymer brush
tethered to a surface, is brought into contact with a parallel plate a small
distance above it that contains reaction sites where photodissociation of bound
polymer and light can occur. Under the appropriate conditions, the collective
effect of these polymers is to apply a force parallel to the plates, converting
incoming light into mechanical work. Numerical work is carried out to
understand this effect, a three dimensional Langevin simulation, solution to
the Fokker Planck equation, and a one dimensional Monte Carlo simulation.
Theoretical analysis of the Fokker Planck equation is used to study a model
where equilibration of the unbound state occurs and equilibration to a
metastable equilibrium is achieved in the bound state. It is shown that the
work per cycle can be made much larger than the thermal energy but at the
expense of requiring a greatly diminished photodissociation rate. Parameters
are discussed in order optimize mechanical energy conversion.


  Turing models have been proposed to explain the emergence of digits during
limb development. However, so far the molecular components that would give rise
to Turing patterns are elusive. We have recently shown that a particular type
of receptor-ligand interaction can give rise to Schnakenberg-type Turing
patterns, which reproduce patterning during lung and kidney branching
morphogenesis. Recent knock-out experiments have identified Smad4 as a key
protein in digit patterning. We show here that the BMP-receptor interaction
meets the conditions for a Schnakenberg-type Turing pattern, and that the
resulting model reproduces available wildtype and mutant data on the expression
patterns of BMP, its receptor, and Fgfs in the apical ectodermal ridge (AER)
when solved on a realistic 2D domain that we extracted from limb bud images of
E11.5 mouse embryos. We propose that receptor-ligand-based mechanisms serve as
a molecular basis for the emergence of Turing patterns in many developing
tissues.


  High-throughput spectrometers are capable of producing data sets containing
thousands of spectra for a single biological sample. These data sets contain a
substantial amount of redundancy from peptides that may get selected multiple
times in a LC-MS/MS experiment. In this paper, we present an efficient
algorithm, CAMS (Clustering Algorithm for Mass Spectra) for clustering mass
spectrometry data which increases both the sensitivity and confidence of
spectral assignment. CAMS utilizes a novel metric, called F-set, that allows
accurate identification of the spectra that are similar. A graph theoretic
framework is defined that allows the use of F-set metric efficiently for
accurate cluster identifications. The accuracy of the algorithm is tested on
real HCD and CID data sets with varying amounts of peptides. Our experiments
show that the proposed algorithm is able to cluster spectra with very high
accuracy in a reasonable amount of time for large spectral data sets. Thus, the
algorithm is able to decrease the computational time by compressing the data
sets while increasing the throughput of the data by interpreting low S/N
spectra.


  In order to understand how animals evolved over time, biologists must learn
how their body parts form during their development. The following is a detailed
description of how one species of sea spider transforms from a hatchling to an
8-legged adult. Pycnogonids, or sea spiders, comprise a primitive lineage of
arthropods. As such, they hold potential to reveal insights into arthropod
evolution. Recent phylogenetic analyses have supported their position as either
basal chelicerates or as a separate, fifth major lineage of extant arthropods.
Disagreements concerning pycnogonid relations to other arthropods is partly a
product of too few primary observations of pycnogonid anatomy and development.
This investigation of post-embryonic development of the pycnogonid
Anoplodactylus eroticus employs multiple techniques of anatomical observation,
including Nomarski optics, scanning electron microscopy and fluorescence
microscopy, in order to thoroughly document the life cycle.
  After the second post-embryonic stage, larvae of A. eroticus burrow within a
hydroid and undergo morphogenesis. Larvae emerge from the hydroid and
simultaneously molt into the juvenile stage. Over the course of post-embryonic
development there are eight stages preceding the mature adult. All structures,
except for the anteriormost appendages, the chelifores, undergo some degree of
transformation. Chelifores are present prior to hatching and remain mobile over
the course of development. Some larger issues important to arthropod evolution
are addressed, such as the equivalent of a germband and labrum in pycnogonids.
Post-embryonic development of A. eroticus provides an example counteracting
previous reports of anamorphic development and a four-segmented head in the
pycnogonid ground pattern, findings that were extrapolated to fit the ground
pattern of Arthropoda.


  Markov state models (MSMs) have become a popular approach for investigating
the conformational dynamics of proteins and other biomolecules. MSMs are
typically built from numerous molecular dynamics simulations by dividing the
sampled configurations into a large number of microstates based on geometric
criteria. The resulting microstate model can then be coarse-grained into a more
understandable macro state model by lumping together rapidly mixing microstates
into larger, metastable aggregates. However, finite sampling often results in
the creation of many poorly sampled microstates. During coarse-graining, these
states are mistakenly identified as being kinetically important because
transitions to/from them appear to be slow. In this paper we propose a
formalism based on an algebraic principle for matrix approximation, i.e. the
Nystrom method, to deal with such poorly sampled microstates. Our scheme builds
a hierarchy of microstates from high to low populations and progressively
applies spectral clustering on sets of microstates within each level of the
hierarchy. It helps spectral clustering identify metastable aggregates with
highly populated microstates rather than being distracted by lowly populated
states. We demonstrate the ability of this algorithm to discover the major
metastable states on two model systems, the alanine dipeptide and TrpZip2.


  Superoxide dismutase-1 (SOD1) is a ubiquitous, Cu and Zn binding, free
radical defense enzyme whose misfolding and aggregation play a potential key
role in amyotrophic lateral sclerosis, an invariably fatal neurodegenerative
disease. Over 150 mutations in SOD1 have been identified with a familial form
of the disease, but it is presently not clear what unifying features, if any,
these mutants share to make them pathogenic. Here, we develop several new
computational assays for probing the thermo-mechanical properties of both
ALS-associated and rationally-designed SOD1 variants. Allosteric interaction
free energies between residues and metals are calculated, and a series of
atomic force microscopy experiments are simulated with variable tether
positions, to quantify mechanical rigidity "fingerprints" for SOD1 variants.
Mechanical fingerprinting studies of a series of C-terminally truncated
mutants, along with an analysis of equilibrium dynamic fluctuations while
varying native constraints, potential energy change upon mutation,
frustratometer analysis, and analysis of the coupling between local frustration
and metal binding interactions for a glycine scan of 90 residues together
reveal that the apo protein is internally frustrated, that these internal
stresses are partially relieved by mutation but at the expense of metal-binding
affinity, and that the frustration of a residue is directly related to its role
in binding metals. This evidence points to apo SOD1 as a strained intermediate
with "self-allostery" for high metal-binding affinity. Thus, the prerequisites
for the function of SOD1 as an antioxidant compete with apo state
thermo-mechanical stability, increasing the susceptibility of the protein to
misfold in the apo state.


  The Health State Function theory is applied to find a quantitative estimate
of the Human Development Stages by defining and calculating the specific age
groups and subgroups. Early and late adolescence stages, first, second and
third stages of adult development are estimated along with the early, middle
and old age groups and subgroups. We briefly present the first exit time theory
used to find the health state function of a population and then we give the
details of the new theoretical approach with the appropriate applications to
support and validate the theoretical assumptions. Our approach is useful for
people working in several scientific fields and especially in medicine,
biology, anthropology, psychology, gerontology, probability and statistics. The
results are connected with the speed and acceleration of the deterioration of
the human organism during age as a consequence of the changes in the first,
second and third differences of the Health State Function and of the
Deterioration Function.
  Keywords: Human development stages, Deterioration, Deterioration function,
Human Mortality Database, HMD, World Health Organization, WHO, Quantitative
methods, Health State Function, Erikson's stages of psychosocial development,
Piaget method, Sullivan method, Disability stages, Light disability, Moderate
disability, Severe disability stage, Old ages, Critical ages.


  Inference of gene regulatory network from expression data is a challenging
task. Many methods have been developed to this purpose but a comprehensive
evaluation that covers unsupervised, semi-supervised and supervised methods,
and provides guidelines for their practical application, is lacking.
  We performed an extensive evaluation of inference methods on simulated
expression data. The results reveal very low prediction accuracies for
unsupervised techniques with the notable exception of the z-score method on
knock-out data. In all other cases the supervised approach achieved the highest
accuracies and even in a semi-supervised setting with small numbers of only
positive samples, outperformed the unsupervised techniques.


  Cellular senescence is thought to play a major role in age-related diseases,
which cause nearly 67% of all human deaths worldwide. Recent research in mice
showed that exercising mice had higher levels of telomerase, an enzyme that
helps maintain telomere length, than non-exercising mice. A commonly used model
for biological aging was proposed by Penna. I propose two modifications of the
Penna model that incorporate senescence and find analytical steady state
solutions following Coe, Mao and Cates. I find that models corresponding to
delayed senescence have younger populations that live longer.


  Simulating protein-membrane interactions is an important and dynamic area of
research. A proper definition of electrostatic forces on membrane surfaces is
necessary for developing electromechanical models of protein-membrane
interactions. Here we modeled the bilayer membrane as a continuum with general
continuous distributions of lipids charges on membrane surfaces. A new
electrostatic potential energy functional was then defined for this solvated
protein-membrane system. We investigated the geometrical transformation
properties of the membrane surfaces under a smooth velocity field. These
properties allows us to apply the Hadamard-Zolesio structure theorem, and the
electrostatic forces on membrane surfaces can be computed as the shape
derivative of the electrostatic energy functional.


  This document is another installment in a series of near real-time weekly
influenza forecasts made during the 2012-2013 influenza season. Here we present
some of the results of forecasts initiated following assimilation of
observations for Week 52 (i.e. the forecast begins December 30, 2012) for
municipalities in the United States. The forecasts were made on January 4,
2013. Results from forecasts initiated the five previous weeks (Weeks 47-51)
are also presented.


  With exquisite precision and reproducibility, cells orchestrate the
cooperative action of thousands of nanometer-sized molecular motors to carry
out mechanical tasks at much larger length scales, such as cell motility,
division and replication. Besides their biological importance, such inherently
non-equilibrium processes are an inspiration for developing biomimetic active
materials from microscopic components that consume energy to generate
continuous motion. Being actively driven, these materials are not constrained
by the laws of equilibrium statistical mechanics and can thus exhibit highly
sought-after properties such as autonomous motility, internally generated flows
and self-organized beating. Starting from extensile microtubule bundles, we
hierarchically assemble active analogs of conventional polymer gels, liquid
crystals and emulsions. At high enough concentration, microtubules form a
percolating active network characterized by internally driven chaotic flows,
hydrodynamic instabilities, enhanced transport and fluid mixing. When confined
to emulsion droplets, 3D networks spontaneously adsorb onto the droplet
surfaces to produce highly active 2D nematic liquid crystals whose streaming
flows are controlled by internally generated fractures and self-healing, as
well as unbinding and annihilation of oppositely charged disclination defects.
The resulting active emulsions exhibit unexpected properties, such as
autonomous motility, which are not observed in their passive analogues. Taken
together, these observations exemplify how assemblages of animate microscopic
objects exhibit collective biomimetic properties that are starkly different
from those found in materials assembled from inanimate building blocks,
challenging us to develop a theoretical framework that would allow for a
systematic engineering of their far-from-equilibrium material properties.


  We investigate the cooperative dynamics of an ensemble of N microtubules
growing against an elastic barrier. Microtubules undergo so-called
catastrophes, which are abrupt stochastic transitions from a growing to a
shrinking state, and rescues, which are transitions back to the growing state.
Microtubules can exert pushing or polymerization forces on an obstacle, such as
an elastic barrier if the growing end is in contact with the obstacle. We use
dynamical mean-field theory and stochastic simulations to analyze a model where
each microtubule undergoes catastrophes and rescues and where microtubules
interact by force sharing. For zero rescue rate, cooperative growth terminates
in a collective catastrophe. The maximal polymerization force before
catastrophes grows linearly with N for small N or a stiff elastic barrier, in
agreement with available experimental results, whereas it crosses over to a
logarithmic dependence for larger N or a soft elastic barrier. For a nonzero
rescue rate and a soft elastic barrier, the dynamics becomes oscillatory with
both collective catastrophe and rescue events, which are part of a robust limit
cycle. Both the average and maximal polymerization forces then grow linearly
with N, and we investigate their dependence on tubulin on-rates and rescue
rates, which can be involved in cellular regulation mechanisms. We further
investigate the robustness of the collective catastrophe and rescue
oscillations with respect to different catastrophe models.


  This paper considers adiabatic reduction in both discrete and continuous
models of stochastic gene expression. In gene expression models, the concept of
bursting is a production of several molecules simultaneously and is generally
represented as a compound Poisson process of random size. In a general
two-dimensional birth and death discrete model, we prove that under specific
assumptions and scaling (that are characteristics of the mRNA-protein system)
an adiabatic reduction leads to a one-dimensional discrete-state space model
with bursting production. The burst term appears through the reduction of the
first variable. In a two-dimensional continuous model, we also prove that an
adiabatic reduction can be performed in a stochastic slow/fast system. In this
gene expression model, the production of mRNA (the fast variable) is assumed to
be bursty and the production of protein (the slow variable) is linear as a
function of mRNA. When the dynamics of mRNA is assumed to be faster than the
protein dynamics (due to a mRNA degradation rate larger than for the protein)
we prove that, with the appropriate scaling, the bursting phenomena can be
transmitted to the slow variable. We show that the reduced equation is either a
stochastic differential equation with a jump Markov process or a deterministic
ordinary differential equation depending on the scaling that is appropriate.
  These results are significant because adiabatic reduction techniques seem to
have not been applied to a stochastic differential system containing a jump
Markov process. Last but not least, for our particular system, the adiabatic
reduction allows us to understand what are the necessary conditions for the
bursting production-like of protein to occur.


  Many important stochastic counting models can be written as general
birth-death processes (BDPs). BDPs are continuous-time Markov chains on the
non-negative integers and can be used to easily parameterize a rich variety of
probability distributions. Although the theoretical properties of general BDPs
are well understood, traditionally statistical work on BDPs has been limited to
the simple linear (Kendall) process, which arises in ecology and evolutionary
applications. Aside from a few simple cases, it remains impossible to find
analytic expressions for the likelihood of a discretely-observed BDP, and
computational difficulties have hindered development of tools for statistical
inference. But the gap between BDP theory and practical methods for estimation
has narrowed in recent years. There are now robust methods for evaluating
likelihoods for realizations of BDPs: finite-time transition, first passage,
equilibrium probabilities, and distributions of summary statistics that arise
commonly in applications. Recent work has also exploited the connection between
continuously- and discretely-observed BDPs to derive EM algorithms for maximum
likelihood estimation. Likelihood-based inference for previously intractable
BDPs is much easier than previously thought and regression approaches analogous
to Poisson regression are straightforward to derive. In this review, we outline
the basic mathematical theory for BDPs and demonstrate new tools for
statistical inference using data from BDPs. We give six examples of BDPs and
derive EM algorithms to fit their parameters by maximum likelihood. We show how
to compute the distribution of integral summary statistics and give an example
application to the total cost of an epidemic. Finally, we suggest future
directions for innovation in this important class of stochastic processes.


  The interplay between global constraints and local material properties of
chain molecules is a subject of emerging interest. Studies of molecules that
are intrinsically chiral, such as double-stranded DNA, is one example. Their
properties generally depend on the local geometry, i.e. on curvature and
torsion, yet the paths of closed molecules are globally restricted by topology.
Molecules that fulfill a twist neutrality condition, a zero sum rule for the
incremental change in the rate of winding along the curve, will behave
neutrally to strain. This has implications for plasmids. For small circular
microDNAs it follows that there must exist a minimum length for these to be
double-stranded. It also follows that all microDNAs longer than the minimum
length must be concave. This counterintuitive result is consistent with the
kink-like appearance which has been observed for circular DNA. A prediction for
the total negative curvature of a circular microDNA is given as a function of
its length.


  An ab-initio numerical study of the density-dependent, evolutionary stable
dispersal strategy is presented. The simulations are based on a simple
discretei generation island model with four processes: reproduction, dispersal,
competition and local catastrophe. We do not impose any a priori constraints on
the dispersal schedule, allowing the entire schedule to evolve. We find that
the system converges at long times to a unique nontrivial dispersal schedule
such that the dispersal probability is a monotonically increasing function of
the density. We have explored the dependence of the selected dispersal strategy
on the various system parameters: mean number of offspring, site carrying
capacity, dispersal cost and system size. A few general scaling laws are seen
to emerge from the data.


  We study the adaptation dynamics of a maladapted asexual population on rugged
fitness landscapes with many local fitness peaks. The distribution of
beneficial fitness effects is assumed to belong to one of the three extreme
value domains, viz. Weibull, Gumbel and Fr{\'e}chet. We work in the strong
selection-weak mutation regime in which beneficial mutations fix sequentially,
and the population performs an uphill walk on the fitness landscape until a
local fitness peak is reached. A striking prediction of our analysis is that
the fitness difference between successive steps follows a pattern of
diminishing returns in the Weibull domain and accelerating returns in the
Fr{\'e}chet domain, as the initial fitness of the population is increased.
These trends are found to be robust with respect to fitness correlations. We
believe that this result can be exploited in experiments to determine the
extreme value domain of the distribution of beneficial fitness effects. Our
work here differs significantly from the previous ones that assume the
selection coefficient to be small. On taking large effect mutations into
account, we find that the length of the walk shows different qualitative trends
from those derived using small selection coefficient approximation.


Reply to Comment on 'Length Scale Dependence of DNA Mechanical Properties'


  During embryonic development, a complex organism is formed from a single
starting cell. These processes of growth and differentiation are driven by
large transcriptional changes, which are following the expression and activity
of transcription factors (TFs). This study sought to compare TF expression
during embryonic development in a diverse group of metazoan animals:
representatives of vertebrates (Danio rerio, Xenopus tropicalis), a chordate
(Ciona intestinalis) and invertebrate phyla such as insects (Drosophila
melanogaster, Anopheles gambiae) and nematodes (Caenorhabditis elegans) were
sampled, The different species showed overall very similar TF expression
patterns, with TF expression increasing during the initial stages of
development. C2H2 zinc finger TFs were over-represented and Homeobox TFs were
under-represented in the early stages in all species. We further clustered TFs
for each species based on their quantitative temporal expression profiles. This
showed very similar TF expression trends in development in vertebrate and
insect species. However, analysis of the expression of orthologous pairs
between more closely related species showed that expression of most individual
TFs is not conserved, following the general model of duplication and
diversification. The degree of similarity between TF expression between Xenopus
tropicalis and Danio rerio followed the hourglass model, with the greatest
similarity occuring during the early tailbud stage in Xenopus tropicalis and
the late segmentation stage in Danio rerio. However, for Drosophila
melanogaster and Anopheles gambiae there were two periods of high TF
transcriptome similarity, one during the Arthropod phylotypic stage at 8-10
hours into Drosophila development and the other later at 16-18 hours into
Drosophila development.


  In the last years, a remarkable theoretical effort has been made in order to
understand stability and complexity in ecological communities. The non-random
structures of real ecological interaction networks has been recognized as one
key ingredient contributing to the coexistence between high complexity and
stability in real ecosystems. However most of the theoretical studies have
considered communities with only one interaction type (either antagonistic,
competitive, or mutualistic). Recently it has been proposed a theoretical
analysis on multiple interaction types in ecological systems, concluding that:
a) Mixture of antagonistic and mutualistic interactions stabilize the system
with respect to the less realistic case of only one interaction type; b)
Complexity, measured in terms of the number of interactions and the number of
species, increases stability in systems with different types of interactions.
By introducing new theoretical investigations and analyzing 21 empirical data
sets representing mutualistic plant-pollinator networks, we show that that
conclusions are incorrect. We will prove that the positive complexity-stability
effect observed in systems with different kind of interactions is a mere
consequence of a rescaling of the interaction strengths, and thus unrelated to
the mixing of interaction types.


  It has been shown that minimum free energy structure for RNAs and RNA-RNA
interaction is often incorrect due to inaccuracies in the energy parameters and
inherent limitations of the energy model. In contrast, ensemble based
quantities such as melting temperature and equilibrium concentrations can be
more reliably predicted. Even structure prediction by sampling from the
ensemble and clustering those structures by Sfold [7] has proven to be more
reliable than minimum free energy structure prediction. The main obstacle for
ensemble based approaches is the computational complexity of the partition
function and base pairing probabilities. For instance, the space complexity of
the partition function for RNA-RNA interaction is $O(n^4)$ and the time
complexity is $O(n^6)$ which are prohibitively large [4,12]. Our goal in this
paper is to give a fast algorithm, based on sparse folding, to calculate an
upper bound on the partition function. Our work is based on the recent
algorithm of Hazan and Jaakkola [10]. The space complexity of our algorithm is
the same as that of sparse folding algorithms, and the time complexity of our
algorithm is $O(MFE(n)\ell)$ for single RNA and $O(MFE(m, n)\ell)$ for RNA-RNA
interaction in practice, in which $MFE$ is the running time of sparse folding
and $\ell \leq n$ ($\ell \leq n + m$) is a sequence dependent parameter.


  Zinc is found saturated in the deposited Amyloid-beta (AB) peptide plaques in
brains of patients subjected to Alzheimer disease (AD). Zinc binding to AB
promotes aggregations, including the toxic soluble AB species. Up to now, only
the region 1-16 of AB complexed with Zinc (AB16-Zn) is defined structurally in
experiment, requiring an efficient theoretical method to present the
interaction between zinc and AB peptide. In order to explore the induced
polarization effect on the global conformation fluctuations and the
experimentally observed coordination mode of AB16-Zn, in this work we consider
an all-atom molecular dynamics (MD) of AB16-Zn solvated in implicit water. In
our model the polarization effect affects the whole peptide is applied. The
induced dipoles are divided into three distinct scales according to their
distances from zinc. Besides, the atomistic polarizability on the coordinating
sidechains is rescaled to describe the electron redistribution effect. As a
comparison, another model which exactly follows the method of Sakharov and Lim
(J. Am. Chem. Soc., 127, 13, 2005) has been discussed also. We show that,
associated with proper van der Waals (vdW) parameters, our model not only
obtains the reasonable coordinating configuration of zinc binding site, but
also retains the global stabilization, especially the N-terminal region, of the
AB16-Zn. We suggest that it is the induced polarization effect that promotes
reasonable solvent exposures of hydrophobic/hydrophilic residues regarding
zinc-induced AB aggregation.


  Despite nearly two scores of research on RNA secondary structure and RNA-RNA
interaction prediction, the accuracy of the state-of-the-art algorithms are
still far from satisfactory. Researchers have proposed increasingly complex
energy models and improved parameter estimation methods in anticipation of
endowing their methods with enough power to solve the problem. The output has
disappointingly been only modest improvements, not matching the expectations.
Even recent massively featured machine learning approaches were not able to
break the barrier. In this paper, we introduce the notion of learnability of
the parameters of an energy model as a measure of its inherent capability. We
say that the parameters of an energy model are learnable iff there exists at
least one set of such parameters that renders every known RNA structure to date
the minimum free energy structure. We derive a necessary condition for the
learnability and give a dynamic programming algorithm to assess it. Our
algorithm computes the convex hull of the feature vectors of all feasible
structures in the ensemble of a given input sequence. Interestingly, that
convex hull coincides with the Newton polytope of the partition function as a
polynomial in energy parameters. We demonstrated the application of our theory
to a simple energy model consisting of a weighted count of A-U and C-G base
pairs. Our results show that this simple energy model satisfies the necessary
condition for less than one third of the input unpseudoknotted
sequence-structure pairs chosen from the RNA STRAND v2.0 database. For another
one third, the necessary condition is barely violated, which suggests that
augmenting this simple energy model with more features such as the Turner loops
may solve the problem. The necessary condition is severely violated for 8%,
which provides a small set of hard cases that require further investigation.


  For DNA sequences of various species we construct the Google matrix G of
Markov transitions between nearby words composed of several letters. The
statistical distribution of matrix elements of this matrix is shown to be
described by a power law with the exponent being close to those of outgoing
links in such scale-free networks as the World Wide Web (WWW). At the same time
the sum of ingoing matrix elements is characterized by the exponent being
significantly larger than those typical for WWW networks. This results in a
slow algebraic decay of the PageRank probability determined by the distribution
of ingoing elements. The spectrum of G is characterized by a large gap leading
to a rapid relaxation process on the DNA sequence networks. We introduce the
PageRank proximity correlator between different species which determines their
statistical similarity from the view point of Markov chains. The properties of
other eigenstates of the Google matrix are also discussed. Our results
establish scale-free features of DNA sequence networks showing their
similarities and distinctions with the WWW and linguistic networks.


  I present a review of the theoretical and computational methodologies that
have been used to model the assembly of viral capsids. I discuss the
capabilities and limitations of approaches ranging from equilibrium continuum
theories to molecular dynamics simulations, and I give an overview of some of
the important conclusions about virus assembly that have resulted from these
modeling efforts. Topics include the assembly of empty viral shells, assembly
around single-stranded nucleic acids to form viral particles, and assembly
around synthetic polymers or charged nanoparticles for nanotechnology or
biomedical applications. I present some examples in which modeling efforts have
promoted experimental breakthroughs, as well as directions in which the
connection between modeling and experiment can be strengthened.


  The finitely repeated Prisoners' Dilemma is a good illustration of the
discrepancy between the strategic behaviour suggested by a game-theoretic
analysis and the behaviour often observed among human players, where
cooperation is maintained through most of the game. A game-theoretic reasoning
based on backward induction eliminates strategies step by step until defection
from the first round is the only remaining choice, reflecting the Nash
equilibrium of the game. We investigate the Nash equilibrium solution for two
different sets of strategies in an evolutionary context, using
replicator-mutation dynamics. The first set consists of conditional
cooperators, up to a certain round, while the second set in addition to these
contains two strategy types that react differently on the first round action:
The "Convincer" strategies insist with two rounds of initial cooperation,
trying to establish more cooperative play in the game, while the "Follower"
strategies, although being first round defectors, have the capability to
respond to an invite in the first round. For both of these strategy sets,
iterated elimination of strategies shows that the only Nash equilibria are
given by defection from the first round. We show that the evolutionary dynamics
of the first set is always characterised by a stable fixed point, corresponding
to the Nash equilibrium, if the mutation rate is sufficiently small (but still
positive). The second strategy set is numerically investigated, and we find
that there are regions of parameter space where fixed points become unstable
and the dynamics exhibits cycles of different strategy compositions. The
results indicate that, even in the limit of very small mutation rate, the
replicator-mutation dynamics does not necessarily bring the system with
Convincers and Followers to the fixed point corresponding to the Nash
equilibrium of the game.


  All non-human great apes are endangered in the wild, and it is therefore
important to gain an understanding of their demography and genetic diversity.
To date, however, genetic studies within these species have largely been
confined to mitochondrial DNA and a small number of other loci. Here, we
present a genome-wide survey of genetic variation in gorillas using a reduced
representation sequencing approach, focusing on the two lowland subspecies. We
identify 3,274,491 polymorphic sites in 14 individuals: 12 western lowland
gorillas (Gorilla gorilla gorilla) and 2 eastern lowland gorillas (Gorilla
beringei graueri). We find that the two species are genetically distinct, based
on levels of heterozygosity and patterns of allele sharing. Focusing on the
western lowland population, we observe evidence for population substructure,
and a deficit of rare genetic variants suggesting a recent episode of
population contraction. In western lowland gorillas, there is an elevation of
variation towards telomeres and centromeres on the chromosomal scale. On a
finer scale, we find substantial variation in genetic diversity, including a
marked reduction close to the major histocompatibility locus, perhaps
indicative of recent strong selection there. These findings suggest that
despite their maintaining an overall level of genetic diversity equal to or
greater than that of humans, population decline, perhaps associated with
disease, has been a significant factor in recent and long-term pressures on
wild gorilla populations.


  The ongoing functional annotation of proteins relies upon the work of
curators to capture experimental findings from scientific literature and apply
them to protein sequence and structure data. However, with the increasing use
of high-throughput experimental assays, a small number of experimental studies
dominate the functional protein annotations collected in databases. Here we
investigate just how prevalent is the "few articles -- many proteins"
phenomenon. We examine the experimentally validated annotation of proteins
provided by several groups in the GO Consortium, and show that the distribution
of proteins per published study is exponential, with 0.14% of articles
providing the source of annotations for 25% of the proteins in the UniProt-GOA
compilation. Since each of the dominant articles describes the use of an assay
that can find only one function or a small group of functions, this leads to
substantial biases in what we know about the function of many proteins.
Mass-spectrometry, microscopy and RNAi experiments dominate high throughput
experiments. Consequently, the functional information derived from these
experiments is mostly of the subcellular location of proteins, and of the
participation of proteins in embryonic developmental pathways. For some
organisms, the information provided by different studies overlap by a large
amount. We also show that the information provided by high throughput
experiments is less specific than those provided by low throughput experiments.
Given the experimental techniques available, certain biases in protein function
annotation due to high-throughput experiments are unavoidable. Knowing that
these biases exist and understanding their characteristics and extent is
important for database curators, developers of function annotation programs,
and anyone who uses protein function annotation data to plan experiments.


  This paper considers the behavior of discrete and continuous mathematical
models for gene expression in the presence of transcriptional/translational
bursting. We treat this problem in generality with respect to the distribution
of the burst size as well as the frequency of bursting, and our results are
applicable to both inducible and repressible expression patterns in prokaryotes
and eukaryotes. We have given numerous examples of the applicability of our
results, especially in the experimentally observed situation that burst size is
geometrically or exponentially distributed.


  Molecular dynamics simulation is utilized to investigate the ionic transport
of NaCl in solution through a graphene nanopore under an applied electric
field. Results show the formation of concentration polarization layers in the
vicinity of the graphene sheet. The non-uniformity of the ion distribution
gives rise to an electric pressure which drives vortical motions in the fluid
if the electric field is sufficiently strong to overcome the influence of
viscosity and thermal fluctuations. The relative importance of hydrodynamic
transport and thermal fluctuations in determining the pore conductivity is
investigated. A second important effect that is observed is the mass transport
of water through the nanopore, with an average velocity proportional to the
applied voltage and independent of the pore diameter. The flux arises as a
consequence of the asymmetry in the ion distribution with respect to reflection
about the plane of the graphene sheet. The accumulation of liquid molecules in
the vicinity of the nanopore due to reorientation of the water dipoles by the
local electric field is seen to result in a local increasein the liquid
density. Results confirm that the electric conductance is proportional to the
nanopore diameter for the parameter regimes that we simulated. The occurrence
of fluid vortices is found to result in an increase in the effective electrical
conductance.


  Bacteriophages infect cells by attaching to the outer membrane and injecting
their DNA into the cell.The phage DNA is then transcribed by the cell's
transcription machinery.A number of physical mechanisms by which DNA can be
translocated from the phage capsid into the cell have been identified. A fast
ejection driven by the elastic and electrostatic potential energy of the
compacted DNA within the viral capsid appears to be used by most phages, at
least to initiate infection.In recent in vitro experiments, the speed of DNA
translocation from a lambda phage capsid has been measured as a function of
ejected length over the entire duration of the event.Here a mechanical model is
proposed that is able to explain the observed dependence of exit velocity on
ejected length, and that is also consistent with the accepted picture of the
geometric arrangement of DNA within the viral capsid.