  Recent experiments elucidated the structure and function of the
cyanobacterial circadian oscillator, which is driven by sunlight intensity
variation and therefore by Earth's rotation. It is known that cyanobacteria
appeared about 3.5 billion years ago and that Earth's rotational speed is
continuously decreasing because of tidal friction. What is the effect of the
continuous slowdown of Earth's rotation on the operation of the cyanobacterial
oscillator? To answer this question we derived the oscillator's equation of
motion directly from experimental data, coupled it with Earth's rotation and
computed its natural periods and its resonance curve. The results show that
there are two resonance peaks of the "cyanobacterial oscillator-rotating Earth"
system, indicating that cyanobacteria used more efficiently the solar energy
during the geological period in which the day length varied from about 11 to 15
hours and make more efficient use of solar energy at the geological period
which started with a day length of 21 hours and will end at a day length of 28
hours.


  We review the mathematical formalism underlying the modelling of
stochasticity in biological systems. Beginning with a description of the system
in terms of its basic constituents, we derive the mesoscopic equations
governing the dynamics which generalise the more familiar macroscopic
equations. We apply this formalism to the analysis of two specific
noise-induced phenomena observed in biologically-inspired models. In the first
example, we show how the stochastic amplification of a Turing instability gives
rise to spatial and temporal patterns which may be understood within the linear
noise approximation. The second example concerns the spontaneous emergence of
cell polarity, where we make analytic progress by exploiting a separation of
time-scales.


  Social network analysis is now widely used to investigate the dynamics of
infectious disease spread from person to person. Vaccination dramatically
disrupts the disease transmission process on a contact network, and indeed,
sufficiently high vaccination rates can disrupt the process to such an extent
that disease transmission on the network is effectively halted. Here, we build
on mounting evidence that health behaviors - such as vaccination, and refusal
thereof - can spread through social networks through a process of complex
contagion that requires social reinforcement. Using network simulations that
model both the health behavior and the infectious disease spread, we find that
under otherwise identical conditions, the process by which the health behavior
spreads has a very strong effect on disease outbreak dynamics. This variability
in dynamics results from differences in the topology within susceptible
communities that arise during the health behavior spreading process, which in
turn depends on the topology of the overall social network. Our findings point
to the importance of health behavior spread in predicting and controlling
disease outbreaks.


  A plausible consequence of rugged energy landscapes of biomolecules is that
functionally competent folded states may not be unique, as is generally
assumed. Indeed, molecule-to-molecule variations in the dynamics of enzymes and
ribozymes under folding conditions have recently been identified in single
molecule experiments. However, systematic quantification and the structural
origin of the observed complex behavior remain elusive. Even for a relatively
simple case of isomerization dynamics in Holliday Junctions (HJs), molecular
heterogeneities persist over a long observation time (Tobs ~ 40 sec). Here,
using concepts in glass physics and complementary clustering analysis, we
provide a quantitative method to analyze the smFRET data probing the
isomerization in HJ dynamics. We show that ergodicity of HJ dynamics is
effectively broken; as a result, the conformational space of HJs is partitioned
into a folding network of kinetically disconnected clusters. While
isomerization dynamics in each cluster occurs rapidly as if the associated
conformational space is fully sampled, distinct patterns of time series
belonging to different clusters do not interconvert on Tobs. Theory suggests
that persistent heterogeneity of HJ dynamics is a consequence of internal
multiloops with varying sizes and flexibilities frozen by Mg2+ ions. An
annealing experiment using Mg2+ pulse that changes the Mg2+ cocentration from
high to low to high values lends support to this idea by explicitly showing
that interconversions can be driven among trajectories with different patterns.


  Upon entering physiological environments, nanoparticles readily assume the
form of a nanoparticle-protein corona that dictates their biological identity.
Understanding the structure and dynamics of nanoparticle-protein corona is
essential for predicting the fate, transport, and toxicity of nanomaterials in
living systems and for enabling the vast applications of nanomedicine. We
combined multiscale molecular dynamics simulations and complementary
experiments to characterize the silver nanoparticle-ubiquitin corona formation.
Specifically, ubiquitins competed with citrates for the nanoparticle surface
and bound to the particle in a specific manner. Under a high
protein/nanoparticle stoichiometry, ubiquitions formed a multi-layer corona on
the particle surface. The binding exhibited an unusual stretched-exponential
behavior, suggesting a rich kinetics originated from protein-protein,
protein-citrate, and protein-nanoparticle interactions. Furthermore, the
binding destabilized the {\alpha}-helices while increasing the {\beta}-sheets
of the proteins. Our results revealed the structural and dynamic complexities
of nanoparticle-protein corona formation and shed light on the origin of
nanotoxicity.


  In a short article submitted to ArXiv [1], Maddamsetti et al. argue that the
variation in the neutral mutation rate among genes in Escherichia coli that we
recently reported [2] might be explained by horizontal gene transfer (HGT). To
support their argument they present a reanalysis of synonymous diversity in 10
E.coli strains together with an analysis of a collection of 1,069 synonymous
mutations found in repair-deficient strains in a long-term in vitro evolution
experiment. Here we respond to this communication. Briefly, we explain that HGT
was carefully accounted for in our study by multiple independent phylogenetic
and population genetic approaches, and we show that there is no new evidence of
HGT affecting our results. We also argue that caution must be exercised when
comparing mutations from repair deficient strains to data from wild-type
strains, as these conditions are dominated by different mutational processes.
Finally, we reanalyse Maddamsetti's collection of mutations from a long-term in
vitro experiment and we report preliminary evidence of non-random variation of
the mutation rate in these repair deficient strains.


  In metapopulations, genetic variation of local populations is influenced by
the genetic content of the founders, and of migrants following establishment.
We analyse the effect of multiple paternity on genetic diversity using a model
in which the highly promiscuous marine snail Littorina saxatilis expands from a
mainland to colonise initially empty islands of an archipelago. Migrant females
carry a large number of eggs fertilised by 1 - 10 mates. We quantify the
genetic diversity of the population in terms of its heterozygosity: initially
during the transient colonisation process, and at long times when the
population has reached an equilibrium state with migration. During
colonisation, multiple paternity increases the heterozygosity by 10 - 300 % in
comparison with the case of single paternity. The equilibrium state, by
contrast, is less strongly affected: multiple paternity gives rise to 10 - 50 %
higher heterozygosity compared with single paternity. Further we find that far
from the mainland, new mutations spreading from the mainland cause bursts of
high genetic diversity separated by long periods of low diversity. This effect
is boosted by multiple paternity. We conclude that multiple paternity
facilitates colonisation and maintenance of small populations, whether or not
this is the main cause for the evolution of extreme promiscuity in Littorina
saxatilis.


  When mutation rates are low, natural selection remains effective, and
increasing the mutation rate can give rise to an increase in adaptation rate.
When mutation rates are high to begin with, however, increasing the mutation
rate may have a detrimental effect because of the overwhelming presence of
deleterious mutations. Indeed, if mutation rates are high enough: 1) adaptation
rate can become negative despite the continued availability of adaptive and/or
compensatory mutations, or 2) natural selection may be disabled because
adaptive and/or compensatory mutations -- whether established or newly-arising
-- are eroded by excessive mutation and decline in frequency. We apply these
two criteria to a standard model of asexual adaptive evolution and derive
mathematical expressions -- some new, some old in new guise -- delineating the
mutation rates under which either adaptive evolution or natural selection is
neutralized. The expressions are simple and require no \emph{a priori}
knowledge of organism- and/or environment-specific parameters. Our discussion
connects these results to each other and to previous theory, showing
convergence or equivalence of the different results in most cases.


  Characterization of protein energy landscape and conformational ensembles is
important for understanding mechanisms of protein folding and function. We
studied ensembles of bound and unbound conformations of six proteins to explore
their binding mechanisms and characterize the energy landscapes in implicit
solvent. First, results show that bound and unbound spectra often significantly
overlap. Moreover, the larger the overlap the smaller the RMSD between bound
and unbound conformational ensembles. Second, the analysis of the
unbound-to-bound changes points to conformational selection as the binding
mechanism for four of the proteins. Third, the center of the unbound spectrum
has a higher energy than the center of the corresponding bound spectrum of the
dimeric and multimeric states for most of the proteins. This suggests that the
unbound states often have larger entropy than the bound states considered
outside of the complex. Fourth, the exhaustively long minimization, making
small intra-rotamer adjustments, dramatically reduces the distance between the
centers of the bound and unbound spectra as well as the spectra extent. It
condenses unbound and bound energy levels into a thin layer at the bottom of
the energy landscape with the energy spacing that varies between 0.8-4.6 and
3.5-10.5 kcal/mol for the unbound and bound states correspondingly. At the same
time, the energy gap between the two lowest states in the full ensemble varies
between 0.9 and 12.1 kcal/mol. Finally, the analysis of protein energy
fluctuations showed that protein vibrations itself can excite the inter-state
transitions, thus supporting the conformational selection theory.


  The cargo motion in living cells transported by two species of motor protein
with different intrinsic directionality is discussed in this study. Similar to
single motor movement, cargo steps forward and backward along microtubule
stochastically. Recent experiments found that, cargo transportation by two
motor species has a memory, it does not change its direction as frequently as
expected, which means that its forward and backward step rates depends on its
previous motion trajectory. By assuming cargo has only the least memory, i.e.
its step direction depends only on the direction of its last step, two cases of
cargo motion are detailed analyzed in this study: {\bf (I)} cargo motion under
constant external load; and {\bf (II)} cargo motion in one fixed optical trap.
Due to the existence of memory, for the first case, cargo can keep moving in
the same direction for a long distance. For the second case, the cargo will
oscillate in the trap. The oscillation period decreases and the oscillation
amplitude increases with the motor forward step rates, but both of them
decrease with the trap stiffness. The most likely location of cargo, where the
probability of finding the oscillated cargo is maximum, may be the same as or
may be different with the trap center, which depends on the step rates of the
two motor species. Meanwhile, if motors are robust, i.e. their forward to
backward step rate ratios are high, there may be two such most likely
locations, located on the two sides of the trap center respectively. The
probability of finding cargo in given location, the probability of cargo in
forward/backward motion state, and various mean first passage times of cargo to
give location or given state are also analyzed.


  Experimental methods based on single particle tracking (SPT) are being
