  A consensus is emerging that continuous (or metric) measures can be useful in
phylogenetic systematics. Many of the methods for coding such characters, how-
ever, employ elements that are arbitrary and therefore should be excluded from
use in cladistic analysis. The continued use of such potentially inappropriate
methods can be attributed to either their simplicity or the availability of
computer programs specifically designed to produce data matrices using these
methods. Conversely, one of the most suitable methods, homogeneous subset
coding (HSC), is underused, probably due to the lack of a suitable software
implementation for this somewhat complex procedure.
  This paper describes TaxMan, a Web-based facility for the coding of
continuous data using HSC. Data are entered using a form accessible via any
internet browser and are automatically converted to a matrix suitable for input
into tree-searching software. This implementation of the HSC technique provides
an uncomplicated procedure for the incorporation of metric data in phylogenetic
systematics. The algorithmic implementation of the HSC procedure, and
interpolation of the Studentised range and maximum modulus distributions
required by it, is described in detail in appendices.


  Evolutionary models on graphs, as an extension of the Moran process, have two
major implementations: birth-death (BD) models (or the invasion process) and
death- birth (DB) models (or voter models). The isothermal theorem states that
the fixation probability of mutants in a large group of graph structures (known
as isothermal graphs, which include regular graphs) coincides with that for the
mixed population. This result has been proven by Lieberman et al (Nature 433:
312-316, 2005) in the case of BD processes, where mutants differ from the wild
types by their birth rate (and not by their death rate). In this paper we
discuss to what extent the isothermal theorem can be formulated for DB
processes, proving that it only holds for mutants that differ from the wild
type by their death rate (and not by their birth rate). For more general BD and
DB processes with arbitrary birth and death rates of mutants, we show that the
fixation probabilities of mutants are different from those obtained in the
mass-action populations. We focus on spatial lattices and show that the
difference between BD and DB processes on 1D and 2D lattices are non-small even
for large population sizes. We support these results with a generating function
approach that can be generalized to arbitrary graph structures. Finally, we
discuss several biological applications of the results.


  The molecular motor protein kinesin plays a key role in fundamental cellular
processes such as intracellular transport, mitotic spindle formation, and
cytokinesis, with important implications for neurodegenerative and cancer
disease pathways. Recently, kinesin has been studied as a paradigm for the
tailored design of nano-bio sensor and other nanoscale systems. As it processes
along a microtubule within the cell, kinesin undergoes a cycle of chemical
state and physical conformation transitions that enable it to take ~100 regular
8.2-nm steps before ending its processive walk. Despite an extensive body of
experimental and theoretical work, a unified microscopic model of kinesin
mechanochemistry does not yet exist. Here we present a methodology that
optimizes a kinetic model for kinesin constructed with a minimum of a priori
assumptions about the underlying processive mechanism. Kinetic models are
preferred for numerical calculations since information about the kinesin
stepping mechanism at all levels, from the atomic to the microscopic scale, is
fully contained within the particular states of the cycle: how states
transition, and the rate constants associated with each transition. We combine
Markov chain calculations and simulated annealing optimization to determine the
rate constants that best fit experimental data on kinesin speed and
processivity.


  In the mammalian suprachiasmatic nucleus (SCN), a population of noisy
cell-autonomous oscillators synchronizes to generate robust circadian rhythms
at the organism-level. Within these cells two isoforms of Cryptochrome, Cry1
and Cry2, participate in a negative feedback loop driving circadian
rhythmicity. Previous work has shown that single, dissociated SCN neurons
respond differently to Cry1 and Cry2 knockouts: Cry1 knockouts are arrhythmic
while Cry2 knockouts display more regular rhythms. These differences have led
to speculation that CRY1 and CRY2 may play different functional roles in the
oscillator. To address this proposition, we have developed a new coupled,
stochastic model focused on the Period (Per) and Cry feedback loop, and
incorporating intercellular coupling via vasoactive intestinal peptide (VIP).
Due to the stochastic nature of molecular oscillations, we demonstrate that
single-cell Cry1 knockout oscillations display partially rhythmic behavior, and
cannot be classified as simply rhythmic or arrhythmic. Our model demonstrates
that intrinsic molecular noise and differences in relative abundance, rather
than differing functions, are sufficient to explain the range of rhythmicity
encountered in Cry knockouts in the SCN. Our results further highlight the
essential role of stochastic behavior in understanding and accurately modeling
the circadian network and its response to perturbation.


  We identify a new natural coalescent structure, which we call the seed-bank
coalescent, that describes the gene genealogy of populations under the
influence of a strong seed-bank effect, where "dormant forms" of individuals
(such as seeds or spores) may jump a significant number of generations before
joining the "active" population. Mathematically, our seed-bank coalescent
appears as scaling limit in a Wright-Fisher model with geometric seed-bank age
structure if the average time of seed dormancy scales with the order of the
total population size $N$. This extends earlier results of Kaj, Krone and
Lascoux [J. Appl. Probab. 38 (2011) 285-300] who show that the genealogy of a
Wright-Fisher model in the presence of a "weak" seed-bank effect is given by a
suitably time-changed Kingman coalescent. The qualitatively new feature of the
seed-bank coalescent is that ancestral lineages are independently blocked at a
certain rate from taking part in coalescence events, thus strongly altering the
predictions of classical coalescent models. In particular, the seed-bank
coalescent "does not come down from infinity," and the time to the most recent
common ancestor of a sample of size $n$ grows like $\log\log n$. This is in
line with the empirical observation that seed-banks drastically increase
genetic variability in a population and indicates how they may serve as a
buffer against other evolutionary forces such as genetic drift and selection.


  Empirical evidence shows that the rate of irregular usage of English verbs
exhibits discontinuity as a function of their frequency: the most frequent
verbs tend to be totally irregular. We aim to qualitatively understand the
origin of this feature by studying simple agent--based models of language
dynamics, where each agent adopts an inflectional state for a verb and may
change it upon interaction with other agents. At the same time, agents are
replaced at some rate by new agents adopting the regular form. In models with
only two inflectional states (regular and irregular), we observe that either
all verbs regularize irrespective of their frequency, or a continuous
transition occurs between a low frequency state where the lemma becomes fully
regular, and a high frequency one where both forms coexist. Introducing a third
(mixed) state, wherein agents may use either form, we find that a third,
qualitatively different behavior may emerge, namely, a discontinuous transition
in frequency. We introduce and solve analytically a very general class of
three--state models that allows us to fully understand these behaviors in a
unified framework. Realistic sets of interaction rules, including the
well-known Naming Game (NG) model, result in a discontinuous transition, in
agreement with recent empirical findings. We also point out that the
distinction between speaker and hearer in the interaction has no effect on the
collective behavior. The results for the general three--state model, although
discussed in terms of language dynamics, are widely applicable.


  So far mutations analysis was performed in terms of transitions and
trasversions, so on the basis of the molecule, or in terms of GC-content and
isochors, through the quantification of GC->AT mutations over AT->GC mutations.
We tried a different approach hypothesizing that probably we are partially
protected (by the genetic code) from amyloid fibrils formation.


  Fundamental biological processes such as transcription and translation, where
a genetic sequence is sequentially read by a macromolecule, have been well
described by a classical model of non-equilibrium statistical physics, the
totally asymmetric exclusion principle (TASEP). This model describes particles
hopping between sites of a one-dimensional lattice, with the particle current
determining the transcription or translation rate. An open problem is how to
analyze a TASEP where particles can pause randomly, as has been observed during
transcription. In this work, we report that surprisingly, a simple mean-field
model predicts well the particle current for all values of the average pause
duration, using a simple description of blocking behind paused particles.


  In explaining altruistic cooperation and punishment, the challenging riddle
is how transcendental rules can emerge within the empirical world.
Game-theoretical studies showed that pool punishment, particularly second-order
punishment, plays a key role in understanding the evolution of cooperation.
Second-order pool punishment, however, is tautological in nature; the
punishment system itself is caused by its own effects. The emergence of pool
punishment poses a logical conundrum that, to date, has been overlooked in the
study of the evolution of social norms and institutions. Here, we tackle the
issue by considering the interplay of (a) cognitive biases in reasoning and (b)
Agamben's notion of homo sacer (Agamben, G. 1998, Homo Sacer: Sovereign Power
and Bare Life, Stanford University Press); that is, a person who may be killed
without legal consequences. Based on the cognitive disposition of reversing the
cause-and-effect relationship, we propose a new system: the preemptive
punishment of homo sacer. This action can lead to retrospectively forming moral
assessments, particularly for second-order pool punishment. This study
potentially provides new insight into the evolutionary process behind the
emergence of social norms and institutions.


  Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately
regulated by the type of ligands that bind to the orthosteric or allosteric
binding sites. To glean the structural and dynamical origin of ligand-dependent
modulation of GPCR activity, we performed total $\sim$ 5 $\mu$sec molecular
dynamics simulations of A$_{2A}$ adenosine receptor (A$_{2A}$AR) in its apo,
antagonist-bound, and agonist-bound forms in an explicit water and membrane
environment, and examined the corresponding dynamics and correlation between
the 10 key structural motifs that serve as the allosteric hotspots in
intramolecular signaling network. We dubbed these 10 structural motifs "binary
switches" as they display molecular interactions that switch between two
distinct states. By projecting the receptor dynamics on these binary switches
that yield $2^{10}$ microstates, we show that (i) the receptors in apo,
antagonist-bound, and agonist-bound states explore vastly different
conformational space; (ii) among the three receptor states the apo state
explores the broadest range of microstates; (iii) in the presence of the
agonist, the active conformation is maintained through coherent couplings among
the binary switches; and (iv) to be most specific, our analysis shows that
W246, located deep inside the binding cleft, can serve as both an agonist
sensor and actuator of ensuing intramolecular signaling for the receptor
activation.Finally, our analysis of multiple trajectories generated by
inserting an agonist to the apo state underscores that the transition of the
receptor from inactive to active form requires the disruption of ionic-lock in
the DRY motif.


  The self-alignment concentrations, $c(x)$, as functions of the length, $x$,
