  Molecular biology knowledge can be systematically represented in a
computer-readable form as a comprehensive map of molecular interactions. There
exist a number of maps of molecular interactions containing detailed
description of various cell mechanisms. It is difficult to explore these large
maps, to comment their content and to maintain them. Though there exist several
tools addressing these problems individually, the scientific community still
lacks an environment that combines these three capabilities together. NaviCell
is a web-based environment for exploiting large maps of molecular interactions,
created in CellDesigner, allowing their easy exploration, curation and
maintenance. NaviCell combines three features: (1) efficient map browsing based
on Google Maps engine; (2) semantic zooming for viewing different levels of
details or of abstraction of the map and (3) integrated web-based blog for
collecting the community feedback. NaviCell can be easily used by experts in
the field of molecular biology for studying molecular entities of their
interest in the context of signaling pathways and cross-talks between pathways
within a global signaling network. NaviCell allows both exploration of detailed
molecular mechanisms represented on the map and a more abstract view of the map
up to a top-level modular representation. NaviCell facilitates curation,
maintenance and updating the comprehensive maps of molecular interactions in an
interactive fashion due to an imbedded blogging system. NaviCell provides an
easy way to explore large-scale maps of molecular interactions, thanks to the
Google Maps and WordPress interfaces, already familiar to many users. Semantic
zooming used for navigating geographical maps is adopted for molecular maps in
NaviCell, making any level of visualization meaningful to the user. In
addition, NaviCell provides a framework for community-based map curation.


  The spatial arrangement of individuals is thought to overcome the dilemma of
cooperation: When cooperators engage in clusters they might share the benefit
of cooperation while being more protected against non-cooperating individuals,
which benefit from cooperation but save the cost of cooperation. This is
paradigmatically shown by the spatial prisoner's dilemma model. Here, we study
this model in one and two spatial dimensions, but explicitly take into account
that in biological setups fitness collection and selection are separated
processes occurring mostly on vastly different time scales. This separation is
particularly important to understand the impact of mobility on the evolution of
cooperation. We find that even small diffusive mobility strongly restricts
cooperation since it enables non-cooperative individuals to invade cooperative
clusters. Thus, in most biological scenarios, where the mobility of competing
individuals is an irrefutable fact, the spatial prisoner's dilemma alone cannot
explain stable cooperation but additional mechanisms are necessary for spatial
structure to promote the evolution of cooperation. The breakdown of cooperation
is analyzed in detail. We confirm the existence of a phase transition, here
controlled by mobility and costs, which distinguishes between purely
cooperative and non-cooperative absorbing states. While in one dimension the
model is in the class of the Voter Model, it belongs to the Directed
Percolation (DP) universality class in two dimensions.


  The theory of slow manifolds is an important tool in the study of
deterministic dynamical systems, giving a practical method by which to reduce
the number of relevant degrees of freedom in a model, thereby often resulting
in a considerable simplification. In this article we demonstrate how the same
basic methodology may also be applied to stochastic dynamical systems, by
examining the behaviour of trajectories conditioned on the event that they do
not depart the slow manifold. We apply the method to two models: one from
ecology and one from epidemiology, achieving a reduction in model dimension and
illustrating the high quality of the analytical approximations.


  Chemical mapping is a widespread technique for structural analysis of nucleic
acids in which a molecule's reactivity to different probes is quantified at
single-nucleotide resolution and used to constrain structural modeling. This
experimental framework has been extensively revisited in the past decade with
new strategies for high-throughput read-outs, chemical modification, and rapid
data analysis. Recently, we have coupled the technique to high-throughput
mutagenesis. Point mutations of a base-paired nucleotide can lead to exposure
of not only that nucleotide but also its interaction partner. Carrying out the
mutation and mapping for the entire system gives an experimental approximation
of the molecules contact map. Here, we give our in-house protocol for this
mutate-and-map strategy, based on 96-well capillary electrophoresis, and we
provide practical tips on interpreting the data to infer nucleic acid
structure.


  The Rosetta molecular modeling software package provides experimentally
tested and rapidly evolving tools for the 3D structure prediction and
high-resolution design of proteins, nucleic acids, and a growing number of
non-natural polymers. Despite its free availability to academic users and
improving documentation, use of Rosetta has largely remained confined to
developers and their immediate collaborators due to the code's difficulty of
use, the requirement for large computational resources, and the unavailability
of servers for most of the Rosetta applications. Here, we present a unified web
framework for Rosetta applications called ROSIE (Rosetta Online Server that
Includes Everyone). ROSIE provides (a) a common user interface for Rosetta
protocols, (b) a stable application programming interface for developers to add
additional protocols, (c) a flexible back-end to allow leveraging of computer
cluster resources shared by RosettaCommons member institutions, and (d)
centralized administration by the RosettaCommons to ensure continuous
maintenance. This paper describes the ROSIE server infrastructure, a
step-by-step 'serverification' protocol for use by Rosetta developers, and the
deployment of the first nine ROSIE applications by six separate developer
teams: Docking, RNA de novo, ERRASER, Antibody, Sequence Tolerance,
Supercharge, Beta peptide design, NCBB design, and VIP redesign. As illustrated
by the number and diversity of these applications, ROSIE offers a general and
speedy paradigm for serverification of Rosetta applications that incurs
negligible cost to developers and lowers barriers to Rosetta use for the
broader biological community. ROSIE is available at
http://rosie.rosettacommons.org.


  We study non-adaptive pooling strategies for detection of rare faulty items.
Given a binary sparse N-dimensional signal x, how to construct a sparse binary
MxN pooling matrix F such that the signal can be reconstructed from the
smallest possible number M of measurements y=Fx? We show that a very low number
of measurements is possible for random spatially coupled design of pools F. Our
design might find application in genetic screening or compressed genotyping. We
show that our results are robust with respect to the uncertainty in the matrix
F when some elements are mistaken.


  Understanding models which represent the invasion of network-based systems by
infectious agents can give important insights into many real-world situations,
including the prevention and control of infectious diseases and computer
viruses. Here we consider Markovian susceptible-infectious-susceptible (SIS)
dynamics on finite strongly connected networks, applicable to several sexually
transmitted diseases and computer viruses. In this context, a theoretical
definition of endemic prevalence is easily obtained via the quasi-stationary
distribution (QSD). By representing the model as a percolation process and
utilising the property of duality, we also provide a theoretical definition of
invasion probability. We then show that, for undirected networks, the
probability of invasion from any given individual is equal to the
(probabilistic) endemic prevalence, following successful invasion, at the
individual (we also provide a relationship for the directed case). The total
(fractional) endemic prevalence in the population is thus equal to the average
invasion probability (across all individuals). Consequently, for such systems,
the regions or individuals already supporting a high level of infection are
likely to be the source of a successful invasion by another infectious agent.
This could be used to inform targeted interventions when there is a threat from
an emerging infectious disease.


  A formalism is developed which describes the extent to which stochastic
oscillations in biochemical models are synchronised. It is based on the
calculation of the complex coherence function within the linear noise
approximation. The method is illustrated on a simple example and then applied
to study the synchronisation of chemical concentrations in social amoeba. The
degree to which variation of rate constants in different cells and the volume
of the cells affects synchronisation of the oscillations is explored, and the
phase lag calculated. In all cases the analytical results are shown to be in
good agreement with those obtained through numerical simulations.


  Both genetic drift and natural selection cause the frequencies of alleles in
a population to vary over time. Discriminating between these two evolutionary
forces, based on a time series of samples from a population, remains an
outstanding problem with increasing relevance to modern data sets. Even in the
idealized situation when the sampled locus is independent of all other loci
this problem is difficult to solve, especially when the size of the population
from which the samples are drawn is unknown. A standard $\chi^2$-based
likelihood ratio test was previously proposed to address this problem. Here we
show that the $\chi^2$ test of selection substantially underestimates the
probability of Type I error, leading to more false positives than indicated by
its $P$-value, especially at stringent $P$-values. We introduce two methods to
correct this bias. The empirical likelihood ratio test (ELRT) rejects
neutrality when the likelihood ratio statistic falls in the tail of the
empirical distribution obtained under the most likely neutral population size.
The frequency increment test (FIT) rejects neutrality if the distribution of
normalized allele frequency increments exhibits a mean that deviates
significantly from zero. We characterize the statistical power of these two
tests for selection, and we apply them to three experimental data sets. We
demonstrate that both ELRT and FIT have power to detect selection in practical
parameter regimes, such as those encountered in microbial evolution
experiments. Our analysis applies to a single diallelic locus, assumed
independent of all other loci, which is most relevant to full-genome selection
scans in sexual organisms, and also to evolution experiments in asexual
organisms as long as clonal interference is weak. Different techniques will be
required to detect selection in time series of co-segregating linked loci.


  We are interested in modeling the Darwinian evolution resulting from the
interplay of phenotypic variation and natural selection through ecological
interactions, in the specific scales of the biological framework of adaptive
dynamics. Adaptive dynamics so far has been put on a rigorous footing only for
direct competition models (Lotka-Volterra models) involving a competition
kernel which describes the competition pressure from one individual to another
one. We extend this to a multi-resources chemostat model, where the competition
between individuals results from the sharing of several resources which have
their own dynamics. Starting from a stochastic birth and death process model,
we prove that, when advantageous mutations are rare, the population behaves on
the mutational time scale as a jump process moving between equilibrium states
(the polymorphic evolution sequence of the adaptive dynamics literature). An
essential technical ingredient is the study of the long time behavior of a
chemostat multi-resources dynamical system. In the small mutational steps limit
this process in turn gives rise to a differential equation in phenotype space
called canonical equation of adaptive dynamics. From this canonical equation
and still assuming small mutation steps, we prove a rigorous characterization
of the evolutionary branching points.


  Network epidemiology's most important assumption is that the contact