increasingly employed in the physical and biological sciences, where nanoscale
objects are visualized with high temporal and spatial resolution. SPT can probe
interactions between a particle and its environment but the price to be paid is
the absence of ensemble averaging and a consequent lack of statistics. Here we
address the benchmark question of how to accurately extract the diffusion
constant of one single Brownian trajectory. We analyze a class of estimators
based on weighted functionals of the square displacement. For a certain choice
of the weight function these functionals provide the true ensemble averaged
diffusion coefficient, with a precision that increases with the trajectory
resolution.


  In this article it is shown how optimized and dedicated microarray
experiments can be used to study the thermodynamics of DNA hybridization for a
large number of different conformations in a highly parallel fashion. In
particular, free energy penalties for mismatches are obtained in two
independent ways and are shown to be correlated with values from melting
experiments in solution reported in the literature. The additivity principle,
which is at the basis of the nearest-neighbor model, and according to which the
penalty for two isolated mismatches is equal to the sum of the independent
penalties, is thoroughly tested. Additivity is shown to break down for a
mismatch distance below 5 nt. The behavior of mismatches in the vicinity of the
helix edges, and the behavior of tandem mismatches are also investigated.
Finally, some thermodynamic outlying sequences are observed and highlighted.
These sequences contain combinations of GA mismatches. The analysis of the
microarray data reported in this article provides new insights on the DNA
hybridization parameters and can help to increase the accuracy of
hybridization-based technologies.


  Patterns of nestedness and specialization asymmetry, where specialist species
interact mainly with generalists while generalists interact with both
generalists and specialists, are often observed in mutualistic and antagonistic
bi-partite ecological networks. These have been explained in terms of the
relative abundance of species, using a null model that assigns links in
proportion to abundance, but doubts have been raised as to whether this offers
a complete explanation. In particular, host-parasite networks offer a variety
of examples in which the reverse patterns are observed.
  We propose that the link between specificity and species-richness may also be
driven by the co-evolution of hosts and parasites, as hosts allocate resources
to optimize defence against parasites, and parasites to optimize attack on
hosts. In this hypothesis, species interactions are a result of resource
allocations. This novel concept, linking together many different arguments for
network structures, is introduced through the adaptive dynamics of a simple
ecological toy system of two hosts and two parasites.
  We analyse the toy model and its functionality, demonstrating that
co-evolution leads to specialization asymmetry in networks with closely related
parasites or fast host mutation rates, but not in networks with more distantly
related species. Having constructed the toy model and tested its applicability,
our model can now be expanded to the full problem of a larger system.


  The UCSC Genome Browser is a popular tool for the exploration and analysis of
reference genomes. Mirrors of the UCSC Genome Browser and its contents exist at
multiple geographic locations, and this mirror procedure has been modified to
support genome sequences not maintained by UCSC and generated by individual
researchers. While straightforward, this procedure is lengthy and tedious and
would benefit from automation, especially when processing many genome
sequences. We present a Unix shell program that facilitates the creation of
custom instances of the UCSC Genome Browser for genome sequences not being
maintained by UCSC. It automates many steps of the browser creation process,
provides password protection for each browser instance, and automates the
creation of basic annotation tracks. As an example we generate a custom UCSC
Genome Browser for a bacterial genome obtained from a massively parallel
sequencing platform.


  High-throughput sequencing (HTS) is revolutionizing biological research by
enabling scientists to quickly and cheaply query variation at a genomic scale.
Despite the increasing ease of obtaining such data, using these data
effectively still poses notable challenges, especially for those working with
organisms without a high-quality reference genome. For every stage of analysis
- from assembly to annotation to variant discovery - researchers have to
distinguish technical artifacts from the biological realities of their data
before they can make inference. In this work, I explore these challenges by
generating a large de novo comparative transcriptomic dataset data for a clade
of lizards and constructing a pipeline to analyze these data. Then, using a
combination of novel metrics and an externally validated variant data set, I
test the efficacy of my approach, identify areas of improvement, and propose
ways to minimize these errors. I find that with careful data curation, HTS can
be a powerful tool for generating genomic data for non-model organisms.


  Biodiversity conservation becoming increasingly urgent. It is important to
find mechanisms of competitive coexistence of species with different fitness in
especially difficult circumstances - on one limiting resource, in isolated
stable uniform habitat, without any trade-offs and cooperative interactions.
Here we show a mechanism of competitive coexistence based on a soliton-like
behaviour of population waves. We have modelled it by the logical axiomatic
deterministic individual-based cellular automata method. Our mechanistic models
of population and ecosystem dynamics are of white-box type and so they provide
direct insight into mechanisms under study. The mechanism provides indefinite
coexistence of two, three and four competing species. This mechanism violates
the known formulations of the competitive exclusion principle. As a
consequence, we have proposed a fully mechanistic and most stringent
formulation of the principle.


  A coarse-grained computational model is used to investigate how the bending
rigidity of a polymer under tension affects the formation of a trefoil knot.
Thermodynamic integration techniques are applied to demonstrate that the
free-energy cost of forming a knot has a minimum at non-zero bending rigidity.
The position of the minimum exhibits a power-law dependence on the applied
tension. For knotted polymers with non-uniform bending rigidity, the knots
preferentially localize in the region with a bending rigidity that minimizes
the free-energy.


  Spatial structure has a profound effect on the outcome of evolutionary games.
In the ultimatum game, it leads to the dominance of much fairer players than
those predicted to evolve in well-mixed settings. Here we show that spatiality
leads to fair ultimatums only if the intervals from which the players are able
to choose how much to offer and how little to accept are sufficiently
fine-grained. Small sets of discrete strategies lead to the stable coexistence
of the two most rational strategies in the set, while larger sets lead to the
dominance of a single yet not necessarily the fairest strategy. The fairest
outcome is obtained for the most accurate strategy imitation, that is in the
limit of a continuous strategy set. Having a multitude of choices is thus
crucial for the evolution of fairness, but not necessary for the evolution of
empathy.


  Accurate identification of haplotypes in sequenced human genomes can provide
invaluable information about population demography and fine-scale correlations
along the genome, thus empowering both population genomic and medical
association studies. Yet phasing unrelated individuals remains a challenging
problem. Incorporating available data from high throughput sequencing into
traditional statistical phasing approaches is a promising avenue to alleviate
these issues. We present a novel statistical method that expands on an existing
graphical haplotype reconstruction method (shapeIT) to incorporate phasing
information from paired-end read data. The algorithm harnesses the haplotype
graph information estimated by shapeIT from genotypes across the population and
refines haplotype likelihoods for a given individual to be compatible with the
sequencing data. Applying the method to HapMap individuals genotyped on the
Affymetrix Axiom chip at 7,745,081 SNPs and on a trio sequenced by Complete
Genomics, we found that the inclusion of paired end read data significantly
improved phasing, with reductions in switch error on the order of 4-15% against
shapeIT across all panels. As expected, the improvements were found to be most
significant at sites harboring rare variants; furthermore, we found that longer
read sizes and higher throughput translated to greater decreases in switching
error, as did higher variance in the size of the insert separating the two
reads--suggesting that multi-platform next generation sequencing may be
exploited to yield particularly accurate haplotypes. Overall, the phasing
improvements afforded by this new method highlight the power of integrating
sequencing read information and population genotype data for reconstructing
haplotypes in unrelated individuals.


  We review and introduce a generalized reaction-diffusion approach to epidemic
spreading in a metapopulation modeled as a complex network. The metapopulation
consists of susceptible and infected individuals that are grouped in
subpopulations symbolising cities and villages that are coupled by human travel
in a transportation network. By analytic methods and numerical simulations we
calculate the fraction of infected people in the metaopoluation in the long
time limit, as well as the relevant parameters characterising the epidemic
threshold that separates an epidemic from a non-epidemic phase. Within this
model, we investigate the effect of a heterogeneous network topology and a
heterogeneous subpopulation size distribution. Such a system is suited for
epidemic modeling where small villages and big cities exist simultaneously in
the metapopulation. We find that the heterogeneous conditions cause the
epidemic threshold to be a non-trivial function of the reaction rates (local
parameters), the network's topology (global parameters) and the cross-over
population size that separates "village dynamics" from "city dynamics".


  The high-throughput data generated by microarray experiments provides
complete set of genes being expressed in a given cell or in an organism under
particular conditions. The analysis of these enormous data has opened a new
dimension for the researchers. In this paper we describe a novel algorithm to
microarray data analysis focusing on the identification of genes that are
differentially expressed in particular internal or external conditions and
which could be potential drug targets. The algorithm uses the time-series gene
expression data as an input and recognizes genes which are expressed
differentially. This algorithm implements standard statistics-based gene
functional investigations, such as the log transformation, mean, log-sigmoid
function, coefficient of variations, etc. It does not use clustering analysis.
The proposed algorithm has been implemented in Perl. The time-series gene
expression data on yeast Saccharomyces cerevisiae from the Stanford Microarray
Database (SMD)consisting of 6154 genes have been taken for the validation of
the algorithm. The developed method extracted 48 genes out of total 6154 genes.
These genes are mostly responsible for the yeast's resistants at a high
temperature.


  The last decade has seen the advent and consolidation of ontology based tools