of the identically matching maximal segments in the genomes of a variety of
species, typically present power-law tails extending to the largest scales,
i.e., $c(x) \propto x^{\alpha}$, with similar or apparently different negative
$\alpha$s ($<-2$). The relevant fundamental processes of molecular evolution
are segmental duplication and point mutation, and that recently the stick
fragmentation phenomenology has been used to account the neutral evolution.
However, disorder is intrinsic to the evolution system and, by freezing it in
time (quenching) for the setup of a simple fragmentation model, we obtain
decaying, steady-state and the general full time-dependent solutions, all
$\propto x^{\alpha}$ for $x\to \infty$, which is in contrast to the only
power-law solution, $x^{-3}$ for $x\to 0$ of the pure model (without disorder).
%Other algebraic terms may dominate at intermediate scales, which seems to be
confirmed by some species, such as rice. We also present self-alignment results
showing more than one scaling regimes, consistent with the theoretical results
of the existence of more than one algebraic terms which dominate at different
regimes.


  A crucial problem in genome assembly is the discovery and correction of
misassembly errors in draft genomes. We develop a method that will enhance the
quality of draft genomes by identifying and removing misassembly errors using
paired short read sequence data and optical mapping data. We apply our method
to various assemblies of the loblolly pine and Francisella tularensis genomes.
Our results demonstrate that we detect more than 54% of extensively
misassembled contigs and more than 60% of locally misassembed contigs in an
assembly of Francisella tularensis, and between 31% and 100% of extensively
misassembled contigs and between 57% and 73% of locally misassembed contigs in
the assemblies of loblolly pine. MISSEQUEL can be downloaded at
http://www.cs.colostate.edu/seq/.


  Human serum transferrin (hTf) transports ferric ions in the blood stream and
inflamed mucosal surfaces with high affinity and delivers them to cells via
receptor mediated endocytosis. A typical hTf is folded into two homologous
lobes; each lobe is further divided into two similar sized domains. Three
different crystal structures of hTf delineate large conformational changes
involved in iron binding/dissociation. However, whether the release process
follows the same trend at serum (~7.4) and endosomal (~5.6) pH remains
unanswered. The specialized role of the two lobes and if communication between
them leads to efficient and controlled release is also debated. Here, we study
the dynamics of the full structure as well as the separate lobes in different
closed, partially open, and open conformations under the nearly neutral pH
conditions in the blood serum and the more acidic one in the endosome. The
results corroborate experimental observations and underscore the distinguishing
effect of pH on the dynamics of hTf. Furthermore, in a total of 2 {\mu}s
molecular dynamics simulation of different forms of hTf, residue fluctuations
elucidate the cross talk between the two lobes correlated by the peptide linker
bridging the two lobes at serum pH, while their correlated motions is lost
under endosomal conditions. At serum pH, the presence of even a single iron on
either lobe leads C-lobe fluctuations to subside, making it the target for
recognition by human cells or hostile bacteria seeking iron sequestration. The
N-lobe, on the other hand, has a propensity to open, making iron readily
available when needed at regions of serum pH. At endosomal pH, both lobes
readily open, making irons available for delivery. The interplay between the
relative mobility of the lobes renders efficient mechanism for the recognition
and release of hTf at the cell surface, and therefore its recycling in the
organism.


  Traditionally, phylogeny and sequence alignment are estimated separately:
first estimate a multiple sequence alignment and then infer a phylogeny based
on the sequence alignment estimated in the previous step. However, uncertainty
in the alignment estimation is ignored, resulting, possibly, in overstated
certainty in phylogeny estimates. We develop a joint model for co-estimating
phylogeny and sequence alignment which improves estimates from the traditional
approach by accounting for uncertainty in the alignment in phylogenetic
inferences. Our insertion and deletion (indel) model allows arbitrary-length
overlapping indel events and a general distribution for indel fragment size. We
employ a Bayesian approach using MCMC to estimate the joint posterior
distribution of a phylogenetic tree and a multiple sequence alignment. Our
approach has a tree and a complete history of indel events mapped onto the tree
as the state space of the Markov Chain while alternative previous approaches
have a tree and an alignment. A large state space containing a complete history
of indel events makes our MCMC approach more challenging, but it enables us to
infer more information about the indel process. The performances of this joint
method and traditional sequential methods are compared using simulated data as
well as real data. Software named BayesCAT (Bayesian Co-estimation of Alignment
and Tree) is available at https://github.com/heejungshim/BayesCAT.


  According to Cobanoglu et al and Murphy, it is now widely acknowledged that
the single target paradigm (one protein or target, one disease, one drug) that
has been the dominant premise in drug development in the recent past is
untenable. More often than not, a drug-like compound (ligand) can be
promiscuous - that is, it can interact with more than one target protein. In
recent years, in in silico target prediction methods the promiscuity issue has
been approached computationally in different ways. In this study we confine
attention to the so-called ligand-based target prediction machine learning
approaches, commonly referred to as target-fishing. With a few exceptions, the
target-fishing approaches that are currently ubiquitous in cheminformatics
literature can be essentially viewed as single-label multi-classification
schemes; these approaches inherently bank on the single target paradigm
assumption that a ligand can home in on one specific target. In order to
address the ligand promiscuity issue, one might be able to cast target-fishing
as a multi-label multi-class classification problem. For illustrative and
comparison purposes, single-label and multi-label Naive Bayes classification
models (denoted here by SMM and MMM, respectively) for target-fishing were
implemented. The models were constructed and tested on 65,587 compounds and 308
targets retrieved from the ChEMBL17 database. SMM and MMM performed
differently: for 16,344 test compounds, the MMM model returned recall and
precision values of 0.8058 and 0.6622, respectively; the corresponding recall
and precision values yielded by the SMM model were 0.7805 and 0.7596,
respectively. However, at a significance level of 0.05 and one degree of
freedom McNemar test performed on the target prediction results returned by SMM
and MMM for the 16,344 test ligands gave a chi-squared value of 15.656, in
favour of the MMM approach.


  A key question in biological systems is whether genetic diversity persists in
the long run under evolutionary competition or whether a single dominant
genotype emerges. Classic work by Kalmus in 1945 has established that even in
simple diploid species (species with two chromosomes) diversity can be
guaranteed as long as the heterozygote individuals enjoy a selective advantage.
Despite the classic nature of the problem, as we move towards increasingly
polymorphic traits (e.g. human blood types) predicting diversity and
understanding its implications is still not fully understood. Our key
contribution is to establish complexity theoretic hardness results implying
that even in the textbook case of single locus diploid models predicting
whether diversity survives or not given its fitness landscape is
algorithmically intractable. We complement our results by establishing that
under randomly chosen fitness landscapes diversity survives with significant
probability. Our results are structurally robust along several dimensions
(e.g., choice of parameter distribution, different definitions of
stability/persistence, restriction to typical subclasses of fitness
landscapes). Technically, our results exploit connections between game theory,
nonlinear dynamical systems, complexity theory and biology and establish
hardness results for predicting the evolution of a deterministic variant of the
well known multiplicative weights update algorithm in symmetric coordination
games which could be of independent interest.


  Epistasis refers to the phenomenon in which phenotypic consequences caused by
mutation of one gene depend on one or more mutations at another gene. Epistasis
is critical for understanding many genetic and evolutionary processes,
including pathway organization, evolution of sexual reproduction, mutational
load, ploidy, genomic complexity, speciation, and the origin of life.
Nevertheless, current understandings for the genome-wide distribution of
epistasis are mostly inferred from interactions among one mutant type per gene,
whereas how epistatic interaction partners change dynamically for different
mutant alleles of the same gene is largely unknown. Here we address this issue
by combining predictions from flux balance analysis and data from a recently
published high-throughput experiment. Our results show that different alleles
can epistatically interact with very different gene sets. Furthermore, between
two random mutant alleles of the same gene, the chance for the allele with more
severe mutational consequence to develop a higher percentage of negative
epistasis than the other allele is 50-70% in eukaryotic organisms, but only
20-30% in bacteria and archaea. We developed a population genetics model that
predicts that the observed distribution for the sign of epistasis can speed up
the process of purging deleterious mutations in eukaryotic organisms. Our
results indicate that epistasis among genes can be dynamically rewired at the
genome level, and call on future efforts to revisit theories that can integrate
epistatic dynamics among genes in biological systems.


  Epistasis describes the phenomenon that mutations at different loci do not
have independent effects with regard to certain phenotypes. Understanding the
global epistatic landscape is vital for many genetic and evolutionary theories.
Current knowledge for epistatic dynamics under multiple conditions is limited
by the technological difficulties in experimentally screening epistatic
relations among genes. We explored this issue by applying flux balance analysis
to simulate epistatic landscapes under various environmental perturbations.
Specifically, we looked at gene-gene epistatic interactions, where the
mutations were assumed to occur in different genes. We predicted that epistasis
tends to become more positive from glucose-abundant to nutrient-limiting
conditions, indicating that selection might be less effective in removing
deleterious mutations in the latter. We also observed a stable core of
epistatic interactions in all tested conditions, as well as many epistatic
interactions unique to each condition. Interestingly, genes in the stable
epistatic interaction network are directly linked to most other genes whereas
genes with condition-specific epistasis form a scale-free network. Furthermore,
genes with stable epistasis tend to have similar evolutionary rates, whereas
this co-evolving relationship does not hold for genes with condition-specific
epistasis. Our findings provide a novel genome-wide picture about epistatic
dynamics under environmental perturbations.


  Bioinformatics, as an emerging and rapidly developing interdisciplinary, has