structure over which infectious diseases propagate can be represented as a
static network. However, contacts are highly dynamic, changing at many time
scales. In this paper, we investigate conceptually simple methods to construct
static graphs for network epidemiology from temporal contact data. We evaluate
these methods on empirical and synthetic model data. For almost all our cases,
the network representation that captures most relevant information is a
so-called exponential-threshold network. In these, each contact contributes
with a weight decreasing exponentially with time, and there is an edge between
a pair of vertices if the weight between them exceeds a threshold. Networks of
aggregated contacts over an optimally chosen time window perform almost as good
as the exponential-threshold networks. On the other hand, networks of
accumulated contacts over the entire sampling time, and networks of concurrent
partnerships, perform worse. We discuss these observations in the context of
the temporal and topological structure of the data sets.


  A new approach to estimate the Shannon entropy of a long-range correlated
sequence is proposed. The entropy is written as the sum of two terms
corresponding respectively to power-law (\emph{ordered}) and exponentially
(\emph{disordered}) distributed blocks (clusters). The approach is illustrated
on the 24 human chromosome sequences by taking the nucleotide composition as
the relevant information to be encoded/decoded. Interestingly, the nucleotide
composition of the \emph{ordered} clusters is found, on the average, comparable
to the one of the whole analyzed sequence, while that of the \emph{disordered}
clusters fluctuates. From the information theory standpoint, this means that
the power-law correlated clusters carry the same information of the whole
analysed sequence. Furthermore, the fluctuations of the nucleotide composition
of the disordered clusters are linked to relevant biological properties, such
as segmental duplications and gene density.


  Cell type (e.g. pluripotent cell, fibroblast) is the end result of many
complex processes that unfold due to evolutionary, developmental, and
transformational stimuli. A cell's phenotype and the discrete, a priori states
that define various cell subtypes (e.g. skin fibroblast, embryonic stem cell)
are ultimately part of a continuum that may predict changes and systematic
variation in cell subtypes. These features can be both observable in existing
cellular states and hypothetical (e.g. unobserved). In this paper, a series of
approaches will be used to approximate the continuous diversity of gene
expression across a series of pluripotent, totipotent, and fibroblast cellular
subtypes. We will use a series of previously-collected datasets and analyze
them using three complementary approaches: the computation of distances based
on the subsampling of diversity, assessing the separability of individual genes
for a specific cell line both within and between cell types, and a hierarchical
soft classification technique that will assign a membership value for specific
genes in specific cell types given a number of different criteria. These
approaches will allow us to assess the observed gene-expression diversity in
these datasets, as well as assess how well a priori cell types characterize
their constituent populations. In conclusion, the application of these findings
to a broader biological context will be discussed.


  The rates of escape and reversion in response to selection pressure arising
from the host immune system, notably the cytotoxic T-lymphocyte (CTL) response,
are key factors determining the evolution of HIV. Existing methods for
estimating these parameters from cross-sectional population data using ordinary
differential equations (ODE) ignore information about the genealogy of sampled
HIV sequences, which has the potential to cause systematic bias and
over-estimate certainty. Here, we describe an integrated approach, validated
through extensive simulations, which combines genealogical inference and
epidemiological modelling, to estimate rates of CTL escape and reversion in HIV
epitopes. We show that there is substantial uncertainty about rates of viral
escape and reversion from cross-sectional data, which arises from the inherent
stochasticity in the evolutionary process. By application to empirical data, we
find that point estimates of rates from a previously published ODE model and
the integrated approach presented here are often similar, but can also differ
several-fold depending on the structure of the genealogy. The model-based
approach we apply provides a framework for the statistical analysis of escape
and reversion in population data and highlights the need for longitudinal and
denser cross-sectional sampling to enable accurate estimate of these key
parameters.


  Why evolvability appears to have increased over evolutionary time is an
important unresolved biological question. Unlike most candidate explanations,
this paper proposes that increasing evolvability can result without any
pressure to adapt. The insight is that if evolvability is heritable, then an
unbiased drifting process across genotypes can still create a distribution of
phenotypes biased towards evolvability, because evolvable organisms diffuse
more quickly through the space of possible phenotypes. Furthermore, because
phenotypic divergence often correlates with founding niches, niche founders may
on average be more evolvable, which through population growth provides a
genotypic bias towards evolvability. Interestingly, the combination of these
two mechanisms can lead to increasing evolvability without any pressure to
out-compete other organisms, as demonstrated through experiments with a series
of simulated models. Thus rather than from pressure to adapt, evolvability may
inevitably result from any drift through genotypic space combined with
evolution's passive tendency to accumulate niches.


  To learn about the past from a sample of genomic sequences, one needs to
understand how evolutionary processes shape genetic diversity. Most population
genetic inference is based on frameworks assuming adaptive evolution is rare.
But if positive selection operates on many loci simultaneously, as has recently
been suggested for many species including animals such as flies, a different
approach is necessary. In this review, I discuss recent progress in
characterizing and understanding evolution in rapidly adapting populations
where random associations of mutations with genetic backgrounds of different
fitness, i.e., genetic draft, dominate over genetic drift. As a result, neutral
genetic diversity depends weakly on population size, but strongly on the rate
of adaptation or more generally the variance in fitness. Coalescent processes
with multiple mergers, rather than Kingman's coalescent, are appropriate
genealogical models for rapidly adapting populations with important
implications for population genetic inference.


  The model considers the human population, the adult mosquito population and
the population of immature stages, which includes eggs, larvae and pupae. The
model also considers the vertical transmission of dengue in the mosquitoes and
the seasonal variation in the mosquito population. From this basic model
describing the dynamics of dengue infection, we deduce thresholds for avoiding
the introduction of the disease and for the elimination of the disease. In
particular, we deduce a Basic Reproduction Number for dengue that includes
parameters related to the immature stages of the mosquito. By neglecting
seasonal variation, we calculate the equilibrium values of the models
variables. We also present a sensitivity analysis of the impact of four
vector-control strategies on the Basic Reproduction Number and on the Force of
Infection of dengue. Each of the strategies was studied separately from the
others. The analysis presented allows us to conclude that of the available
vector control strategies, adulticide application is the most effective,
followed by the reduction of the exposure to mosquito bites, locating and
destroying breeding places and, finally, larvicides. Current vector-control
methods are concentrated in mechanical destruction of mosquitoes breeding
places. Our results suggest that reducing the contact between vector and hosts
(biting rates) are as efficient as the logistically difficult but very
efficient adult mosquitos control.


  In this work, we model the zero-bias conductance for the four different DNA
strands that were used in conductance measurement experiment [A. K. Mahapatro,
K. J. Jeong, G. U. Lee, and D. B. Janes, Nanotechnology 18, 195202 (2007)]. Our
approach consists of three elements: (i) ab initio calculations of DNA, (ii)
Green's function approach for transport calculations, and (iii) the use of two
parameters to determine the decoherence rates. We first study the role of the
backbone. We find that the backbone can alter the coherent transmission
significantly at some energy points by interacting with the bases, though the
overall shape of the transmission stays similar for the two cases. More
importantly, we find that the coherent electrical conductance is tremendously
smaller than what the experiments measure. We consider DNA strands under a
variety of different experimental conditions and show that even in the most
ideal cases, the calculated coherent conductance is much smaller than the
experimental conductance. To understand the reasons for this, we carefully look
at the effect of decoherence. By including decoherence, we show that our model
can rationalize the measured conductance of the four strands, both
qualitatively and quantitatively. We find that the effect of decoherence on G:C
base pairs is crucial in getting agreement with the experiments. However, the
decoherence on G:C base pairs alone does not explain the experimental
conductance in strands containing a number of A:T base pairs. Including
decoherence on A:T base pairs is also essential. By fitting the experimental
trends and magnitudes in the conductance of the four different DNA molecules,
we estimate for the first time that the deocherence rate is 6 meV for G:C and
1.5 meV for A:T base pairs.


  Smolin (1992) proposed that the fine-tuning problem for parameters of the
Standard Model might be accounted for by a Darwinian process of universe
reproduction - Cosmological Natural Selection (CNS) - in which black holes give
rise to offspring universes with slightly altered parameters. The laws for
variation and inheritance of the parameters are also subject to CNS if
variation in transmission laws occurs. This is the strategy introduced by Nei
(1967) to understand genetic transmission, through the evolutionary theory of
modifier genes, whose methodology is adopted here. When mechanisms of variation
themselves vary, they are subject to Feldman's (1972) evolutionary Reduction
Principle that selection favors greater faithfulness of replication. A theorem
of Karlin (1982) allows one to generalize this principle beyond biological
genetics to the unknown inheritance laws that would operate in CNS. The
reduction principle for CNS is illustrated with a general multitype branching
process model of universe creation containing competing inheritance laws. The
most faithful inheritance law dominates the ensemble of universes. The
Reduction Principle thus provides a mechanism to account for high fidelity of
inheritance between universes. Moreover, it reveals that natural selection in
the presence of variation in inheritance mechanisms has two distinct objects:
maximization of both fitness and faithful inheritance. Tradeoffs between
fitness and faithfulness open the possibility that evolved fundamental
parameters are compromises, and not local optima to maximize universe
production, in which case their local non-optimality may point to their
involvement in the universe inheritance mechanisms.


  Accessibility percolation is a new type of percolation problem inspired by
evolutionary biology. To each vertex of a graph a random number is assigned and
a path through the graph is called accessible if all numbers along the path are
in ascending order. For the case when the random variables are independent and
identically distributed, we derive an asymptotically exact expression for the
probability that there is at least one accessible path from the root to the
leaves in an $n$-tree. This probability tends to 1 (0) if the branching number
is increased with the height of the tree faster (slower) than linearly. When
the random variables are biased such that the mean value increases linearly
with the distance from the root, a percolation threshold emerges at a finite
value of the bias.


  Stochastic modeling of phylogenies raises five questions that have received
