  Aptamers are single stranded DNA, RNA or peptide sequences having the ability
to bind a variety of specific targets (proteins, molecules as well as ions).
Therefore, aptamer production and selection for therapeutic and diagnostic
applications is very challenging. Usually they are in vitro generated, but,
recently, computational approaches have been developed for the in silico
selection, with a higher affinity for the specific target. Anyway, the
mechanism of aptamer-ligand formation is not completely clear, and not obvious
to predict. This paper aims to develop a computational model able to describe
aptamer-ligand affinity performance by using the topological structure of the
corresponding graphs, assessed by means of numerical tools such as the
conventional degree distribution, but also the rank-degree distribution
(hierarchy) and the node assortativity. Calculations are applied to the
thrombin binding aptamer (TBA), and the TBA-thrombin complex, produced in the
presence of Na+ or K+. The topological analysis reveals different affinity
performances between the macromolecules in the presence of the two cations, as
expected by previous investigations in literature. These results nominate the
graph topological analysis as a novel theoretical tool for testing affinity.
Otherwise, starting from the graphs, an electrical network can be obtained by
using the specific electrical properties of amino acids and nucleobases.
Therefore, a further analysis concerns with the electrical response, which
reveals that the resistance sensitively depends on the presence of sodium or
potassium thus posing resistance as a crucial physical parameter for testing
affinity.


  Biological sensory systems generally operate out of equilibrium, which often
leads to their improved performance. Here, we study the sensitivity of ligand
concentration for a general receptor model, which is generally in the
non-equilibrium stationary state, in the framework of a stochastic diffusion
equation. We derived a general formula of the maximum sensitivity.
Specifically, the sensitivity is limited universally by the Berg-Purcell limit
[Biophys. J ., 1977], regardless of whether the receptor is in an equilibrium
or non-equilibrium state.


  Stochastic chemical reaction networks (CRNs) are complex systems which
combine the features of concurrent transformation of multiple variables in each
elementary reaction event, and nonlinear relations between states and their
rates of change. Most general results concerning CRNs are limited to restricted
cases where a topological characteristic known as deficiency takes value 0 or
1. Here we derive equations of motion for fluctuation moments at all orders for
stochastic CRNs at general deficiency. We show, for the case of the mass-action
rate law, that the generator of the stochastic process acts on the hierarchy of
factorial moments with a finite representation. Whereas simulation of
high-order moments for many-particle systems is costly, this representation
reduces solution of moment hierarchies to a complexity comparable to solving a
heat equation. At steady states, moment hierarchies for finite CRNs interpolate
between low-order and high-order scaling regimes, which may be approximated
separately by distributions similar to those for deficiency-0 networks, and
connected through matched asymptotic expansions. In CRNs with multiple stable
or metastable steady states, boundedness of high-order moments provides the
starting condition for recursive solution downward to low-order moments,
reversing the order usually used to solve moment hierarchies. A basis for a
subset of network flows defined by having the same mean-regressing property as
the flows in deficiency-0 networks gives the leading contribution to low-order
moments in CRNs at general deficiency, in a $1/n$-expansion in large particle
numbers. Our results give a physical picture of the different informational
roles of mean-regressing and non-mean-regressing flows, and clarify the
dynamical meaning of deficiency not only for first-moment conditions but for
all orders in fluctuations.


  In ecology it is widely recognised that many landscapes comprise a network of
discrete patches of habitat. The species that inhabit the patches interact with
each other through a foodweb, the network of feeding interactions. The
meta-foodweb model proposed by Pillai et al. combines the feeding relationships
at each patch with the dispersal of species between patches, such that the
whole system is represented by a network of networks. Previous work on
meta-foodwebs has focussed on landscape networks that do not have an explicit
spatial embedding, but in real landscapes the patches are usually distributed
in space. Here we compare the dispersal of a meta-foodweb on \ER networks, that
do not have a spatial embedding, and random geometric networks, that do have a
spatial embedding. We found that local structure and large network distances in
spatially embedded networks, lead to meso-scale patterns of patch occupation by
both specialist and omnivorous species. In particular, we found that spatial
separations make the coexistence of competing species more likely. Our results
highlight the effects of spatial embeddings for meta-foodweb models, and the
need for new analytical approaches to them.


  Genomics has revolutionized biology, enabling the interrogation of whole
transcriptomes, genome-wide binding sites for proteins, and many other
molecular processes. However, individual genomic assays measure elements that
interact in vivo as components of larger molecular machines. Understanding how
these high-order interactions drive gene expression presents a substantial
statistical challenge. Building on Random Forests (RF), Random Intersection
Trees (RITs), and through extensive, biologically inspired simulations, we
developed the iterative Random Forest algorithm (iRF). iRF trains a
feature-weighted ensemble of decision trees to detect stable, high-order
interactions with same order of computational cost as RF. We demonstrate the
utility of iRF for high-order interaction discovery in two prediction problems:
enhancer activity in the early Drosophila embryo and alternative splicing of
primary transcripts in human derived cell lines. In Drosophila, among the 20
pairwise transcription factor interactions iRF identifies as stable (returned
in more than half of bootstrap replicates), 80% have been previously reported
as physical interactions. Moreover, novel third-order interactions, e.g.
between Zelda (Zld), Giant (Gt), and Twist (Twi), suggest high-order
relationships that are candidates for follow-up experiments. In human-derived
cells, iRF re-discovered a central role of H3K36me3 in chromatin-mediated
splicing regulation, and identified novel 5th and 6th order interactions,
indicative of multi-valent nucleosomes with specific roles in splicing
regulation. By decoupling the order of interactions from the computational cost
of identification, iRF opens new avenues of inquiry into the molecular
mechanisms underlying genome biology.


  We introduce a minimal model for the evolution of functional
protein-interaction networks using a sequence-based mutational algorithm, and
apply the model to study neutral drift in networks that yield oscillatory
dynamics. Starting with a functional core module, random evolutionary drift
increases network complexity even in the absence of specific selective
pressures. Surprisingly, we uncover a hidden order in sequence space that gives
rise to long-term evolutionary memory, implying strong constraints on network
evolution due to the topology of accessible sequence space.


  The number of complete chloroplastic genomes increases day after day, making
it possible to rethink plants phylogeny at the biomolecular era. Given a set of
close plants sharing in the order of one hundred of core chloroplastic genes,
this article focuses on how to extract the largest subset of sequences in order
to obtain the most supported species tree. Due to computational complexity, a
discrete and distributed Particle Swarm Optimization (DPSO) is proposed. It is
finally applied to the core genes of Rosales order.


  In this paper, we study the dynamics of contagious spreading processes taking
place in complex contact networks. We specifically present a lower-bound on the
decay rate of the number of nodes infected by a
susceptible-infected-susceptible (SIS) stochastic spreading process. A precise
quantification of this decay rate is crucial for designing efficient strategies
to contain epidemic outbreaks. However, existing lower-bounds on the decay rate
based on first-order mean-field approximations are often accompanied by a large
error resulting in inefficient containment strategies. To overcome this
deficiency, we derive a lower-bound based on a second-order moment-closure of
the stochastic SIS processes. The proposed second-order bound is theoretically
guaranteed to be tighter than existing first-order bounds. We also present
various numerical simulations to illustrate how our lower-bound drastically
improves the performance of existing first-order lower-bounds in practical
scenarios, resulting in more efficient strategies for epidemic containment.


  Functional responses are widely used to describe interactions and resources
exchange between individuals in ecology. The form given to functional responses
dramatically affects the dynamics and stability of populations and communities.
Despite their importance, functional responses are generally considered with a
phenomenological approach, without clear mechanistic justifications from
individual traits and behaviors. Here, we develop a bottom-up stochastic
framework grounded in Renewal Theory showing how functional responses emerge
from the level of the individuals through the decomposition of interactions
into different activities. Our framework has many applications for conceptual,
theoretical and empirical purposes. First, we show how the mean and variance of
classical functional responses are obtained with explicit ecological
assumptions, for instance regarding foraging behaviors. Second, we give
examples in specific ecological contexts, such as in nuptial-feeding species or
size dependent handling times. Finally, we demonstrate how to analyze data with
our framework, especially highlighting that observed variability in the number
of interactions can be used to infer parameters and compare functional response
models.


  Photosynthetic materials form the basis of quantum biology. An important
attribute of quantum biology is correlation and coherence of spin states. Such
correlated spin states are targets of static magnetic field. In this paper, we
report magnetic properties and spectroscopically realizable static magnetic
field effect in photosynthetic materials. Two classes of nano-scale assembly of
chlorophyll (NC) are used for such a study. Magnetic measurements are made
using a superconducting quantum interference device (SQUID). Both ferromagnetic
and superparamagnetic states are observed in NC along with a blocking
temperature around 250 K. Low temperature quantum (liquid nitrogen)
spectroscopy is employed to see how optical transitions are affected in
presence of static magnetic field. Plausible practical application aspects of
magnetic properties of this optically active material are discussed in the
text.


  The two-step contagion model is a simple toy model for understanding pandemic
outbreaks that occur in the real world. The model takes into account that a
susceptible person either gets immediately infected or weakened when getting
into contact with an infectious one. As the number of weakened people
increases, they eventually can become infected in a short time period and a
pandemic outbreak occurs. The time required to reach such a pandemic outbreak
allows for intervention and is often called golden time. Understanding the
size-dependence of the golden time is useful for controlling pandemic outbreak.
Here we find that there exist two types of golden times in the two-step
contagion model, which scale as $O(N^{1/3})$ and $O(N^{\zeta})$ with the system
size $N$ on Erd\H{o}s-R\'enyi networks, where the measured $\zeta$ is slightly
larger than $1/4$. They are distinguished by the initial number of infected
nodes, $o(N)$ and $O(N)$, respectively. While the exponent $1/3$ of the
$N$-dependence of the golden time is universal even in other models showing
discontinuous transitions induced by cascading dynamics, the measured $\zeta$
exponents are all close to $1/4$ but show model-dependence. It remains open
whether or not $\zeta$ reduces to $1/4$ in the asymptotically large-$N$ limit.


  Immunogenicity is a major problem during the development of biotherapeutics