become a promising and popular research field in 21st century. Extracting and
explaining useful biological information from huge amount of genetic data is an
urgent issue in post-genome era. In eukaryotic DNA sequences, gene consists of
exons and introns. To predict the location of exons which carry most genetic
information accurately has become one of the most essential issues in
bioinformatics. Here, we have used biological characteristics of introns to
find the candidate initial and final exon sections. Then we select candidate
exon sections by using Support Vector Machine (SVM). Next, we predict exon
sections accurately based on Discrete Fourier Transform (DFT) and using
three-base periodicity of DNA sequence signals. This paper provides a gene
prediction method based on statistics and signal processing and also, the
improvement and prospect for this method in the future are discussed.


  Biodiversity and extinction are central issues in evolution. Dynamical
balance among different species in ecosystems is often described by
deterministic replicator equations with moderate success. However, fluctuations
are inevitable, either caused by external environment or intrinsic random
competitions in finite populations, and the evolutionary dynamics is stochastic
in nature. Here we show that, after appropriate coarse-graining, random
fluctuations generate dissipation towards extinction because the evolution
trajectories in the phase space of all competing species possess positive
curvature. As a demonstrating example, we compare the fluctuation-induced
dissipative dynamics in Lotka-Volterra model with numerical simulations and
find impressive agreement. Our finding is closely related to the
fluctuation-dissipation theorem in statistical mechanics but the marked
difference is the non-equilibrium essence of the generic evolutionary dynamics.
As the evolving ecosystems are far from equilibrium, the relation between
fluctuations and dissipations is often complicated and dependent on microscopic
details. It is thus remarkable that the generic positivity of the trajectory
curvature warrants dissipation arisen from the seemingly harmless fluctuations.
The unexpected dissipative dynamics is beyond the reach of conventional
replicator equations and plays a crucial role in investigating the biodiversity
in ecosystems.


  Viruses can amplify their genomes following different replication modes (RMs)
ranging from the stamping machine replication (SMR) model to the geometric
replication (GR) model. Different RMs are expected to produce different
evolutionary and dynamical outcomes in viral quasispecies due to differences in
the mutations accumulation rate. Theoretical and computational models revealed
that while SMR may provide RNA viruses with mutational robustness, GR may
confer a dynamical advantage against genomes degradation. Here, recent advances
in the investigation of the RM in positive-sense single-stranded RNA viruses
are reviewed. Dynamical experimental quantification of Turnip mosaic virus RNA
strands, together with a nonlinear mathematical model, indicated the SMR model
for this pathogen. The same mathematical model for natural infections is here
further analyzed, and we prove that the interior equilibrium involving
coexistence of both positive and negative viral strands is globally
asymptotically stable.


  Cancer is a complex disease and thus is complicated to model. However, simple
models that describe the main processes involved in tumoral dynamics, e.g.,
competition and mutation, can give us clues about cancer behaviour, at least
qualitatively, also allowing us to make predictions. Here we analyze a
simplified quasispecies mathematical model given by differential equations
describing the time behaviour of tumor cells populations with different levels
of genomic instability. We find the equilibrium points, also characterizing
their stability and bifurcations focusing on replication and mutation rates. We
identify a transcritical bifurcation at increasing mutation rates of the tumor
cells population. Such a bifurcation involves an scenario with dominance of
healthy cells and impairment of tumor populations. Finally, we characterize the
transient times for this scenario, showing that a slight increase beyond the
critical mutation rate may be enough to have a fast response towards the
desired state (i.e., low tumor populations) during directed mutagenic
therapies.


  Acquired resistance is one of the major barriers to successful cancer
therapy. The development of resistance is commonly attributed to genetic
heterogeneity. However, heterogeneity of drug penetration of the tumor
microenvironment both on the microscopic level within solid tumors as well as
on the macroscopic level across metastases may also contribute to acquired drug
resistance. Here we use mathematical models to investigate the effect of drug
heterogeneity on the probability of escape from treatment and time to
resistance. Specifically we address scenarios with sufficiently efficient
therapies that suppress growth of all preexisting genetic variants in the
compartment with highest drug concentration. To study the joint effect of drug
heterogeneity, growth rate, and evolution of resistance we analyze a multitype
stochastic branching process describing growth of cancer cells in two
compartments with different drug concentration and limited migration between
compartments. We show that resistance is more likely to arise first in the low
drug compartment and from there populate the high drug compartment. Moreover,
we show that only below a threshold rate of cell migration does spatial
heterogeneity accelerate resistance evolution, otherwise deterring drug
resistance with excessively high migration rates. Our results provide new
insights into understanding why cancers tend to quickly become resistant, and
that cell migration and the presence of sanctuary sites with little drug
exposure are essential to this end.


  Ticks are important vectors of emerging zoonotic diseases. While adults of
many tick species parasitize mammals, immature ticks are often found on wild
birds. In the tropics, difficulties in species-level identification of immature
ticks hinder studies of tick ecology and tick-borne disease transmission,
including any potential role for birds. In Panama, we found immature ticks on
227 out of 3,498 birds representing 93 host species, about 1/8th of the entire
Panamanian terrestrial avifauna. Tick parasitism rates did not vary with
temperature or rainfall, but parasitism rates did vary with host ecological
traits: non-migratory residents, forest dwelling birds, bark insectivores,
terrestrial foragers and lowland species were most likely to be infested with
ticks. Using a molecular library developed from adult ticks specifically for
this study, we identified 130 immature ticks obtained from wild birds,
corresponding to eleven tick species, indicating that a substantial portion of
the Panamanian avifauna is parasitized by a variety of tick species.
Furthermore, we found evidence that immature ticks show taxonomic or ecological
specificity to avian hosts. Finally, our data indicate that Panamanian birds
are not parasitized regularly by the tick species responsible for most known
tick-borne diseases. However, they are frequent hosts of other tick species
known to carry a variety of rickettsial parasites of unknown pathogenicity.
Given the broad interaction between tick and avian biodiversity in the
Neotropics, future work on emerging tropical tick-borne disease should
explicitly consider wild birds as vertebrate hosts.


  We consider a generic proto-cell model consisting of any conservative
chemical reaction network embedded within a membrane. The membrane results from
the self-assembly of one of the chemical species (membrane precursor) and is
semi-permeable to some other chemical species (nutrients) diffusing from an
outside growth medium into the proto-cell. Inside the proto-cell, nutrients are
metabolized into all other chemical species including the membrane precursor,
and the membrane grows in area and the proto-cell in volume. Investigating the
conditions under which such a proto-cell may reach stationary growth, we prove
that a simple necessary condition is that each moiety be fed with some nutrient
flux; and that a sufficient condition for the existence of a stationary growth
regime is that every siphon containing any species participating in the
membrane precursor incorporation kinetics also contains the support of a moiety
that is fed with some nutrient flux. These necessary and sufficient conditions
hold regardless of chemical reaction kinetics, membrane parameters or nutrient
flux diffusion characteristics.


  Binets and trinets are phylogenetic networks with two and three leaves,
respectively. Here we consider the problem of deciding if there exists a binary
level-1 phylogenetic network displaying a given set $\mathcal{T}$ of binary
binets or trinets over a set $X$ of taxa, and constructing such a network
whenever it exists. We show that this is NP-hard for trinets but
polynomial-time solvable for binets. Moreover, we show that the problem is
still polynomial-time solvable for inputs consisting of binets and trinets as
long as the cycles in the trinets have size three. Finally, we present an
$O(3^{|X|} poly(|X|))$ time algorithm for general sets of binets and trinets.
The latter two algorithms generalise to instances containing level-1 networks
with arbitrarily many leaves, and thus provide some of the first supernetwork
algorithms for computing networks from a set of rooted phylogenetic networks.


  An increasing number of dissident voices claim that the standard
neo-Darwinian view of genes as 'leaders' and phenotypes as 'followers' during
the process of adaptive evolution should be turned on its head. This idea is
older than the rediscovery of Mendel's laws of inheritance and has been given
several names before its final 'Baldwin effect' label. A condition for this
effect is that environmentally induced variation such as phenotypic plasticity
or learning is crucial for the initial establishment of a population. This
gives the necessary time for natural selection to act on genetic variation and
the adaptive trait can be eventually encoded in the genotype. An influential
paper published in the late 1980s showed the Baldwin effect to happen in
computer simulations, and claimed that it was crucial to solve a difficult
adaptive task. This generated much excitement among scholars in various
disciplines that regard neo-Darwinian accounts to explain the evolutionary
emergence of high-order phenotypic traits such as consciousness or language
almost hopeless. Here, we use analytical and computational approaches to show
that a standard population genetics treatment can easily crack what the
scientific community has granted as an unsolvable adaptive problem without
learning. The Baldwin effect is once again in need of convincing theoretical
foundations.


  The molecular chaperone Hsp70 plays a central role in the import of
cytoplasmic proteins into organelles, driving their translocation by binding
them from the organellar interior. Starting from the experimentally-determined
structure of the \textit{E. coli} Hsp70, we computed, by means of molecular
simulations, the effective free-energy profile for substrate translocation upon
chaperone binding. We then used the resulting free energy to quantitatively
characterize the kinetics of the import process, whose comparison with
unassisted translocation highlights the essential role played by Hsp70 in
importing cytoplasmic proteins.


  Both multi-cell biochemical assays and single cell fluorescence measurements