varying levels of attention from quantitatively inclined biologists. 1) How
large do we expect (from the model) the ration of maximum historical diversity
to current diversity to be? 2) From a correct phylogeny of the extant species
of a clade, what can we deduce about past speciation and extinction rates? 3)
What proportion of extant species are in fact descendants of still-extant
ancestral species, and how does this compare with predictions od models? 4)
When one moves from trees on species to trees on sets of species (whether
traditional higher order taxa or clades from PhyloCode), does one expect trees
to become more unbiased as a purely logical consequence of tree structure,
without signifying any real biological phenomenon? 5) How do we expect that
fluctuation rates for counts of higher order taxa should compare with
fluctuation rates for number of species? WE present a mathematician's view
based on an oversimplified modeling framework in which all these questions can
be studied coherently.


  We study a stochastic two-species chemical reaction system with two
mechanisms. One mechanism consists of chemical interactions which govern the
overall drift of species amounts in the system; the other mechanism consists of
resampling, branching or splitting which makes unbiased perturbative changes to
species amounts. Our results show that in a system with a large but bounded
capacity, certain combinations of these two types of interactions can lead to
stochastically-induced bistability. Depending on the relative magnitudes of the
rates of these two sets of interactions, bistability can occur in two distinct
ways with different dynamical signatures.


  In this paper, we study a weak law of large numbers for the total internal
length of the Bolthausen-Szmitman coalescent. As a consequence, we obtain the
weak limit law of the centered and rescaled total external length. The latter
extends results obtained by Dhersin & M\"ohle \cite{DM12}. An application to
population genetics dealing with the total number of mutations in the
genealogical tree is also given.


  We examine possible sources of a substantial increase in tree ring 14C
measurements for the years AD 774-775. Contrary to claims regarding a coronal
mass ejection (CME), the required CME energy is not several orders of magnitude
greater than known solar events. We consider solar proton events (SPEs) with
three different fluences and two different spectra. The data may be explained
by an event with fluence about one order of magnitude beyond the October 1989
SPE. Two hard spectrum cases considered here result in moderate ozone
depletion, so no mass extinction is implied, though we do predict increases in
erythema and damage to plants from enhanced solar UV. We are able to rule out
an event with a very soft spectrum that causes severe ozone depletion and
subsequent biological impacts. Nitrate enhancements are consistent with their
apparent absence in ice core data. The modern technological implications of
such an event may be extreme, and considering recent confirmation of
superflares on solar-type stars, this issue merits attention.


  The shapes of cell membranes are largely regulated by membrane associated,
curvature active, proteins. We use a numerical model of the membrane with
elongated membrane inclusions, recently developed by us, which posses
spontaneous directional curvatures that could be different along and
perpendicular to its long axis. We show that, due to membrane mediated
interactions these curvature inducing membrane nematogens can oligomerize
spontaneously, even at low concentrations, and change the local shape of the
membrane. We demonstrate that for a large group of such inclusions, where the
two spontaneous curvatures have equal sign, the tubular conformation and
sometime the sheet conformation of the membrane are the common equilibrium
shapes. We elucidate the factors necessary for the formation of these {\it
protein lattices}. Furthermore, the elastic properties of the tubes, like their
compressional stiffness and persistence length are calculated. Finally, we
discuss the possible role of nematic disclination in capping and branching of
the tubular membranes.


  It is the first time, when genetic diversity and the common kilka population
structure were investigated throughout the areal. Data about the species
condition at the Upper Volga basins were updated accordingly to modern
condition. Physiological and ecological adaptations to northern water basins
were evaluated. Significance of interaction between some loci and the most
important abiotic environmental factors at the selection was demonstrated. The
results suggest that the common kilka Clupeonella cultriventris (Nordmann,
1840) is presented by uniform set of population throughout the areal. The
successive expansion of the species through the cascade of the Volga's water
reservoirs can be explained by complex of genetic and biochemical adaptation
for sweet water habitat. It was supposed that the sweet water populations
originated from Saratov's backwaters sweet water population. Seasonal
fluctuations of abiotic and biotic environmental factors have significant
influence on genotype distribution at the newly formed population. This book is
intended for use by ichthyologists, ecologists, environmental protection and
management of natural resources specialists. In Russian.


  Randomness is an unavoidable feature of the intracellular environment due to
chemical reactants being present in low copy number. That phenomenon, predicted
by Delbr\"uck long ago \cite{delbruck40}, has been detected in both prokaryotic
\cite{elowitz02,cai06} and eukaryotic \cite{blake03} cells after the
development of the fluorescence techniques. On the other hand, developing
organisms, e.g. {\em D. melanogaster}, exhibit strikingly precise
spatio-temporal patterns of protein/mRNA concentrations
\cite{gregor07b,manu09a,manu09b,boettiger09}. Those two characteristics of
living organisms are in apparent contradiction: the precise patterns of protein
concentrations are the result of multiple mutually interacting random chemical
reactions. The main question is to establish biochemical mechanisms for
coupling random reactions so that canalization, or fluctuations reduction
instead of amplification, takes place. Here we explore a model for coupling two
stochastic processes where the noise of the combined process can be smaller
than that of the isolated ones. Such a canalization occurs if, and only if,
there is negative covariance between the random variables of the model. Our
results are obtained in the framework of a master equation for a negatively
self-regulated -- or externally regulated -- binary gene and show that the
precise control due to negative self regulation \cite{becskei00} is because it
may generate negative covariance. Our results suggest that negative covariance,
in the coupling of random chemical reactions, is a theoretical mechanism
underlying the precision of developmental processes.


  The conformations available to polypeptides are determined by the interatomic
forces acting on the peptide units, whereby backbone torsion angles are
restricted as described by the Ramachandran plot. Although typical proteins are
composed predominantly from {\alpha}-helices and {\beta}-sheets, they
nevertheless adopt diverse tertiary structure, each folded as dictated by its
unique amino-acid sequence. Despite such uniqueness, however, the functioning
of many proteins involves changes between quite different conformations. The
study of large-scale conformational changes, particularly in large systems, is
facilitated by a coarse-grained representation such as provided by virtually
bonded C{\alpha} atoms. We have developed a virtual atom molecular mechanics
(VAMM) force field to describe conformational dynamics in proteins and a
VAMM-based algorithm for computing conformational transition pathways. Here we
describe the stereochemical analysis of proteins in this coarse-grained
representation, comparing the relevant plots in coarse-grained conformational
space to the corresponding Ramachandran plots, having contoured each at levels
determined statistically from residues in a large database. The distributions
shown for an all-{\alpha} protein, two all-{\beta} proteins and one
{\alpha}+{\beta} protein serve to relate the coarse-grained distributions to
the familiar Ramachandran plot.


  We report on atomistic simulation of the folding of a natively-knotted
protein, MJ0366, based on a realistic force field. To the best of our knowledge
this is the first reported effort where a realistic force field is used to
investigate the folding pathways of a protein with complex native topology. By
using the dominant-reaction pathway scheme we collected about 30 successful
folding trajectories for the 82-amino acid long trefoil-knotted protein.
Despite the dissimilarity of their initial unfolded configuration, these
trajectories reach the natively-knotted state through a remarkably similar
succession of steps. In particular it is found that knotting occurs essentially
through a threading mechanism, involving the passage of the C-terminal through
an open region created by the formation of the native beta-sheet at an earlier
stage. The dominance of the knotting by threading mechanism is not observed in
MJ0366 folding simulations using simplified, native-centric models. This points
to a previously underappreciated role of concerted amino acid interactions,
including non-native ones, in aiding the appropriate order of contact formation
to achieve knotting.


  Pattern analysis of tandem repeats in gene is an indispensable computational
approach to the understanding of the gene expression and pathogenesis of
diseases. This research applied a computational motif model and database
techniques to study the distribution of tandem repeats in Nlrp1 gene, which is
a critical gene to detect the invading pathogens in the immunologic mechanisms.
The frequency of tandem repeats in Nlrp1 gene was studied for mono-, di-, tri-,
and tetranucleotides. Mutations of Nlrp1 gene were analyzed to identify the
insertion,deletion, and substitution of nucleotides. The results of this
research provide a basis for future work in computational drug design and
biomedical engineering in tackling diseases associated with immune system.


  Tandem repeats are ubiquitous in the genome of organisms and their mutated
forms play a vital role in pathogenesis. In this study, tandem repeats in
Gastric Intrinsic Factor (GIF) of gastric parietal cells have been investigated
using an in silico approach. Six types of the nucleotide tandem repeat motifs
have been investigated, including mono-, di-, tri-, tetra-, penta- and
hexanucleotide. The distribution of the repeat motifs in the gene was analyzed.
The results of this study provide an insight into the biomolecular mechanisms
and pathogenesis implicated by the GIF of gastric parietal cells. Based on the
findings of the tandem repeats in GIF of gastric parietal cells, therapeutic
strategies and disease markers may be developed accordingly by the biomedical
scientists.


  Growing evidence suggests that many vertebrate lineages are evolving at
significantly different rates. As a first approximation of evolutionary rates,
we assessed the amount of neutral (dS) and non-neutral (dN) substitutions that
have accumulated within and across sister clades since the time of their
divergence. We found that in fish, tetraodontiformes (pufferfish) are evolving
at faster rates than cypriniformes (fresh water teleosts), while cypriniformes
are evolving faster than elasmobranchs (sharks, skates and rays). A similar
rate variation was observed in salamanders: plethodontidae were found to evolve
at a rate nearly two fold faster than the hydromantes lineage. We discuss
possible explanations for this striking variation in substitution rates among
different vertebrate lineages that occupy widely diverse habitats and niches.


  Genetic redundancy is ubiquitous and can be found in any organism. However,