for the identification and biological interpretation of classes of genes, such
as the Gene Ontology. The information accumulated time-by-time and included in
the GO is encoded in the definition of terms and in the setting up of semantic
relations amongst terms. This approach might be usefully complemented by a
bottom-up approach based on the knowledge of relationships amongst genes. To
this end, we investigate the Gene Ontology from a complex network perspective.
We consider the semantic network of terms naturally associated with the
semantic relationships provided by the Gene Ontology consortium and a
gene-based weighted network in which the nodes are the terms and a link between
any two terms is set up whenever genes are annotated in both terms. One aim of
the present paper is to understand whether the semantic and the gene-based
network share the same structural properties or not. We then consider network
communities. The identification of communities in the SVNs network can
therefore be the basis of a simple protocol aiming at fully exploiting the
possible relationships amongst terms, thus improving the knowledge of the
semantic structure of GO. This is also important from a biomedical point of
view, as it might reveal how genes over-expressed in a certain term also affect
other biological functions not directly linked by the GO semantics. As a
by-product, we present a simple methodology that allows to have a first glance
insight about the biological characterization of groups of GO terms.


  Breakage-Fusion-Bridge cycles in cancer arise when a broken segment of DNA is
duplicated and an end from each copy joined together. This structure then
`unfolds' into a new piece of palindromic DNA. This is one mechanism
responsible for the localised amplicons observed in cancer genome data. The
process has parallels with paper folding sequences that arise when a piece of
paper is folded several times and then unfolded. Here we adapt such methods to
study the breakage-fusion-bridge structures in detail. We firstly consider
discrete representations of this space with 2-d trees to demonstrate that there
are 2^(n(n-1)/2) qualitatively distinct evolutions involving n
breakage-fusion-bridge cycles. Secondly we consider the stochastic nature of
the fold positions, to determine evolution likelihoods, and also describe how
amplicons become localised. Finally we highlight these methods by inferring the
evolution of breakage-fusion-bridge cycles with data from primary tissue cancer
samples.


  Many organisms live in populations structured by space and by class, exhibit
plastic responses to their social partners, and are subject to non-additive
ecological and fitness effects. Social evolution theory has long recognized
that all of these factors can lead to different selection pressures but has
only recently attempted to synthesize how these factors interact. Using models
for both discrete and continuous phenotypes, we show that analyzing these
factors in a consistent framework reveals that they interact with one another
in ways previously overlooked. Specifically, behavioral responses
(reciprocity), genetic relatedness, and synergy interact in non-trivial ways
that cannot be easily captured by simple summary indices of assortment. We
demonstrate the importance of these interactions by showing how they have been
neglected in previous synthetic models of social behavior both within and
between species. These interactions also affect the level of behavioral
responses that can evolve in the long run; proximate biological mechanisms are
evolutionarily stable when they generate enough responsiveness relative to the
level of responsiveness that exactly balances the ecological costs and
benefits. Given the richness of social behavior across taxa, these interactions
should be a boon for empirical research as they are likely crucial for
describing the complex relationship linking ecology, demography, and social
behavior.


  Natural supports of information are given by random copolymers such as DNA or
RNA where information is coded in the sequence of covalent bonds. At the
molecular scale, the stochastic growth of a single copolymer with or without a
template proceeds by successive random attachments or detachments of monomers
continuously supplied by the surrounding solution. The thermodynamics of
copolymerization shows that fundamental links already exist between information
and thermodynamics at the molecular scale, which opens new perspectives to
understand the dynamical aspects of information in biology.


  Well-Tempered Metadynamics (WTmetaD) is an efficient method to enhance the
reconstruction of the free-energy surface of proteins. WTmetaD guarantees a
faster convergence in the long time limit in comparison with the standard
metadynamics. It still suffers however from the same limitation, i.e. the non
trivial choice of pertinent collective variables (CVs). To circumvent this
problem, we couple WTmetaD with a set of CVs generated from a dihedral
Principal Component Analysis (dPCA) on the Ramachadran dihedral angles
describing the backbone structure of the protein. The dPCA provides a generic
method to extract relevant CVs built from internal coordinates. We illustrate
the robustness of this method in the case of the small and very diffusive
Metenkephalin pentapeptide, and highlight a criterion to limit the number of
CVs necessary to biased the metadynamics simulation. The free-energy landscape
(FEL) of Met-enkephalin built on CVs generated from dPCA is found rugged
compared with the FEL built on CVs extracted from PCA of the Cartesian
coordinates of the atoms.


  Many species reproduce both sexually and asexually producing genetically
unique and clonal recruits. The relative contribution of sexual and asexual
reproduction is mostly considered a consequence of natural selection, for
example, disfavouring sexual propagules during invasions. Here, we present a
novel model for the invasion of species into a new habitat. The model is fully
neutral with respect to the dispersal and the survival of sexual and asexual
recruits. In addition to local dispersal through sexual and asexual
reproduction, long-distance dispersal is implemented through a non-zero
probability for long-range relocations. We parameterized our model using
empirical data on the distribution of asexually and sexually recruited
individuals of the recently established macroalga Fucus radicans over its
current area of distribution, the 8000 year old Baltic Sea. As a solely
stochastic alternative mechanism to natural selection, our simulation results
suggest local abundance of one sex, that during a persistent phase of
establishment will result in a heterogeneous distribution of sexual and asexual
recruits. Specifically, the model suggests an initial general dominance of one
or a few clones, and this pattern is observed in Fucus radicans as well as in
several other species with recent colonization histories.


  Error correction of sequenced reads remains a difficult task, especially in
single-cell sequencing projects with extremely non-uniform coverage. While
existing error correction tools designed for standard (multi-cell) sequencing
data usually come up short in single-cell sequencing projects, algorithms
actually used for single-cell error correction have been so far very
simplistic.
  We introduce several novel algorithms based on Hamming graphs and Bayesian
subclustering in our new error correction tool BayesHammer. While BayesHammer
was designed for single-cell sequencing, we demonstrate that it also improves
on existing error correction tools for multi-cell sequencing data while working
much faster on real-life datasets. We benchmark BayesHammer on both $k$-mer
counts and actual assembly results with the SPAdes genome assembler.


  One of the most fundamental concepts of evolutionary dynamics is the
"fixation" probability, i.e. the probability that a mutant spreads through the
whole population. Most natural communities are geographically structured into
habitats exchanging individuals among each other and can be modeled by an
evolutionary graph (EG), where directed links weight the probability for the
offspring of one individual to replace another individual in the community.
Very few exact analytical results are known for EGs. We show here how by using
the techniques of the fixed point of Probability Generating Function, we can
uncover a large class of of graphs, which we term bithermal, for which the
exact fixation probability can be simply computed.


  MicroRNAs are small noncoding RNAs that regulate genes post-transciptionally
by binding and degrading target eukaryotic mRNAs. We use a quantitative model
to study gene regulation by inhibitory microRNAs and compare it to gene
regulation by prokaryotic small non-coding RNAs (sRNAs). Our model uses a
combination of analytic techniques as well as computational simulations to
calculate the mean-expression and noise profiles of genes regulated by both
microRNAs and sRNAs. We find that despite very different molecular machinery
and modes of action (catalytic vs stoichiometric), the mean expression levels
and noise profiles of microRNA-regulated genes are almost identical to genes
regulated by prokaryotic sRNAs. This behavior is extremely robust and persists
across a wide range of biologically relevant parameters. We extend our model to
study crosstalk between multiple mRNAs that are regulated by a single microRNA
and show that noise is a sensitive measure of microRNA-mediated interaction
between mRNAs. We conclude by discussing possible experimental strategies for
uncovering the microRNA-mRNA interactions and testing the competing endogenous
RNA (ceRNA) hypothesis.


  The exclusion of freeriders from common privileges or public acceptance is
widely found in the real world. Current models on the evolution of cooperation
with incentives mostly assume peer sanctioning, whereby a punisher imposes
penalties on freeriders at a cost to itself. It is well known that such costly
punishment has two substantial difficulties. First, a rare punishing cooperator
barely subverts the asocial society of freeriders, and second, natural
selection often eliminates punishing cooperators in the presence of
non-punishing cooperators (namely, "second-order" freeriders). We present a
game-theoretical model of social exclusion in which a punishing cooperator can
exclude freeriders from benefit sharing. We show that such social exclusion can
overcome the above-mentioned difficulties even if it is costly and stochastic.
The results do not require a genetic relationship, repeated interaction,
reputation, or group selection. Instead, only a limited number of freeriders
are required to prevent the second-order freeriders from eroding the social
immune system.


  One strategy for winning a coevolutionary struggle is to evolve rapidly. Most
of the literature on host-pathogen coevolution focuses on this phenomenon, and
looks for consequent evidence of coevolutionary arms races. An alternative
strategy, less often considered in the literature, is to deter rapid
evolutionary change by the opponent. To study how this can be done, we
construct an evolutionary game between a controller that must process
information, and an adversary that can tamper with this information processing.
In this game, a species can foil its antagonist by processing information in a
way that is hard for the antagonist to manipulate. We show that the structure
of the information processing system induces a fitness landscape on which the
adversary population evolves. Complex processing logic can carve long, deep
fitness valleys that slow adaptive evolution in the adversary population. We
suggest that this type of defensive complexity on the part of the vertebrate
adaptive immune system may be an important element of coevolutionary dynamics
between pathogens and their vertebrate hosts. Furthermore, we cite evidence
that the immune control logic is phylogenetically conserved in mammalian
lineages. Thus our model of defensive complexity suggests a new hypothesis for
the lower rates of evolution for immune control logic compared to other immune
structures.


  Although hypoxia is a major stress on physiological processes, several human
populations have survived for millennia at high altitudes, suggesting that they
have adapted to hypoxic conditions. This hypothesis was recently corroborated
by studies of Tibetan highlanders, which showed that polymorphisms in candidate
genes show signatures of natural selection as well as well-replicated
association signals for variation in hemoglobin levels. We extended genomic
analysis to two Ethiopian ethnic groups: Amhara and Oromo. For each ethnic
group, we sampled low and high altitude residents, thus allowing genetic and
phenotypic comparisons across altitudes and across ethnic groups. Genome-wide
SNP genotype data were collected in these samples by using Illumina arrays. We
find that variants associated with hemoglobin variation among Tibetans or other
variants at the same loci do not influence the trait in Ethiopians. However, in
the Amhara, SNP rs10803083 is associated with hemoglobin levels at genome-wide
levels of significance. No significant genotype association was observed for
oxygen saturation levels in either ethnic group. Approaches based on allele
frequency divergence did not detect outliers in candidate hypoxia genes, but
the most differentiated variants between high- and lowlanders have a clear role
in pathogen defense. Interestingly, a significant excess of allele frequency
divergence was consistently detected for genes involved in cell cycle control,
DNA damage and repair, thus pointing to new pathways for high altitude
adaptations. Finally, a comparison of CpG methylation levels between high- and
lowlanders found several significant signals at individual genes in the Oromo.


  A common metaphor for describing development is a rugged "epigenetic
landscape" where cell fates are represented as attracting valleys resulting
from a complex regulatory network. Here, we introduce a framework for
explicitly constructing epigenetic landscapes that combines genomic data with
techniques from spin-glass physics. Each cell fate is a dynamic attractor, yet
cells can change fate in response to external signals. Our model suggests that
partially reprogrammed cells are a natural consequence of high-dimensional
landscapes, and predicts that partially reprogrammed cells should be hybrids
that co-express genes from multiple cell fates. We verify this prediction by
reanalyzing existing datasets. Our model reproduces known reprogramming
protocols and identifies candidate transcription factors for reprogramming to
novel cell fates, suggesting epigenetic landscapes are a powerful paradigm for
understanding cellular identity.


  The fruit fly Drosophila is a classic model organism to study adaptation as
well as the relationship between genetic variation and phenotypes. Although
associated bacterial communities might be important for many aspects of
Drosophila biology, knowledge about their diversity, composition, and factors
shaping them is limited. We used 454-based sequencing of a variable region of
the bacterial 16S ribosomal RNA gene to characterize the bacterial communities
associated with wild and laboratory Drosophila isolates. In order to
specifically investigate effects of food source and host species on bacterial
communities, we analyzed samples from wild Drosophila melanogaster and D.
simulans collected from a variety of natural substrates, as well as from adults
and larvae of nine laboratory reared Drosophila species. We find no evidence
for host species effects in lab reared flies, instead lab of origin and
stochastic effects, which could influence studies of Drosophila phenotypes, are
pronounced. In contrast, the natural Drosophila associated microbiota appears
to be predominantly shaped by food substrate with an additional but smaller
effect of host species identity. We identify a core member of this natural
microbiota that belongs to the genus Gluconobacter and is common to all wild
caught flies in this study, but absent from the laboratory. This makes it a
strong candidate for being part of what could be a natural D. melanogaster and
D. simulans core microbiome. Furthermore we were able to identify candidate
pathogens in natural fly isolates.


  Stable mixtures of cooperators and defectors are often seen in nature. This
fact is at odds with predictions based on linear public goods games under weak
selection. That model implies fixation either of cooperators or of defectors,
and the former scenario requires a level of group relatedness larger than the
cost/benefit ratio, being therefore expected only if there is either kin
recognition or a very low cost/benefit ratio, or else under stringent
conditions with low gene flow. This motivated us to study here social evolution
in a large class of group structured populations, with arbitrary
multi-individual interactions among group members and random migration among
groups. Under the assumption of weak selection, we analyze the equilibria and
their stability. For some significant models of social evolution with
non-linear fitness functions, including contingent behavior in iterated public
goods games and threshold models, we show that three regimes occur, depending
on the migration rate among groups. For sufficiently high migration rates, a
rare cooperative allele A cannot invade a monomorphic population of asocial
alleles N. For sufficiently low values of the migration rate, allele A can
invade the population, when rare, and then fixate, eliminating N. For
intermediate values of the migration rate, allele A can invade the population,
when rare, producing a polymorphic equilibrium, in which it coexists with N.
The equilibria and their stability do not depend on the details of the
population structure. The levels of migration (gene flow) and group relatedness
that allow for invasion of the cooperative allele leading to polymorphic
equilibria with the non-cooperative allele are common in nature.


  Though algebraic geometry over $\mathbb C$ is often used to describe the
closure of the tensors of a given size and complex rank, this variety includes
tensors of both smaller and larger rank. Here we focus on the $n\times n\times
n$ tensors of rank $n$ over $\mathbb C$, which has as a dense subset the orbit
of a single tensor under a natural group action. We construct polynomial
invariants under this group action whose non-vanishing distinguishes this orbit
from points only in its closure. Together with an explicit subset of the
defining polynomials of the variety, this gives a semialgebraic description of
the tensors of rank $n$ and multilinear rank $(n,n,n)$. The polynomials we
construct coincide with Cayley's hyperdeterminant in the case $n=2$, and thus
generalize it. Though our construction is direct and explicit, we also recast
our functions in the language of representation theory for additional insights.
  We give three applications in different directions: First, we develop basic
topological understanding of how the real tensors of complex rank $n$ and
multilinear rank $(n,n,n)$ form a collection of path-connected subsets, one of
which contains tensors of real rank $n$. Second, we use the invariants to
develop a semialgebraic description of the set of probability distributions
that can arise from a simple stochastic model with a hidden variable, a model
that is important in phylogenetics and other fields. Third, we construct simple
examples of tensors of rank $2n-1$ which lie in the closure of those of rank
$n$.


  Background: The idea that autocatalytic sets played an important role in the
origin of life is not new. However, the likelihood of autocatalytic sets
emerging spontaneously has long been debated. Recently, progress has been made
along two different lines. Experimental results have shown that autocatalytic
sets can indeed emerge in real chemical systems, and theoretical work has shown
that the existence of such self-sustaining sets is highly likely in formal
models of chemical systems. Here, we take a first step towards merging these
two lines of work by constructing and investigating a formal model of a real
chemical system of RNA replicators exhibiting autocatalytic sets.
  Results: We show that the formal model accurately reproduces recent
experimental results on an RNA replicator system, in particular how the system
goes through a sequence of larger and larger autocatalytic sets, and how a
cooperative (autocatalytic) system can outcompete an equivalent selfish system.
Moreover, the model provides additional insights that could not be obtained
from experiments alone, and it suggests several experimentally testable
hypotheses.
  Conclusions: Given these additional insights and predictions, the modeling
framework provides a better and more detailed understanding of the nature of
chemical systems in general and the emergence of autocatalytic sets in
particular. This provides an important first step in combining experimental and
theoretical work on autocatalytic sets in the context of the orgin of life.


  In living cells, proteins combine 3D bulk diffusion and 1D sliding along the
DNA to reach a target faster. This process is known as facilitated diffusion,
and we investigate its dynamics in the physiologically relevant case of
confined DNA. The confining geometry and DNA elasticity are key parameters: we
find that facilitated diffusion is most efficient inside an isotropic volume,
and on a flexible polymer. By considering the typical copy numbers of proteins
in vivo, we show that the speedup due to sliding becomes insensitive to fine
tuning of parameters, rendering facilitated diffusion a robust mechanism to
speed up intracellular diffusion-limited reactions. The parameter range we
focus on is relevant for in vitro systems and for facilitated diffusion on
yeast chromatin.


  Hubbell's neutral theory of biodiversity has successfully explained the
observed composition of many ecological communities but it relies on strict
demographic equivalence among species and provides no room for evolutionary
processes like selection, adaptation and speciation. Here we show how to embed
the neutral theory within the Darwinian framework. In a fitness landscape with
a quadratic maximum, typical of quantitative traits, selection restricts the
extant species to have traits close to optimal, so that the fitness differences
between surviving species are small. For sufficiently small mutation steps, the
community structure fits perfectly to the Fisher log-series species abundance
distribution. The theory is relatively insensitive to moderate amounts of
environmental noise, wherein the location of the fitness maximum changes by
amounts of order the width of the noise-free distribution. Adding very large
environmental noise to the model qualitatively changes the abundance
distributions, converting the exponential fall-off of large species to a
power-law decay, typical of a neutral model with environmental noise.


  The Influenza A virus belongs to the Orthomyxoviridae family. Influenza virus
infection occurs yearly in all countries of the world. It usually kills between
250,000 and 500,000 people and causes severe illness in millions more. Over the
last century alone we have seen 3 global influenza pandemics. The great human
and financial cost of this disease has made it the second most studied virus
today, behind HIV. Recently, several genome-wide RNA interference studies have
focused on identifying host molecules that participate in Influenza infection.
We used nine of these studies for this meta-analysis. Even though the overlap
among genes identified in multiple screens was small, network analysis
indicates that similar protein complexes and biological functions of the host
were present. As a result, several host gene complexes important for the
Influenza virus life cycle were identified. The biological function and the
relevance of each identified protein complex in the Influenza virus life cycle
is further detailed in this paper.


  In this paper, we study a simple one-dimensional model of reaction-diffusion
with bistable non-linearity and in heterogeneous media. The bistable term
accounts for the so-called Allee effect, and the heterogeneity in the media is
localized. We recall the definition of a transition wave used in similar
situations by well-known authors, and propose an alternative definition for
discussion. We call it generalized traveling wave. As a consequence, we give
new results of pinning in such media, due to Allee effect.


  "The investigation into the possible effects of cosmic rays on living
organisms will also offer great interest." - Victor F. Hess, Nobel Lecture,
December 12, 1936
  High-energy radiation bursts are commonplace in our Universe. From nearby
solar flares to distant gamma ray bursts, a variety of physical processes
accelerate charged particles to a wide range of energies, which subsequently
reach the Earth. Such particles contribute to a number of physical processes
occurring in the Earth system. A large fraction of the energy of charged
particles gets deposited in the atmosphere, ionizing the atmosphere, causing
changes in its chemistry and affecting the global electric circuit. Remaining
secondary particles contribute to the background dose of cosmic rays on the
surface and parts of the subsurface region. Life has evolved over the past ~ 3
billion years in presence of this background radiation, which itself has varied
considerably during the period. As demonstrated by the Miller-Urey experiment,
lightning plays a very important role in the formation of complex organic
molecules, which are the building blocks of more complex structures forming
life. There is growing evidence of increase in the lightning rate with
increasing flux of charged particles. Is there a connection between enhanced
rate of cosmic rays and the origin of life? Cosmic ray secondaries are also
known to damage DNA and cause mutations, leading to cancer and other diseases.
It is now possible to compute radiation doses from secondary particles, in
particular muons and neutrons. Have the variations in cosmic ray flux affected
the evolution of life on earth? We describe the mechanisms of cosmic rays
affecting terrestrial life and review the potential implications of the
variation of high-energy astrophysical radiation on the history of life on
earth.


  Sexual selection theory models evolution of sexual signals and preferences
using simple life histories. However, life-history models predict that males
benefit from increasing sexual investment approaching old age, producing
age-dependent sexual traits. Age-dependent traits require time and energy to
grow, and will not fully mature before individuals enter mating competition.
Early evolutionary stages pose several problems for these traits. Age-dependent
traits suffer from strong viability selection and gain little benefit from mate
choice when rare. Few males will grow large traits, and they will rarely
encounter choosy females. The evolutionary origins of age-dependent traits
therefore remain unclear. I used numerical simulations to analyze evolution of
preferences, condition (viability) and traits in an age-structured population.
Traits in the model depended on age and condition ("good genes") in a
population with no genetic drift. I asked (1) if age-dependent indicator traits
and their preferences can originate depending on the strength of selection and
the size of the trait; (2) which mode of development (age-dependent versus
age-independent) eventually predominates when both modes occur in the
population; and (3) if age-independent traits can invade a population with
age-dependent traits. Age-dependent traits evolve under weaker selection and at
smaller sizes than age-independent traits. Evolution of age-independent traits
depends only on trait size, whereas evolution of age-dependent traits depends
on both strength of selection and growth rate. Invasion of age-independence
into populations with established traits followed a similar pattern with
age-dependence predominating at small trait sizes. I suggest that reduced adult
mortality facilitates sexual selection by favoring the evolution of
age-dependent sexual signals under weak selection.


  The equations of evolutionary change by natural selection are commonly
expressed in statistical terms. Fisher's fundamental theorem emphasizes the
variance in fitness. Quantitative genetics expresses selection with covariances
and regressions. Population genetic equations depend on genetic variances. How
can we read those statistical expressions with respect to the meaning of
natural selection? One possibility is to relate the statistical expressions to
the amount of information that populations accumulate by selection. However,
the connection between selection and information theory has never been
compelling. Here, I show the correct relations between statistical expressions
for selection and information theory expressions for selection. Those relations
link selection to the fundamental concepts of entropy and information in the
theories of physics, statistics, and communication. We can now read the
equations of selection in terms of their natural meaning. Selection causes
populations to accumulate information about the environment.


  We study the internal resonance, energy transfer, activation mechanism, and
control of a model of DNA division via parametric resonance. While the system
is robust to noise, this study shows that it is sensitive to specific fine
scale modes and frequencies that could be targeted by low intensity
electro-magnetic fields for triggering and controlling the division. The DNA
model is a chain of pendula in a Morse potential. While the (possibly
parametrically excited) system has a large number of degrees of freedom and a
large number of intrinsic time scales, global and slow variables can be
identified by (i) first reducing its dynamic to two modes exchanging energy
between each other and (ii) averaging the dynamic of the reduced system with
respect to the phase of the fastest mode. Surprisingly the global and slow
dynamic of the system remains Hamiltonian (despite the parametric excitation)
and the study of its associated effective potential shows how parametric
excitation can turn the unstable open state into a stable one. Numerical
experiments support the accuracy of the time-averaged reduced Hamiltonian in
capturing the global and slow dynamic of the full system.


  We propose PALPS, a Process Algebra with Locations for Population Systems.
PALPS allows us to produce spatially-explicit, individual-based models and to
reason about their behavior. Our calculus has two levels: at the first level we
may define the behavior of an individual of a population while, at the second
level, we may specify a system as the collection of individuals of various
species located in space, moving through their life cycle while changing their
location, if they so wish, and interacting with each other in various ways such
as preying on each other. Furthermore, we propose a probabilistic temporal
logic for reasoning about the behavior of PALPS processes. We illustrate our
framework via models of dispersal in metapopulations.


  RNA interference (RNAi) is a mechanism whereby small RNAs (siRNAs) directly
control gene expression without assistance from proteins. This mechanism
consists of interactions between RNAs and small RNAs both of which may be
single or double stranded. The target of the mechanism is mRNA to be degraded
or aberrated, while the initiator is double stranded RNA (dsRNA) to be cleaved
into siRNAs. Observing the digital nature of RNAi, we represent RNAi as a
Minsky register machine such that (i) The two registers hold single and double
stranded RNAs respectively, and (ii) Machine's instructions are interpreted by
interactions of enzyme (Dicer), siRNA (with RISC com- plex) and polymerization
(RdRp) to the appropriate registers. Interpreting RNAi as a computational
structure, we can investigate the computational meaning of RNAi, especially its
complexity. Initially, the machine is configured as a Chemical Ground Form
(CGF), which generates incorrect jumps. To remedy this problem, the system is
remodeled as recursive RNAi, in which siRNA targets not only mRNA but also the
machine instructional analogues of Dicer and RISC. Finally, probabilistic
termination is investigated in the recursive RNAi system.


  In this paper we study the effect of rare mutations, driven by a marked point
process, on the evolutionary behavior of a population. We derive a Kolmogorov
equation describing the expected values of the different frequencies and prove
some rigorous analytical results about their behavior. Finally, in a simple
case of two different quasispecies, we are able to prove that the rarity of
mutations increases the survival opportunity of the low fitness species.


  Trypanosoma cruzi and Giardia intestinalis are two human pathogens and
protozoan parasites responsible for the diseases Chagas disease and giardiasis,
respectively. Both diseases cause suffering and illness in several million
individuals. The former disease occurs primarily in South America and Central
America, and the latter disease occurs worldwide. Current therapeutics are
toxic and lack efficacy, and potential vaccines are far from the market.
Increased knowledge about the biology of these parasites is essential for drug
and vaccine development, and new diagnostic tests. In this thesis,
high-throughput sequencing was applied together with extensive bioinformatic
analyses to yield insights into the biology and evolution of Trypanosoma cruzi
and Giardia intestinalis. Bioinformatics analysis of DNA and RNA sequences was
performed to identify features that may be of importance for parasite biology
and functional characterization. This thesis is based on five papers (i-v).
Paper i and ii describe comparative genome studies of three distinct genotypes
of Giardia intestinalis (A, B and E). Paper iii describes a genome comparison
of the human infecting Trypanosoma cruzi with the bat-restricted subspecies
Trypanosoma cruzi marinkellei. Paper iv describes the repertoire of small
non-coding RNAs in Trypanosoma cruzi epimastigotes. Paper v describes
transcriptome analysis using paired-end RNA-Seq of three distinct genotypes of
Giardia intestinalis (A, B and E).


  We introduce a system for experimental evolution consisting of populations of
short oligonucleotides (Oli populations) evolving in a modified quantitative
Polymerase Chain Reaction (qPCR). It is tractable at the genetic, genomic,
phenotypic and fitness levels. The Oli system uses DNA hairpins designed to
form structures that self-prime under defined conditions. Selection acts on the
phenotype of self-priming, after which differences in fitness are amplified and
quantified using qPCR. We outline the methodological and bioinformatics tools
for the Oli system here, and demonstrate that it can be used as a conventional
experimental evolution model system by test-driving it in an experiment
investigating adaptive evolution under different rates of environmental change.


  The relative solvent accessibility (RSA) of a residue in a protein measures
the extent of burial or exposure of that residue in the 3D structure. RSA is
frequently used to describe a protein's biophysical or evolutionary properties.
To calculate RSA, a residue's solvent accessibility (ASA) needs to be
normalized by a suitable reference value for the given amino acid; several
normalization scales have previously been proposed. However, these scales do
not provide tight upper bounds on ASA values frequently observed in empirical
crystal structures. Instead, they underestimate the largest allowed ASA values,
by up to 20%. As a result, many empirical crystal structures contain residues
that seem to have RSA values in excess of one. Here, we derive a new
normalization scale that does provide a tight upper bound on observed ASA
values. We pursue two complementary strategies, one based on extensive analysis
of empirical structures and one based on systematic enumeration of
biophysically allowed tripeptides. Both approaches yield congruent results that
consistently exceed published values. We conclude that previously published ASA
normalization values were too small, primarily because the conformations that
maximize ASA had not been correctly identified. As an application of our
results, we show that empirically derived hydrophobicity scales are sensitive
to accurate RSA calculation, and we derive new hydrophobicity scales that show
increased correlation with experimentally measured scales.


  When dealing with large scale gene expression studies, observations are
commonly contaminated by unwanted variation factors such as platforms or
batches. Not taking this unwanted variation into account when analyzing the
data can lead to spurious associations and to missing important signals. When
the analysis is unsupervised, e.g., when the goal is to cluster the samples or
to build a corrected version of the dataset - as opposed to the study of an
observed factor of interest - taking unwanted variation into account can become
a difficult task. The unwanted variation factors may be correlated with the
unobserved factor of interest, so that correcting for the former can remove the
latter if not done carefully. We show how negative control genes and replicate
samples can be used to estimate unwanted variation in gene expression, and
discuss how this information can be used to correct the expression data or
build estimators for unsupervised problems. The proposed methods are then
evaluated on three gene expression datasets. They generally manage to remove
unwanted variation without losing the signal of interest and compare favorably
to state of the art corrections.


  In aiming to explain the establishment, maintenance and stability of
methylation pattern in gene body of Arabidopsis we propose here a theoretical
framework for understanding how the methylated and unmethylated states of
cytosine residues are maintained and transmitted during DNA replication. Routed
in statistical mechanics, the framework built herein is used to explore minimal
models of epigenetic inheritance and identify the necessary conditions for
stability of methylated/unmethylated states of cytosine over rounds of DNA
replication. The models are flexible enough to allow adding new biological
concepts and information.


  In the beeswax combs of honey bees, the cells of brood, pollen, and honey
have a consistent spatial pattern that is sustained throughout the life of a
colony. This spatial pattern is believed to emerge from simple behavioral rules
that specify how the queen moves, where foragers deposit honey/pollen and how
honey/pollen is consumed from cells. Prior work has shown that a set of such
rules can explain the formation of the allocation pattern starting from an
empty comb. We show that these rules cannot maintain the pattern once the brood
start to vacate their cells, and we propose new, biologically realistic rules
that better sustain the observed allocation pattern. We analyze the three
resulting models by performing hundreds of simulation runs over many
gestational periods and a wide range of parameter values. We develop new
metrics for pattern assessment and employ them in analyzing pattern retention
over each simulation run. Applied to our simulation results, these metrics show
alteration of an accepted model for honey/pollen consumption based on local
information can stabilize the cell allocation pattern over time. We also show
that adding global information, by biasing the queen's movements towards the
center of the comb, expands the parameter regime over which pattern retention
occurs.


  The esteemed physicist Erwin Schroedinger, whose name is associated with the
most notorious equation of quantum mechanics, also wrote a brief essay entitled
'What is Life?', asking: 'How can the events in space and time which take place
within the spatial boundary of a living organism be accounted for by physics
and chemistry?' The 60+ years following this seminal work have seen enormous
developments in our understanding of biology on the molecular scale, physics
playing a key role in solving many central problems through the development and
application of new physical science techniques, biophysical analysis and
rigorous intellectual insight. The early days of single molecule biophysics
research was centred around molecular motors and biopolymers, largely divorced
from a real physiological context. The new generation of single molecule
bioscience investigations has much greater scope, involving robust methods for
understanding molecular level details of the most fundamental biological
processes in far more realistic, and technically challenging, physiological
contexts, emerging into a new field of 'single molecule cellular biophysics'.
Here, I outline how this new field has evolved, discuss the key active areas of
current research, and speculate on where this may all lead in the near future.


  An important challenge in cancer systems biology is to uncover the complex
network of interactions between genes (tumor suppressor genes and oncogenes)
implicated in cancer. Next generation sequencing provides unparalleled ability
to probe the expression levels of the entire set of cancer genes and their
transcript isoforms. However, there are onerous statistical and computational
issues in interpreting high-dimensional sequencing data and inferring the
underlying genetic network. In this study, we analyzed RNA-Seq data from
lymphoblastoid cell lines derived from a population of 69 human individuals and
implemented a probabilistic framework to construct biologically-relevant
genetic networks. In particular, we employed a graphical lasso analysis,
motivated by considerations of the maximum entropy formalism, to estimate the
sparse inverse covariance matrix of RNA-Seq data. Gene ontology, pathway
enrichment and protein-protein path length analysis were all carried out to
validate the biological context of the predicted network of interacting cancer
gene isoforms.


  In phylogenetic studies, the evolution of molecular sequences is assumed to
have taken place along the phylogeny traced by the ancestors of extant species.
In the presence of lateral gene transfer (LGT), however, this may not be the
case, because the species lineage from which a gene was transferred may have
gone extinct or not have been sampled. Because it is not feasible to specify or
reconstruct the complete phylogeny of all species, we must describe the
evolution of genes outside the represented phylogeny by modelling the
speciation dynamics that gave rise to the complete phylogeny. We demonstrate
that if the number of sampled species is small compared to the total number of
existing species, the overwhelming majority of gene transfers involve
speciation to, and evolution along extinct or unsampled lineages. We show that
the evolution of genes along extinct or unsampled lineages can to good
approximation be treated as those of independently evolving lineages described
by a few global parameters. Using this result, we derive an algorithm to
calculate the probability of a gene tree and recover the maximum likelihood
reconciliation given the phylogeny of the sampled species. Examining 473 near
universal gene families from 36 cyanobacteria, we find that nearly a third of
transfer events -- 28% -- appear to have topological signatures of evolution
along extinct species, but only approximately 6% of transfers trace their
ancestry to before the common ancestor of the sampled cyanobacteria.


  We explore a generic mechanism whereby a droplet of active matter acquires
motility by the spontaneous breakdown of a discrete symmetry. The model we
study offers a simple representation of a "cell extract" comprising, e.g., a
droplet of actomyosin solution. (Such extracts are used experimentally to model
the cytoskeleton.) Actomyosin is an active gel whose polarity describes the
mean sense of alignment of actin fibres. In the absence of polymerization and
depolymerization processes ('treadmilling'), the gel's dynamics arises solely
from the contractile motion of myosin motors; this should be unchanged when
polarity is inverted. Our results suggest that motility can arise in the
absence of treadmilling, by spontaneous symmetry breaking (SSB) of polarity
inversion symmetry. Adapting our model to wall-bound cells in two dimensions,
we find that as wall friction is reduced, treadmilling-induced motility falls
but SSB-mediated motility rises. The latter might therefore be crucial in three
dimensions where frictional forces are likely to be modest. At a supra-cellular
level, the same generic mechanism can impart motility to aggregates of
non-motile but active bacteria; we show that SSB in this (extensile) case leads
generically to rotational as well as translational motion.


  Species interactions shape predator-prey networks, impacting community
structure and, potentially, ecological dynamics. It is likely that global
climatic perturbations that occur over long periods of time have a significant
impact on species interactions patterns. However, observations of how these
patterns change over time are typically limited to extant communities, which is
particularly problematic for communities with long-lived species. Here we
integrate stable isotope analysis and network theory to reconstruct patterns of
trophic interactions for six independent mammalian communities that inhabited
mammoth steppe environments spanning western Europe to eastern Alaska during
the Pleistocene. We use a Bayesian mixing model to quantify the proportional
contribution of prey to the diets of local predators, and assess how the
structure of trophic interactions changed across space and the Last Glacial
Maximum (LGM), a global climatic event that severely impacted mammoth steppe
communities. We find that large felids had diets that were more constrained
than other co-occurring predators, and largely influenced by an increase in
{\it Rangifer} abundance after the LGM. Moreover, the structural organization
of Beringian and European communities strongly differed: compared to Europe,
species interactions in Beringian communities before the LGM were highly
compartmentalized, or modular. This modularity was lost during the LGM, and
partially recovered after the glacial retreat, and we suggest that changes in
modularity among predators and prey may have been driven by geographic
insularity.


  Foraging is constrained by the energy within resources and the mechanics of
acquisition and assimilation. Thick molar enamel, a character trait
differentiating hominins from African apes, is predicted to mitigate the
mechanical costs of chewing obdurate foods. The classic expression of
hyperthick enamel together with relatively massive molars, termed megadontia,
is most evident in {\it Paranthropus}, a lineage of hominins that lived ca. 2.7
to 1.2 million years ago. Among contemporary primates, thicker molar enamel
corresponds with the consumption of stiffer, deformation-resistant foods,
possibly because thicker enamel can better resist cracking under high
compressive loads. Accordingly, plant underground storage organs (USOs) are
thought to be a central food resource for hominins such as {\it Paranthropus}
due to their abundance, isotopic composition, and mechanical properties. Here,
we present a process-based model to investigate foraging constraints as a
function of energetic demands and enamel wear among human ancestors. Our
framework allows us to determine the fitness benefits of megadontia, and to
explore under what conditions stiff foods such as USOs are predicted to be
chosen as fallback, rather than preferred, resources. Our model predictions
bring consilience to the noted disparity between functional interpretations of
megadontia and microwear evidence, particularly with respect to {\it
Paranthropus boisei}.


  Long-distance bidirectional transport of organelles depends on the motor
proteins kinesin and dynein. Using quantitative data obtained from a fungal
model system, we previously developed ASEP-models of bidirectional motion of
motors along unipolar microtubules (MTs) near the cell ends of the elongated
hyphal cells (herein referred as "unipolar section"). However, recent
quantitative live cell imaging in this system has demonstrated that long-range
motility of motors and their endosomal cargo mainly occurs along extended
antipolar microtubule bundles within the central part of the cell (herein
referred to as "bipolar section"). Dynein and kinesin-3 motors coordinate their
activity to move early endosomes (EEs) in a bidirectional fashion, with dynein
mediating retrograde motility along the unipolar section near the cell poles,
whereas kinesin-3 is responsible for bidirectional motions along the antipolar
section. Here we extend our modelling approach to simulate bidirectional
motility along an antipolar microtubule bundle. In our model, cargos
(particles) change direction on each MT with a turning rate $\Omega$ and the
MTs are linked to each other at the minus ends where particles can hop between
MTs with a rate $q_1$ (obstacle-induced switching rate) or $q_2$ (end-induced
switching rate). By numerical simulations and mean-field approximations, we
investigate the distribution of particles along the MTs for different overall
densities $\Theta$. We find that even if $\Theta$ is low, the system can
exhibit shocks in the density profiles near plus and minus ends caused by
queueing of particles. We also discuss how the switching rates $q_{1,2}$
influence the type of motor that dominates the active transport in the bundle.


  The protein p53 has a well established role in protecting genomic integrity
in human cells. When DNA is damaged p53 induces the cell cycle arrest to
prevent the transmission of the damage to cell progeny, triggers the production
of proteins for DNA repair and ultimately calls for apoptosis. In particular,
the p53-Mdm2 feedback loop seems to be the key circuit in this response of
cells to damage. For many years, based on measurements over populations of
cells it was believed that the p53-Mdm2 feedback loop was the responsible for
the existence of damped oscillations in the levels of p53 and Mdm2 after DNA
damage. However, recent measurements in individual human cells have shown that
p53 and its regulator Mdm2 develop sustained oscillations over long periods of
time even in the absence of stress. These results have attracted a lot of
interest, first because they open a new experimental framework to study the p53
and its interactions and second because they challenge years of mathematical
models with new and accurate data on single cells. Inspired by these
experiments standard models of the p53-Mdm2 circuit were modified introducing
ad-hoc some biologically motivated noise that becomes responsible for the
stability of the oscillations. Here, we follow an alternative approach
proposing that the noise that stabilizes the fluctuations is the intrinsic
noise due to the finite nature of the populations of p53 and Mdm2 in a single
cell.


  Non-equilibrium processes which convert chemical energy into mechanical
motion enable the motility of organisms. Bundles of inextensible filaments
driven by energy transduction of molecular motors form essential components of
micron-scale motility engines like cilia and flagella. The mimicry of
cilia-like motion in recent experiments on synthetic active filaments supports
the idea that generic physical mechanisms may be sufficient to generate such
motion. Here we show, theoretically, that the competition between the
destabilising effect of hydrodynamic interactions induced by force-free and
torque-free chemomechanically active flows, and the stabilising effect of
nonlinear elasticity, provides a generic route to spontaneous oscillations in
active filaments. These oscillations, reminiscent of prokaryotic and eukaryotic
flagellar motion, are obtained without having to invoke structural complexity
or biochemical regulation. This minimality implies that biomimetic
oscillations, previously observed only in complex bundles of active filaments,
can be replicated in simple chains of generic chemomechanically active beads.


  Nucleic acids are negatively charged macromolecules and their structure
properties are strongly coupled to metal ions in solutions. In this paper, the
salt effects on the flexibility of single stranded (ss) nucleic acid chain
ranging from 12 to 120 nucleotides are investigated systematically by the
coarse grained Monte Carlo simulations where the salt ions are considered
explicitly and the ss chain is modeled with the virtual bond structural model.
Our calculations show that, the increase of ion concentration causes the
structural collapse of ss chain and multivalent ions are much more efficient in
causing such collapse, and trivalent and small divalent ions can both induce
more compact state than a random relaxation state. We found that monovalent,
divalent and trivalent ions can all overcharge ss chain, and the dominating
source for such overcharging changes from ion exclusion volume effect to ion
Coulomb correlations. In addition, the predicted Na and Mg dependent
persistence length lp of ss nucleic acid are in accordance with the available
experimental data, and through systematic calculations, we obtained the
empirical formulas for lp as a function of Na, Mg and chain length.


  Here we have systematically studied the cooperative binding of substrate
molecules on the active sites of a single oligomeric enzyme in a chemiostatic
condition. The average number of bound substrate and the net velocity of the
enzyme catalyzed reaction are studied by the formulation of stochastic master
equation for the cooperative binding classified here as spatial and temporal.
We have estimated the entropy production for the cooperative binding schemes
based on single trajectory analysis using a kinetic Monte Carlo technique. It
is found that the total as well as the medium entropy production show the same
generic diagnostic signature for detecting the cooperativity, usually
characterized in terms of the net velocity of the reaction. This feature is
also found to be valid for the total entropy production rate at the
nonequilibrium steady state. We have introduced an index of cooperativity, C,
defined in terms of the ratio of the surprisals or equivalently, the stochastic
system entropy associated with the fully bound state of the cooperative and
non-cooperative cases. The criteria of cooperativity in terms of C is compared
with that of the Hill coefficient and gives a microscopic insight on the
cooperative binding of substrate on a single oligomeric enzyme which is usually
characterized by macroscopic reaction rate.


  In multi-domain proteins, the domains are connected by a flexible
unstructured region called as protein domain linker. The accurate demarcation
of these linkers holds a key to understanding of their biochemical and
evolutionary attributes. This knowledge helps in designing a suitable linker
for engineering stable multi-domain chimeric proteins. Here we propose a novel
method for the demarcation of the linker based on a three-dimensional protein
structure and a domain definition. The proposed method is based on biological
knowledge about structural flexibility of the linkers. We performed structural
analysis on a linker probable region (LPR) around domain boundary points of
known SCOP domains. The LPR was described using a set of overlapping peptide
fragments of fixed size. Each peptide fragment was then described by geometric
invariants (GIs) and subjected to clustering process where the fragments
corresponding to actual linker come up as outliers. We then discover the actual
linkers by finding the longest continuous stretch of outlier fragments from
LPRs. This method was evaluated on a benchmark dataset of 51 continuous
multi-domain proteins, where it achieves F1 score of 0.745 (0.83 precision and
0.66 recall). When the method was applied on 725 continuous multi-domain
proteins, it was able to identify novel linkers that were not reported
previously. This method can be used in combination with supervised / sequence
based linker prediction methods for accurate linker demarcation.


  Although the polymorphism of calcium carbonate is well known, and its
polymorphs-calcite, aragonite, and vaterite-have been highly studied in the
context of biomineralization, polyamorphism is a much more recently discovered
phenomenon, and the existence of more than one amorphous phase of calcium
carbonate in biominerals has only very recently been understood. Here we
summarize what is known about polyamorphism in calcium carbonate as well as
what is under- stood about the role of amorphous calcium carbonate in
biominerals. We show that consideration of the amorphous forms of calcium
carbonate within the physical notion of polyamorphism leads to new insights
when it comes to the mechanisms by which polymorphic structures can evolve in
the first place. This not only has implications for our understanding of
biomineralization, but also of the means by which crystallization may be
controlled in medical, pharmaceutical, and industrial contexts.


  Complexes of physically interacting proteins are one of the fundamental
functional units responsible for driving key biological mechanisms within the
cell. Their identification is therefore necessary not only to understand
complex formation but also the higher level organization of the cell. With the
advent of high-throughput techniques in molecular biology, significant amount
of physical interaction data has been cataloged from organisms such as yeast,
which has in turn fueled computational approaches to systematically mine
complexes from the network of physical interactions among proteins (PPI
network). In this survey, we review, classify and evaluate some of the key
computational methods developed till date for the identification of protein
complexes from PPI networks. We present two insightful taxonomies that reflect
how these methods have evolved over the years towards improving automated
complex prediction. We also discuss some open challenges facing accurate
reconstruction of complexes, the crucial ones being presence of high proportion
of errors and noise in current high-throughput datasets and some key aspects
overlooked by current complex detection methods. We hope this review will not
only help to condense the history of computational complex detection for easy
reference, but also provide valuable insights to drive further research in this
area.


  Finite dynamical systems (e.g. Boolean networks and logical models) have been
used in modeling biological systems to focus attention on the qualitative
features of the system, such as the wiring diagram. Since the analysis of such
systems is hard, it is necessary to focus on subclasses that have the
properties of being general enough for modeling and simple enough for
theoretical analysis. In this paper we propose the class of AND-NOT networks
for modeling biological systems and show that it provides several advantages.
Some of the advantages include: Any finite dynamical system can be written as
an AND-NOT network with similar dynamical properties. There is a one-to-one
correspondence between AND-NOT networks, their wiring diagrams, and their
dynamics. Results about AND-NOT networks can be stated at the wiring diagram
level without losing any information. Results about AND-NOT networks are
applicable to any Boolean network. We apply our results to a Boolean model of
Th-cell differentiation.


  Biological systems are modular, and this modularity evolves over time and in
different environments. A number of observations have been made of increased
modularity in biological systems under increased environmental pressure. We
here develop a quasispecies theory for the dynamics of modularity in
populations of these systems. We show how the steady-state fitness in a
randomly changing environment can be computed. We derive a fluctuation
dissipation relation for the rate of change of modularity and use it to derive
a relationship between rate of environmental changes and rate of growth of
modularity. We also find a principle of least action for the evolved modularity
at steady state. Finally, we compare our predictions to simulations of protein
evolution and find them to be consistent.


  We have investigated the influence of the plasma membrane environment on the
molecular evolution of G protein-coupled receptors (GPCRs), the largest
receptor family in Metazoa. In particular, we have analyzed the site-specific
rate variation across the two primary structural partitions, transmembrane (TM)
and extramembrane (EM), of these membrane proteins. We find that transmembrane
domains evolve more slowly than do extramembrane domains, though TM domains
display increased rate heterogeneity relative to their EM counterparts.
Although the majority of residues across GPCRs experience strong to weak
purifying selection, many GPCRs experience positive selection at both TM and EM
residues, albeit with a slight bias towards the EM. Further, a subset of GPCRs,
chemosensory receptors (including olfactory and taste receptors), exhibit
increased rates of evolution relative to other GPCRs, an effect which is more
pronounced in their TM spans. Although it has been previously suggested that
the TM's low evolutionary rate is caused by their high percentage of buried
residues, we show that their attenuated rate seems to stem from the strong
biophysical constraints of the membrane itself, or by functional requirements.
In spite of the strong evolutionary constraints acting on the transmembrane
spans of GPCRs, positive selection and high levels of evolutionary rate
variability are common. Thus, biophysical constraints should not be presumed to
preclude a protein's ability to evolve.


  Access to combination antiretroviral treatment (ART) has improved greatly
over recent years. At the end of 2011, more than eight million HIV infected
people were receiving antiretroviral therapy in low-income and middle-income
countries. ART generally works well in keeping the virus suppressed and the
patient healthy. However, treatment only works as long as the virus is not
resistant against the drugs used. In the last decades, HIV treatments have
become better and better at slowing down the evolution of drug resistance, so
that some patients are treated for many years without having any resistance
problems. However, for some patients, especially in low-income countries, drug
resistance is still a serious threat to their health. This essay will review
what is known about transmitted and acquired drug resistance, multi-class drug
resistance, resistance to newer drugs, resistance due to treatment for the
prevention of mother-to-child transmission, the role of minority variants
(low-frequency drug-resistance mutations), and resistance due to pre-exposure
prophylaxis.


  This paper presents a study of the two-predators-two-preys discrete-time
Lotka-Volterra model with self- inhibition terms for preys with direct
applications to ecological problems. Parameters in the model are modified so
that each of them has its own biological meaning, enabling more intuitive
interpretation of biological conditions to the model . Moreover, the modified
version is applicable to simulate a part of a large ecosystem, not only a
closed predator-prey system with four species. An easy graphical method of
analysis of the conditions on parameters ensuring long persistence under
coexistence is presented. As an application, it is explained that a predator
specie who feed on relatively small number of preys compared to the other
predator species must be selective on the available preys in order for long
persistence of the ecosystem. This may be regarded as a theoretical explanation
of the existence of flush-pursuer birds, those who uses highly specialized
hunting strategy and cross-adapts to the ecosystem relative to the ordinary
bird species.


  The opportunistic fungal pathogen Candida albicans causes serious hematogenic
hospital acquired candidiasis with worldwide impact on public health. Because
of its importance as a nosocomial etiologic agent, C. albicans genome has been
largely studied to identify intraspecific variation and several typing methods
have been developed to distinguish closely related strains. Mitochondrial DNA
can be useful for this purpose because, as compared to nuclear DNA, its higher
mutational load and evolutionary rate readily reveals microvariants.
Accordingly, we sequenced and assembled, with 8 fold coverage, the
mitochondrial genomes of two C. albicans clinical isolates (L296 and L757) and
compared these sequences with the genome sequence of reference strain SC5314.
The genome alignment of 33,928 positions revealed 372 polymorphic sites being
230 in coding and 142 in non-coding regions. Three intergenic regions located
between genes tRNAGly/COX1, NAD3/COB and ssurRNA/NAD4L, named IG1, IG2 and IG3
respectively, which showed high number of neutral substitutions, were amplified
and sequenced from 18 clinical isolates from different locations in Latin
America and 2 ATCC standard C. albicans strains. High variability of sequence
and size were observed, ranging up to 56bp size difference and phylogenies
based on IG1, IG2 and IG3 revealed three groups. Insertions of up to 49bp were
observed exclusively in Argentinean strains relative to the other sequences
which could suggest clustering by geographical polymorphism. Because of neutral
evolution, high variability, easy isolation by PCR and full length sequencing
these mitochondrial intergenic regions can contribute with a novel perspective
in molecular studies of C. albicans isolates, complementing well established
multilocus sequence typing methods.


  Mutations and oxidative modification in the protein Cu,Zn superoxide
dismutase (SOD1) have been implicated in the death of motor neurons in
amyotrophic lateral sclerosis (ALS), a presently incurable, invariably fatal
neurodegenerative disease. Here we employ steered, all-atom molecular dynamics
simulations in implicit solvent to investigate the significance of either
mutations or post-translational modifications (PTMs) to SOD1 on metal affinity,
dimer stability, and mechanical malleability. The work required to induce
moderate structural deformations as a function of sequence index constitutes a
"mechanical fingerprint" measuring structural rigidity in the native basin,
from which we are able to unambiguously distinguish wild-type (WT) SOD1 from
PTM variants, and measure the severity of a given PTM on structural integrity.
The cumulative distribution of work values provided a way to cleanly
discriminate between SOD1 variants. Disulfide reduction destabilizes dimer
stability more than the removal of either metal, but not moreso than the
removal of both metals. Intriguingly, we found that disulfide reduction
mechanically stabilizes apo SOD1 monomer, underscoring the differences between
native basin mechanical properties and equilibrium thermodynamic stabilities,
and elucidating the presence of internal stress in the apo state. All PTMs and
ALS-associated mutants studied showed an increased tendency to lose either Cu
or Zn, and to monomerize- processes known to be critical in the progression of
ALS. The valence of Cu strongly modulates its binding free energy. As well,
several mutants were more susceptible to loss of metals and monomerization than
the disulfide-reduced or apo forms of SOD1. Distance constraints are required
to calculate free energies for metal binding and dimer separation, which are
validated using thermodynamic cycles.


  The Allen Brain Atlas project (ABA) generated a genome-scale collection of
gene-expression profiles using in-situ hybridization. These profiles were
co-registered to the three-dimensional Allen Reference Atlas (ARA) of the adult
mouse brain. A set of more than 4,000 such volumetric data are available for
the full brain, at a resolution of 200 microns. These data are presented in a
voxel-by-gene matrix. The ARA comes with several systems of annotation,
hierarchical (40 cortical regions, 209 sub-cortical regions in the whole
brain), or non-hierarchical (12 regions in the left hemisphere, with refinement
into 94 regions, and cortical layers). The high-dimensional nature of this
dataset and the possible connection between anatomy and gene expression pose
challenges to data analysis. We developed the Brain Gene Expression Analysis
Toolbox, whose functionalities include: determination of marker genes for brain
regions, statistical analysis of brain-wide co-expression patterns, and the
computation of brain-wide correlation maps with cell-type specific microarray
data.


  In a shotgun proteomics experiment, proteins are the most biologically
meaningful output. The success of proteomics studies depends on the ability to
accurately and efficiently identify proteins. Many methods have been proposed
to facilitate the identification of proteins from the results of peptide
identification. However, the relationship between protein identification and
peptide identification has not been thoroughly explained before.
  In this paper, we are devoted to a combinatorial perspective of the protein
inference problem. We employ combinatorial mathematics to calculate the
conditional protein probabilities (Protein probability means the probability
that a protein is correctly identified) under three assumptions, which lead to
a lower bound, an upper bound and an empirical estimation of protein
probabilities, respectively. The combinatorial perspective enables us to obtain
a closed-form formulation for protein inference.
  Based on our model, we study the impact of unique peptides and degenerate
peptides on protein probabilities. Here, degenerate peptides are peptides
shared by at least two proteins. Meanwhile, we also study the relationship of
our model with other methods such as ProteinProphet. A probability confidence
interval can be calculated and used together with probability to filter the
protein identification result. Our method achieves competitive results with
ProteinProphet in a more efficient manner in the experiment based on two
datasets of standard protein mixtures and two datasets of real samples.
  We name our program ProteinInfer. Its Java source code is available at
http://bioinformatics.ust.hk/proteininfer


  Limited spectrum coverage is a problem in shotgun proteomics. Replicates are
generated to improve the spectrum coverage. When integrating peptide
identification results obtained from replicates, the state-of-the-art algorithm
PeptideProphet combines Peptide-Spectrum Matches (PSMs) before building the
statistical model to calculate peptide probabilities.
  In this paper, we find the connection between merging results of replicates
and Bagging, which is a standard routine to improve the power of statistical
methods. Following Bagging's philosophy, we propose to run PeptideProphet
separately on each replicate and combine the outputs to obtain the final
peptide probabilities. In our experiments, we show that the proposed routine
can improve PeptideProphet consistently on a standard protein dataset, a Human
dataset and a Yeast dataset.


  The position of an invasion front, propagating into an unstable state,
fluctuates because of the shot noise coming from the discreteness of reacting
particles and stochastic character of the reactions and diffusion. A recent
macroscopic theory [Meerson and Sasorov, Phys. Rev. E 84, 030101(R) (2011)]
yields the probability of observing, during a long time, an unusually slow
front. The theory is formulated as an effective classical Hamiltonian field
theory which operates with the density field and the conjugate "momentum"
field. Further, the theory assumes that the most probable density field history
of an unusually slow front represents, up to small corrections, a traveling
front solution of the Hamilton equations. Here we verify this assumption by
solving the Hamilton equations numerically for models belonging to the directed
percolation universality class.


  Non-coding RNA sequences play a great role in controlling a number of
cellular functions, thus raising the need to understand their complex
conformational dynamics in quantitative detail. In this perspective, we first
show that single molecule pulling experiments when combined with with theory
and simulations can be used to quantitatively explore the folding landscape of
nucleic acid hairpins, and riboswitches with tertiary interactions.
Applications to riboswitches, which are non-coding RNA elements that control
gene expression by undergoing dynamical conformational changes in response to
binding of metabolites, lead to an organization principle that assembly of RNA
is determined by the stability of isolated helices. We also point out the
limitations of single molecule pulling experiments, with molecular extension as
the only accessible parameter, in extracting key parameters of the folding
landscapes of RNA molecules.


  Replicator equations (RE) are among the basic tools in mathematical theory of
selection and evolution. The method of Hidden keystone variables (HKV) for
reducing a wide class of RE of high or even infinite dimensionality to
low-dimensional escort systems of ODEs, which can be explored analytically in
many cases was developed. The method has potential for different applications
to models of biological populations and communities