have revealed that the elongation rate of Polymerase II (PolII) in eukaryotes
varies largely across different cell types and genes. However, there is not yet
a consensus whether intrinsic factors such as the position, local mobility or
the engagement by an active molecular mechanism of a genetic locus could be the
determinants of the observed heterogeneity. Employing high-speed 3D
fluorescence nanoimaging we resolve here at the single cell level multiple,
distinct regions of mRNA synthesis within a labeled transgene array. By
employing phasor analysis, a fluorescence fluctuation spectroscopy technique,
we demonstrate that these regions are active transcription sites that release
mRNA molecules in the nucleoplasm, and we extract the local PolII elongation
rate. While we detect a range of 10-100 bp/s for PolII elongation from cell to
cell, we are now also able to measure up to a four-fold variation in the
average elongation between identical copies of the same gene measured
simultaneously within the same cell. Furthermore, we are able to visualize
changes of PolII elongation as a function of time. We observe a correlation
between the average elongation rate measured in a locus and its local mobility.
Finally, by cross-correlating the transcriptional activity with the nm-sized
movements of the active loci, we provide evidence of an active molecular
mechanism determining displacements of the transcription sites concomitant to
increases in transcriptional activity. Together these observations demonstrate
that local factors, such as chromatin local mobility and the microenvironment
of the transcription site, are an important source of transcription kinetics
variability.


  Stochastic encounter-mating (SEM) models describe monogamous permanent pair
formation in finite zoological populations of multitype females and males. In
this article, we study SEM with Poisson firing times. First, we prove that the
model enjoys a fluid limit as the population size diverges, i.e., the
stochastic dynamics converges to a deterministic system governed by coupled
ODEs. Then, we convert these ODEs to the well-known Lotka-Volterra and
replicator equations from population dynamics. Next, under the so-called fine
balance condition which characterizes panmixia, we solve the corresponding
replicator equations and give an exact expression for the fluid limit. Finally,
we consider the case with two types of females and males. Without the fine
balance assumption, but under certain symmetry conditions, we give an explicit
formula for the limiting mating pattern, and then use it to characterize
assortative mating.


  We present new empirical evidence, based on millions of interactions on
Twitter, confirming that human contacts scale with population sizes. We
integrate such observations into a reaction-diffusion metapopulation framework
providing an analytical expression for the global invasion threshold of a
contagion process. Remarkably, the scaling of human contacts is found to
facilitate the spreading dynamics. Our results show that the scaling properties
of human interactions can significantly affect dynamical processes mediated by
human contacts such as the spread of diseases, and ideas.


  - Genome-scan methods are used for screening genome-wide patterns of DNA
polymorphism to detect signatures of positive selection. There are two main
types of methods: (i) 'outlier' detection methods based on Fst that detect loci
with high differentiation compared to the rest of the genomes, and (ii)
environmental association methods that test the association between allele
frequencies and environmental variables.
  - We present a new Fst-based genome-scan method, BayeScEnv, which
incorporates environmental information in the form of 'environmental
differentiation'. It is based on the F-model, but, as opposed to existing
approaches, it considers two locus-specific effects; one due to divergent
selection, and another one due to various other processes different from local
adaptation (e.g. range expansions, differences in mutation rates across loci or
background selection). The method was developped in C++ and is avaible at
http://github.com/devillemereuil/bayescenv.
  - Simulation studies shows that our method has a much lower false positive
rate than an existing Fst-based method, BayeScan, under a wide range of
demographic scenarios. Although it has lower power, it leads to a better
compromise between power and false positive rate.
  - We apply our method to human and salmon datasets and show that it can be
used successfully to study local adaptation. We discuss its scope and compare
its mechanics to other existing methods.


  In this paper, we focus on the influence of heterogeneity and stochasticity
of the population on the dynamical structure of a basic
susceptible-infected-susceptible (SIS) model. First we prove that, upon a
suitable mathematical reformulation of the basic reproduction number, the
homogeneous system and the heterogeneous system exhibit a completely analogous
global behaviour. Then we consider noise terms to incorporate the fluctuation
effects and the random import of the disease into the population and analyse
the influence of heterogeneity on warning signs for critical transitions (or
tipping points). This theory shows that one may be able to anticipate whether a
bifurcation point is close before it happens. We use numerical simulations of a
stochastic fast-slow heterogeneous population SIS model and show various
aspects of heterogeneity have crucial influences on the scaling laws that are
used as early-warning signs for the homogeneous system. Thus, although the
basic structural qualitative dynamical properties are the same for both
systems, the quantitative features for epidemic prediction are expected to
change and care has to be taken to interpret potential warning signs for
disease outbreaks correctly.


  We prove polynomial upper and lower bounds on the expected size of the
maximum agreement subtree of two random binary phylogenetic trees under both
the uniform distribution and Yule-Harding distribution. This positively answers
a question posed in earlier work. Determining tight upper and lower bounds
remains an open problem.


  Experiments show that cellular uptake of nanoparticles, via receptor-mediated
endocytosis, strongly depends on nanoparticle size. There is an optimal size,
approximately 50 nm in diameter, at which cellular uptake is the highest. In
addition, there is a maximum size, approximately 200 nm, beyond which uptake
via receptor-mediated endocytosis does not occur. By comparing results from
different experiments, we found that these sizes weakly depend on the type of
cells, nanoparticles, and ligands used in the experiments. Here, we argue that
these observations are consequences of the energetics and assembly dynamics of
the protein coat that forms on the cytoplasmic side of the outer cell membrane
during receptor-mediated endocytosis. Specifically, we show that the energetics
of coat formation imposes an upper bound on the size of the nanoparticles that
can be internalized, whereas the nanoparticle-size-dependent dynamics of coat
assembly results in the optimal nanoparticle size. The weak dependence of the
optimal and maximum sizes on cell-nanoparticle-ligand type also follows
naturally from our analysis.


  Background: While many infectious disease epidemics are initially
characterized by an exponential growth in time, we show that district-level
Ebola virus disease (EVD) outbreaks in West Africa follow slower
polynomial-based growth kinetics over several generations of the disease.
Methods: We analyzed epidemic growth patterns at three different spatial scales
(regional, national, and subnational) of the Ebola virus disease epidemic in
Guinea, Sierra Leone and Liberia by compiling publicly available weekly time
series of reported EVD case numbers from the patient database available from
the World Health Organization website for the period 05-Jan to 17-Dec 2014.
Results: We found significant differences in the growth patterns of EVD cases
at the scale of the country, district, and other subnational administrative
divisions. The national cumulative curves of EVD cases in Guinea, Sierra Leone,
and Liberia show periods of approximate exponential growth. In contrast, local
epidemics are asynchronous and exhibit slow growth patterns during 3 or more
EVD generations, which can be better approximated by a polynomial than an
exponential. Conclusions: The slower than expected growth pattern of local EVD
outbreaks could result from a variety of factors, including behavior changes,
success of control interventions, or intrinsic features of the disease such as
a high level of clustering. Quantifying the contribution of each of these
factors could help refine estimates of final epidemic size and the relative
impact of different mitigation efforts in current and future EVD outbreaks.


  It is well-known that the formation of amyloid fiber may cause invertible
damage to cells, while the underlying mechanism has not been fully uncovered.
In this paper, we construct a mathematical model, consisting of infinite ODEs
in the form of mass-action equations together with two reaction-convection
PDEs, and then simplify it to a system of 5 ODEs by using the maximum entropy
principle. This model is based on four simple assumptions, one of which is that
cell damage is raised by oligomers rather than mature fibrils. With the
simplified model, the effects of nucleation and elongation, fragmentation,
protein and seeds concentrations on amyloid formation and cell damage are
extensively explored and compared with experiments. We hope that our results
can provide a valuable insight into the processes of amyloid formation and cell
damage thus raised.


  This chapter reviews the differential geometry-based solvation and
electrolyte transport for biomolecular solvation that have been developed over
the past decade. A key component of these methods is the differential geometry
of surfaces theory, as applied to the solvent-solute boundary. In these
approaches, the solvent-solute boundary is determined by a variational
principle that determines the major physical observables of interest, for
example, biomolecular surface area, enclosed volume, electrostatic potential,
ion density, electron density, etc. Recently, differential geometry theory has
been used to define the surfaces that separate the microscopic (solute) domains
for biomolecules from the macroscopic (solvent) domains. In these approaches,
the microscopic domains are modeled with atomistic or quantum mechanical
descriptions, while continuum mechanics models (including fluid mechanics,
elastic mechanics, and continuum electrostatics) are applied to the macroscopic
domains. This multiphysics description is integrated through an energy
functional formalism and the resulting Euler-Lagrange equation is employed to
derive a variety of governing partial differential equations for different
solvation and transport processes; e.g., the Laplace-Beltrami equation for the
solvent-solute interface, Poisson or Poisson-Boltzmann equations for
electrostatic potentials, the Nernst-Planck equation for ion densities, and the
Kohn-Sham equation for solute electron density. Extensive validation of these
models has been carried out over hundreds of molecules, including proteins and
ion channels, and the experimental data have been compared in terms of
solvation energies, voltage-current curves, and density distributions. We also
propose a new quantum model for electrolyte transport.


  The case of fisheries management illustrates how the inherent structural