it has been argued that genetic redundancy reduces total population fitness,
and therefore, redundancy is unlikely to evolve. In this letter, we study an
evolutionary model with high-dimensional genotype-phenotype mapping (GPM) to
investigate the relevance of complexity in GPM to the evolution of genetic
redundancy. By applying the replica method to deal with quenched randomness,
the redundancy dependence of the fitness is analytically obtained, which
demonstrates that genetic redundancy can indeed evolve, provided that the GPM
is complex. Our result provides a novel insight into how genetic redundancy
evolves.


  Complementary strands in DNA double helix show temporary fluctuational
openings which are essential to biological functions such as transcription and
replication of the genetic information. Such large amplitude fluctuations,
known as the breathing of DNA, are generally localized and, microscopically,
are due to the breaking of the hydrogen bonds linking the base pairs
(\emph{bps}). I apply imaginary time path integral techniques to a mesoscopic
Hamiltonian which accounts for the helicoidal geometry of a short circular DNA
molecule. The \emph{bps} displacements with respect to the ground state are
interpreted as time dependent paths whose amplitudes are consistent with the
model potential for the hydrogen bonds. The portion of the paths configuration
space contributing to the partition function is determined by selecting the
ensemble of paths which fulfill the second law of thermodynamics. Computations
of the thermodynamics in the denaturation range show the energetic advantage
for the equilibrium helicoidal geometry peculiar of B-DNA. I discuss the
interplay between twisting of the double helix and anharmonic stacking along
the molecule backbone suggesting an interesting relation between intrinsic
nonlinear character of the microscopic interactions and molecular topology.


  Consider a birth and death chain to model the number of types of a given
virus. Each type gives birth to a new type at rate $\lambda$ and dies at rate
1. Each type is also assigned a fitness. When a death occurs either the least
fit type dies (with probability $1-r$) or we kill a type at random (with
probability $r$). We show that this random killing has a large effect (for any
$r>0$) on the behavior of the model when $\lambda<1$. The behavior of the model
with $r>0$ and $\lambda<1$ is consistent with features of the phylogenetic tree
of influenza.


  DNA is the most important biological molecule, and its hydration contributes
essentially to the structure and functions of the double helix. We analyze the
microhydration of the individual bases of nucleic acids and their methyl
derivatives using methods of molecular mechanics (MM) with the Poltev-Malenkov
(PM), AMBER and OPLS force fields, as well as ab initio Quantum Mechanics (QM)
calculations at MP2/6-31G(d,p) level of theory. A comparison is made between
the calculated interaction energies and the experimental enthalpies of
microhydration of bases, obtained from mass spectrometry at low temperatures.
Each local water-base interaction energy minimum obtained with MM corresponds
to the minimum obtained with QM. General qualitative agreement was observed in
the geometrical characteristics of the local minima obtained via the two groups
of methods. MM minima correspond to slightly more coplanar structures than
those obtained via QM methods, and the absolute MM energy values overestimate
corresponding values obtained with QM. For Adenine and Thymine the QM local
minima energy values are closer to those obtained by the PM potential (average
of 0.72 kcal/mol) than by the AMBER force field (1.86 kcal/mol). The
differences in energy between MM and QM results are more pronounced for Guanine
and Cytosine, especially for minima with the water molecule forming H-bonds
with two proton-acceptor centers of the base. Such minima are the deepest ones
obtained via MM methods while QM calculations result in the global minima
corresponding to water molecule H-bonded to one acceptor and one donor site of
the base. Calculations for trimethylated bases with a water molecule
corroborate the MM results. The energy profiles were obtained with some degrees
of freedom of the water molecule being frozen. This data will contribute to the
improvement of the molecular mechanics force fields.


  Phylogenetic networks are used to display the relationship of different
species whose evolution is not treelike, which is the case, for instance, in
the presence of hybridization events or horizontal gene transfers. Tree
inference methods such as Maximum Parsimony need to be modified in order to be
applicable to networks. In this paper, we discuss two different definitions of
Maximum Parsimony on networks, "hardwired" and "softwired", and examine the
complexity of computing them given a network topology and a character. By
exploiting a link with the problem Multicut, we show that computing the
hardwired parsimony score for 2-state characters is polynomial-time solvable,
while for characters with more states this problem becomes NP-hard but is still
approximable and fixed parameter tractable in the parsimony score. On the other
hand we show that, for the softwired definition, obtaining even weak
approximation guarantees is already difficult for binary characters and
restricted network topologies, and fixed-parameter tractable algorithms in the
parsimony score are unlikely. On the positive side we show that computing the
softwired parsimony score is fixed-parameter tractable in the level of the
network, a natural parameter describing how tangled reticulate activity is in
the network. Finally, we show that both the hardwired and softwired parsimony
score can be computed efficiently using Integer Linear Programming. The
software has been made freely available.


  Traveling wavetrains in generalized two-species predator-prey models and
two-component reaction-diffusion equations are considered. The stability of the
fixed points of the traveling wave ODEs (in the usual "spatial" variable) is
considered. For general functional forms of the nonlinear prey birthrate/prey
deathrate or reaction terms, a Hopf bifurcation is shown to occur at two
different critical values of the traveling wave speed. The post-bifurcation
dynamics is investigated for five different functional forms of the
nonlinearities. In cases where the bifurcation is supercritical, the
post-bifurcation behaviour yields stable periodic orbits of the traveling-wave
ODEs in the spatial variable. These correspond to stable periodic wavetrains of
the full PDEs. Subcritical Hopf bifurcations yield more complex
post-bifurcation dynamics in the PDE wavetrains. In special cases where the
subcritical bifurcation marks the end of the regime of stability, the
post-bifurcation behavior in the spatial ODEs is chaotic, corresponding to
wavetrains of the original PDEs which are spatially coherent, but have chaotic
temporal dynamics. All the models are integrated numerically to investigate the
post-bifurcation dynamics and chaotic regimes are characterized by computing
power spectra, autocorrelation functions, and fractal dimensions.


  Intracellular transport based on molecular motors and its regulation are
crucial to the functioning of cells. Filamentary tracks of the cells are
abundantly decorated with non-motile microtubule-associated proteins, such as
tau. Motivated by experiments on kinesin-tau interactions [Dixit et al. Science
319, 1086 (2008)] we developed a stochastic model of interacting single-headed
motor proteins KIF1A that also takes into account the interactions between
motor proteins and tau molecules. Our model reproduce experimental observations
and predicts significant effects of tau on bound time and run length which
suggest an important role of tau in regulation of kinesin-based transport.


  Hybrid models of chemotaxis combine agent-based models of cells with partial
differential equation models of extracellular chemical signals. In this paper,
travelling wave properties of hybrid models of bacterial chemotaxis are
investigated. Bacteria are modelled using an agent-based (individual-based)
approach with internal dynamics describing signal transduction. In addition to
the chemotactic behaviour of the bacteria, the individual-based model also
includes cell proliferation and death. Cells consume the extracellular nutrient
field (chemoattractant) which is modelled using a partial differential
equation. Mesoscopic and macroscopic equations representing the behaviour of
the hybrid model are derived and the existence of travelling wave solutions for
these models is established. It is shown that cell proliferation is necessary
for the existence of non-transient (stationary) travelling waves in hybrid
models. Additionally, a numerical comparison between the wave speeds of the
continuum models and the hybrid models shows good agreement in the case of weak
chemotaxis and qualitative agreement for the strong chemotaxis case. In the
case of slow cell adaptation, we detect oscillating behaviour of the wave,
which cannot be explained by mean-field approximations.


  FoF1-ATP synthase is the membrane protein catalyzing the synthesis of the
'biological energy currency' adenosine triphosphate (ATP). The enzyme uses
internal subunit rotation for the mechanochemical conversion of a proton motive
force to the chemical bond. We apply single-molecule F\"orster resonance energy
transfer (FRET) to monitor subunit rotation in the two coupled motors F1 and
Fo. Therefore, enzymes have to be isolated from the plasma membranes of
Escherichia coli, fluorescently labeled and reconstituted into 120-nm sized
lipid vesicles to yield proteoliposomes. These freely diffusing proteoliposomes
occasionally traverse the confocal detection volume resulting in a burst of
photons. Conformational dynamics of the enzyme are identified by sequential
changes of FRET efficiencies within a single photon burst. The observation
times can be extended by capturing single proteoliposomes in an anti-Brownian
electrokinetic trap (ABELtrap, invented by A. E. Cohen and W. E. Moerner). Here
we describe the preparation procedures of FoF1-ATP synthase and simulate FRET
efficiency trajectories for 'trapped' proteoliposomes. Hidden Markov Models are
applied at signal-to-background ratio limits for identifying the dwells and
substeps of the rotary enzyme when running at low ATP concentrations, excited
by low laser power, and confined by the ABELtrap.


  Identifying drivers of complex traits from the noisy signals of genetic
variation obtained from high throughput genome sequencing technologies is a
central challenge faced by human geneticists today. We hypothesize that the
variants involved in complex diseases are likely to exhibit non-neutral
evolutionary signatures. Uncovering the evolutionary history of all variants is
therefore of intrinsic interest for complex disease research. However, doing so
necessitates the simultaneous elucidation of the targets of natural selection
and population-specific demographic history. Here we characterize the action of
natural selection operating across complex disease categories, and use
population genetic simulations to evaluate the expected patterns of genetic
variation in large samples. We focus on populations that have experienced
historical bottlenecks followed by explosive growth (consistent with most human
populations), and describe the differences between evolutionarily deleterious
mutations and those that are neutral. Genes associated with several complex
disease categories exhibit stronger signatures of purifying selection than
non-disease genes. In addition, loci identified through genome-wide association
studies of complex traits also exhibit signatures consistent with being in
regions recurrently targeted by purifying selection. Through simulations, we
show that population bottlenecks and rapid growth enables deleterious rare
variants to persist at low frequencies just as long as neutral variants, but
low frequency and common variants tend to be much younger than neutral
variants. This has resulted in a large proportion of modern-day rare alleles
that have a deleterious effect on function, and that potentially contribute to
disease susceptibility.


  Due to inherent noise in intracellular networks cellular decisions can be