since it can lead to rapid clearance of the drug and adverse reactions. The
challenge for biotherapeutic design is therefore to identify mutants of the
protein sequence that minimize immunogenicity in a target population whilst
retaining pharmaceutical activity and protein function. Current approaches are
moderately successful in designing sequences with reduced immunogenicity, but
do not account for the varying frequencies of different human leucocyte antigen
alleles in a specific population and in addition, since many designs are
non-functional, require costly experimental post-screening. Here we report a
new method for de-immunization design using multi-objective combinatorial
optimization that simultaneously optimizes the likelihood of a functional
protein sequence at the same time as minimizing its immunogenicity tailored to
a target population. We bypass the need for three-dimensional protein structure
or molecular simulations to identify functional designs by automatically
generating sequences using probabilistic models that have been used previously
for mutation effect prediction and structure prediction. As proof-of-principle
we designed sequences of the C2 domain of Factor VIII and tested them
experimentally, resulting in a good correlation with the predicted
immunogenicity of our model.


  Disease, opinions, ideas, gossip, etc. all spread on social networks. How
these networks are connected (the network structure) influences the dynamics of
the spreading processes. By investigating these relationships one gains
understanding both of the spreading itself and the structure and function of
the contact network. In this chapter, we will summarize the recent literature
using simulation of spreading processes on top of empirical contact data. We
will mostly focus on disease simulations on temporal proximity networks --
networks recording who is close to whom, at what time -- but also cover other
types of networks and spreading processes. We analyze 29 empirical networks to
illustrate the methods.


  A birth-death-sampling model gives rise to phylogenetic trees with samples
from the past and the present. Interpreting "birth" as branching speciation,
"death" as extinction, and "sampling" as fossil preservation and recovery, this
model -- also referred to as the fossilized birth-death (FBD) model -- gives
rise to phylogenetic trees on extant and fossil samples. The model has been
mathematically analyzed and successfully applied to a range of datasets on
different taxonomic levels, such as penguins, plants, and insects. However, the
current mathematical treatment of this model does not allow for a group of
temporally distinct fossil specimens to be assigned to the same species. In
this paper, we provide a general mathematical FBD modeling framework that
explicitly takes "stratigraphic ranges" into account, with a stratigraphic
range being defined as the lineage interval associated with a single species,
ranging through time from the first to the last fossil appearance of the
species. To assign a sequence of fossil samples in the phylogenetic tree to the
same species, i.e., to specify a stratigraphic range, we need to define the
mode of speciation. We provide expressions to account for three common
speciation modes: budding (or asymmetric) speciation, bifurcating (or
symmetric) speciation, and anagenetic speciation. Our equations allow for
flexible joint Bayesian analysis of paleontological and neontological data.
Furthermore, our framework is directly applicable to epidemiology, where a
stratigraphic range is the observed duration of infection of a single patient,
"birth" via budding is transmission, "death" is recovery, and "sampling" is
sequencing the pathogen of a patient. Thus, we present a model that allows for
incorporation of multiple observations through time from a single patient.


  DNA methylation is an epigenetic mechanism whose important role in
development has been widely recognized. This epigenetic modification results in
heritable changes in gene expression not encoded by the DNA sequence. The
underlying mechanisms controlling DNA methylation are only partly understood
and recently different mechanistic models of enzyme activities responsible for
DNA methylation have been proposed. Here we extend existing Hidden Markov
Models (HMMs) for DNA methylation by describing the occurrence of spatial
methylation patterns over time and propose several models with different
neighborhood dependencies. We perform numerical analysis of the HMMs applied to
bisulfite sequencing measurements and accurately predict wild-type data. In
addition, we find evidence that the enzymes' activities depend on the left 5'
neighborhood but not on the right 3' neighborhood.


  The cellular cytoskeleton is an active material, driven out of equilibrium by
molecular motor proteins. It is not understood how the collective behaviors of
cytoskeletal networks emerge from the properties of the network's constituent
motor proteins and filaments. Here we present experimental results on networks
of stabilized microtubules in Xenopus oocyte extracts, which undergo
spontaneous bulk contraction driven by the motor protein dynein, and
investigate the effects of varying the initial microtubule density and length
distribution. We find that networks contract to a similar final density,
irrespective of the length of microtubules or their initial density, but that
the contraction timescale varies with the average microtubule length. To gain
insight into why this microscopic property influences the macroscopic network
contraction time, we developed simulations where microtubules and motors are
explicitly represented. The simulations qualitatively recapitulate the
variation of contraction timescale with microtubule length, and allowed stress
contributions from different sources to be estimated and decoupled.


  A new approach to computing the equilibria and steady-states of biomolecular
systems modelled by bond graphs is presented. The approach is illustrated using
a model of a biomolecular cycle representing a membrane transporter and a model
of the mitochondrial electron transport chain.


  Game theory has been widely applied to many areas including economics,
biology and social sciences. However, it is still challenging to quantify the
global stability and global dynamics of the game theory. We developed a
landscape and flux framework to quantify the global stability and global
dynamics of the game theory. As an example, we investigated the models of
three-strategy games: a special replicator-mutator game, the repeated prison
dilemma model. In this model, one stable state, two stable states and limit
cycle can emerge under different parameters. The repeated Prisoner's Dilemma
system has Hopf bifurcation transitions from one stable state to limit cycle
state, and then to another one stable state or two stable states, or vice
versa. We explored the global stability of the repeated Prisoner's Dilemma
system and the kinetic paths between the basins of attractor. The paths are
irreversible due to the non-zero flux. One can explain the game for $Peace$ and
$War$.


  We extend a technique of approximation of the long-term behavior of a
supercritical stochastic epidemic model, using the WKB approximation and a
Hamiltonian phase space, to the subcritical case. The limiting behavior of the
model and approximation are qualitatively different in the subcritical case,
requiring a novel analysis of the limiting behavior of the Hamiltonian system
away from its deterministic subsystem. This yields a novel, general technique
of approximation of the quasistationary distribution of stochastic epidemic and
birth-death models, and may lead to techniques for analysis of these models
beyond the quasistationary distribution. For a classic SIS model, the
approximation found for the quasistationary distribution is very similar to
published approximations but not identical. For a birth-death process without
depletion of susceptibles, the approximation is exact. Dynamics on the phase
plane similar to those predicted by the Hamiltonian analysis are demonstrated
in cross-sectional data from trachoma treatment trials in Ethiopia, in which
declining prevalences are consistent with subcritical epidemic dynamics.


  The Adaptive Poisson-Boltzmann Solver (APBS) software was developed to solve
the equations of continuum electrostatics for large biomolecular assemblages
that has provided impact in the study of a broad range of chemical, biological,
and biomedical applications. APBS addresses three key technology challenges for
understanding solvation and electrostatics in biomedical applications: accurate
and efficient models for biomolecular solvation and electrostatics, robust and
scalable software for applying those theories to biomolecular systems, and
mechanisms for sharing and analyzing biomolecular electrostatics data in the
scientific community. To address new research applications and advancing
computational capabilities, we have continually updated APBS and its suite of
accompanying software since its release in 2001. In this manuscript, we discuss
the models and capabilities that have recently been implemented within the APBS
software package including: a Poisson-Boltzmann analytical and a
semi-analytical solver, an optimized boundary element solver, a geometry-based
geometric flow solvation model, a graph theory based algorithm for determining
p$K_a$ values, and an improved web-based visualization tool for viewing
electrostatics.


  It is generally difficult to predict the positions of mutations in genomic
DNA at the nucleotide level. Retroviral DNA insertion is one mode of mutation,
resulting in host infections that are difficult to treat. This mutation process
involves the integration of retroviral DNA into the host-infected cellular
genomic DNA following the interaction between host DNA and a pre-integration
complex consisting of retroviral DNA and integrase. Here, we report that
retroviral insertion sites around a hotspot within the Zfp521 and N-myc genes
can be predicted by a periodic function that is deduced using the diffraction
lattice model. In conclusion, the mutagenesis process is described by a
biophysical model for DNA-DNA interactions.


  In this work maximum entropy distributions in the space of steady states of
metabolic networks are defined upon constraining the first and second moment of
the growth rate. Inherent bistability of fast and slow phenotypes, akin to a
Van-Der Waals picture, emerges upon considering control on the average growth
(optimization/repression) and its fluctuations (heterogeneity). This is applied
to the carbon catabolic core of E.coli where it agrees with some stylized facts
on the persisters phenotype and it provides a quantitative map with metabolic
fluxes, opening for the possibility to detect coexistence from flux data.
Preliminary analysis on data for E.Coli cultures in standard conditions shows,
on the other hand, degeneracy for the inferred parameters that extend in the
coexistence region.


  In evolutionary biology, it is common to study how various entities evolve
together, for example, how parasites coevolve with their host, or genes with
their species. Coevolution is commonly modelled by considering certain maps or
reconciliations from one evolutionary tree $P$ to another $H$, all of which
induce the same map $\phi$ between the leaf-sets of $P$ and $H$ (corresponding
to present-day associations). Recently, there has been much interest in
studying spaces of reconciliations, which arise by defining some metric $d$ on
the set $Rec(P,H,\phi)$ of all possible reconciliations between $P$ and $H$.
  In this paper, we study the following question: How do we compute a geometric