instability of ecosystems can have deep-running policy implications. We
contrast ten types of management plans to achieve maximum sustainable yields
(MSY) from multiple stocks and compare their effectiveness based on a
management strategy evaluation (MSE) that uses complex food webs in its
operating model. Plans that target specific stock sizes ($B_{\text{MSY}}$)
consistently led to higher yields than plans targeting specific fishing
pressures ($F_{\text{MSY}}$). A new self-optimising control rule, introduced
here for its robustness to structural instability, led to intermediate yields.
Most plans outperformed single-species management plans with pressure targets
set without considering multispecies interactions. However, more refined plans
to "maximise the yield from each stock separately", in the sense of a Nash
equilibrium, produced total yields comparable to plans aiming to maximise total
harvested biomass, and were more robust to structural instability. Our analyses
highlight trade-offs between yields, amenability to negotiations, pressures on
biodiversity, and continuity with current approaches in the European context.
Based on these results, we recommend directions for developments of EU
fisheries policy.


  Mathematical modelling can help to explain the nature and dynamics of
infection transmissions, as well as support a policy for implementing those
strategies that are most likely to bring public health and economic benefits.
The paper addresses the application of optimal control strategies in a
tuberculosis model. The model consists of a system of ordinary differential
equations, which considers reinfection and post-exposure interventions. We
propose a multiobjective optimization approach to find optimal control
strategies for the minimization of active infectious and persistent latent
individuals, as well as the cost associated to the implementation of the
control strategies. Optimal control strategies are investigated for different
values of the model parameters. The obtained numerical results cover a whole
range of the optimal control strategies, providing valuable information about
the tuberculosis dynamics and showing the usefulness of the proposed approach.


  The relationship between network topology and system dynamics has significant
implications for unifying our understanding of the interplay among metabolic,
gene-regulatory, and ecosystem network architecures. Here we analyze the
stability and robustness of a large class of dynamics on such networks. We
determine the probability distribution of robustness as a function of network
topology and show that robustness is classified by the number of links between
modules of the network. We also demonstrate that permutation of these modules
is a fundamental symmetry of dynamical robustness. Analysis of these findings
leads to the conclusion that the most robust systems have the most hierarchical
structure. This relationship provides a means by which evolutionary selection
for a purely dynamical phenomenon may shape network architectures across scales
of the biological hierarchy.


  We study a variant of Hanski's incidence function model that allows habitat
patch characteristics to vary over time following a Markov process. The widely
studied case where patches are classified as either suitable or unsuitable is
included as a special case. For large metapopulations, we determine a recursion
for the probability that a given habitat patch is occupied. This recursion
enables us to clarify the role of landscape dynamics in the survival of a
metapopulation. In particular, we show that landscape dynamics affects the
persistence and equilibrium level of the metapopulation primarily through its
effect on the distribution of a local population's life span.


  This paper addresses the decomposition of biochemical networks into
functional modules that preserve their dynamic properties upon interconnection
with other modules, which permits the inference of network behavior from the
properties of its constituent modules. The modular decomposition method
developed here also has the property that any changes in the parameters of a
chemical reaction only affect the dynamics of a single module. To illustrate
our results, we define and analyze a few key biological modules that arise in
gene regulation, enzymatic networks, and signaling pathways. We also provide a
collection of examples that demonstrate how the behavior of a biological
network can be deduced from the properties of its constituent modules, based on
results from control systems theory.


  Practically, all chemotherapeutic agents lead to drug resistance. Clinically,
it is a challenge to determine whether resistance arises prior to, or as a
result of, cancer therapy. Further, a number of different intracellular and
microenvironmental factors have been correlated with the emergence of drug
resistance. With the goal of better understanding drug resistance and its
connection with the tumor microenvironment, we have developed a hybrid
discrete-continuous mathematical model. In this model, cancer cells described
through a particle-spring approach respond to dynamically changing oxygen and
DNA damaging drug concentrations described through partial differential
equations. We thoroughly explored the behavior of our self-calibrated model
under the following common conditions: a fixed layout of the vasculature, an
identical initial configuration of cancer cells, the same mechanism of drug
action, and one mechanism of cellular response to the drug. We considered one
set of simulations in which drug resistance existed prior to the start of
treatment, and another set in which drug resistance is acquired in response to
treatment. This allows us to compare how both kinds of resistance influence the
spatial and temporal dynamics of the developing tumor, and its clonal
diversity. We show that both pre-existing and acquired resistance can give rise
to three biologically distinct parameter regimes: successful tumor eradication,
reduced effectiveness of drug during the course of treatment (resistance), and
complete treatment failure.


  Protein oligomers have been implicated as toxic agents in a wide range of
amyloid-related diseases. Yet it has remained unsolved whether the oligomers
are a necessary step in the formation of amyloid fibrils, or just a dangerous
by-product. Analogously, it has not been resolved if the amyloid nucleation
process is a classical one-step nucleation process, or a two-step process
involving pre-nucleation clusters. We use coarse-grained computer simulations
to study the effect of non-specific attractions between peptides on the primary
nucleation process underlying amyloid fibrillization. We find that for peptides
that do not attract, the classical one-step nucleation mechanism is possible,
but only at non-physiologically high peptide concentrations. At low peptide
concentrations, which mimic the physiologically relevant regime, attractive
inter-peptide interactions are essential for fibril formation. Nucleation then
inevitably takes place through a two-step mechanism involving prefibrillar
oligomers. We show that oligomers not only help peptides meet each other, but
create an environment that facilitates the conversion of monomers into the
$\beta$-sheet rich form characteristic of fibrils. Nucleation typically does
not proceed via the most prevalent oligomers, but via an oligomer size that is
only observed in rare fluctuations, which is why such aggregates might be hard
to capture experimentally. Finally, we find that the nucleation of amyloid
fibrils cannot be described by classical nucleation theory: in the two-step
mechanism the critical nucleus size increases both with an increase in
concentration and in the inter-peptide interactions, in direct contrast with
predictions from classical nucleation theory.


  As an emergent infectious disease outbreak unfolds, public health response is
reliant on information on key epidemiological quantities, such as transmission
potential and serial interval. Increasingly, transmission models fit to
incidence data are used to estimate these parameters and guide policy. Some
widely-used modeling practices lead to potentially large errors in parameter
estimates and, consequently, errors in model-based forecasts. Even more
worryingly, in such situations, confidence in parameter estimates and forecasts
can itself be far over-estimated, leading to the potential for large errors
that mask their own presence. Fortunately, straightforward and computationally
inexpensive alternatives exist that avoid these problems. Here, we first use a
simulation study to demonstrate potential pitfalls of the standard practice of
fitting deterministic models to cumulative incidence data. Next, we demonstrate
an alternative based on stochastic models fit to raw data from an early phase
of 2014 West Africa Ebola Virus Disease outbreak. We show not only that bias is
thereby reduced, but that uncertainty in estimates and forecasts is better
quantified and that, critically, lack of model fit is more readily diagnosed.
We conclude with a short list of principles to guide the modeling response to
future infectious disease outbreaks.


  Mixed feedback loops combining transcriptional and post-transcriptional
regulations are common in cellular regulatory networks. They consist of two
genes, encoding a transcription factor and a small non-coding RNA (sRNA), which
mutually regulate each other's expression. We present a theoretical and
numerical study of coherent mixed feedback loops of this type, in which both
regulations are negative. Under suitable conditions, these feedback loops are
expected to exhibit bistability, namely two stable states, one dominated by the
transcriptional repressor and the other dominated by the sRNA. We use
deterministic methods based on rate equation models, in order to identify the
range of parameters in which bistability takes place. However, the
deterministic models do not account for the finite lifetimes of the bistable
states and the spontaneous, fluctuation-driven transitions between them.
Therefore, we use stochastic methods to calculate the average lifetimes of the
two states. It is found that these lifetimes strongly depend on rate
coefficients such as the transcription rates of the transcriptional repressor
and the sRNA. In particular, we show that the fraction of time the system
spends in the sRNA dominated state follows a monotonically decreasing sigmoid
function of the transcriptional repressor transcription rate. The biological
relevance of these results is discussed in the context of such mixed feedback
loops in {\it Escherichia coli}.


  We propose a model of multispecies populations surviving on distributed
resources. System dynamics are investigated under changes in abiotic factors
such as the climate, as parameterized through environmental temperature. In
particular, we introduce a feedback between species abundances and resources
via abiotic factors. This model is apparently the first of its kind to include
a feedback mechanism coupling climate and population dynamics. Moreover, we
take into account self-limitation effects. The model explains the coexistence
of many species, yet also displays the possibility of catastrophic
bifurcations, where all species become extinct under the influence of abiotic
factors. We show that as these factors change there are different regimes of
ecosystem behavior, including a possibly chaotic regime when abiotic influences
are sufficiently strong.


  Despite the large quantity of information available, thorough researches in
various biological databases are still needed in order to reconstruct and
understand the steps that lead to known or new phenomena. By using
protein-protein interaction networks and algorithms to extract relevant
interconnections among proteins of interest, it is possible to assemble
subnetworks from global interactomes. Using these extracted networks it is
possible to use algorithms to predict signal directions while activation and
inhibition effects can be predicted using RNA interference screenings. The
result of this approach is the automatic generation of boolean networks. This
way of modelling dynamical systems allows the discovery of steady states and
the prediction of stimuli response.


  Classification standards such as the Mammal Species of the World (MSW) aim to