random, so genetically identical cells can display different phenotypic
behavior even in identical environments. Most previous work in understanding
the decision-making process has focused on the role of intrinsic noise in these
systems. Yet, especially in the high copy-number regime, extrinsic noise has
been shown to be much more significant. Here, using a prototypical example of a
bistable self-regulating gene model, we develop a theoretical framework
describing the combined effect of intrinsic and extrinsic noise on the dynamics
of stochastic genetic switches. Employing our theory and Monte Carlo
simulations, we show that extrinsic noise not only significantly alters the
lifetimes of the phenotypic states, but can induce bistability in unexpected
regions of parameter space, and may fundamentally change the escape mechanism.
These results have implications for interpreting experimentally observed
heterogeneity in cellular populations and for stochastic modeling of cellular
decision processes.


  In order to analyze data from cancer genome sequencing projects, we need to
be able to distinguish causative, or "driver," mutations from "passenger"
mutations that have no selective effect. Toward this end, we prove results
concerning the frequency of neutural mutations in exponentially growing
multitype branching processes that have been widely used in cancer modeling.
Our results yield a simple new population genetics result for the site
frequency spectrum of a sample from an exponentially growing population.


  Living systems often exhibit internal driving: active, molecular processes
drive nonequilibrium phenomena such as metabolism or migration. Active gels
constitute a fascinating class of internally driven matter, where molecular
motors exert localized stresses inside polymer networks. There is evidence that
network crosslinking is required to allow motors to induce macroscopic
contraction. Yet a quantitative understanding of how network connectivity
enables contraction is lacking. Here we show experimentally that myosin motors
contract crosslinked actin polymer networks to clusters with a scale-free size
distribution. This critical behavior occurs over an unexpectedly broad range of
crosslink concentrations. To understand this robustness, we develop a
quantitative model of contractile networks that takes into account network
restructuring: motors reduce connectivity by forcing crosslinks to unbind.
Paradoxically, to coordinate global contractions, motor activity should be low.
Otherwise, motors drive initially well-connected networks to a critical state
where ruptures form across the entire network.


  The tempo and mode of an adaptive process is strongly determined by the
structure of the fitness landscape that underlies it. In order to be able to
predict evolutionary outcomes (even on the short term), we must know more about
the nature of realistic fitness landscapes than we do today. For example, in
order to know whether evolution is predominantly taking paths that move upwards
in fitness and along neutral ridges, or else entails a significant number of
valley crossings, we need to be able to visualize these landscapes: we must
determine whether there are peaks in the landscape, where these peaks are
located with respect to one another, and whether evolutionary paths can connect
them. This is a difficult task because genetic fitness landscapes (as opposed
to those based on traits) are high-dimensional, and tools for visualizing such
landscapes are lacking. In this contribution, we focus on the predictability of
evolution on rugged genetic fitness landscapes, and determine that peaks in
such landscapes are highly clustered: high peaks are predominantly close to
other high peaks. As a consequence, the valleys separating such peaks are
shallow and narrow, such that evolutionary trajectories towards the highest
peak in the landscape can be achieved via a series of valley crossings


  In contrast with animal communication systems, diversity is characteristic of
almost every aspect of human language. Languages variously employ tones,
clicks, or manual signs to signal differences in meaning; some languages lack
the noun-verb distinction (e.g., Straits Salish), whereas others have a
proliferation of fine-grained syntactic categories (e.g., Tzeltal); and some
languages do without morphology (e.g., Mandarin), while others pack a whole
sentence into a single word (e.g., Cayuga). A challenge for evolutionary
biology is to reconcile the diversity of languages with the high degree of
biological uniformity of their speakers. Here, we model processes of language
change and geographical dispersion and find a consistent pressure for flexible
learning, irrespective of the language being spoken. This pressure arises
because flexible learners can best cope with the observed high rates of
linguistic change associated with divergent cultural evolution following human
migration. Thus, rather than genetic adaptations for specific aspects of
language, such as recursion, the coevolution of genes and fast-changing
linguistic structure provides the biological basis for linguistic diversity.
Only biological adaptations for flexible learning combined with cultural
evolution can explain how each child has the potential to learn any human
language.


  We developed a novel detection method for osteopontin (OPN), a new biomarker
for prostate cancer, by attaching a genetically engineered single chain
variable fragment (scFv) protein with high binding affinity for OPN to a carbon
nanotube field-effect transistor (NTFET). Chemical functionalization using
diazonium salts is used to covalently attach scFv to NT-FETs, as confirmed by
atomic force microscopy, while preserving the activity of the biological
binding site for OPN. Electron transport measurements indicate that
functionalized NT-FET may be used to detect the binding of OPN to the
complementary scFv protein. A concentration-dependent increase in the
source-drain current is observed in the regime of clinical significance, with a
detection limit of approximately 30 fM. The scFv-NT hybrid devices exhibit
selectivity for OPN over other control proteins. These devices respond to the
presence of OPN in a background of concentrated bovine serum albumin, without
loss of signal. Based on these observations, the detection mechanism is
attributed to changes in scattering at scFv protein-occupied defect sites on
the carbon nanotube sidewall. The functionalization procedure described here is
expected to be generalizable to any antibody containing an accessible amine
group, and to result in biosensors appropriate for detection of corresponding
complementary proteins at fM concentrations.


  We examined the potential of antibody-functionalized single-walled carbon
nanotube (SWNT) field-effect transistors (FETs) for use as a fast and accurate
sensor for a Lyme disease antigen. Biosensors were fabricated on oxidized
silicon wafers using chemical vapor deposition grown carbon nanotubes that were
functionalized using diazonium salts. Attachment of Borrelia burgdorferi (Lyme)
flagellar antibodies to the nanotubes was verified by Atomic Force Microscopy
and electronic measurements. A reproducible shift in the turn-off voltage of
the semiconducting SWNT FETs was seen upon incubation with Borrelia burgdorferi
flagellar antigen, indicative of the nanotube FET being locally gated by the
residues of flagellar protein bound to the antibody. This sensor effectively
detected antigen in buffer at concentrations as low as 1 ng/ml, and the
response varied strongly over a concentration range coinciding with levels of
clinical interest. Generalizable binding chemistry gives this biosensing
platform the potential to be expanded to monitor other relevant antigens,
enabling a multiple vector sensor for Lyme disease. The speed and sensitivity
of this biosensor make it an ideal candidate for development as a medical
diagnostic test.


  Phylogenetic analyses of gene expression have great potential for addressing
a wide range of questions. These analyses will, for example, identify genes
that have evolutionary shifts in expression that are correlated with
evolutionary changes in morphological, physiological, and developmental
characters of interest. This will provide entirely new opportunities to
identify genes related to particular phenotypes. There are, however, three key
challenges that must be addressed for such studies to realize their potential.
First, gene expression data must be measured from multiple species, some of
which may be field collected, and parameterized in such a way that they can be
compared across species. Second, it will be necessary to develop phylogenetic
comparative methods suitable for large multidimensional datasets. In most
phylogenetic comparative studies to date, the number n of independent
observations (independent contrasts) has been greater than the number p of
variables (characters). The behavior of comparative methods for these classic
n>p problems are now well understood under a wide variety of conditions. In
gene expression studies, and studies based on other high-throughput tools, the
number n of samples is dwarfed by the number p of variables. The estimated
covariance matrices will be singular, complicating their analysis and
interpretation, and prone to spurious results. Third, new approaches are needed
to investigate the expression of the many genes whose phylogenies are not
congruent with species phylogenies due to gene loss, gene duplication, and
incomplete lineage sorting. Here we outline general project design
considerations for phylogenetic analyses of gene expression, and suggest
solutions to these three categories of challenges. These topics are relevant to
high-throughput phenotypic data well beyond gene expression.


  Flocks of starlings exhibit a remarkable ability to maintain cohesion as a
group in highly uncertain environments and with limited, noisy information.
Recent work demonstrated that individual starlings within large flocks respond
to a fixed number of nearest neighbors, but until now it was not understood why
this number is seven. We analyze robustness to uncertainty of consensus in
empirical data from multiple starling flocks and show that the flock
interaction networks with six or seven neighbors optimize the trade-off between
group cohesion and individual effort. We can distinguish these numbers of
neighbors from fewer or greater numbers using our systems-theoretic approach to
measuring robustness of interaction networks as a function of the network
structure, i.e., who is sensing whom. The metric quantifies the disagreement
within the network due to disturbances and noise during consensus behavior and
can be evaluated over a parameterized family of hypothesized sensing strategies
(here the parameter is number of neighbors). We use this approach to further
show that for the range of flocks studied the optimal number of neighbors does
not depend on the number of birds within a flock; rather, it depends on the
shape, notably the thickness, of the flock. The results suggest that robustness
to uncertainty may have been a factor in the evolution of flocking for
starlings. More generally, our results elucidate the role of the interaction
network on uncertainty management in collective behavior, and motivate the
application of our approach to other biological networks.


  Regulating physical size is an essential problem that biological organisms