median for a given subset $\Psi$ of $Rec(P,H,\phi)$ relative to $d$, i.e. an
element $\psi_{med} \in Rec(P,H,\phi)$ such that $$ \sum_{\psi' \in \Psi}
d(\psi_{med},\psi') \le \sum_{\psi' \in \Psi} d(\psi,\psi') $$ holds for all
$\psi \in Rec(P,H,\phi)$? For a model where so-called host-switches or
transfers are not allowed, and for a commonly used metric $d$ called the
edit-distance, we show that although the cardinality of $Rec(P,H,\phi)$ can be
super-exponential, it is still possible to compute a geometric median for a set
$\Psi$ in $Rec(P,H,\phi)$ in polynomial time. We expect that this result could
be useful for computing a summary or consensus for a set of reconciliations
(e.g. for a set of suboptimal reconciliations).


  Analysis of DNA samples is an important step in forensics, and the speed of
analysis can impact investigations. Comparison of DNA sequences is based on the
analysis of short tandem repeats (STRs), which are short DNA sequences of 2-5
base pairs. Current forensics approaches use 20 STR loci for analysis. The use
of single nucleotide polymorphisms (SNPs) has utility for analysis of complex
DNA mixtures. The use of tens of thousands of SNPs loci for analysis poses
significant computational challenges because the forensic analysis scales by
the product of the loci count and number of DNA samples to be analyzed. In this
paper, we discuss the implementation of a DNA sequence comparison algorithm by
re-casting the algorithm in terms of linear algebra primitives. By developing
an overloaded matrix multiplication approach to DNA comparisons, we can
leverage advances in GPU hardware and algoithms for Dense Generalized
Matrix-Multiply (DGEMM) to speed up DNA sample comparisons. We show that it is
possible to compare 2048 unknown DNA samples with 20 million known samples in
under 6 seconds using a NVIDIA K80 GPU.


  We present a parameterized analytical model of alchemical molecular binding.
The model describes accurately the free energy profiles of linear
single-decoupling alchemical binding free energy calculations. The parameters
of the model, which are physically motivated, are obtained by fitting model
predictions to numerical simulations. The validity of the model has been
assessed on a set of host-guest complexes. The model faithfully reproduces both
the binding free energy profiles and the probability densities of the
perturbation energy as a function of the alchemical progress parameter
$\lambda$. The model offers a rationalization for the characteristic shape of
the free energy profiles. The parameters obtained from the model are
potentially useful descriptors of the association equilibrium of molecular
complexes.


  Decisions in the cell that lead to its ultimate fate are important for
cellular functions such as proliferation, growth, differentiation, development
and death. Understanding this decision process is imperative for advancements
in the treatment of diseases such as cancer. It is clear that underlying gene
regulatory networks and surrounding environments of the cells are crucial for
function. The self-repressor is a very abundant gene regulatory motif, and is
often believed to have only one cell fate. In this study, we elucidate the
effects of microenvironments mimicking the epigenetic effects on cell fates
through the introduction of inducers capable of binding to a self-repressing
gene product (protein), thus regulating the associated gene. This alters the
effective regulatory binding speed of the self-repressor regulatory protein to
its destination DNA without changing the gene itself. The steady state
observations and real time monitoring of the self-repressor expression dynamics
reveal the emergence of the two cell fates, The simulations are consistent with
the experimental findings. We provide physical and quantitative explanations
for the origin of the two phenotypic cell fates. We find that two cell fates,
rather than a single fate, and their associated switching dynamics emerge from
a change in effective gene regulation strengths. The switching time scale is
quantified. Our results reveal a new mechanism for the emergence of multiple
cell fates. This provides an origin for the heterogeneity often observed among
cell states, while illustrating the influence of microenvironments on cell
fates and their decision-making processes without genetic changes


  We present a deterministic mathematical model for malaria transmission with
waning immunity. The model consists of five non-linear system of differential
equations. We used next generation matrix to derive the basic reproduction
number $R_0$. The disease free equilibrium was computed and its local stability
has been shown by the virtue of the Jacobean matrix. Moreover, using Lyapunov
function theory and LaSalle Invariance Principle we have proved that the
disease free equilibrium is globally asymptotically stable. Conditions for
existence of endemic equilibrium point have been established. A qualitative
study based on bifurcation theory reveals that backward bifurcation occur in
the model. The stable disease free equilibrium of the model coexists with the
stable endemic equilibrium when $R_0<1$. Furthermore, we have shown that
bringing the number of disease (malaria) induced death rate below some
threshold is sufficient enough to eliminate backward bifurcation in the model.


  We prove local existence for classical solutions of a free boundary problem
which arises in one of the biological selection models proposed by Brunet and
Derrida, [3]. The problem we consider describes the limit evolution of
branching brownian particles on the line with death of the leftmost particle at
each creation time as studied in [7]. We use extensively results in [1] and
[2].


  This paper studies mathematical properties of reaction systems that was
introduced by Enrenfeucht and Rozenberg as computational models inspired by
biochemical reaction in the living cells. In particular, we continue the study
on the generative power of functions specified by minimal reaction systems
under composition initiated by Salomaa. Allowing degenerate reaction systems,
functions specified by minimal reaction systems over a quarternary alphabet
that are permutations generate the alternating group on the power set of the
background set.


  The bacterial genome is organized in a structure called the nucleoid by a
variety of associated proteins. These proteins can form complexes on DNA that
play a central role in various biological processes, including chromosome
segregation. A prominent example is the large ParB-DNA complex, which forms an
essential component of the segregation machinery in many bacteria. ChIP-Seq
experiments show that ParB proteins localize around centromere-like parS sites
on the DNA to which ParB binds specifically, and spreads from there over large
sections of the chromosome. Recent theoretical and experimental studies suggest
that DNA-bound ParB proteins can interact with each other to condense into a
coherent 3D complex on the DNA. However, the structural organization of this
protein-DNA complex remains unclear, and a predictive quantitative theory for
the distribution of ParB proteins on DNA is lacking. Here, we propose the
Looping and Clustering (LC) model, which employs a statistical physics approach
to describe protein-DNA complexes. The LC model accounts for the extrusion of
DNA loops from a cluster of interacting DNA-bound proteins. Conceptually, the
structure of the protein-DNA complex is determined by a competition between
attractive protein interactions and the configurational and loop entropy of
this protein-DNA cluster. Indeed, we show that the protein interaction strength
determines the "tightness" of the loopy protein-DNA complex. With this approach
we consider the genomic organization of such a protein-DNA cluster around a
single high-affinity binding site. Thus, our model provides a theoretical
framework to quantitatively compute the binding profiles of ParB-like proteins
around a cognate (parS) binding site.


  Signal transduction networks can form highly interconnected systems within
cells due to network crosstalk, the sharing of input signals between multiple
downstream responses. To better understand the evolutionary design principles
underlying such networks, we study the evolution of crosstalk and the emergence
of specificity for two parallel signaling pathways that arise via gene
duplication and are subsequently allowed to diverge. We focus on a sequence
based evolutionary algorithm and evolve the network based on two physically
motivated fitness functions related to information transmission. Surprisingly,
we find that the two fitness functions lead to very different evolutionary
outcomes, one with a high degree of crosstalk and the other without.


  In the era of high throughput DNA sequencing (HTS) technologies, calculating
the edit distance (i.e., the minimum number of substitutions, insertions, and
deletions between a pair of sequences) for billions of genomic sequences is the
computational bottleneck in todays read mappers. The shifted Hamming distance
(SHD) algorithm proposes a fast filtering strategy that can rapidly filter out
invalid mappings that have more edits than allowed. However, SHD shows high
inaccuracy in its filtering by admitting invalid mappings to be marked as
correct ones. This wastes the execution time and imposes a large computational
burden. In this work, we comprehensively investigate four sources that lead to
the filtering inaccuracy. We propose MAGNET, a new filtering strategy that
maintains high accuracy across different edit distance thresholds and data
sets. It significantly improves the accuracy of pre-alignment filtering by one
to two orders of magnitude. The MATLAB implementations of MAGNET and SHD are
open source and available at: https://github.com/BilkentCompGen/MAGNET.


  Aging remains a fundamental open problem in modern biology. Although there
exist a number of theories on aging on the cellular scale, nearly nothing is
known about how microscopic failures cascade to macroscopic failures of
tissues, organs and ultimately the organism. The goal of this work is to bridge
microscopic cell failure to macroscopic manifestations of aging. We use tissue
engineered constructs to control the cellular-level damage and cell-cell
distance in individual tissues to establish the role of complex interdependence
and interactions between cells in aging tissues. We found that while
microscopic mechanisms drive aging, the interdependency between cells plays a
major role in tissue death, providing evidence on how cellular aging is
connected to its higher systemic consequences.


  S-systems are simple examples of power-law dynamical systems (polynomial
systems with real exponents). For planar S-systems, we study global stability
of the unique positive equilibrium and solve the center problem. Further, we
construct a planar S-system with two limit cycles.


  Meaningful laws of nature must be independent of the units employed to
measure the variables. The principle of similitude (Rayleigh 1915) or
dimensional homogeneity, states that only commensurable quantities (ones having
the same dimension) may be compared, therefore, meaningful laws of nature must
be homogeneous equations in their various units of measurement, a result which
was formalized in the $\rm \Pi$ theorem (Vaschy 1892; Buckingham 1914).
However, most relations in allometry do not satisfy this basic requirement,
including the `3/4 Law' (Kleiber 1932) that relates the basal metabolic rate
and body mass, which it is sometimes claimed to be the most fundamental
biological rate (Brown et al. 2004) and the closest to a law in life sciences
(West \& Brown 2004). Using the $\rm \Pi$ theorem, here we show that it is
possible to construct a unique homogeneous equation for the metabolic rates, in
agreement with data in the literature. We find that the variations in the
dependence of the metabolic rates on body mass are secondary, coming from
variations in the allometric dependence of the heart frequencies. This includes
not only different classes of animals (mammals, birds, invertebrates) but also
different exercise conditions (basal and maximal). Our results demonstrate that
most of the differences found in the allometric exponents (White et al. 2007)
are due to compare incommensurable quantities and that our dimensionally
homogenous formula, unify these differences into a single formulation. We
discuss the ecological implications of this new formulation in the context of
the Malthusian's, Fenchel's and the total energy consumed in a lifespan
relations.


  It is known from the wave-like motion of microtubules in motility assays that
the piconewton forces that motors produce can be sufficient to bend the
filaments. In cellular phenomena such as cytosplasmic streaming, molecular
motors translocate along cytoskeletal filaments, carrying cargo which entrains
fluid. When large numbers of such forced filaments interact through the
surrounding fluid, as in particular stages of oocyte development in
$Drosophila~melanogaster$, complex dynamics are observed, but the detailed
mechanics underlying them has remained unclear. Motivated by these
observations, we study here perhaps the simplest model for these phenomena: an
elastic filament, pinned at one end, acted on by a molecular motor treated as a
point force. Because the force acts tangential to the filament, no matter what
its shape, this "follower-force" problem is intrinsically non-variational, and
thereby differs fundamentally from Euler buckling, where the force has a fixed
direction, and which, in the low Reynolds number regime, ultimately leads to a
stationary, energy-minimizing shape. Through a combination of linear stability
theory, analytical study of a solvable simplified "two-link" model, and
numerical studies of the full elastohydrodynamic equations of motion we
elucidate the Hopf bifurcation that occurs with increasing forcing of a
filament, leading to flapping motion analogous to the high Reynolds number
oscillations of a garden hose with a free end.


  Background: There is uncertainty about the burden of severe outcomes
associated with RSV and influenza in children.
  Methods: We have applied previously developed methodology to estimate the
rates of influenza and RSV-associated hospitalizations with different diagnoses
in different age subgroups of US children between the 2003-04 through the
2009-10 seasons.
  Results: The average annual rates (per 100,000) of influenza and
RSV-associated hospitalizations for any cause excluding asthma (ICD-9 code 493)
in the principal diagnosis with a respiratory cause (ICD-9 codes 460-519) in
the diagnosis (either principal or secondary) were 173.1(95%CI (134.3,212.6))
(influenza) vs. 2336(2209,2463) (RSV) (age <1y); 77.9(51.1,105.2) vs.
625.7(540,711.5) (age 1y); 56.3(40.9,72.2) vs. 324.2(275.1,375) (age 2y);
44.6(32.9,56.6) vs. 172.2(134.2,211) (age 3y); 36.4(27.7,45.3) vs.
89.7(60.6,117.8) (age 4y); 36.6(29.3,43.9) vs. 62.3(38.8,86.4) (ages 5-6y);
20.5(16.3,24.6) vs. 30.2(16.1,43.9) (ages 7-11y); and 16.1(10.6,21.7) vs.
36.9(17.8,55.9) (ages 12-17y). The rates of RSV-associated hospitalizations
with asthma as a secondary (non-principal) diagnosis were 162.6(143.6,181.6)
(age <1y), 144.7(120.8,168.2) (age 1y), 99.5(82.3,116.8) (age 2y),
60.9(46.4,75.3) (age 3y), 23.5(10.7,36.4) (age 4y), 16.1(5.1,27) (ages 5-6y),
11.4(4.7,18.3) (ages 7-11y), and 15.8(4.9,26.7) (ages 12-17y).
  Conclusions: The estimated rates of RSV-associated hospitalizations in young
children are high, with those rates declining rapidly with age. Additionally,
the estimated rates of RSV-associated hospitalizations with asthma as a
secondary diagnosis in young children, combined with data on the prevalence of
asthma in those age groups suggest very high rates of RSV-associated
hospitalizations in young children with asthma. Those children may be
considered as potential target groups for RSV-related mitigation efforts.


  Spirochetes belonging to the Borrelia burgdoferi sensu lato (sl) group cause
Lyme Borreliosis (LB), which is the most commonly reported vector-borne
zoonosis in Europe. B. burgdorferi sl is maintained in nature in a complex
cycle involving Ixodes ricinus ticks and several species of vertebrate hosts.
The transmission dynamics of B. burgdorferi sl is complicated by the varying
competence of animals for different genospecies of spirochetes that, in turn,
vary in their capability of causing disease. In this study, a set of difference
equations simplifying the complex interaction between vectors and their hosts
(competent and not for Borrelia) is built to gain insights into conditions
underlying the dominance of B. lusitaniae (transmitted by lizards to
susceptible ticks) and the maintenance of B. afzelii (transmitted by wild
rodents) observed in a study area in Tuscany, Italy. Findings, in agreement
with field observations, highlight the existence of a threshold for the
fraction of larvae feeding on rodents below which the persistence of B. afzelii
is not possible. Furthermore, thresholds change as nonlinear functions of the
expected number of nymph bites on mice, and the transmission and recovery
probabilities. In conclusion, our model provided an insight into mechanisms
underlying the relative frequency of different Borrelia genospecies, as
observed in field studies.


  Many biologically important ligands of proteins are large, flexible, and
often charged molecules that bind to extended regions on the protein surface.
It is infeasible or expensive to locate such ligands on proteins with standard
methods such as docking or molecular dynamics (MD) simulation. The alternative
approach proposed here is the scanning of a spatial and angular grid around the
protein with smaller fragments of the large ligand. Energy values for complete
grids can be computed efficiently with a well-known Fast Fourier Transform
accelerated algorithm and a physically meaningful interaction model. We show
that the approach can readily incorporate flexibility of protein and ligand.
The energy grids (EGs) resulting from the ligand fragment scans can be
transformed into probability distributions, and then directly compared to
probability distributions estimated from MD simulations and experimental
structural data. We test the approach on a diverse set of complexes between
proteins and large, flexible ligands, including a complex of Sonic Hedgehog
protein and heparin, three heparin sulfate substrates or non-substrates of an
epimerase, a multi-branched supramolecular ligand that stabilizes a
protein-peptide complex, and a flexible zwitterionic ligand that binds to a
surface basin of a Kringle domain. In all cases the EG approach gives results
that are in good agreement with experimental data or MD simulations.


  The time for a DNA sequence to find its homologous depends on a long random
search process inside the cell nucleus. Using polymer models, we model and
compute here the mean first encounter time (MFET) between two sites located on
two different polymer chains and confined by potential wells. We find that
reducing the potential (tethering) forces results in a local polymer
decondensation near the loci and numerical simulations of the polymer model
show that these changes are associated with a reduction of the MFET by several
orders of magnitude. We derive here new asymptotic formula for the MFET,
confirmed by Brownian simulations. We conclude that the acceleration of the
search process after local chromatin decondensation can be used to analyze the
local search step during homology search.


  Cell cycle is an indispensable process in the proliferation and development.
Despite significant efforts, global quantification and physical understanding
are still challenging. In this study, we explored the mechanisms of Xenopus
laevis embryonic cell cycle by quantifying the underlying landscape and flux.
We uncovered the irregular Mexican hat landscape of the Xenopus laevis
embryonic cell cycle with several local basins and barriers on the oscillation
path. The local basins characterize the different phases of Xenopus laevis
embryonic cell cycle and the local barriers represent the checkpoints. The
checkpoint mechanism of cell cycle is revealed by the landscape basins and
barriers. While landscape shape determines the stabilities of the states on the
oscillation path, the curl flux force determines the stability of the cell
cycle flow. Replication is fundamental for biology of living. From our
quantitative study here, we see that replication can not proceed without energy
input. In fact, the curl flux originated from energy or nutrition supply
determines the speed of the cell cycle and guarantees the progression. Speed of
cell cycle is a hallmark of cancer. Through landscape and flux analysis, one
can identify the key elements for controlling the speed. This can help to
design effective strategy for drug discovery against cancer.


  Direct reciprocity is a mechanism for sustaining mutual cooperation in
repeated social dilemma games, where a player would keep cooperation to avoid
being retaliated by a co-player in the future. So-called zero-determinant (ZD)
strategies enable a player to unilaterally set a linear relationship between
the player's own payoff and the co-player's payoff regardless of the strategy
of the co-player. In the present study, we analytically study zero-determinant
strategies in finitely repeated (two-person) prisoner's dilemma games with a
general payoff matrix. Our results are as follows. First, we present the forms
of solutions that extend the known results for infinitely repeated games (with
a discount factor w of unity) to the case of finitely repeated games (0 < w <
1). Second, for the three most prominent ZD strategies, the equalizers,
extortioners, and generous strategies, we derive the threshold value of w above
which the ZD strategies exist. Third, we show that the only strategies that
enforce a linear relationship between the two players' payoffs are either the
ZD strategies or unconditional strategies, where the latter independently
cooperates with a fixed probability in each round of the game, proving a
conjecture previously made for infinitely repeated games.


  We have uncovered that the non-equilibrium network dynamics and global
properties are determined by two essential features: the potential landscape
and the flux landscape. We have found that the funneled potential landscape is
crucial for the stability of the states on the cell cycle, however, the
stabilities of the oscillation states cannot guarantee the stable directional
flows. We have uncovered that the funneled flux landscape is important for the
emergence and maintenance of the stable limit cycle oscillation flow. This work
will allow us to identify the key factors and structure elements of the
networks in determining the stability, speed and robustness of the fission
yeast cell cycle oscillations. We see that the nonequilibriumness characterized
by the degree of detailed balance breaking from the energy pump quantified by
the flux is the cause of the energy dissipation for initiating and sustaining
the replications essential for the origin and evolution of life. Regulating the
cell cycle speed is crucial for designing the prevention and curing strategy of
cancer.


  In abstract terms, ecosystem ecology is about determining when two
ecosystems, superficially different, are alike in some deeper way. An external
observer can choose any ecosystem property as being important. In contrast, two
ecosystems are equivalent from the point of view of the organisms they contain
if and only if for each species, the proportional population growth rate does
not differ between the ecosystems. Comparative studies of ecosystems should
therefore focus on patterns in proportional population growth rates, rather
than patterns in other properties such as relative abundances. Popular
activities such as measuring dissimilarity, and representing dissimilarity via
ordination, can then be done from the point of view of the organisms in
ecosystems. Summarizing the state of an ecosystem under this approach remains
challenging. In general, the dynamics on equivalence classes of ecosystems
defined in this way are structurally different from the dynamics of ecosystems
as seen by an external observer. This may limit the extent to which natural
selection can act on ecosystem structure.


  We use a coarse-grained model to study the conformational changes in two
barley proteins, LTP1 and its ligand adduct isoform LTP1b, that result from
their adsorption to the air-water interface. The model introduces the interface
through hydropathy indices. We justify the model by all-atom simulations. The
choice of the proteins is motivated by making attempts to understand formation
and stability of foam in beer. We demonstrate that both proteins flatten out at
the interface and can make a continuous stabilizing and denser film. We show
that the degree of the flattening depends on the protein -- the layers of LTP1b
should be denser than those of LTP1 -- and on the presence of glycation. It
also depends on the number ($\le 4$) of the disulfide bonds in the proteins.
The geometry of the proteins is sensitive to the specificity of the absent
bonds. We provide estimates of the volume of cavities of the proteins when away
from the interface.


  Using numerical simulations we investigate the space-time properties of a
system in which spirals emerge within coarsening domains, thus giving rise to
non-trivial internal dynamics. Initially proposed in the context of population
dynamics, the studied six-species model exhibits growing domains composed of
three species in a rock-paper-scissors relationship. Through the investigation
of different quantities, such as space-time correlations and the derived
characteristic length, autocorrelation, density of empty sites, and interface
width, we demonstrate that the non-trivial dynamics inside the domains affects
the coarsening process as well as the properties of the interfaces separating
different domains. Domain growth, aging, and interface fluctuations are shown
to be governed by exponents whose values differ from those expected in systems
with curvature driven coarsening.


  We investigate the global dynamics of a general Kermack-McKendrick-type
epidemic model formulated in terms of a system of renewal equations.
Specifically, we consider a renewal model for which both the force of infection
and the infected removal rates are arbitrary functions of the infection age,
$\tau$, and use the direct Lyapunov method to establish the global asymptotic
stability of the equilibrium solutions. In particular, we show that the basic
reproduction number, $R_0$, represents a sharp threshold parameter such that
for $R_0\leq 1$, the infection-free equilibrium is globally asymptotically
stable; whereas the endemic equilibrium becomes globally asymptotically stable
when $R_0 > 1$, i.e. when it exists.


  Biomolecular structure comparison not only reveals evolutionary
relationships, but also sheds light on biological functional properties.
However, traditional definitions of structure or sequence similarity always
involve superposition or alignment and are computationally inefficient. In this
paper, I propose a new method called persistent similarity, which is based on a
newly-invented method in algebraic topology, known as persistent homology.
Different from all previous topological methods, persistent homology is able to
embed a geometric measurement into topological invariants, thus provides a
bridge between geometry and topology. Further, with the proposed persistent
Betti function (PBF), topological information derived from the persistent
homology analysis can be uniquely represented by a series of continuous
one-dimensional (1D) functions. In this way, any complicated biomolecular
structure can be reduced to several simple 1D PBFs for comparison. Persistent
similarity is then defined as the quotient of sizes of intersect areas and
union areas between two correspondingly PBFs. If structures have no significant
topological properties, a pseudo-barcode is introduced to insure a better
comparison. Moreover, a multiscale biomolecular representation is introduced
through the multiscale rigidity function. It naturally induces a multiscale
persistent similarity. The multiscale persistent similarity enables an
objective-oriented comparison. State differently, it facilitates the comparison
of structures in any particular scale of interest. Finally, the proposed method
is validated by four different cases. It is found that the persistent
similarity can be used to describe the intrinsic similarities and differences
between the structures very well.


  Variability on external conditions has important consequences for the
dynamics and the organization of biological systems. In many cases, the
characteristic timescale of environmental changes as well as their correlations
play a fundamental role in the way living systems adapt and respond to it. A
proper mathematical approach to understand population dynamics, thus, requires
approaches more refined than, e.g., simple white-noise approximations. To shed
further light onto this problem, in this paper we propose a unifying framework
based on different analytical and numerical tools available to deal with
"colored" environmental noise. In particular, we employ a "unified colored
noise approximation" to map the original problem into an effective one with
white noise, and then we apply a standard path integral approach to gain
analytical understanding. For the sake of specificity, we present our approach
using as a guideline a variation of the contact process--which can also be seen
as a birth-death process of the Malthus-Verhulst class--where the propagation
or birth rate varies stochastically in time. Our approach allows us to tackle
in a systematic manner some of the relevant questions concerning population
dynamics under environmental variability, such as determining the stationary
population density, establishing the conditions under which a population may
become extinct, and estimating extinction times. We focus on the emerging phase
diagram and its possible phase transitions, underlying how these are affected
by the presence of environmental noise time-correlations.


  We consider multi-chain protein native structures and propose a criterion
that determines whether two chains in the system are entangled or not. The
criterion is based on the behavior observed by pulling at both temini of each
chain simultaneously in the two chains. We have identified about 900 entangled
systems in the Protein Data Bank and provided a more detailed analysis for
several of them. We argue that entanglement enhances the thermodynamic
stability of the system but it may have other functions: burying the
hydrophobic residues at the interface, and increasing the DNA or RNA binding
area. We also study the folding and stretching properties of the knotted
dimeric proteins MJ0366, YibK and bacteriophytochrome. These proteins have been
studied theoretically in their monomeric versions so far. The dimers are seen
to separate on stretching through the tensile mechanism and the characteristic
unraveling force depends on the pulling direction.


  Phylogenetic networks are often constructed by merging multiple conflicting
phylogenetic signals into a directed acyclic graph. It is interesting to
explore whether a network constructed in this way induces biologically-relevant
phylogenetic signals that were not present in the input. Here we show that,
given a multiple alignment A for a set of taxa X and a rooted phylogenetic
network N whose leaves are labelled by X, it is NP-hard to locate the most
parsimonious phylogenetic tree displayed by N (with respect to A) even when the
level of N - the maximum number of reticulation nodes within a biconnected
component - is 1 and A contains only 2 distinct states. (If, additionally, gaps
are allowed the problem becomes APX-hard.) We also show that under the same
conditions, and assuming a simple binary symmetric model of character
evolution, finding the most likely tree displayed by the network is NP-hard.
These negative results contrast with earlier work on parsimony in which it is
shown that if A consists of a single column the problem is fixed parameter
tractable in the level. We conclude with a discussion of why, despite the
NP-hardness, both the parsimony and likelihood problem can likely be
well-solved in practice.


  In these lecture notes, the basic principles of stochastic thermodynamics are
developed starting with a closed system in contact with a heat bath. A
trajectory undergoes Markovian transitions between observable meso-states that
correspond to a coarse-grained description of, e.g., a biomolecule or a
biochemical network. By separating the closed system into a core system and
into reservoirs for ligands and reactants that bind to, and react with the core
system, a description as an open system controlled by chemical potentials and
possibly an external force is achieved. Entropy production and further
thermodynamic quantities defined along a trajectory obey various fluctuation
theorems. For describing fluctuations in a non-equilibrium steady state in the
long-time limit, the concept of a rate function for large deviations from the
mean behaviour is derived from the weight of a trajectory. Universal bounds on
this rate function follow which prove and generalize the thermodynamic
uncertainty relation that quantifies the inevitable trade-off between cost and
precision of any biomolecular process. Specific illustrations are given for
molecular motors, Brownian clocks and enzymatic networks that show how these
tools can be used for thermodynamic inference of hidden properties of a system.


  DDX3X is a human DEAD-box RNA helicase implicated in many cancers and in
viral progression. In addition to the RecA-like catalytic core, DDX3X contains
N- and C-terminal domains. Here, we investigate the substrate and protein
requirements to support the ATPase activity of a DDX3X construct lacking 80
residues from its C-terminal domain. Our data confirmed previous results that
for an RNA molecule to support the ATPase activity of DDX3X it must contain a
single-stranded-double-stranded region. We investigated protein and RNA
structural reasons for this requirement. First, the RNA substrates consisting
only of a double-helix were unable to support DDX3X binding. A single-stranded
RNA substrate supported DDX3X binding, while an RNA substrate consisting of a
single-stranded-double-stranded region not only supported the binding of DDX3X
to RNA, but also promoted DDX3X trimer formation. Thus, the
single-stranded-double-stranded RNA region is needed for DDX3X trimer
formation, and trimer formation is required for ATPase activity. Interestingly,
the dependence of ATP hydrolysis on the protein concentration suggests that the
DDX3X trimer hydrolyzes only two molecules of ATP. Lastly, a DNA substrate that
contains single-stranded-double-stranded regions does not support the ATPase
activity of DDX3X.


  Generally, genotypes and phenotypes are expected to be spatially congruent,
however, in widespread species complexes with few barriers to dispersal,
multiple contact zones, and limited reproductive isolation, discordance between
phenotypes and phylogeographic groups is more probable. Wagtails (Aves:
Motacilla) are a genus of birds with striking plumage pattern variation across
Eurasia. Up to 13 subspecies are recognized within a single species, yet
previous studies using mitochondrial DNA have supported phylogeographic groups
that are inconsistent with subspecies plumage characteristics. In this study,
we investigate the link between phenotypes and genotype by comparing
populations thought to be at different stages along the speciation continuum.
We take a phylogeographic approach by estimating population structure, testing
for isolation by distance, conducting demographic modeling, and estimating the
first time-calibrated species tree for the genus. Our study provides strong
evidence for species-level patterns of differentiation in wagtails, however
population-level differentiation is less pronounced. We find evidence that
three of four widespread Eurasian species exhibit an east-west divide that
contradicts both subspecies taxonomy and phenotypic variation. Both the
geographic location of this divide and time estimates from demographic models
are overlapping in two sympatric species, indicating that coincident
Pleistocene events shaped their histories.


  A theoretical analysis of the unfolding pathway of simple modular proteins in
length- controlled pulling experiments is put forward. Within this framework,
we predict the first module to unfold in a chain of identical units,
emphasizing the ranges of pulling speeds in which we expect our theory to hold.
These theoretical predictions are checked by means of steered molecular
dynamics of a simple construct, specifically a chain composed of two
coiled-coils motives, where anisotropic features are revealed. These
simulations also allow us to give an estimate for the range of pulling
velocities in which our theoretical approach is valid.


  RNA secondary structure folding kinetics is known to be important for the
biological function of certain processes, such as the hok/sok system in E.
coli. Although linear algebra provides an exact computational solution of
secondary structure folding kinetics with respect to the Turner energy model
for tiny (~ 20 nt) RNA sequences, the folding kinetics for larger sequences can
only be approximated by binning structures into macrostates in a coarse-grained
model, or by repeatedly simulating secondary structure folding with either the
Monte Carlo algorithm or the Gillespie algorithm. Here we investigate the
relation between the Monte Carlo algorithm and the Gillespie algorithm. We
prove that asymptotically, the expected time for a K-step trajectory of the
Monte Carlo algorithm is equal to <N> times that of the Gillespie algorithm,
where <N> denotes the Boltzmann expected network degree. If the network is
regular (i.e. every node has the same degree), then the mean first passage time
(MFPT) computed by the Monte Carlo algorithm is equal to MFPT computed by the
Gillespie algorithm multiplied by <N>; however, this is not true for
non-regular networks. In particular, RNA secondary structure folding kinetics,
as computed by the Monte Carlo algorithm, is not equal to the folding kinetics,
as computed by the Gillespie algorithm, although the mean first passage times
are roughly correlated. Simulation software for RNA secondary structure folding
according to the Monte Carlo and Gille- spie algorithms is publicly available,
as is our software to compute the expected degree of the net- work of secondary
structures of a given RNA sequence { see http://bioinformatics.bc.edu/clote/
RNAexpNumNbors.


  A central question in ecology is to understand the ecological processes that
shape community structure. Niche-based theories have emphasized the important
role played by competition for maintaining species diversity. Many of these
insights have been derived using MacArthur's consumer resource model (MCRM) or
its generalizations. Most theoretical work on the MCRM has focused on small
ecosystems with a few species and resources. However theoretical insights
derived from small ecosystems many not scale up large ecosystems with many
resources and species because large systems with many interacting components
often display new emergent behaviors that cannot be understood or deduced from
analyzing smaller systems. To address this shortcoming, we develop a
sophisticated statistical physics inspired cavity method to analyze MCRM when
both the number of species and the number of resources is large. We find that
in this limit, species generically and consistently perturb their environments
and significantly modify available ecological niches. We show how our cavity
approach naturally generalizes niche theory to large ecosystems by accounting
for the effect of this emergent environmental engineering on species invasion
and ecological stability. Our work suggests that environmental engineering is a
generic feature of large, natural ecosystems and must be taken into account
when analyzing and interpreting community structure. It also highlights the
important role that statistical-physics inspired approaches can play in
furthering our understanding of ecology.


  It is a long-standing question whether human sexual and reproductive cycles
are affected predominantly by biology or culture. The literature is mixed with
respect to whether biological or cultural factors best explain the reproduction
cycle phenomenon, with biological explanations dominating the argument. The
biological hypothesis proposes that human reproductive cycles are an adaptation
to the seasonal cycles caused by hemisphere positioning, while the cultural
hypothesis proposes that conception dates vary mostly due to cultural factors,
such as vacation schedule or religious holidays. However, for many countries,
common records used to investigate these hypotheses are incomplete or
unavailable, biasing existing analysis towards primarily Christian countries in
the Northern Hemisphere. Here we show that interest in sex peaks sharply online
during major cultural and religious celebrations, regardless of hemisphere
location. This online interest, when shifted by nine months, corresponds to
documented human birth cycles, even after adjusting for numerous factors such
as language, season, and amount of free time due to holidays. We further show
that mood, measured independently on Twitter, contains distinct collective
emotions associated with those cultural celebrations, and these collective
moods correlate with sex search volume outside of these holidays as well. Our
results provide converging evidence that the cyclic sexual and reproductive
behavior of human populations is mostly driven by culture and that this
interest in sex is associated with specific emotions, characteristic of, but
not limited to, major cultural and religious celebrations.


  Sensors often serve at least two purposes: predicting their input and
minimizing dissipated heat. However, determining whether or not a particular
sensor is evolved or designed to be accurate and efficient is difficult. This
arises partly from the functional constraints being at cross purposes and
partly since quantifying the predictive performance of even in silico sensors
can require prohibitively long simulations. To circumvent these difficulties,
we develop expressions for the predictive accuracy and thermodynamic costs of
the broad class of conditionally Markovian sensors subject to unifilar hidden
semi-Markov (memoryful) environmental inputs. Predictive metrics include the
instantaneous memory and the mutual information between present sensor state
and input future, while dissipative metrics include power consumption and the
nonpredictive information rate. Success in deriving these formulae relies
heavily on identifying the environment's causal states, the input's minimal
sufficient statistics for prediction. Using these formulae, we study the
simplest nontrivial biological sensor model---that of a Hill molecule,
characterized by the number of ligands that bind simultaneously, the sensor's
cooperativity. When energetic rewards are proportional to total predictable
information, the closest cooperativity that optimizes the total energy budget
generally depends on the environment's past hysteretically. In this way, the
sensor gains robustness to environmental fluctuations. Given the simplicity of
the Hill molecule, such hysteresis will likely be found in more complex
predictive sensors as well. That is, adaptations that only locally optimize
biochemical parameters for prediction and dissipation can lead to sensors that
"remember" the past environment.


  Hematopoietic stem cells in mammals are known to reside mostly in the bone
marrow, but also transitively passage in small numbers in the blood.
Experimental findings have suggested that they exist in a dynamic equilibrium,
continuously migrating between these two compartments. Here we construct an
individual-based mathematical model of this process, which is parametrised
using existing empirical findings from mice. This approach allows us to
quantify the amount of migration between the bone marrow niches and the
peripheral blood. We use this model to investigate clonal hematopoiesis, which
is a significant risk factor for hematologic cancers. We also analyse the
engraftment of donor stem cells into non-conditioned and conditioned hosts,
quantifying the impact of different treatment scenarios. The simplicity of the
model permits a thorough mathematical analysis, providing deeper insights into
the dynamics of both the model and of the real-world system. We predict the
time taken for mutant clones to expand within a host, as well as chimerism
levels that can be expected following transplantation therapy, and the
probability that a preconditioned host is reconstituted by donor cells.


  The reconstruction of a species phylogeny from genomic data faces two
significant hurdles: 1) the trees describing the evolution of each individual
gene--i.e., the gene trees--may differ from the species phylogeny and 2) the
molecular sequences corresponding to each gene often provide limited
information about the gene trees themselves. In this paper we consider an
approach to species tree reconstruction that addresses both these hurdles.
Specifically, we propose an algorithm for phylogeny reconstruction under the
multispecies coalescent model with a standard model of site substitution. The
multispecies coalescent is commonly used to model gene tree discordance due to
incomplete lineage sorting, a well-studied population-genetic effect.
  In previous work, an information-theoretic trade-off was derived in this
context between the number of loci, $m$, needed for an accurate reconstruction
and the length of the locus sequences, $k$. It was shown that to reconstruct an
internal branch of length $f$, one needs $m$ to be of the order of $1/[f^{2}
\sqrt{k}]$. That previous result was obtained under the molecular clock
assumption, i.e., under the assumption that mutation rates (as well as
population sizes) are constant across the species phylogeny.
  Here we generalize this result beyond the restrictive molecular clock
assumption, and obtain a new reconstruction algorithm that has the same data
requirement (up to log factors). Our main contribution is a novel reduction to
the molecular clock case under the multispecies coalescent. As a corollary, we
also obtain a new identifiability result of independent interest: for any
species tree with $n \geq 3$ species, the rooted species tree can be identified
from the distribution of its unrooted weighted gene trees even in the absence
of a molecular clock.


  We chemically characterize the symmetries underlying the exact solutions of a
stochastic negatively self-regulating gene. The breaking of symmetry at low
molecular number causes three effects. Two branches of the solution exist,
having high and low switching rates, such that the low switching rate branch
approaches deterministic behavior and the high switching rate branch exhibits
sub-Fano behavior. Average protein number differs from the deterministically
expected value. Bimodal probability distributions appear as the protein number
becomes a readout of the ON/OFF state of the gene.


  Hepatocellular carcinoma (HCC) is the most common primary liver malignancy
and is a leading cause of cancer-related death worldwide. However, cure is not
possible with currently used therapies, and there is not so much approved
targeted therapy for HCC despite numerous attempts and clinical trials. So, it
is essential to identify additional therapeutic strategies to block the growth
of HCC tumors. As a cancer disease, it is associated with aberrant genomic and
transcriptional landscapes. We sought to use a systematic drug repositioning
bioinformatics approach to identify novel candidate drugs to treat HCC, which
considers not only aberrant genomic information, but also the changes of
transcriptional landscapes. First, we screen the collection of HCC feature
genes that frequently mutated in most samples of HCC based on human mutation
data. Then, the gene expression data of HCC in TCGA are combined to classify
the kernel genes of HCC. Finally, the therapeutic score (TS) of each drug is
calculated based on the kolmogorov-smirnov statistical method. Using this
strategy, we identified five drugs that associated with HCC, including three
drugs that could treat HCC and two drugs that might have side-effect on HCC. In
addition, we also make Connectivity Map (CMap) profiles similarity analysis and
KEGG enrichment analysis on drug targets. All these findings suggest that our
approach is effective for accurate discovering novel therapeutic options for
HCC and easily to be extended to other tumors.


  Sucralose is a commonly employed artificial sweetener that appears to
destabilize protein native structures. This is in direct contrast to the
bio-preservative nature of its natural counterpart, sucrose, which enhances the
stability of biomolecules against environmental stress. We have further
explored the molecular interactions of sucralose as compared to sucrose to
illuminate the origin of the differences in their bio-preservative efficacy. We
show that the mode of interactions of sucralose and sucrose in bulk solution
differ subtly using hydration dynamics measurement and computational
simulation. Sucralose does not appear to disturb the native state of proteins
for moderate concentrations (<0.2 M) at room temperature. However, as the
concentration increases, or in the thermally stressed state, sucralose appears
to differ in its interactions with protein leading to the reduction of native
state stability. This difference in interaction appears weak. We explored the
difference in the preferential exclusion model using time-resolved
spectroscopic techniques and observed that both molecules appear to be
effective reducers of bulk hydration dynamics. However, the chlorination of
sucralose appears to slightly enhance the hydrophobicity of the molecule, which
reduces the preferential exclusion of sucralose from the protein-water
interface. The weak interaction of sucralose with hydrophobic pockets on the
protein surface differs from the behavior of sucrose. We experimentally
followed up upon the extent of this weak interaction using isothermal titration
calorimetry (ITC) measurements. We propose this as a possible origin for the
difference in their bio-preservative properties.


  Motivation: Understanding functions of proteins in specific human tissues is
essential for insights into disease diagnostics and therapeutics, yet
prediction of tissue-specific cellular function remains a critical challenge
for biomedicine.
  Results: Here we present OhmNet, a hierarchy-aware unsupervised node feature
learning approach for multi-layer networks. We build a multi-layer network,
where each layer represents molecular interactions in a different human tissue.
OhmNet then automatically learns a mapping of proteins, represented as nodes,
to a neural embedding based low-dimensional space of features. OhmNet
encourages sharing of similar features among proteins with similar network
neighborhoods and among proteins activated in similar tissues. The algorithm
generalizes prior work, which generally ignores relationships between tissues,
by modeling tissue organization with a rich multiscale tissue hierarchy. We use
OhmNet to study multicellular function in a multi-layer protein interaction
network of 107 human tissues. In 48 tissues with known tissue-specific cellular
functions, OhmNet provides more accurate predictions of cellular function than
alternative approaches, and also generates more accurate hypotheses about
tissue-specific protein actions. We show that taking into account the tissue
hierarchy leads to improved predictive power. Remarkably, we also demonstrate
that it is possible to leverage the tissue hierarchy in order to effectively
transfer cellular functions to a functionally uncharacterized tissue. Overall,
OhmNet moves from flat networks to multiscale models able to predict a range of
phenotypes spanning cellular subsystems


  A central problem in population ecology is understanding the consequences of
stochastic fluctuations. Analytically tractable models with Gaussian driving
noise have led to important, general insights, but they fail to capture rare,
catastrophic events, which are increasingly observed at scales ranging from
global fisheries to intestinal microbiota. Due to mathematical challenges,
growth processes with random catastrophes are less well characterized and it
remains unclear how their consequences differ from those of Gaussian processes.
In the face of a changing climate and predicted increases in ecological
catastrophes, as well as increased interest in harnessing microbes for
therapeutics, these processes have never been more relevant. To better
understand them, I revisit here a differential equation model of logistic
growth coupled to density-independent catastrophes that arrive as a Poisson
process, and derive new analytic results that reveal its statistical structure.
First, I derive exact expressions for the model's stationary moments, revealing
a single effective catastrophe parameter that largely controls low order
statistics. Then, I use weak convergence theorems to construct its Gaussian
analog in a limit of frequent, small catastrophes, keeping the stationary
population mean constant for normalization. Numerically computing statistics
along this limit shows how they transform as the dynamics shifts from
catastrophes to diffusions, enabling quantitative comparisons. For example, the
mean time to extinction increases monotonically by orders of magnitude,
demonstrating significantly higher extinction risk under catastrophes than
under diffusions. Together, these results provide insight into a wide range of
stochastic dynamical systems important for ecology and conservation.


  The most frequent form of pairwise synthetic lethality (SL) in metabolic
networks is known as plasticity synthetic lethality (PSL). It occurs when the
simultaneous inhibition of paired functional and silent metabolic reactions or
genes is lethal, while the default of the functional reaction or gene in the
pair is backed up by the activation of the silent one. Based on a complex
systems approach and by using computational techniques on bacterial
genome-scale metabolic reconstructions, we found that the failure of the
functional PSL partner triggers a critical reorganization of fluxes to ensure
viability in the mutant which not only affects the SL pair but a significant
fraction of other interconnected reactions forming what we call a SL cluster.
Interestingly, SL clusters show a strong entanglement both in terms of silent
coessential reactions, which band together to form backup systems, and of other
functional and silent reactions in the metabolic network. This strong overlap,
also detected at the level of genes, mitigates the acquired vulnerabilities and
increased structural and functional costs that pay for the robustness provided
by essential plasticity. Finally, the participation of coessential reactions
and genes in different SL clusters is very heterogeneous and those at the
intersection of many SL clusters could serve as supertargets for more efficient
drug action in the treatment of complex diseases and to elucidate improved
strategies directed to reduce undesired resistance to chemicals in pathogens.


  A population genetics model based on a multitype branching process, or
equivalently a Galton-Watson branching process for multiple alleles, is pre-
sented. The diffusion limit forward Kolmogorov equation is derived for the case
of neutral mutations. The asymptotic stationary solution is obtained and has
the property that the extant population partitions into subpopulations whose
relative sizes are determined by mutation rates. An approximate time-dependent
solution is obtained in the limit of low mutation rates. This solution has the
property that the system undergoes a rapid transition from a drift-dominated
phase to a mutation-dominated phase in which the distribution collapses onto
the asymptotic stationary distribution. The changeover point of the transition
is determined by the per-generation growth factor and mutation rate. The
approximate solution is confirmed using numerical simulations.


  In a metapopulation network, infectious diseases spread widely because of the
travel of individuals. In the present study, we consider a modified
metapopulation Susceptible-Infected-Removed (SIR) model with a latent period,
which we call the SHIR model. In the SHIR model, an infectious period is
divided into two stages. In the first stage, which corresponds to the latent
period, infectious individuals can travel. However, in the second stage, the
same individuals cannot travel since they are seriously ill. Final size
distributions of the metapopulation SIR and SHIR models are simulated with two
different methods and compared. In Monte Carlo simulations, in which the
population is treated as an integer, the distributions show similar behavior.
However, in reaction-diffusion systems, in which the population is treated as a
real number, the final size distribution of the SHIR model has a discontinuous
jump, and that of the SIR model shows a continuous transition. The
discontinuous jump is found to be an artifact that occurs owing to an
inappropriate termination condition.


  Protein misfolding is implicated in many diseases, including the
serpinopathies. For the canonical inhibitory serpin {\alpha}1-antitrypsin
(A1AT), mutations can result in protein deficiencies leading to lung disease,
and misfolded mutants can accumulate in hepatocytes leading to liver disease.
Using all-atom simulations based on the recently developed Bias Functional
algorithm we elucidate how wild-type A1AT folds and how the disease-associated
S (Glu264Val) and Z (Glu342Lys) mutations lead to misfolding. The deleterious Z
mutation disrupts folding at an early stage, while the relatively benign S
mutant shows late stage minor misfolding. A number of suppressor mutations
ameliorate the effects of the Z mutation and simulations on these mutants help
to elucidate the relative roles of steric clashes and electrostatic
interactions in Z misfolding. These results demonstrate a striking correlation
between atomistic events and disease severity and shine light on the mechanisms
driving chains away from their correct folding routes.


  Biological populations are subject to two types of noise: demographic
stochasticity due to fluctuations in the reproductive success of individuals,
and environmental variations that affect coherently the relative fitness of
entire populations. The rate in which the average fitness of a community
increases has been considered so far using models with pure demographic
stochasticity; here we present some theoretical considerations and numerical
results for the general case where environmental variations are taken into
account. When the competition is pairwise, fitness fluctuations are shown to
reduce the speed of evolution, while under global competition the speed
increases due to environmental stochasticity.


  Analogies between evolutionary dynamics and statistical mechanics, such as
Fisher's second-law-like "fundamental theorem of natural selection" and
Wright's "fitness landscapes", have had a deep and fruitful influence on the
development of evolutionary theory. Here I discuss a new conceptual link
between evolution and statistical physics. I argue that natural selection can
be viewed as a coarsening phenomenon, similar to the growth of domain size in
quenched magnets or to Ostwald ripening in alloys and emulsions. In particular,
I show that the most remarkable features of coarsening---scaling and
self-similarity---have strict equivalents in evolutionary dynamics. This
analogy has three main virtues: it brings a set of well-developed mathematical
tools to bear on evolutionary dynamics; it suggests new problems in theoretical
evolution; and it provides coarsening physics with a new exactly soluble model.


  We establish necessary and sufficient conditions for consistent root
reconstruction in continuous-time Markov models with countable state space on
bounded-height trees. Here a root state estimator is said to be consistent if
the probability that it returns to the true root state converges to 1 as the
number of leaves tends to infinity. We also derive quantitative bounds on the
error of reconstruction. Our results answer a question of Gascuel and Steel and
have implications for ancestral sequence reconstruction in a classical
evolutionary model of nucleotide insertion and deletion.


  In evolutionary biology, the speciation history of living organisms is
represented graphically by a phylogeny, that is, a rooted tree whose leaves
correspond to current species and branchings indicate past speciation events.
Phylogenies are commonly estimated from molecular sequences, such as DNA
sequences, collected from the species of interest. At a high level, the idea
behind this inference is simple: the further apart in the Tree of Life are two
species, the greater is the number of mutations to have accumulated in their
genomes since their most recent common ancestor. In order to obtain accurate
estimates in phylogenetic analyses, it is standard practice to employ
statistical approaches based on stochastic models of sequence evolution on a
tree. For tractability, such models necessarily make simplifying assumptions
about the evolutionary mechanisms involved. In particular, commonly omitted are
insertions and deletions of nucleotides -- also known as indels.
  Properly accounting for indels in statistical phylogenetic analyses remains a
major challenge in computational evolutionary biology. Here we consider the
problem of reconstructing ancestral sequences on a known phylogeny in a model
of sequence evolution incorporating nucleotide substitutions, insertions and
deletions, specifically the classical TKF91 process. We focus on the case of
dense phylogenies of bounded height, which we refer to as the taxon-rich
setting, where statistical consistency is achievable. We give the first
polynomial-time ancestral reconstruction algorithm with provable guarantees
under constant rates of mutation. Our algorithm succeeds when the phylogeny
satisfies the "big bang" condition, a necessary and sufficient condition for
statistical consistency in this context.


  Trees have long been used as a graphical representation of species
relationships. However complex evolutionary events, such as genetic
reassortments or hybrid speciations which occur commonly in viruses, bacteria
and plants, do not fit into this elementary framework. Alternatively, various
network representations have been developed. Circular networks are a natural
generalization of leaf-labeled trees interpreted as split systems, that is,
collections of bipartitions over leaf labels corresponding to current species.
Although such networks do not explicitly model specific evolutionary events of
interest, their straightforward visualization and fast reconstruction have made
them a popular exploratory tool to detect network-like evolution in genetic
datasets.
  Standard reconstruction methods for circular networks, such as Neighbor-Net,
rely on an associated metric on the species set. Such a metric is first
estimated from DNA sequences, which leads to a key difficulty: distantly
related sequences produce statistically unreliable estimates. This is
problematic for Neighbor-Net as it is based on the popular tree reconstruction
method Neighbor-Joining, whose sensitivity to distance estimation errors is
well established theoretically. In the tree case, more robust reconstruction
methods have been developed using the notion of a distorted metric, which
captures the dependence of the error in the distance through a radius of
accuracy. Here we design the first circular network reconstruction method based
on distorted metrics. Our method is computationally efficient. Moreover, the
analysis of its radius of accuracy highlights the important role played by the
maximum incompatibility, a measure of the extent to which the network differs
from a tree.


  The effect of stochasticity, in the form of Gaussian white noise, in a
predator-prey model with two distinct time-scales is presented. A supercritical
singular Hopf bifurcation yields a Type II excitability in the deterministic
model. We explore the effect of stochasticity in the excitable regime, leading
to dynamics that are not anticipated by its deterministic counterpart. The
stochastic model admits several kinds of noise-driven mixed-mode oscillations
which capture the intermediate dynamics between two cycles of population
outbreaks. Depending on the strength of noise, the prey population exhibits
intermediate to high-amplitude fluctuations (related to moderate or severe
outbreaks respectively). We classify these fluctuations as isolated or
intermittent or as clusters depending on their recurrences. We study the
distribution of the random variable $N$, representing the number of small
oscillations between successive spikes, as a function of the noise intensity
and the distance to the Hopf bifurcation. The distribution of $N$ is
"asymptotically geometric" with the corresponding parameter related to the
principal eigenvalue of a substochastic Markov chain. Finally, the stochastic
model is transformed into its "normal form" which is used to obtain an estimate
of the probability of repeated outbreaks.


  We develop a mesoscopic modeling framework for diffusion in a crowded
environment, particularly targeting applications in the modeling of living
cells. Through homogenization techniques we effectively coarse-grain a detailed
microscopic description into a previously developed internal state diffusive
framework. The observables in the mesoscopic model correspond to solutions of
macroscopic partial differential equations driven by stochastically varying
diffusion fields in space and time. Analytical solutions and numerical
experiments illustrate the framework.


  How smart can a micron-sized bag of chemicals be? How can an artificial or
real cell make inferences about its environment? From which kinds of
probability distributions can chemical reaction networks sample? We begin
tackling these questions by showing four ways in which a stochastic chemical
reaction network can implement a Boltzmann machine, a stochastic neural network
model that can generate a wide range of probability distributions and compute
conditional probabilities. The resulting models, and the associated theorems,
provide a road map for constructing chemical reaction networks that exploit
their native stochasticity as a computational resource. Finally, to show the
potential of our models, we simulate a chemical Boltzmann machine to classify
and generate MNIST digits in-silico.


  Complex mathematical models of interaction networks are routinely used for
prediction in systems biology. However, it is difficult to reconcile network
complexities with a formal understanding of their behavior. Here, we propose a
simple procedure (called $\bar \varphi$) to reduce biological models to
functional submodules, using statistical mechanics of complex systems combined
with a fitness-based approach inspired by $\textit{in silico}$ evolution. $\bar
\varphi$ works by putting parameters or combination of parameters to some
asymptotic limit, while keeping (or slightly improving) the model performance,
and requires parameter symmetry breaking for more complex models. We illustrate
$\bar \varphi$ on biochemical adaptation and on different models of immune
recognition by T cells. An intractable model of immune recognition with close
to a hundred individual transition rates is reduced to a simple two-parameter
model. $\bar \varphi$ extracts three different mechanisms for early immune
recognition, and automatically discovers similar functional modules in
different models of the same process, allowing for model classification and
comparison. Our procedure can be applied to biological networks based on rate
equations using a fitness function that quantifies phenotypic performance.


  Reversibility of state transitions is intensively studied topic in many
scientific disciplines over many years. In cell biology, it plays an important
role in epigenetic variation of phenotypes, known as phenotypic plasticity.
More interestingly, the cell state reversibility is probably crucial in the
adaptation of population phenotypic heterogeneity to environmental fluctuations
by evolving bet-hedging strategy, which might confer to cancer cells resistance
to therapy. In this article, we propose a formalization of the evolution of
highly reversible states in the environments of periodic variability. Two
interrelated models of heterogeneous cell populations are proposed and their
behavior is studied. The first model captures selection dynamics of the cell
clones for the respective levels of phenotypic reversibility. The second model
focuses on the interplay between reversibility and drug resistance in the
particular case of cancer. Overall, our results show that the threshold
dependencies are emergent features of the investigated model with eventual
therapeutic relevance. Presented examples demonstrate importance of taking into
account cell to cell heterogeneity within a system of clones with different
reversibility quantified by appropriately chosen genetic and epigenetic entropy
measures.


  The size structure of fish-communities is an emergent high-level property of
marine food webs responsive to changes in structure and function. To measure
this food web property using data arising from routine fisheries surveys, a
simple metric known as Typical Length has been proposed as more suitable than
the Large Fish Indicator, which has been highly engineered to be responsive to
fishing pressure. Typical Length avoids the inherent dependence of the Large
Fish Indicator on a parameter that requires case-by-case adjustments. Using
IBTS survey time series for five spatial subdivisions of the Greater North Sea,
we show that the Typical Length can provide information equivalent to the Large
Fish Indicator when fishing is likely the strongest driver, but differences can
also arise. In this example, Typical Length exhibits smaller random
fluctuations ("noise") than the Large Fish Indicator. Typical Length is also
more adaptable than the Large Fish Indicator and can be easily applied to
monitor pelagic fish in addition to demersal fish, and together with
information on the potential growth of the fish community, a proxy of which can
be derived from the Mean Maximum Length indicator, it is possible to partition
change in community composition from change in size structure. This suggests
that Typical Length is an improvement over the Large Fish Indicator as a food
web indicator with the potential to offer further insight when considered in
conjunction with indicators of community composition.


  Humanity faces numerous problems of common-pool resource appropriation. This
class of multi-agent social dilemma includes the problems of ensuring
sustainable use of fresh water, common fisheries, grazing pastures, and
irrigation systems. Abstract models of common-pool resource appropriation based
on non-cooperative game theory predict that self-interested agents will
generally fail to find socially positive equilibria---a phenomenon called the
tragedy of the commons. However, in reality, human societies are sometimes able
to discover and implement stable cooperative solutions. Decades of behavioral
game theory research have sought to uncover aspects of human behavior that make
this possible. Most of that work was based on laboratory experiments where
participants only make a single choice: how much to appropriate. Recognizing
the importance of spatial and temporal resource dynamics, a recent trend has
been toward experiments in more complex real-time video game-like environments.
However, standard methods of non-cooperative game theory can no longer be used
to generate predictions for this case. Here we show that deep reinforcement
learning can be used instead. To that end, we study the emergent behavior of
groups of independently learning agents in a partially observed Markov game
modeling common-pool resource appropriation. Our experiments highlight the
importance of trial-and-error learning in common-pool resource appropriation
and shed light on the relationship between exclusion, sustainability, and
inequality.


  Network-based computational approaches to predict unknown genes associated
with certain diseases are of considerable significance for uncovering the
molecular basis of human diseases. In this paper, we proposed a kind of new
disease-gene-prediction methods by combining the path-based similarity with the
community structure in the human protein-protein interaction network. Firstly,
we introduced a set of path-based similarity indices, a novel community-based
similarity index, and a new similarity combining the path-based similarity
index. Then we assessed the statistical significance of the measures in
distinguishing the disease genes from non-disease genes, to confirm their
availability in predicting disease genes. Finally, we applied these measures to
the disease-gene prediction of single disease-gene family, and analyzed the
performance of these measures in disease-gene prediction, especially the effect
of the community structure on the prediction performance in detail. The results
indicated that genes associated with the same or similar diseases commonly
reside in the same community of the protein-protein interaction network, and
the community structure is greatly helpful for the disease-gene prediction.


  In the article I study the evolutionary adaptivity of two simple population
models, based on either altruistic or egoistic law of energy exchange. The
computational experiments show the convincing advantage of the altruists, which
brings us to a small discussion about genetic algorithms and extraterrestrial
life.