unify name usages at the global scale, but may nevertheless experience
significant levels of taxonomic change from one edition to the next. This
circumstance challenges the biodiversity and phylogenetic data communities to
develop more granular identifiers to track taxonomic congruence and
incongruence in ways that both humans and machines can process, i.e., to
logically represent taxonomic provenance across multiple classification
hierarchies. Here we show that reasoning over taxonomic provenance is feasible
for two classifications of primates corresponding to the second and third MSW
editions. Our approach entails three main components: (1) individuation of name
usages as taxonomic concepts, (2) articulation of concepts via human-asserted
Region Connection Calculus (RCC-5) relationships, and (3) the use of an Answer
Set Programming toolkit to infer and visualize logically consistent alignments
of these taxonomic input constraints. Our use case entails the Primates sec.
Groves (1993; MSW2 - 317 taxonomic concepts; 233 at the species level) and
Primates sec. Groves (2005; MSW3 - 483 taxonomic concepts; 376 at the species
level). Using 402 concept-to-concept input articulations, the reasoning process
yields a single, consistent alignment, and infers 153,111 Maximally Informative
Relations that constitute a comprehensive provenance resolution map for every
concept pair in the Primates sec. MSW2/MSW3. The entire alignment and various
partitions facilitate quantitative analyses of name/meaning dissociation,
revealing that approximately one in three paired name usages across treatments
is not reliable - in the sense of the same name identifying congruent taxonomic
meanings. We conclude with an optimistic outlook for logic-based provenance
tools in next-generation biodiversity and phylogeny data platforms.


  In biology and ecology, individuals or communities of individuals living in
unpredictable environments often alternate between different evolutionary
strategies to spread and reduce risks. Such behavior is commonly referred to as
"bet-hedging." Long-term survival probabilities and population sizes can be
much enhanced by exploiting such hybrid strategies. Here, we study the simplest
possible birth-death stochastic model in which individuals can choose among a
poor but safe strategy, a better but risky alternative, or a combination of
both. We show analytically and computationally that the benefits derived from
bet-hedging strategies are much stronger for higher environmental variabilities
(large external noise) and/or for small spatial dimensions (large intrinsic
noise). These circumstances are typically encountered by living systems, thus
providing us with a possible justification for the ubiquitousness of
bet-hedging in nature.


  The increased affordability of whole genome sequencing has motivated its use
for phenotypic studies. We address the problem of learning interpretable models
for discrete phenotypes from whole genomes. We propose a general approach that
relies on the Set Covering Machine and a k-mer representation of the genomes.
We show results for the problem of predicting the resistance of Pseudomonas
Aeruginosa, an important human pathogen, against 4 antibiotics. Our results
demonstrate that extremely sparse models which are biologically relevant can be
learnt using this approach.


  For more than 30 years, the only enzymatic system known to catalyze the
elimination of superoxide was superoxide dismutase, SOD. SOD has been found in
almost all organisms living in the presence of oxygen, including some anaerobic
bacteria, supporting the notion that superoxide is a key and general component
of oxidative stress. Recently, a new concept in the field of the mechanisms of
cellular defense against superoxide has emerged. It was discovered that
elimination of superoxide in some anaerobic and microaerophilic bacteria could
occur by reduction, a reaction catalyzed by a small metalloenzyme thus named
superoxide reductase, SOR. Having played a major role in this discovery, we
describe here how the concept of superoxide reduction emerged and how it was
experimentally substantiated independently in our laboratory.


  The rapid development of sequencing technologies represents new opportunities
for population genetics research. It is expected that genomic data will
increase our ability to reconstruct the history of populations. While this
increase in genetic information will likely help biologists and anthropologists
to reconstruct the demographic history of populations, it also represents new
challenges. Recent work has shown that structured populations generate signals
of population size change. As a consequence it is often difficult to determine
whether demographic events such as expansions or contractions (bottlenecks)
inferred from genetic data are real or due to the fact that populations are
structured in nature. Given that few inferential methods allow us to account
for that structure, and that genomic data will necessarily increase the
precision of parameter estimates, it is important to develop new approaches. In
the present study we analyse two demographic models. The first is a model of
instantaneous population size change whereas the second is the classical
symmetric island model. We (i) re-derive the distribution of coalescence times
under the two models for a sample of size two, (ii) use a maximum likelihood
approach to estimate the parameters of these models (iii) validate this
estimation procedure under a wide array of parameter combinations, (iv)
implement and validate a model choice procedure by using a Kolmogorov-Smirnov
test. Altogether we show that it is possible to estimate parameters under
several models and perform efficient model choice using genetic data from a
single diploid individual.


  The theory of natural selection has two forms. Deductive theory describes how
populations change over time. One starts with an initial population and some
rules for change. From those assumptions, one calculates the future state of
the population. Deductive theory predicts how populations adapt to
environmental challenge. Inductive theory describes the causes of change in
populations. One starts with a given amount of change. One then assigns
different parts of the total change to particular causes. Inductive theory
analyzes alternative causal models for how populations have adapted to
environmental challenge. This chapter emphasizes the inductive analysis of
cause.


  Homologous gene shuffling between DNA promotes genetic diversity and is an
important pathway for DNA repair. For this to occur, homologous genes need to
find and recognize each other. However, despite its central role in homologous
recombination, the mechanism of homology recognition is still an unsolved
puzzle. While specific proteins are known to play a role at later stages of
recombination, an initial coarse grained recognition step has been proposed.
This relies on the sequence dependence of the DNA structural parameters, such
as twist and rise, mediated by intermolecular interactions, in particular
electrostatic ones. In this proposed mechanism, sequences having the same base
pair text, or are homologous, have lower interaction energy than those
sequences with uncorrelated base pair texts; the difference termed the
recognition energy. Here, we probe how the recognition energy changes when one
DNA fragment slides past another, and consider, for the first time, homologous
sequences in antiparallel alignment. This dependence on sliding was termed the
recognition well. We find that there is recognition well for anti-parallel,
homologous DNA tracts, but only a very shallow one, so that their interaction
will differ little from the interaction between two nonhomologous tracts. This
fact may be utilized in single molecule experiments specially targeted to test
the theory. As well as this, we test previous theoretical approximations in
calculating the recognition well for parallel molecules against MC simulations,
and consider more rigorously the optimization of the orientations of the
fragments about their long axes. The more rigorous treatment affects the
recognition energy a little, when the molecules are considered rigid. However
when torsional flexibility of the DNA molecules is introduced, we find
excellent agreement between analytical approximation and simulation.


  Persistent homology captures the evolution of topological features of a model
as a parameter changes. The most commonly used summary statistics of persistent
homology are the barcode and the persistence diagram. Another summary
statistic, the persistence landscape, was recently introduced by Bubenik. It is
a functional summary, so it is easy to calculate sample means and variances,
and it is straightforward to construct various test statistics. Implementing a
permutation test we detect conformational changes between closed and open forms
of the maltose-binding protein, a large biomolecule consisting of 370 amino
acid residues. Furthermore, persistence landscapes can be applied to machine
learning methods. A hyperplane from a support vector machine shows the clear
separation between the closed and open proteins conformations. Moreover,
because our approach captures dynamical properties of the protein our results
may help in identifying residues susceptible to ligand binding; we show that
the majority of active site residues and allosteric pathway residues are
located in the vicinity of the most persistent loop in the corresponding
filtered Vietoris-Rips complex. This finding was not observed in the classical
anisotropic network model.


  Endocytosis underlies many cellular functions including signaling and
nutrient uptake. The endocytosed cargo gets redistributed across a dynamic
network of endosomes undergoing fusion and fission. Here, a theoretical
approach is reviewed which can explain how the microscopic properties of
endosome interactions cause the emergent macroscopic properties of cargo
trafficking in the endosomal network. Predictions by the theory have been
tested experimentally and include the inference of dependencies and parameter
values of the microscopic processes. This theory could also be used to infer
mechanisms of signal-trafficking crosstalk. It is applicable to in vivo systems
since fixed samples at few time points suffice as input data.


  Cancer cells co-cultured in vitro reveal unexpected differential growth rates
that classical exponential growth models cannot account for. Two
non-interacting cell lines were grown in the same culture, and counts of each
species were recorded at periodic times. The relative growth of population
ratios was found to depend on the initial proportion, in contradiction with the
traditional exponential growth model. The proposed explanation is the
variability of growth rates for clones inside the same cell line. This leads to
a log-quadratic growth model that provides both a theoretical explanation to
the phenomenon that was observed, and a better fit to our growth data.


  Although the structure, function, conformational dynamics, and controlled
thermodynamics of proteins are manifested by their corresponding amino acid
sequences, the natural rules for molecular design and their corresponding
interplay remain obscure. In this study, we focused on the role of internal
cavities of proteins in conformational dynamics. We investigated the
pressure-induced responses from the cavity-enlarged L99A mutant of T4 lysozyme,
using high-pressure NMR spectroscopy. The signal intensities of the methyl
groups in the 1H/13C HSQC spectra, particularly those around the enlarged
cavity, decreased with the increasing pressure, and disappeared at 200 MPa,
without the appearance of new resonances, thus indicating the presence of
heterogeneous conformations around the cavity within the ground state ensemble.
Above 200 MPa, the signal intensities of more than 20 methyl groups gradually
decreased with the increasing pressure, without the appearance of new
resonances. Interestingly, these residues closely matched those sensing a large
conformational change between the ground- and high-energy states, at
atmospheric pressure. 13C and 1H NMR line-shape simulations showed that the
pressure-induced loss in the peak intensity could be explained by the increase
in the high-energy state population. In this high-energy state, the aromatic
side chain of F114 gets flipped into the enlarged cavity. The accommodation of
the phenylalanine ring into the efficiently packed cavity may decrease the
partial molar volume of the high-energy state, relative to the ground state. We
suggest that the enlarged cavity is involved in the conformational transition
to high-energy states and in the volume fluctuation of the ground state.


  When the process underlying DNA substitutions varies across evolutionary