must solve from the subcellular to the organismal scales, but it is not well
understood what physical principles and mechanisms organisms use to sense and
regulate their size. Any biophysical size-regulation scheme operates in a noisy
environment and must be robust to other cellular dynamics and fluctuations.
This work develops theory of filament length regulation inspired by recent
experiments on kinesin-8 motor proteins, which move with directional bias on
microtubule filaments and alter microtubule dynamics. Purified kinesin-8 motors
can depolymerize chemically-stabilized microtubules. In the length-dependent
depolymerization model, the rate of depolymerization tends to increase with
filament length, because long filaments accumulate more motors at their tips
and therefore shorten more quickly. When balanced with a constant filament
growth rate, this mechanism can lead to a fixed polymer length. However, the
mechanism by which kinesin-8 motors affect the length of dynamic microtubules
in cells is less clear. We study the more biologically realistic problem of
microtubule dynamic instability modulated by a motor-dependent increase in the
filament catastrophe frequency. This leads to a significant decrease in the
mean filament length and a narrowing of the filament length distribution. The
results improve our understanding of the biophysics of length regulation in
cells.


  We apply our general method of duality, introduced in [Giardina', Kurchan,
Redig, J. Math. Phys. 48, 033301 (2007)], to models of population dynamics. The
classical dualities between forward and ancestral processes can be viewed as a
change of representation in the classical creation and annihilation operators,
both for diffusions dual to coalescents of Kingman's type, as well as for
models with finite population size. Next, using SU(1,1) raising and lowering
operators, we find new dualities between the Wright-Fisher diffusion with $d$
types and the Moran model, both in presence and absence of mutations. These new
dualities relates two forward evolutions. From our general scheme we also
identify self-duality of the Moran model.


  The last twenty years have witnessed molecular data emerge as a primary
research instrument in most branches of mycology. Fungal systematics, taxonomy,
and ecology have all seen tremendous progress and have undergone rapid,
far-reaching changes as disciplines in the wake of continual improvement in DNA
sequencing technology. A taxonomic study that draws from molecular data
involves a long series of steps, ranging from taxon sampling through the
various laboratory procedures and data analysis to the publication process. All
steps are important and influence the results and the way they are perceived by
the scientific community. The present paper provides a reflective overview of
all major steps in such a project with the purpose to assist research students
about to begin their first study using DNA-based methods. We also take the
opportunity to discuss the role of taxonomy in biology and the life sciences in
general in the light of molecular data. While the best way to learn molecular
methods is to work side by side with someone experienced, we hope that the
present paper will serve to lower the learning threshold for the reader.


  In birds, diverticula of the lungs and air sacs pneumatize specific regions
of the postcranial skeleton. Relationships among pulmonary components and
skeletal regions they pneumatize allow inferences about pulmonary anatomy in
non-avian dinosaurs. Fossae, foramina and chambers in the postcranial skeletons
of pterosaurs and saurischian dinosaurs are diagnostic for pneumaticity. In
basal saurischians only the cervical skeleton is pneumatized, by cervical air
sacs. In more derived saurischians pneumatization of posterior dorsal, sacral,
and caudal vertebrae indicates abdominal air sacs. Abdominal air sacs in
sauropods are also indicated by a pneumatic hiatus (a gap in vertebral
pneumatization) in Haplocanthosaurus. Minimally, saurischians had dorsally
attached diverticular lungs plus anterior and posterior air sacs, and all the
pulmonary prerequisites for flow-through lung ventilation like that of birds.
Pneumaticity reduced skeletal mass in saurischians. I propose the Air Space
Proportion (ASP) as a measure of proportional volume of air in pneumatic bones.
The mean ASP of a sample of sauropod and theropod vertebrae is 0.61, so on
average, air occupied more than half the volume of these vertebrae. In
Diplodocus, pneumatization lightened the living animal by 7-10 percent, and
that does not include extraskeletal diverticula, air sacs, lungs, or trachea.
If all these air reservoirs included, the specific gravity of Diplodocus is
0.80, higher than published values for birds but lower than those for squamates
and crocodilians. Pneumatization of cervical vertebrae facilitated evolution of
long necks in sauropods. Necks longer than nine meters evolved at least four
times, in mamenchisaurs, diplodocids, brachiosaurids, and titanosaurs.
Increases in the number of cervical vertebrae, their proportional lengths, and
their internal complexity occurred in parallel in most of these lineages.


  Populations can be genetically isolated both by geographic distance and by
differences in their ecology or environment that decrease the rate of
successful migration. Empirical studies often seek to investigate the
relationship between genetic differentiation and some ecological variable(s)
while accounting for geographic distance, but common approaches to this problem
(such as the partial Mantel test) have a number of drawbacks. In this article,
we present a Bayesian method that enables users to quantify the relative
contributions of geographic distance and ecological distance to genetic
differentiation between sampled populations or individuals. We model the allele
frequencies in a set of populations at a set of unlinked loci as spatially
correlated Gaussian processes, in which the covariance structure is a
decreasing function of both geographic and ecological distance. Parameters of
the model are estimated using a Markov chain Monte Carlo algorithm. We call
this method Bayesian Estimation of Differentiation in Alleles by Spatial
Structure and Local Ecology (BEDASSLE), and have implemented it in a
user-friendly format in the statistical platform R. We demonstrate its utility
with a simulation study and empirical applications to human and teosinte
datasets.


  Connectivity networks have recently become widely used in biology due to
increasing amounts of information on the physical and functional links between
individual proteins. This connectivity data provides valuable material for
expanding our knowledge far beyond the experimentally validated via
mathematical analysis and theoretical predictions of new functional
interactions. In this paper we demonstrate an application of several algorithms
developed for the ranking of potential gene-expression regulators within the
context of an associated network. We analyze how different types of
connectivity between genes and proteins affect the topology of the integral
C.elegans functional network and thereby validate algorithmic performance. We
demonstrate the possible definition of co-expression gene clusters within a
network context from their specific motif distribution signatures. We also show
that the method based on the shortest path function (SPF) applied to gene
interactions sub-network of the co-expression gene cluster, efficiently
predicts novel regulatory transcription factors (TFs). Simultaneous application
of other methods, including only interactions with neighborhood genes, allows
rapid ranking of potential regulators that could be functionally linked with
the group of co-expressed genes. Predicting functions of regulators for a
cluster of ribosomal/mRNA metabolic genes we highlight a role of mRNA
translation and decay in a longevity of organisms.


  Muller's ratchet is a paradigmatic model for the accumulation of deleterious
mutations in a population of finite size. A click of the ratchet occurs when
all individuals with the least number of deleterious mutations are lost
irreversibly due to a stochastic fluctuation. In spite of the simplicity of the
model, a quantitative understanding of the process remains an open challenge.
In contrast to previous works, we here study a Moran model of the ratchet with
overlapping generations. Employing an approximation which describes the fittest
individuals as one class and the rest as a second class, we obtain closed
analytical expressions of the ratchet rate in the rare clicking regime. As a
click in this regime is caused by a rare large fluctuation from a metastable
state, we do not resort to a diffusion approximation but apply an approximation
scheme which is especially well suited to describe extinction events from
metastable states. This method also allows for a derivation of expressions for
the quasi-stationary distribution of the fittest class. Additionally, we
confirm numerically that the formulation with overlapping generations leads to
the same results as the diffusion approximation and the corresponding
Wright-Fisher model with non-overlapping generations.


  The impact of the presence of risk of destructive event on the silvicultural
practice of a forest stand is investigated. For that, we consider a model of
population dynamics. This model has allowed us to make the comparison without
and with risk, and highlight the influence of the presence of risk of
destructive event on optimal thinning and optimal rotation period.


  Collective actions, from city marathons to labor strikes, are often
mass-driven and subject to the snowball effect. Motivated by this, we study
evolutionary advantages of conditional punishment in the spatial public goods
game. Unlike unconditional punishers who always impose the same fines on
defectors, conditional punishers do so proportionally with the number of other
punishers in the group. Phase diagrams in dependence on the punishment fine and
cost reveal that the two types of punishers cannot coexist. Spontaneous
coarsening of the two strategies leads to an indirect territorial competition
with the defectors, which is won by unconditional punishers only if the
sanctioning is inexpensive. Otherwise conditional punishers are the victors of
the indirect competition, indicating that under more realistic conditions they
are indeed the more effective strategy. Both continuous and discontinuous phase
transitions as well as tricritical points characterize the complex evolutionary
dynamics, which is due to multipoint interactions that are introduced by
conditional punishment. We propose indirect territorial competition as a
generally applicable mechanism relying on pattern formation, by means of which
spatial structure can be utilized by seemingly subordinate strategies to avoid
evolutionary extinction.


  RNA sequencing (RNA-seq) has been rapidly adopted for the profiling of
transcriptomes in many areas of biology, including studies into gene
regulation, development and disease. Of particular interest is the discovery of
differentially expressed genes across different conditions (e.g., tissues,
perturbations), while optionally adjusting for other systematic factors that
affect the data collection process. There are a number of subtle yet critical
aspects of these analyses, such as read counting, appropriate treatment of
biological variability, quality control checks and appropriate setup of
statistical modeling. Several variations have been presented in the literature,
and there is a need for guidance on current best practices. This protocol
presents a "state-of-the-art" computational and statistical RNA-seq
differential expression analysis workflow largely based on the free open-source
R language and Bioconductor software and in particular, two widely-used tools
DESeq and edgeR. Hands-on time for typical small experiments (e.g., 4-10
samples) can be <1 hour, with computation time <1 day using a standard desktop
PC.


  Quantifying diversity is of central importance for the study of structure,