history, the standard Markov models underlying standard phylogenetic methods
are mathematically inconsistent. The most prominent example is the general time
reversible model (GTR) together with some, but not all, of its submodels. To
rectify this deficiency, Lie Markov models have been developed as the class of
models that are consistent in the face of a changing process of DNA
substitutions. Some well-known models in popular use are within this class, but
are either overly simplistic (e.g. the Kimura two-parameter model) or overly
complex (the general Markov model). On a diverse set of biological data sets,
we test a hierarchy of Lie Markov models spanning the full range of parameter
richness. Compared against the benchmark of the ever-popular GTR model, we find
that as a whole the Lie Markov models perform remarkably well, with the best
performing models having eight parameters and the ability to recognise the
distinction between purines and pyrimidines.


  This paper investigates the dynamics of Ebola virus transmission in West
Africa during 2014. The reproduction numbers for the total period of epidemic
and for different consequent time intervals are estimated based on a newly
suggested linear model. It contains one major variable - the average time of
infectiousness (time from onset to hospitalization) that is considered as a
parameter for controlling the future dynamics of epidemics.
  Numerical implementations are carried out on data collected from three
countries Guinea, Sierra Leone and Liberia as well as the total data collected
worldwide. Predictions are provided by considering different scenarios
involving the average times of infectiousness for the next few months and the
end of the current epidemic is estimated according to each scenario.


  We study metanetworks arising in genotype and phenotype spaces, in the
context of a model population of Boolean graphs evolved under selection for
short dynamical attractors. We define the adjacency matrix of a graph as its
genotype, which gets mutated in the course of evolution, while its phenotype is
its set of dynamical attractors. Metanetworks in the genotype and phenotype
spaces are formed, respectively, by genetic proximity and by phenotypic
similarity, the latter weighted by the sizes of the basins of attraction of the
shared attractors. We find that populations of evolved networks form giant
clusters in genotype space, have Poissonian degree distributions but exhibit
hierarchically organized $\kappa$-core decompositions. Nevertheless, at large
scales, they form tree-like expander graphs. Random populations of Boolean
graphs are typically so far removed from each other genetically that they
cannot form a metanetwork. In phenotype space, the metanetworks of evolved
populations are super robust both under the elimination of weak connections and
random removal of nodes.


  Relatedness and synergy affect the selection pressure on cooperation and
altruism. Although early work investigated the effect of these factors
independently of each other, recent efforts have been aimed at exploring their
interplay. Here, we contribute to this ongoing synthesis in two distinct but
complementary ways. First, we integrate models of $n$-player matrix games into
the direct fitness approach of inclusive fitness theory, hence providing a
framework to consider synergistic social interactions between relatives in
family and spatially structured populations. Second, we illustrate the
usefulness of this framework by delineating three distinct types of helping
traits ("whole-group", "nonexpresser-only" and "expresser-only"), which are
characterized by different synergies of kind (arising from differential fitness
effects on individuals expressing or not expressing helping) and can be
subjected to different synergies of scale (arising from economies or
diseconomies of scale). We find that relatedness and synergies of kind and
scale can interact to generate nontrivial evolutionary dynamics, such as cases
of bistable coexistence featuring both a stable equilibrium with a positive
level of helping and an unstable helping threshold. This broadens the
qualitative effects of relatedness (or spatial structure) on the evolution of
helping.


  The newly developed deep-sequencing technologies make it possible to acquire
both quantitative and qualitative information regarding transcript biology. By
measuring messenger RNA levels for all genes in a sample, RNA-seq provides an
attractive option to characterize the global changes in transcription. RNA-seq
is becoming the widely used platform for gene expression profiling. However,
real transcription signals in the RNA-seq data are confounded with measurement
and sequencing errors, and other random biological/technical variation. How to
appropriately take the variability due to errors and sequencing technology
variation into account is essential issue in the RNA-seq data analysis. To
extract biologically useful transcription process from the RNA-seq data, we
propose to use the second ODE for modeling the RNA-seq data. We use
differential principal analysis to develop statistical methods for estimation
of location-varying coefficients of the ODE. We validate the accuracy of the
ODE model to fit the RNA-seq data by prediction analysis and 5 fold cross
validation. We find the accuracy of the second ODE model to predict the gene
expression level across the gene is very high and the second ODE model to fit
the RNA-seq data very well. To further evaluate the performance of the ODE
model for RNA-seq data analysis, we used the location-varying coefficients of
the second ODE as features to classify the normal and tumor cells. We
demonstrate that even using the ODE model for single gene we can achieve high
classification accuracy. We also conduct response analysis to investigate how
the transcription process respond to the perturbation of the external signals
and identify dozens of genes that are related to cancer.


  All living things exhibit adaptations that enable them to survive and
reproduce in the natural environment that they inhabit. From a biological
standpoint, it has long been understood that adaptation comes from natural
selection, whereby maladapted individuals do not pass their traits effectively
to future generations. However, we may also consider the phenomenon of
adaptation from the standpoint of physics, and ask whether it is possible to
delineate what the difference is in terms of physical properties between
something that is well-adapted to its surrounding environment, and something
that is not. In this work, we undertake to address this question from a
theoretical standpoint. Building on past fundamental results in
far-from-equilibrium statistical mechanics, we demonstrate a generalization of
the Helmholtz free energy for the finite-time stochastic evolution of driven
Newtonian matter. By analyzing this expression term by term, we are able to
argue for a general tendency in driven many-particle systems towards
self-organization into states formed through exceptionally reliable absorption
and dissipation of work energy from the surrounding environment. Subsequently,
we illustrate the mechanism of this general tendency towards physical
adaptation by analyzing the process of random hopping in driven energy
landscapes.


  We present a data based statistical study on the effects of seasonal
variations in the growth rates of the gastro-intestinal (GI) parasitic
infection in livestock. The alluded growth rate is estimated through the
variation in the number of eggs per gram (EPG) of faeces in animals. In
accordance with earlier studies, our analysis too shows that rainfall is the
dominant variable in determining EPG infection rates compared to other
macro-parameters like temperature and humidity. Our statistical analysis
clearly indicates an oscillatory dependence of EPG levels on rainfall
fluctuations. Monsoon recorded the highest infection with a comparative
increase of at least 2.5 times compared to the next most infected period
(summer). A least square fit of the EPG versus rainfall data indicates an
approach towards a super diffusive (i. e. root mean square displacement growing
faster than the square root of the elapsed time as obtained for simple
diffusion) infection growth pattern regime for low rainfall regimes
(technically defined as zeroth level dependence) that gets remarkably augmented
for large rainfall zones. Our analysis further indicates that for low
fluctuations in temperature (true on the bulk data), EPG level saturates beyond
a critical value of the rainfall, a threshold that is expected to indicate the
onset of the nonlinear regime. The probability density functions (PDFs) of the
EPG data show oscillatory behavior in the large rainfall regime (greater than
500 mm), the frequency of oscillation, once again, being determined by the
ambient wetness (rainfall, and humidity). Data recorded over three pilot
projects spanning three measures of rainfall and humidity bear testimony to the
universality of this statistical argument.


  In this paper, a two-step strategy for parametric sensitivity analysis for
such systems is proposed, exploiting advantages and synergies between two
recently proposed sensitivity analysis methodologies for stochastic dynamics.
The first method performs sensitivity analysis of the stochastic dynamics by
means of the Fisher Information Matrix on the underlying distribution of the
trajectories; the second method is a reduced-variance, finite-difference,
gradient-type sensitivity approach relying on stochastic coupling techniques
for variance reduction. Here we demonstrate that these two methods can be
combined and deployed together by means of a new sensitivity bound which
incorporates the variance of the quantity of interest as well as the Fisher
Information Matrix estimated from the first method. The first step of the
proposed strategy labels sensitivities using the bound and screens out the
insensitive parameters in a controlled manner based also on the new sensitivity
bound. In the second step of the proposed strategy, the finite-difference
method is applied only for the sensitivity estimation of the (potentially)
sensitive parameters that have not been screened out in the first step. Results
on an epidermal growth factor network with fifty parameters and on a protein
homeostasis with eighty parameters demonstrate that the proposed strategy is
able to quickly discover and discard the insensitive parameters and in the
remaining potentially sensitive parameters it accurately estimates the
sensitivities. The new sensitivity strategy can be several times faster than
current state-of-the-art approaches that test all parameters, especially in
"sloppy" systems. In particular, the computational acceleration is quantified
by the ratio between the total number of parameters over the number of the
sensitive parameters.


  We examined protein residue networks (PRNs) from a local search perspective
to understand why PRNs are highly clustered when having short paths is
important for protein functionality. We found that by adopting a local search
perspective, this conflict between form and function is resolved as increased
clustering actually helps to reduce path length in PRNs. Further, the paths
found via our EDS local search algorithm are more congruent with the
characteristics of intra-protein communication. EDS identifies a subset of PRN
edges called short-cuts that are distinct, have high usage, impacts EDS path
length, diversity and stretch, and are dominated by short-range contacts. The
short-cuts form a network (SCN) that increases in size and transitivity as a
protein folds. The structure of a SCN supports its function and formation, and
the function of a SCN influences its formation. Several significant differences
in terms of SCN structure, function and formation is found between successful
and unsuccessful MD trajectories. By connecting the static and the dynamic
aspects of PRNs, the protein folding process becomes a problem of graph
formation with the purpose of forming suitable pathways within proteins.


  It is well known that cooperation cannot be an evolutionary stable strategy