function and evolution of microbial communities. The estimation of microbial
diversity has received renewed attention with the advent of large-scale
metagenomic studies. Here, we consider what the diversity observed in a sample
tells us about the diversity of the community being sampled. First, we argue
that one cannot reliably estimate the absolute and relative number of microbial
species present in a community without making unsupported assumptions about
species abundance distributions. The reason for this is that sample data do not
contain information about the number of rare species in the tail of species
abundance distributions. We illustrate the difficulty in comparing species
richness estimates by applying Chao's estimator of species richness to a set of
in silico communities: they are ranked incorrectly in the presence of large
numbers of rare species. Next, we extend our analysis to a general family of
diversity metrics ("Hill diversities"), and construct lower and upper estimates
of diversity values consistent with the sample data. The theory generalizes
Chao's estimator, which we retrieve as the lower estimate of species richness.
We show that Shannon and Simpson diversity can be robustly estimated for the in
silico communities. We analyze nine metagenomic data sets from a wide range of
environments, and show that our findings are relevant for empirically-sampled
communities. Hence, we recommend the use of Shannon and Simpson diversity
rather than species richness in efforts to quantify and compare microbial
diversity.


  It is a challenge to predict the response of a large, complex system to a
perturbation. Recent attempts to predict the behaviour of food webs have
revealed that the effort needed to understand a system grows quickly with its
complexity, because increasingly precise information on the elements of the
system is required. Here, we show that not all elements of the system need to
be measured equally well. This suggests that a more efficient allocation of
effort to understand a complex systems is possible. We develop an iterative
technique for determining an efficient measurement strategy. Finally, in our
assessment of model food webs, we find that it is most important to precisely
measure the mortality and predation rates of long-lived, generalist, top
predators. Prioritizing the study of such species will make it easier to
understand the response of complex food webs to perturbations.


  The presence of different transcripts of a gene across samples can be
analysed by whole-transcriptome microarrays. Reproducing results from published
microarray data represents a challenge due to the vast amounts of data and the
large variety of pre-processing and filtering steps employed before the actual
analysis is carried out. To guarantee a firm basis for methodological
development where results with new methods are compared with previous results
it is crucial to ensure that all analyses are completely reproducible for other
researchers. We here give a detailed workflow on how to perform reproducible
analysis of the GeneChip Human Exon 1.0 ST Array at probe and probeset level
solely in R/Bioconductor, choosing packages based on their simplicity of use.
To exemplify the use of the proposed workflow we analyse differential splicing
and differential gene expression in a publicly available dataset using various
statistical methods. We believe this study will provide other researchers with
an easy way of accessing gene expression data at different annotation levels
and with the sufficient details needed for developing their own tools for
reproducible analysis of the GeneChip Human Exon 1.0 ST Array.


  The word-stock of a language is a complex dynamical system in which words can
be created, evolve, and become extinct. Even more dynamic are the short-term
fluctuations in word usage by individuals in a population. Building on the
recent demonstration that word niche is a strong determinant of future rise or
fall in word frequency, here we introduce a model that allows us to distinguish
persistent from temporary increases in frequency. Our model is illustrated
using a 10^8-word database from an online discussion group and a 10^11-word
collection of digitized books. The model reveals a strong relation between
changes in word dissemination and changes in frequency. Aside from their
implications for short-term word frequency dynamics, these observations are
potentially important for language evolution as new words must survive in the
short term in order to survive in the long term.


  This paper deals with the stability properties of a closed market, where
capital and labour force are acting like a predator-prey system in
population-dynamics. The spatial movement of the capital and labour force are
taken into account by cross-diffusion effect. First, we are showing two
possible ways for modeling this system in only one country's market (applying a
simple functional response and a Holling-type ratio-dependent response as
well), examining the conditions of their stability properties. We extend the
ratio-dependent model into two countries common market where two kind of
cross-diffusion effects are present, and find those additional conditions,
whose are necessary for the stability of the global common market besides the
stability of each countries local markets. Our four-dimensional model
highlights that a hectic movement of the capital toward labour force can cause
a Turing instability.


  Background: The development, optimization and validation of protein modeling
methods require efficient tools for structural comparison. Frequently, a large
number of models need to be compared with the target native structure. The main
reason for the development of Clusco software was to create a high-throughput
tool for all-versus-all comparison, because calculating similarity matrix is
the one of the bottlenecks in the protein modeling pipeline. Results: Clusco is
fast and easy-to-use software for high-throughput comparison of protein models
with different similarity measures (cRMSD, dRMSD, GDT_TS, TM-Score, MaxSub,
Contact Map Overlap) and clustering of the comparison results with standard
methods: K-means Clustering or Hierarchical Agglomerative Clustering.
Conclusions: The application was highly optimized and written in C/C++,
including the code for parallel execution on CPU and GPU version of cRMSD,
which resulted in a significant speedup over similar clustering and scoring
computation programs.


  Eukaryotic cells possess motility mechanisms allowing them not only to
self-propel but also to exert forces on obstacles (to push) and to carry
cargoes (to pull). To study the inherent asymmetry between active pushing and
pulling we model a crawling acto-myosin cell extract as a 1D layer of active
gel subjected to external forces. We show that pushing is controlled by
protrusion and that the macroscopic signature of the protrusion dominated
motility mechanism is concavity of the force velocity relation. Instead,
pulling is driven by protrusion only at small values of the pulling force and
it is replaced by contraction when the pulling force is sufficiently large.
This leads to more complex convex-concave structure of the force velocity
relation, in particular, competition between protrusion and contraction can
produce negative mobility in a biologically relevant range. The model
illustrates active readjustment of the force generating machinery in response
to changes in the dipole structure of external forces. The possibility of
switching between complementary active mechanisms implies that if necessary
'pushers' can replace 'pullers' and visa versa.


  Computational approaches to exploring "chemical universes", i.e., very large
sets, potentially infinite sets of compounds that can be constructed by a
prescribed collection of reaction mechanisms, in practice suffer from a
combinatorial explosion. It quickly becomes impossible to test, for all pairs
of compounds in a rapidly growing network, whether they can react with each
other. More sophisticated and efficient strategies are therefore required to
construct very large chemical reaction networks.
  Undirected labeled graphs and graph rewriting are natural models of chemical
compounds and chemical reactions. Borrowing the idea of partial evaluation from
functional programming, we introduce partial applications of rewrite rules.
Binding substrate to rules increases the number of rules but drastically prunes
the substrate sets to which it might match, resulting in dramatically reduced
resource requirements. At the same time, exploration strategies can be guided,
e.g. based on restrictions on the product molecules to avoid the explicit
enumeration of very unlikely compounds. To this end we introduce here a generic
framework for the specification of exploration strategies in graph-rewriting
systems. Using key examples of complex chemical networks from sugar chemistry
and the realm of metabolic networks we demonstrate the feasibility of a
high-level strategy framework.
  The ideas presented here can not only be used for a strategy-based chemical
space exploration that has close correspondence of experimental results, but
are much more general. In particular, the framework can be used to emulate
higher-level transformation models such as illustrated in a small puzzle game.


  Complex networks describe a wide range of systems in nature and society.
Frequently cited examples include Internet, WWW, a network of chemicals linked
by chemical reactions, social relationship networks, citation networks, etc.
The research of complex networks has attracted many scientists' attention.
Physicists have shown that these networks exhibit some surprising characters,
such as high clustering coefficient, small diameter, and the absence of the
thresholds of percolation. Scientists in mathematical epidemiology discovered
that the threshold of infectious disease disappears on contact networks that
following Scale-Free distribution. Researchers in economics and public health
also find that the imitation behavior could lead to cluster phenomena of
vaccination and un-vaccination. In this note, we will review the basic concepts
of complex networks; Basic epidemic models; the development of complex networks
and epidemiology.


  In this communication, direction of evolutionary variability of parameters of
genome size and structurally functional activity of plants in angiosperm taxa
among life forms, are analysed. It is shown that, in the Cretaceous-Cenozoic
era, the nuclear genome of the plants tended to increase. Functional genome
efficiency (intensity of functions per picogram of DNA) decreased as much as
possible from highest, at trees and lianas of rain and monsoonal forests of the
Paleogene, to minimum, at shrubs, perennial and annual grasses of meadow-steppe
vegetation appeared in the Neogene. Environmental changes in temperature,
humidity and carbon dioxide concentration in adverse direction, critical for
vegetation, are discussed as the cause of growth of evolutionary genome size
and loss in its functional efficiency. The growth of the genome in the Cenozoic
did not lead to the intensification of functions, but rather led to the
expansion of the adaptive capacity of species. Growth of nuclear
deoxyribonucleic acid content can be considered as one of the effective tools
of an adaptogenesis.


  A central question is how the conformational changes of proteins affect their
function and the inhibition of this function by drug molecules. Many enzymes
change from an open to a closed conformation upon binding of substrate or
inhibitor molecules. These conformational changes have been suggested to follow
an induced-fit mechanism in which the molecules first bind in the open
conformation in those cases where binding in the closed conformation appears to
be sterically obstructed such as for the HIV-1 protease. In this article, we
present a general model for the catalysis and inhibition of enzymes with
induced-fit binding mechanism. We derive general expressions that specify how
the overall catalytic rate of the enzymes depends on the rates for binding, for
the conformational changes, and for the chemical reaction. Based on these
expressions, we analyze the effect of mutations that mainly shift the
conformational equilibrium on catalysis and inhibition. If the overall
catalytic rate is limited by product unbinding, we find that mutations that
destabilize the closed conformation relative to the open conformation increase
the catalytic rate in the presence of inhibitors by a factor exp(ddG/RT) where
ddG is the mutation-induced shift of the free-energy difference between the
conformations. This increase in the catalytic rate due to changes in the
conformational equilibrium is independent of the inhibitor molecule and, thus,
may help to understand how non-active-site mutations can contribute to the
multi-drug-resistance that has been observed for the HIV-1 protease. A
comparison to experimental data for the non-active-site mutation L90M of the
HIV-1 protease indicates that the mutation slightly destabilizes the closed
conformation of the enzyme.


  Biological systems exhibit two structural features on many levels of
organization: sparseness, in which only a small fraction of possible
interactions between components actually occur; and modularity - the near
decomposability of the system into modules with distinct functionality. Recent
work suggests that modularity can evolve in a variety of circumstances,
including goals that vary in time such that they share the same subgoals
(modularly varying goals), or when connections are costly. Here, we studied the
origin of modularity and sparseness focusing on the nature of the mutation
process, rather than on connection cost or variations in the goal. We use
simulations of evolution with different mutation rules. We found that commonly
used sum-rule mutations, in which interactions are mutated by adding random
numbers, do not lead to modularity or sparseness except for in special
situations. In contrast, product-rule mutations in which interactions are
mutated by multiplying by random numbers - a better model for the effects of
biological mutations - led to sparseness naturally. When the goals of evolution
are modular, in the sense that specific groups of inputs affect specific groups
of outputs, product-rule mutations also lead to modular structure; sum-rule
mutations do not. Product-rule mutations generate sparseness and modularity
because they tend to reduce interactions, and to keep small interaction terms
small.


  Competition between independently arising beneficial mutations is enhanced in
spatial populations due to the linear rather than exponential growth of clones.
Recent theoretical studies have pointed out that the resulting fitness dynamics
is analogous to a surface growth process, where new layers nucleate and spread
stochastically, leading to the build up of scale-invariant roughness. This
scenario differs qualitatively from the standard view of adaptation in that the
speed of adaptation becomes independent of population size while the fitness
variance does not. Here we exploit recent progress in the understanding of
surface growth processes to obtain precise predictions for the universal,
non-Gaussian shape of the fitness distribution for one-dimensional habitats,
which are verified by simulations. When the mutations are deleterious rather
than beneficial the problem becomes a spatial version of Muller's ratchet. In
contrast to the case of well-mixed populations, the rate of fitness decline
remains finite even in the limit of an infinite habitat, provided the ratio
$U_d/s^2$ between the deleterious mutation rate and the square of the
(negative) selection coefficient is sufficiently large. Using again an analogy
to surface growth models we show that the transition between the stationary and
the moving state of the ratchet is governed by directed percolation.


  In the context of the genome-wide association studies (GWAS), one has to
solve long sequences of generalized least-squares problems; such a task has two
limiting factors: execution time --often in the range of days or weeks-- and
data management --data sets in the order of Terabytes. We present an algorithm
that obviates both issues. By pipelining the computation, and thanks to a
sophisticated transfer strategy, we stream data from hard disk to main memory
to GPUs and achieve sustained peak performance; with respect to a
highly-optimized CPU implementation, our algorithm shows a speedup of 2.6x.
Moreover, the approach lends itself to multiple GPUs and attains almost perfect
scalability. When using 4 GPUs, we observe speedups of 9x over the
aforementioned implementation, and 488x over a widespread biology library.


  The identification of transcription factor binding sites (TFBSs) on genomic
DNA is of crucial importance for understanding and predicting regulatory
elements in gene networks. TFBS motifs are commonly described by Position
Weight Matrices (PWMs), in which each DNA base pair independently contributes
to the transcription factor (TF) binding, despite mounting evidence of
interdependence between base pairs positions. The recent availability of
genome-wide data on TF-bound DNA regions offers the possibility to revisit this
question in detail for TF binding {\em in vivo}. Here, we use available fly and
mouse ChIPseq data, and show that the independent model generally does not
reproduce the observed statistics of TFBS, generalizing previous observations.
We further show that TFBS description and predictability can be systematically
improved by taking into account pairwise correlations in the TFBS via the
principle of maximum entropy. The resulting pairwise interaction model is
formally equivalent to the disordered Potts models of statistical mechanics and
it generalizes previous approaches to interdependent positions. Its structure
allows for co-variation of two or more base pairs, as well as secondary motifs.
Although models consisting of mixtures of PWMs also have this last feature, we
show that pairwise interaction models outperform them. The significant pairwise
interactions are found to be sparse and found dominantly between consecutive
base pairs. Finally, the use of a pairwise interaction model for the
identification of TFBSs is shown to give significantly different predictions
than a model based on independent positions.


  Considering a multi-pathway structure in a light-harvesting complex of
photosynthesis, we investigate the role of energy-level mismatches between
antenna molecules in transferring the absorbed energy to a reaction center. We
find a condition in which the antenna molecules faithfully play their roles:
Their effective absorption ratios are larger than those of the receiver
molecule directly coupled to the reaction center. In the absence of
energy-level mismatches and dephasing noise, there arises quantum destructive
interference between multiple paths that restricts the energy transfer. On the
other hand, the destructive interference diminishes as asymmetrically biasing
the energy-level mismatches and/or introducing quantum noise of dephasing for
the antenna molecules, so that the transfer efficiency is greatly enhanced to
near unity. Remarkably, the near-unity efficiency can be achieved at a wide
range of asymmetric energy-level mismatches. Temporal characteristics are also
optimized at the energy-level mismatches where the transfer efficiency is near
unity. We discuss these effects, in particular, for the Fenna-Matthews-Olson
complex.


  Studying all non-redundant proteins in 76 most-commonly found structural
domains, the present work attempts to decipher latent patterns that
characterize acceptable and unacceptable symmetries in residue-residue
interactions in functional proteins. We report that cutting across the
structural classes, a select set of pairwise interactions are universally
favored by geometrical and evolutionary constraints, termed 'acceptable'
structural and evolutionary tunnels, respectively. An equally small subset of
residue-residue interactions, the 'unacceptable' structural and evolutionary
tunnels, is found to be universally disliked by structural and evolutionary
constraints. Non-trivial overlapping is detected among acceptable structural
and evolutionary tunnels, as also among unacceptable structural and
evolutionary tunnels. A subset of tunnels is found to have equal relative
importance, structurally and evolutionarily, in different structural classes.
The MET-MET tunnel is detected to be universally most unacceptable by both
structural and evolutionary constraints, whereas the ASP-LEU tunnel was found
to be the closest approximation to be universally most acceptable. Residual
populations in structural and evolutionary tunnels are found to be independent
of stereochemical properties of individual residues. It is argued with examples
that tunnels are emergent features that connect extent of symmetry in
residue-residue interactions to the level of quaternary structural
organization.


  Complex spatial and temporal patterns of gene expression underlie embryo
differentiation, yet methods do not yet exist for the efficient genome-wide
determination of spatial expression patterns during development. In situ
imaging of transcripts and proteins is the gold-standard, but it is difficult
and time consuming to apply to an entire genome, even when highly automated.
Sequencing, in contrast, is fast and genome-wide, but is generally applied to
homogenized tissues, thereby discarding spatial information. It is likely that
these methods will ultimately converge, and we will be able to sequence RNAs in
situ, simultaneously determining their identity and location. As a step along
this path, we developed methods to cryosection individual blastoderm stage
Drosophila melanogaster embryos along the anterior-posterior axis and sequence
the mRNA isolated from each 25 micron slice. The spatial patterns of gene
expression we infer closely match patterns previously determined by in situ
hybridization and microscopy. We applied this method to generate a genome-wide
timecourse of spatial gene expression from shortly after fertilization through
gastrulation. We identify numerous genes with spatial patterns that have not
yet been described in the several ongoing systematic in situ based projects.
This simple experiment demonstrates the potential for combining careful
anatomical dissection with high-throughput sequencing to obtain spatially
resolved gene expression on a genome-wide scale.


  A detailed study of energy differences between the highest occupied and
lowest unoccupied molecular orbitals (HOMO-LUMO gaps) in protein systems and
water clusters is presented. Recent work questioning the applicability of
Kohn-Sham density-functional theory to proteins and large water clusters (E.
Rudberg, J. Phys.: Condens. Mat. 2012, 24, 072202) has demonstrated vanishing
HOMO-LUMO gaps for these systems, which is generally attributed to the
treatment of exchange in the functional used. The present work shows that the
vanishing gap is, in fact, an electrostatic artefact of the method used to
prepare the system. Practical solutions for ensuring the gap is maintained when
the system size is increased are demonstrated. This work has important
implications for the use of large-scale density-functional theory in
biomolecular systems, particularly in the simulation of photoemission, optical
absorption and electronic transport, all of which depend critically on
differences between energies of molecular orbitals.


  Morphologic diversity is observed across all families of viruses. Yet these
supra-molecular assemblies are produced most of the time in a spontaneous way
through complex molecular self-assembly scenarios. The modeling of these
phenomena remains a challenging problem within the emerging field of Physical
Virology. We present in this work a theoretical analysis aiming at highlighting
the particular role of configuration entropy in the control of viral particle
size distribution. Specializing this model to retroviruses like HIV-1, we
predict a new mechanism of entropic control of both RNA uptake into the viral
particle, and of the particle's size distribution. Evidence of this peculiar
behavior has been recently reported experimentally.


  Many biological regulatory systems process signals out of steady state and
respond with a physiological delay. A simple model of regulation which respects
these features shows how the ability of a delayed output to transmit
information is limited: at short times by the timescale of the dynamic input,
at long times by that of the dynamic output. We find that topologies of
maximally informative networks correspond to commonly occurring biological
circuits linked to stress response and that circuits functioning out of steady
state may exploit absorbing states to transmit information optimally.


  Using active gel theory we study theoretically the properties of the cortical
actin layer of animal cells. The cortical layer is described as a
non-equilibrium wetting film on the cell membrane. The actin density is
approximately constant in the layer and jumps to zero at its edge. The layer
thickness is determined by the ratio of the polymerization velocity and the
depolymerization rate of actin.