for a non-iterative game in a well-mixed population. In contrast, structured
populations favor cooperation since cooperators can benefit each other by
forming local clusters. Previous studies have shown that scale-free networks
strongly promote cooperation. However, little is known about the invasion
mechanism of cooperation in scale-free networks. To study microscopic and
macroscopic behaviors of cooperators' invasion, we conducted computational
experiments of the evolution of cooperation in scale-free networks where,
starting from all defectors, cooperators can spontaneously emerge by mutation.
Since the evolutionary dynamics are influenced by the definition of fitness, we
tested two commonly adopted fitness functions: accumulated payoff and average
payoff. Simulation results show that cooperation is strongly enhanced with the
accumulated payoff fitness compared to the average payoff fitness. However, the
difference between the two functions decreases as the average degree increases.
As the average degree increases, cooperation decreases with the accumulated
payoff fitness, while it increases with the average payoff fitness. Moreover,
with the average payoff fitness, low-degree nodes play a more important role in
spreading cooperative strategies compared to the case of the accumulated payoff
fitness.


  Persistent homology provides a new approach for the topological
simplification of big data via measuring the life time of intrinsic topological
features in a filtration process and has found its success in scientific and
engineering applications. However, such a success is essentially limited to
qualitative data characterization, identification and analysis (CIA). In this
work, we outline a general protocol to construct objective-oriented persistent
homology methods. The minimization of the objective functional leads to a
Laplace-Beltrami operator which generates a multiscale representation of the
initial data and offers an objective oriented filtration process. The resulting
differential geometry based objective-oriented persistent homology is able to
preserve desirable geometric features in the evolutionary filtration and
enhances the corresponding topological persistence. The consistence between
Laplace-Beltrami flow based filtration and Euclidean distance based filtration
is confirmed on the Vietoris-Rips complex for a large amount of numerical
tests. The convergence and reliability of the present Laplace-Beltrami flow
based cubical complex filtration approach are analyzed over various spatial and
temporal mesh sizes. The efficiency and robustness of the present method are
verified by more than 500 fullerene molecules. It is shown that the proposed
persistent homology based quantitative model offers good predictions of total
curvature energies for ten types of fullerene isomers. The present work offers
the first example to design objective-oriented persistent homology to enhance
or preserve desirable features in the original data during the filtration
process and then automatically detect or extract the corresponding topological
traits from the data.


  Persistent homology is a relatively new tool often used for
\emph{qualitative} analysis of intrinsic topological features in images and
data originated from scientific and engineering applications. In this paper, we
report novel \emph{quantitative} predictions of the energy and stability of
fullerene molecules, the very first attempt in employing persistent homology in
this context. The ground-state structures of a series of small fullerene
molecules are first investigated with the standard Vietoris-Rips complex. We
decipher all the barcodes, including both short-lived local bars and long-lived
global bars arising from topological invariants, and associate them with
fullerene structural details. By using accumulated bar lengths, we build
quantitative models to correlate local and global Betti-2 bars respectively
with the heat of formation and total curvature energies of fullerenes. It is
found that the heat of formation energy is related to the local hexagonal
cavities of small fullerenes, while the total curvature energies of fullerene
isomers are associated with their sphericities, which are measured by the
lengths of their long-lived Betti-2 bars. Excellent correlation coefficients
($>0.94$) between persistent homology predictions and those of quantum or
curvature analysis have been observed. A correlation matrix based filtration is
introduced to further verify our findings.


  In this work, we introduce persistent homology for the analysis of
cryo-electron microscopy (cryo-EM) density maps. We identify the topological
fingerprint or topological signature of noise, which is widespread in cryo-EM
data. For low signal to noise ratio (SNR) volumetric data, intrinsic
topological features of biomolecular structures are indistinguishable from
noise. To remove noise, we employ geometric flows which are found to preserve
the intrinsic topological fingerprints of cryo-EM structures and diminish the
topological signature of noise. In particular, persistent homology enables us
to visualize the gradual separation of the topological fingerprints of cryo-EM
structures from those of noise during the denoising process, which gives rise
to a practical procedure for prescribing a noise threshold to extract cryo-EM
structure information from noise contaminated data after certain iterations of
the geometric flow equation. To further demonstrate the utility of persistent
homology for cryo-EM data analysis, we consider a microtubule intermediate
structure (EMD-1129). Three helix models, an alpha-tubulin monomer model, an
alpha- and beta-tubulin model, and an alpha- and beta-tubulin dimer model, are
constructed to fit the cryo-EM data. The least square fitting leads to
similarly high correlation coefficients, which indicates that structure
determination via optimization is an ill-posed inverse problem. By a comparison
of the topological fingerprints of the original data and those of three models,
we found that the third model is topologically favored. The present work offers
persistent homology based new strategies for topological denoising and for
resolving ill-posed inverse problems.


  Many proteins and peptides have an intrinsic capacity to sense and induce
membrane curvature, and play crucial roles for organizing and remodelling cell
membranes. However, the molecular driving forces behind these processes are not
well understood. Here, we describe a new approach to study curvature sensing,
by simulating the direction-dependent interactions of single molecules with a
buckled lipid bilayer. We analyse three amphipathic antimicrobial peptides, a
class of membrane-associated molecules that specifically target and destabilize
bacterial membranes, and find qualitatively different sensing characteristics
that would be difficult to resolve with other methods. These findings provide
new insights into the curvature sensing mechanisms of amphipathic peptides and
challenge existing theories of hydrophobic insertion. Our approach is generally
applicable to a wide range of curvature sensing molecules, and our results
provide strong motivation to develop new experimental methods to track position
and orientation of membrane proteins.


  Pre-mRNA splicing relies on the poorly understood dynamic interplay between
>150 protein components of the spliceosome. The steps at which splicing can be
regulated remain largely unknown. We systematically analyzed the effect of
knocking down the components of the splicing machinery on alternative splicing
events relevant for cell proliferation and apoptosis and used this information
to reconstruct a network of functional interactions. The network accurately
captures known physical and functional associations and identifies new ones,
revealing remarkable regulatory potential of core spliceosomal components,
related to the order and duration of their recruitment during spliceosome
assembly. In contrast with standard models of regulation at early steps of
splice site recognition, factors involved in catalytic activation of the
spliceosome display regulatory properties. The network also sheds light on the
antagonism between hnRNP C and U2AF, and on targets of antitumor drugs, and can
be widely used to identify mechanisms of splicing regulation.


  Motivation: Bisulphite sequencing enables the detection of cytosine
methylation. The sequence of the methylation states of cytosines on any given
read forms a methylation pattern that carries substantially more information
than merely studying the average methylation level at individual positions. In
order to understand better the complexity of DNA methylation landscapes in
biological samples, it is important to study the diversity of these methylation
patterns. However, the accurate quantification of methylation patterns is
subject to sequencing errors and spurious signals due to incomplete bisulphite
conversion of cytosines. Results: A statistical model is developed which
accounts for the distribution of DNA methylation patterns at any given locus.
The model incorporates the effects of sequencing errors and spurious reads, and
enables estimation of the true underlying distribution of methylation patterns.
Conclusions: Calculation of the estimated distribution over methylation
patterns is implemented in the R Bioconductor package MPFE. Source code and
documentation of the package are also available for download at
http://bioconductor.org/packages/3.0/bioc/html/MPFE.html.


  We consider biological individuality in terms of information theoretic and
graphical principles. Our purpose is to extract through an algorithmic
decomposition system-environment boundaries supporting individuality. We infer
or detect evolved individuals rather than assume that they exist. Given a set
of consistent measurements over time, we discover a coarse-grained or quantized
description on a system, inducing partitions (which can be nested). Legitimate
individual partitions will propagate information from the past into the future,
whereas spurious aggregations will not. Individuals are therefore defined in
terms of ongoing, bounded information processing units rather than lists of
static features or conventional replication-based definitions which tend to
fail in the case of cultural change. One virtue of this approach is that it
could expand the scope of what we consider adaptive or biological phenomena,
particularly in the microscopic and macroscopic regimes of molecular and social
phenomena.


  We develop an theoretical approach for predicting biodiversity in
multi-dimensional niche spaces, arising due to ecological drivers such as
competitive exclusion. The novelty of our approach relies on the fact that
ecological niches are described by sequences of strings, which allows us to
describe multiple traits. We define the mathematical framework for analyzing
pattern forming instabilities in these models, showing surprisingly that the
analytic linear theory predicts the asymptotically long time population
distributions of niches in the model. We propose a test for identifying
ecological drivers in biodiversity distributions, based on representing
ecosystem data by means of a certain transform introduced in the theory.


  Superoxide radical (O2.-) is the univalent reduction product of molecular
oxygen and belongs to the group of the so-called toxic oxygen derivatives. For
years the only enzymatic system known to catalyze the elimination of superoxide
was the superoxide dismutase (SOD), which catalyzes dismutation of superoxide
radical anions to hydrogen peroxide and molecular oxygen


  Proteins are the most important biomolecules for living organisms. The
understanding of protein structure, function, dynamics and transport is one of
most challenging tasks in biological science. In the present work, persistent
homology is, for the first time, introduced for extracting molecular
topological fingerprints (MTFs) based on the persistence of molecular
topological invariants. MTFs are utilized for protein characterization,
identification and classification. The method of slicing is proposed to track
the geometric origin of protein topological invariants. Both all-atom and
coarse-grained representations of MTFs are constructed. A new cutoff-like
filtration is proposed to shed light on the optimal cutoff distance in elastic
network models. Based on the correlation between protein compactness, rigidity
and connectivity, we propose an accumulated bar length generated from
persistent topological invariants for the quantitative modeling of protein
flexibility. To this end, a correlation matrix based filtration is developed.
This approach gives rise to an accurate prediction of the optimal
characteristic distance used in protein B-factor analysis. Finally, MTFs are
employed to characterize protein topological evolution during protein folding
and quantitatively predict the protein folding stability. An excellent
consistence between our persistent homology prediction and molecular dynamics
simulation is found. This work reveals the topology-function relationship of
proteins.