  It is highly believed that the individuals' mobility plays an important role
in phase transition in animal collective motion. Here, we propose a model to
study the effects of individuals' mobility in a distributed animal collective
decision-making process, during which each individual faces two options with
equal quality. We implement the quorum response rule, a type of social
interaction rule which is taxonomically recognized in animal collective
decision-making, as the sole interaction rule. After the introduction of
individuals' mobility, we find that the group can reach a consensus decision at
one of the options at some critical points even the interaction is local. This
result is an obvious contrast to the stationary individuals, the population of
which is always equally distributed between the two options with fluctuations.
In order to explore the information dynamics, we introduce an important
information-theoretic measure, mutual information, to study the critical
behaviors. Furthermore, we study the case when individuals interact globally,
and also find some qualitative similar critical behaviors.


  Statistical analysis of evolutionary-related protein sequences provides
insights about their structure, function, and history. We show that Restricted
Boltzmann Machines (RBM), designed to learn complex high-dimensional data and
their statistical features, can efficiently model protein families from
sequence information. We here apply RBM to twenty protein families, and present
detailed results for two short protein domains, Kunitz and WW, one long
chaperone protein, Hsp70, and synthetic lattice proteins for benchmarking. The
features inferred by the RBM are biologically interpretable: they are related
to structure (such as residue-residue tertiary contacts, extended secondary
motifs ($\alpha$-helix and $\beta$-sheet) and intrinsically disordered
regions), to function (such as activity and ligand specificity), or to
phylogenetic identity. In addition, we use RBM to design new protein sequences
with putative properties by composing and turning up or down the different
modes at will. Our work therefore shows that RBM are a versatile and practical
tool to unveil and exploit the genotype-phenotype relationship for protein
families.


  Groups of humans are often able to find ways to cooperate with one another in
complex, temporally extended social dilemmas. Models based on behavioral
economics are only able to explain this phenomenon for unrealistic stateless
matrix games. Recently, multi-agent reinforcement learning has been applied to
generalize social dilemma problems to temporally and spatially extended Markov
games. However, this has not yet generated an agent that learns to cooperate in
social dilemmas as humans do. A key insight is that many, but not all, human
individuals have inequity averse social preferences. This promotes a particular
resolution of the matrix game social dilemma wherein inequity-averse
individuals are personally pro-social and punish defectors. Here we extend this
idea to Markov games and show that it promotes cooperation in several types of
sequential social dilemma, via a profitable interaction with policy
learnability. In particular, we find that inequity aversion improves temporal
credit assignment for the important class of intertemporal social dilemmas.
These results help explain how large-scale cooperation may emerge and persist.


  The replicator equation is one of the fundamental tools to study evolutionary
dynamics in well-mixed populations. This paper contributes to the literature on
evolutionary graph theory, providing a version of the replicator equation for a
family of connected networks with communities, where nodes in the same
community have the same degree. This replicator equation is applied to the
study of different classes of games, exploring the impact of the graph
structure on the equilibria of the evolutionary dynamics.


  Motivation: Word-based or `alignment-free' methods for phylogeny
reconstruction are much faster than traditional approaches, but they are
generally less accurate. Most of these methods calculate pairwise distances for
a set of input sequences, for example from word frequencies, from so-called
spaced-word matches or from the average length of common substrings.
  Results: In this paper, we propose the first word-based approach to tree
reconstruction that is based on multiple sequence comparison and Maximum
Likelihood. Our algorithm first samples small, gap-free alignments involving
four taxa each. For each of these alignments, it then calculates a quartet tree
and, finally, the program Quartet MaxCut is used to infer a super tree topology
for the full set of input taxa from the calculated quartet trees. Experimental
results show that trees calculated with our approach are of high quality.
  Availability: The source code of the program is available at
https://github.com/tdencker/multi-SpaM
  Contact: thomas.dencker@stud.uni-goettingen.de


  Efficient energy transduction is one driver of evolution; and thus
understanding biomolecular energy transduction is crucial to understanding
living organisms. As an energy-orientated modelling methodology, bond graphs
provide a useful approach to describing and modelling the efficiency of living
systems. This paper gives some new results on the efficiency of metabolism
based on bond graph models of the key metabolic processes: glycolysis.


  We introduce KiNetX, a fully automated meta-algorithm for the kinetic
analysis of complex chemical reaction networks derived from semi-accurate but
efficient electronic structure calculations. It is designed to (i) accelerate
the automated exploration of such networks, and (ii) cope with model-inherent
errors in electronic structure calculations on elementary reaction steps. We
developed and implemented KiNetX to possess three features. First, KiNetX
evaluates the kinetic relevance of every species in a (yet incomplete) reaction
network to confine the search for new elementary reaction steps only to those
species that are considered possibly relevant. Second, KiNetX identifies and
eliminates all kinetically irrelevant species and elementary reactions to
reduce a complex network graph to a comprehensible mechanism. Third, KiNetX
estimates the sensitivity of species concentrations toward changes in
individual rate constants (derived from relative free energies), which allows
us to systematically select the most efficient electronic structure model for
each elementary reaction given a predefined accuracy. The novelty of KiNetX
consists in the rigorous propagation of correlated free-energy uncertainty
through all steps of our kinetic analyis. To examine the performance of KiNetX,
we developed AutoNetGen. It semirandomly generates chemistry-mimicking reaction
networks by encoding chemical logic into their underlying graph structure.
AutoNetGen allows us to consider a vast number of distinct chemistry-like
scenarios and, hence, to discuss assess the importance of rigorous uncertainty
propagation in a statistical context. Our results reveal that KiNetX reliably
supports the deduction of product ratios, dominant reaction pathways, and
possibly other network properties from semi-accurate electronic structure data.


  Habitat destruction threatens the viability of many populations, but its full
consequences can take considerable time to unfold. Much of the discourse
surrounding extinction debts--the number of species that persist transiently
following habitat loss, despite being headed for extinction--frames ultimate
population crashes as the means of settling the debt. However, slow population
decline also opens an opportunity to repay the debt by restoring habitat. The
timing necessary for such habitat restoration to rescue a population from
extinction has not been well studied. Here we determine habitat restoration
deadlines for a spatially implicit Levins/Tilman population model modified by
an Allee effect. We find that conditions that hinder detection of an extinction
debt also provide forgiving restoration timeframes. Our results highlight the
importance of transient dynamics in restoration and suggest the beginnings of
an analytic theory to understand outcomes of temporary press perturbations in a
broad class of ecological systems.


  The assembly and maturation of viruses with icosahedral capsids must be
coordinated with icosahedral symmetry. The icosahedral symmetry imposes also
the restrictions on the cooperative specific interactions between genomic
RNA/DNA and coat proteins that should be reflected in quasi-regular
segmentation of viral genomic sequences. Combining discrete direct and double
Fourier transforms, we studied the quasi-regular large-scale segmentation in
genomic sequences of different ssRNA, ssDNA, and dsDNA viruses. The particular
representatives included satellite tobacco mosaic virus and the strains of
satellite tobacco necrosis virus, STNV-C, STNV-1, STNV-2, Escherichia phages
MS2, phiX174, alpha3, and HK97, and Simian virus 40. In all their genomes, we
found the significant quasi-regular segmentation of genomic sequences related
to the virion assembly and the genome packaging within icosahedral capsid. We
also found good correspondence between our results and available cryo-electron
microscopy data on capsid structures and genome packaging in these viruses.
Fourier analysis of genomic sequences provides the additional insight into
mechanisms of hierarchical genome packaging and may be used for verification of
the concepts of 3-fold or 5-fold intermediates in virion assembly. The results
of sequence analysis should be taken into account at the choice of models and
data interpretation. They also may be helpful for the development of antiviral
drugs.


  This work investigates the interactions of H-NS proteins and bacterial
genomic DNA through computer simulations performed with a coarse-grained model.
The model was developed specifically to study the switch of H-NS proteins from
the DNA-stiffening to the DNA-bridging mode, which has been observed repeatedly
upon addition of multivalent cations to the buffer, but is still not
understood. Unravelling the corresponding mechanism is all the more crucial, as
the regulation properties of H-NS proteins, as well as other nucleoid proteins,
are linked to their DNA-binding properties. The simulations reported here
support a mechanism, according to which the primary role of multivalent cations
consists in decreasing the strength of H-NS/DNA interactions compared to
H-NS/H-NS interactions, with the latter ones becoming energetically favored
with respect to the former ones above a certain threshold of the effective
valency of the cations of the buffer. Below the threshold, H-NS dimers form
filaments, which stretch along the DNA molecule but are quite inefficient in
bridging genomically distant DNA sites (DNA-stiffening mode). In contrast, just
above the threshold, H-NS dimers form 3-dimensional clusters, which are able to
connect DNA sites that are distant from the genomic point of view (DNA-bridging
mode). The model provides clear rationales for the experimental observations
that the switch between the two modes is a threshold effect and that the
ability of H-NS dimers to form higher order oligomers is crucial for their
bridging capabilities.


  Genomic data are becoming increasingly valuable as we develop methods to
utilize the information at scale and gain a greater understanding of how
genetic information relates to biological function. Advances in synthetic
biology and the decreased cost of sequencing are increasing the amount of
privately held genomic data. As the quantity and value of private genomic data
grows, so does the incentive to acquire and protect such data, which creates a
need to store and process these data securely. We present an algorithm for the
Secure Interrogation of Genomic DataBases (SIG-DB). The SIG-DB algorithm
enables databases of genomic sequences to be searched with an encrypted query
sequence without revealing the query sequence to the Database Owner or any of
the database sequences to the Querier. SIG-DB is the first application of its
kind to take advantage of locality-sensitive hashing and homomorphic encryption
to allow generalized sequence-to-sequence comparisons of genomic data.


  Autophagy is an intracellular degradation process that employs complex
membrane dynamics to isolate and break down cellular components. However, many
unanswered questions remain concerning remodeling of autophagic membranes.
Here, we focus on the advantages of theoretical modelling to study the
formation of autophagosomes and to understand the origin of autophagosomal
membranes. Starting from the well-defined geometry of final autophagosomes we
ask the question of how these organelles can be formed by combining various
preautophagosomal membranes such as vesicles, membrane tubules or sheets. We
analyse the geometric constraints of autophagosome formation by taking the area
of the precursor membranes and their internal volume into account. Our results
suggest that vesicle fusion contributes little to the formation of
autophagosomes. In the second part we quantify the curvature of the precursors
and report that the formation of autophagosomes is associated with a strong
relaxation of membrane curvature energy. This effect we find for a wide range
of membrane asymmetries. It is especially strong for small distances between
both autophagosomal membranes, as observed in vivo. We quantify the membrane
bending energies of all precursors by considering membrane asymmetries. We
propose that the generation and supply of pre-autophagosomal membranes is one
limiting step for autophagosome formation.


  Genome and metagenome comparisons based on large amounts of next-generation
sequencing (NGS) data pose significant challenges for alignment-based
approaches due to the huge data size and the relatively short length of the
reads. Alignment-free approaches based on the counts of word patterns in NGS
data do not depend on the complete genome and are generally computationally
efficient. Thus, they contribute significantly to genome and metagenome
comparison. Recently, novel statistical approaches have been developed for the
comparison of both long and shotgun sequences. These approaches have been
applied to many problems including the comparison of gene regulatory regions,
genome sequences, metagenomes, binning contigs in metagenomic data,
identification of virus-host interactions, and detection of horizontal gene
transfers. We provide an updated review of these applications and other related
developments of word-count based approaches for alignment-free sequence
analysis.


  Biological systems reach hierarchical complexity that has no counterpart
outside the realm of biology. Undoubtedly, biological entities obey the
fundamental physical laws. Can today's physics provide an explanatory framework
for understanding the evolution of biological complexity? We argue here that
the physical foundation for understanding the origin and evolution of
complexity can be envisaged at the interface between the theory of frustrated
states resulting in pattern formation in glass-like media and the theory of
self-organized criticality (SOC). On the one hand, SOC has been shown to emerge
in spin glass systems of high dimensionality. On the other hand, SOC is often
viewed as the most appropriate physical description of evolutionary transitions
in biology. We unify these two faces of SOC by showing that emergence of
complex features in biological evolution typically if not always is triggered
by frustration that is caused by competing interactions at different
organizational levels. Competing interactions and frustrated states permeate
biology at all organizational levels and are tightly linked to the ubiquitous
competition for limiting resources. This perspective extends from the
comparatively simple phenomena occurring in glasses to large-scale events of
biological evolution, such as major evolutionary transitions. We therefore
submit that frustration caused by competing interactions in multidimensional
systems is the general driving force behind the emergence of complexity, within
and beyond the domain of biology.


  We revisit the model by Wiser, Ribeck, and Lenski (Science \textbf{342}
(2013), 1364--1367), which describes how the mean fitness increases over time
due to beneficial mutations in Lenski's long-term evolution experiment. We
develop the model further both conceptually and mathematically. Conceptually,
we describe the experiment with the help of a Cannings model with mutation and
selection, where the latter includes diminishing returns epistasis. The
analysis sheds light on the growth dynamics within every single day and reveals
a runtime effect, that is, the shortening of the daily growth period with
increasing fitness; and it allows to clarify the contribution of epistasis to
the mean fitness curve. Mathematically, we explain rigorous results in terms of
a law of large numbers (in the limit of infinite population size and for a
certain asymptotic parameter regime), and present approximations based on
heuristics and supported by simulations for finite populations.


  Giant unilamellar vesicles (GUVs), are a convenient tool to study
membrane-bound processes using optical microscopy. An increasing number of
studies highlights the potential of these model membranes when addressing
questions in membrane biophysics and cell biology. Among them, phase
transitions and domain formation, dynamics and stability in raft-like mixtures
are probably some of the most intensively investigated. In doing so, many
research teams rely on standard protocols for GUV preparation and handling
involving the use of sugar solutions. Here, we demonstrate that following such
a standard approach can lead to abnormal formation of micron-sized domains in
GUVs grown from only a single phospholipid. The membrane heterogeneity is
visualized by means of a small fraction (0.1 mol%) of a fluorescent lipid dye.
For dipalmitoylphosphatidylcholine GUVs, different types of membrane
heterogeneities were detected. First, an unexpected formation of micron-sized
dye-depleted domains was observed upon cooling. These domains nucleated about
10 K above the lipid main phase transition temperature, TM. In addition, upon
further cooling of the GUVs down to the immediate vicinity of TM, stripe-like
dye-enriched structures around the domains are detected. The micron-sized
domains in quasi single-component GUVs were observed also when using two other
lipids. Whereas the stripe structures are related to the phase transition of
the lipid, the dye-excluding domains seem to be caused by traces of impurities
present in the glucose. Supplementing glucose solutions with nm-sized liposomes
at millimolar lipid concentration suppresses the formation of the micron-sized
domains, presumably by providing competitive binding of the impurities to the
liposome membrane in excess. It is likely that such traces of impurities can
significantly alter lipid phase diagrams and cause differences among reported
ones.


  Predicting how cooperative behavior arises in the thermodynamic limit is one
of the outstanding problems in evolutionary game theory. For two player games,
cooperation is seldom the Nash equilibrium. However, in the thermodynamic limit
cooperation is the natural recourse regardless of whether we are dealing with
humans or animals. In this work, we use the analogy with the Ising model to
predict how cooperation arises in the thermodynamic limit.


  Macromolecular crowding plays a principal role in a wide range of biological
processes including gene expression, chromosomal compaction, and viral
infection. However, the impact that crowding has on the dynamics of nucleic
acids remains a topic of debate. To address this problem, we use
single-molecule fluorescence microscopy and custom particle-tracking algorithms
to investigate the impact of varying macromolecular crowding conditions on the
transport and conformational dynamics of large DNA molecules. Specifically, we
measure the mean-squared center-of-mass displacements, as well as the
conformational size, shape, and fluctuations, of individual 115 kbp DNA
molecules diffusing through various in vitro solutions of crowding polymers. We
determine the role of crowder structure and concentration, as well as ionic
conditions, on the diffusion and configurational dynamics of DNA. We find that
branched, compact crowders (10 kDa PEG, 420 kDa Ficoll) drive DNA to compact,
whereas linear, flexible crowders (10 kDa, 500 kDa dextran) cause DNA to
elongate. Interestingly, the extent to which DNA mobility is reduced by
increasing crowder concentrations appears largely insensitive to crowder
structure (branched vs linear), despite the highly different configurations DNA
assumes in each case. We also characterize the role of ionic conditions on
crowding-induced DNA dynamics. We show that both DNA diffusion and
conformational size exhibit an emergent non-monotonic dependence on salt
concentration that is not seen in the absence of crowders.


  Double-stranded DNA translocates through sufficiently large nanopores either
in a linear, single-file fashion or in a folded hairpin conformation when
captured somewhere along its length. We show that the folding state of DNA can
be controlled by changing the electrolyte concentration, pH and PEG content of
the measurement buffer. At 1 M LiCl or 0.35 M KCl in neutral pH, single-file
translocations make up more than 90% of the total. We attribute the effect to
the onset of electroosmotic flow from the pore at low ionic strength. Our
hypothesis on the critical role of flows is supported by the preferred
orientation of entry of a strand that has been folded into a multi-helix
structure at one end. Control over DNA folding is critical for nanopore sensing
approaches that use modifications along a DNA strand and the associated
secondary current drops to encode information.


  The Machine Recognition of Crystallization Outcomes (MARCO) initiative has
assembled roughly half a million annotated images of macromolecular
crystallization experiments from various sources and setups. Here,
state-of-the-art machine learning algorithms are trained and tested on
different parts of this data set. We find that more than 94% of the test images
can be correctly labeled, irrespective of their experimental origin. Because
crystal recognition is key to high-density screening and the systematic
analysis of crystallization experiments, this approach opens the door to both
industrial and fundamental research applications.


  A classical example of a mathematical model for oscillations in a biological
system is the Selkov oscillator, which is a simple description of glycolysis.
It is a system of two ordinary differential equations which, when expressed in
dimensionless variables, depends on two parameters. Surprisingly it appears
that no complete rigorous analysis of the dynamics of this model has ever been
given. In this paper several properties of the dynamics of solutions of the
model are established. With a view to studying unbounded solutions a thorough
analysis of the Poincar\'e compactification of the system is given. It is
proved that for any values of the parameters there are solutions which tend to
infinity at late times and are eventually monotone. It is shown that when the
unique steady state is stable any bounded solution converges to the steady
state at late times. When the steady state is unstable it is shown that for
given values of the parameters either there is a unique periodic solution to
which all bounded solutions other than the steady state converge at late times
or there is no periodic solution and all solutions other than the steady state
are unbounded. In the latter case each unbounded solution which tends to
infinity is eventually monotone and each unbounded solution which does not tend
to infinity has the property that each variable takes on arbitrarily large and
small values at arbitrarily late times.


  Assume one has the capability of determining whether a node in a network is
infectious or not by probing them. Then problem of optimizing sentinel
surveillance in networks is to identify the nodes to probe such that an
emerging disease outbreak can be discovered early or reliably. Whether the
emphasis should be on early or reliable detection depends on the scenario in
question. We investigate three objective measures from the literature
quantifying the performance of nodes in sentinel surveillance -- the time to
detection or extinction, the time to detection, and the frequency of detection.
As a basis for the comparison, we use the susceptible-infectious-recovered
model on static and temporal networks of human contacts. We show that, for some
regions of parameter space, the three objective measures can rank the nodes
very differently. This means sentinel surveillance is a class of problems, and
solutions need to chose an objective measure for the particular scenario in
question. As opposed to other problems in network epidemiology, we draw similar
conclusions from the static and temporal networks. Furthermore, we do not find
one type of network structure that predicts the objective measures -- that
depends both on the data set and the SIR parameter values.


  A conflict exists between field biologists and physiologists ("functional
biologists" or "evolutionary ecologists") on the one hand and those working in
molecular evolution ("evolutionary biologists" or "population geneticists") on
the other concerns the relative importance of natural selection and genetic
drift. This paper is concerned with this issue in the case of vertebrates such
as birds, fishes, mammals, and specifically humans, and views the issue in that
context from a multilevel perspective. It proposes that the resolution is that
adaptive selection outcomes occurring at the organism level chain down to
determine outcomes at the genome level. The multiple realizability of higher
level processes at lower levels then causes the adaptive nature of such
processes at the organism level to be largely hidden at the genomic level. The
discussion is further related to the "negative view" of selection, the Evo-Devo
and Extended Evolutionary Synthesis views, and the levels of selection debate,
where processes at the population level can also chain down in a similar way


  A widely used approach to describe the dynamics of gene regulatory networks
is based on the chemical master equation, which considers probability
distributions over all possible combinations of molecular counts. The analysis
of such models is extremely challenging due to their large discrete state
space. We therefore propose a hybrid approximation approach based on a system
of partial differential equations, where we assume a continuous-deterministic
evolution for the protein counts. We discuss efficient analysis methods for
both modeling approaches and compare their performance. We show that the hybrid
approach yields accurate results for sufficiently large molecule counts, while
reducing the computational effort from one ordinary differential equation for
each state to one partial differential equation for each mode of the system.
Furthermore, we give an analytical steady-state solution of the hybrid model
for the case of a self-regulatory gene.


  Several dihedral angles prediction methods were developed for protein
structure prediction and their other applications. However, distribution of
predicted angles would not be similar to that of real angles. To address this
we employed generative adversarial networks (GAN). Generative adversarial
networks are composed of two adversarially trained networks: a discriminator
and a generator. A discriminator distinguishes samples from a dataset and
generated samples while a generator generates realistic samples. Although the
discriminator of GANs is trained to estimate density, GAN model is intractable.
On the other hand, noise-contrastive estimation (NCE) was introduced to
estimate a normalization constant of an unnormalized statistical model and thus
the density function. In this thesis, we introduce noise-contrastive estimation
generative adversarial networks (NCE-GAN) which enables explicit density
estimation of a GAN model. And a new loss for the generator is proposed. We
also propose residue-wise variants of auxiliary classifier GAN (AC-GAN) and
Semi-supervised GAN to handle sequence information in a window. In our
experiment, the conditional generative adversarial network (C-GAN), AC-GAN and
Semi-supervised GAN were compared. And experiments done with improved
conditions were invested. We identified a phenomenon of AC-GAN that
distribution of its predicted angles is composed of unusual clusters. The
distribution of the predicted angles of Semi-supervised GAN was most similar to
the Ramachandran plot. We found that adding the output of the NCE as an
additional input of the discriminator is helpful to stabilize the training of
the GANs and to capture the detailed structures. Adding regression loss and
using predicted angles by regression loss only model could improve the
conditional generation performance of the C-GAN and AC-GAN.


  Background: Network-based interventions are most powerful against epidemics
when the full network structure is known. However, resource constraints often
require decisions based on partial network data. We investigated how the
effectiveness of network-based vaccination schemes varied based on the accuracy
of network data available. Methods: We simulated propagating an SIR process on
static empirical social networks from 75 rural Indian villages. First, we used
regression analysis to predict the percentage of individuals ever infected
(cumulative incidence) based on village-level networks. Second, we simulated
vaccinating 10% of each village at baseline, selecting vaccinees using 1 of 5
network-based approaches: at random (Random); random contacts of random
individuals (Nomination); random high-degree individuals (High Degree); highest
degree (Highest Degree); or most central (Central). The first three approaches
require only sample data; the latter two require full network data. We also
simulated imposing a limit on how many contacts an individual can nominate
(Fixed Choice Design, FCD), which reduces data collection burden but generates
only partially observed networks. Results: We found the mean and standard
deviation of the degree distribution to strongly predict cumulative incidence.
In simulations, Nomination reduced cumulative incidence by one-sixth compared
to Random vaccination; full network methods reduced infection by two-thirds.
High Degree had intermediate effectiveness. Somewhat surprisingly, FCD
truncating individuals' degrees at three was as effective as using complete
networks. Conclusions: Using even partial network information to prioritize
vaccines at either the village or individual level substantially improved
epidemic outcomes. Such approaches may be feasible and effective in outbreak
settings, and full ascertainment of network structure may not be required.


  The two major approaches to studying macroevolution in deep time are the
fossil record and reconstructed relationships among extant taxa from molecular
data. Results based on one approach sometimes conflict with those based on the
other, with inconsistencies often attributed to inherent flaws of one (or the
other) data source. What is unquestionable is that both the molecular and
fossil records are limited reflections of the same evolutionary history, and
any contradiction between them represents a failure of our existing models to
explain the patterns we observe. Fortunately, the different limitations of each
record provide an opportunity to test or calibrate the other, and new
methodological developments leverage both records simultaneously. However, we
must reckon with the distinct relationships between sampling and time in the
fossil record and molecular phylogenies. These differences impact our
recognition of baselines, and the analytical incorporation of age estimate
uncertainty. These differences in perspective also influence how different
practitioners view the past and evolutionary time itself, bearing important
implications for the generality of methodological advancements, and differences
in the philosophical approach to macroevolutionary theory across fields.


  Reliable identification of molecular biomarkers is essential for accurate
patient stratification. While state-of-the-art machine learning approaches for
sample classification continue to push boundaries in terms of performance, most
of these methods are not able to integrate different data types and lack
generalization power, limiting their application in a clinical setting.
Furthermore, many methods behave as black boxes, and we have very little
understanding about the mechanisms that lead to the prediction. While
opaqueness concerning machine behaviour might not be a problem in deterministic
domains, in health care, providing explanations about the molecular factors and
phenotypes that are driving the classification is crucial to build trust in the
performance of the predictive system. We propose Pathway Induced Multiple
Kernel Learning (PIMKL), a novel methodology to reliably classify samples that
can also help gain insights into the molecular mechanisms that underlie the
classification. PIMKL exploits prior knowledge in the form of a molecular
interaction network and annotated gene sets, by optimizing a mixture of
pathway-induced kernels using a Multiple Kernel Learning (MKL) algorithm, an
approach that has demonstrated excellent performance in different machine
learning applications. After optimizing the combination of kernels for
prediction of a specific phenotype, the model provides a stable molecular
signature that can be interpreted in the light of the ingested prior knowledge
and that can be used in transfer learning tasks.


  The folding dynamics of proteins at the single molecule level has been
studied with single-molecule force spectroscopy (SMFS) experiments for twenty
years, but a common standardized method for the analysis of the collected data
and for the sharing among the scientific community members is still not
available. We have developed a new open source tool, Fodis, for the analysis of
the Force-distance curves obtained in SMFS experiments, providing an almost
automatic processing, analysis and classification of the obtained data. Our
method provides also a classification of the possible unfolding pathways and
structural heterogeneity, present during the unfolding of proteins.


  Preserving biodiversity and ecosystem stability is a challenge that can be
pursued through modern statistical mechanics modeling. Here we introduce a
variational maximum entropy-based algorithm to evaluate the entropy in a
minimal ecosystem on a lattice in which two species struggle for survival. The
method quantitatively reproduces the scale-free law of the prey shoals size,
where the simpler mean-field approach fails: the direct near neighbor
correlations are found to be the fundamental ingredient describing the system
self-organized behavior. Furthermore, entropy allows the measurement of
structural ordering, that is found to be a key ingredient in characterizing two
different coexistence behaviors, one where predators form localized patches in
a sea of preys and another where species display more complex patterns. The
general nature of the introduced method paves the way for its application in
many other systems of interest.


  We establish the means by which Mr. Boddy came to transition from a stable
trajectory within the global phase space of Philadelphia, Pennsylvania to a
stable point on the cement floor of an aviary near the west bank of the
Schuylkill River. There exist no documented murder motives, and so the
dynamical interaction leading to the crime must be reconstructed from
circumstantial data. Our investigation proceeds in two stages. First we take an
audio stream recorded within the aviary near the time of death to identify the
local embedding dimension, thereby enumerating the suspects. Second, we
characterize Mr. Boddy's pre- and post-mortem behavior in the phase space in
terms of an attractor that undergoes an abrupt change in stability. A
supercritical Andronov-Hopf bifurcation can explain this transition. Then we
uniquely identify the murderer. Finally, we note long-term plans to construct
an underlying dynamical model capable of predicting the stability of equilibria
in different parameter regimes, in the event that Mr. Boddy is ever murdered
again.


  Individuals in low socioeconomic brackets are considered at-risk for
developing influenza-related complications and often exhibit higher than
average influenza-related hospitalization rates. This disparity has been
attributed to various factors, including restricted access to preventative and
therapeutic health care, limited sick leave, and household structure. Adequate
influenza surveillance in these at-risk populations is a critical precursor to
accurate risk assessments and effective intervention. However, the United
States of America's primary national influenza surveillance system (ILINet)
monitors outpatient healthcare providers, which may be largely inaccessible to
lower socioeconomic populations. Recent initiatives to incorporate
internet-source and hospital electronic medical records data into surveillance
systems seek to improve the timeliness, coverage, and accuracy of outbreak
detection and situational awareness. Here, we use a flexible statistical
framework for integrating multiple surveillance data sources to evaluate the
adequacy of traditional (ILINet) and next generation (BioSense 2.0 and Google
Flu Trends) data for situational awareness of influenza across poverty levels.
We find that zip codes in the highest poverty quartile are a critical
blind-spot for ILINet that the integration of next generation data fails to
ameliorate.


  Nutrients from a flowering plant are shared by its pollinators, giving rise
to competition in the latter. Such exploitative competition of pollinators can
limit their abundance and affect the global organization of the mutualistic
partnership in the plant-pollinator mutualistic community. Here we investigate
the effects of the exploitative competition between pollinators on the
structure and the species abundance of the mutualistic networks which evolve by
changing mutualistic partnership towards higher abundance of species.
Simulations show different emergent network characteristics between plants and
animals; hub plants connected to many pollinators are very rare while a few
super-hub pollinators appear with the exploitative competition included, in
contrast to equally many hubs of both types without the exploitative
competition. More interestingly, the abundance of plant species increases with
increasing the exploitative competition strength. We analyze the inverse of the
generalized interaction matrix in the weak-interaction limit to identify the
leading structural factors relevant to the species abundance, which are shown
to be instrumental in optimizing the network structure to increase the
mutualistic benefit and lower the cost of exploitative competition.


  During the early developmental process of organisms, the formation of the
left-right laterality requires a subtle mechanism, as it is associated with
other principal body axes. Any inherent chiral feature in an egg cell can in
principal trigger this spontaneous breaking of chiral symmetry. Individual
microtubules, major cytoskeletal filaments, are known as chiral objects.
However, to date there lacks convincing evidence of a hierarchical connection
of the molecular nature of microtubules to large-scale chirality, particularly
at the length scale of an entire cell. Here we assemble an in-vitro active
layer, consisting of microtubules and kinesin motor proteins, on a glass
surface. Upon inclusion of methyl cellulose, the layered system exhibits a
long-range active nematic phase, characterized by the global alignment of
gliding MTs. This nematic order spans over the entire system size in the
millimeter range and, remarkably, allows hidden collective chirality to emerge
as counterclockwise global rotation of the active MT layer. The analysis based
on our theoretical model suggests that the emerging global nematic order
results from the local alignment of MTs, stabilized by methyl cellulose. It
also suggests that the global rotation arises from the MTs' intrinsic
curvature, leading to preferential handedness. Given its flexibility, this
layered in-vitro cytoskeletal system enables the study of membranous protein
behavior responsible for important cellular developmental processes.


  Interactions are ubiquitous across biological systems. These interactions can
be abstracted as patterns of connections among distinct units such as genes,
proteins, individual organisms, or species which form a hierarchy of biological
organisation. Connections in this hierarchy are arranged in a nested structure:
gene and protein networks shape phenotypic traits and together constitute
individuals, individuals are embedded within populations, populations within
communities, and communities within ecosystems. This pervasive "nestedness" of
networks can result in propagation of the effects from within-level
interactions at one level to units at higher or lower levels of organization.
The concept of nested biological networks is implicit in a variety of
disciplines ranging from the study of genetic circuits regulating phenotypic
trait expression (Babu et al. 2004), to the study of predator-prey interactions
influencing community composition (Poisot et al. 2016). However, studies
typically only address interactions within and among directly neighbouring
hierarchical levels, such as genotypes and phenotypic traits, or populations
and communities. Here, we formalise nested networks as having nodes that can
contain, or be embedded in, other nodes, and where edges can bridge connections
between sets of embedded nodes. We then argue that explicitly accounting for
network nestedness across levels of organization will encourage integrative
thinking on new interdisciplinary research fronts. We focus on two phenomena in
particular: (i) indirect connections among units can arise from the structure
of connections at higher or lower levels of organisation, (ii) the propagation
of effects across neighbouring hierarchical levels of organization. This
framework of nested interaction networks provides a tool for researchers across
disciplines to conceptualize their work as elements on a common scaffold.


  Multiple-scale mobility is ubiquitous in nature and has become instrumental
for understanding and modeling animal foraging behavior. However, the impact of
individual movements on the long-term stability of populations remains largely
unexplored. We analyze deterministic and stochastic Lotka Volterra systems,
where mobile predators consume scarce resources (prey) confined in patches. In
fragile systems (that is, those unfavorable to species coexistence), the
predator species has a maximized abundance and is resilient to degraded prey
conditions when individual mobility is multiple scaled. Within the L\'evy
flight model, highly superdiffusive foragers rarely encounter prey patches and
go extinct, whereas normally diffusing foragers tend to proliferate within
patches, causing extinctions by overexploitation. L\'evy flights of
intermediate index allow a sustainable balance between patch exploitation and
regeneration over wide ranges of demographic rates. Our analytical and
simulated results can explain field observations and suggest that scale-free
random movements are an important mechanism by which entire populations adapt
to scarcity in fragmented ecosystems.


  In the Ultimatum Game (UG) one player, named "proposer", has to decide how to
allocate a certain amount of money between herself and a "responder". If the
offer is greater than or equal to the responder's minimum acceptable offer
(MAO), then the money is split as proposed, otherwise, neither the proposer nor
the responder get anything. The UG has intrigued generations of behavioral
scientists because people in experiments blatantly violate the equilibrium
predictions that self-interested proposers offer the minimum available non-zero
amount, and self-interested responders accept. Why are these predictions
violated? Previous research has mainly focused on the role of social
preferences. Little is known about the role of general moral preferences for
doing the right thing, preferences that have been shown to play a major role in
other social interactions (e.g., Dictator Game and Prisoner's Dilemma). Here I
develop a theoretical model and an experiment designed to pit social
preferences against moral preferences. I find that, although people recognize
that offering half and rejecting low offers are the morally right things to do,
moral preferences have no causal impact on UG behavior. The experimental data
are indeed well fit by a model according to which: (i) high UG offers are
motivated by inequity aversion and, to a lesser extent, self-interest; (ii)
high MAOs are motivated by inequity aversion.


  Many pressing medical challenges - such as diagnosing disease, enhancing
directed stem cell differentiation, and classifying cancers - have long been
hindered by limitations in our ability to quantify proteins in single cells.
Mass-spectrometry (MS) is poised to transcend these limitations by developing
powerful methods to routinely quantify thousands of proteins and proteoforms
across many thousands of single cells. We outline specific technological
developments and ideas that can increase the sensitivity and throughput of
single cell MS by orders of magnitude and usher in this new age. These advances
will transform medicine and ultimately contribute to understanding biological
systems on an entirely new level.


  We performed a gene co-expression analysis on Lung Squamous Cell Carcinoma
data to find modules (groups) of genes that may highly impact the growth of
these type of tumors. Additionally, we used cancer survival data to relate
modules to prognostic significance in terms of survival time. Analysis on
RNA-seq data revealed modules which are significant in gene enrichment
analysis. Specifically, two genes - RFC4 and ECT2 - have been found to be
significant in terms of survival time. We also performed a second gene
co-expression analysis on a second dataset of microarray data, and many
significant genes found in this analysis could also be found in the RNA-seq
data implying that these genes might indeed play a crucial role in Lung
Squamous Cancer. All the R code for the analysis can be found at:
\url{https://github.com/nafizh/Gene_coexpression_analysis_lung_cancer}


  A new critical effect is predicted in population dispersal. It is based on
the fact that a trade-off between the advantages of mobility and the cost of
mobility breaks with a significant deterioration in living conditions. The
recently developed model of purposeful kinesis (Gorban \& \c{C}abuko\v{g}lu,
Ecological Complexity 33, 2018) is based on the "let well enough alone" idea:
mobility decreases for high reproduction coefficient and, therefore, animals
stay longer in good conditions and leave quicker bad conditions. Mobility has a
cost, which should be measured in the changes of the reproduction coefficient.
Introduction of the cost of mobility into the reproduction coefficient leads to
an equation for mobility. It can be solved in a closed form using Lambert
$W$-function.
  Surprisingly, the "let well enough alone" models with the simple linear cost
of mobility have an intrinsic phase transition: when conditions worsen then the
mobility increases up to some critical value of the reproduction coefficient.
For worse conditions, there is no solution for mobility. We interpret this
critical effect as the complete loss of mobility that is degeneration of
diffusion. Qualitatively, this means that mobility increases with worsening of
conditions up to some limit, and after that, mobility is nullified.


  In this work we consider the computation of Groebner bases of the steady
state ideal of reaction networks equipped with mass-action kinetics.
Specifically, we focus on the role of intermediate species and the relation
between the extended network (with intermediate species) and the core network
(without intermediate species).
  We show that a Groebner basis of the steady state ideal of the core network
always lifts to a Groebner basis of the steady state ideal of the extended
network by means of linear algebra, with a suitable choice of monomial order.
As illustrated with examples, this contributes to a substantial reduction of
the computation time, due mainly to the reduction in the number of variables
and polynomials. We further show that if the steady state ideal of the core
network is binomial, then so is the case for the extended network, as long as
an extra condition is fulfilled. For standard networks, this extra condition
can be visually explored from the network structure alone.


  A topological approach to the study of genetic recombination, based on
persistent homology, was introduced by Chan, Carlsson, and Rabad\'an in 2013.
This associates a sequence of signatures called barcodes to genomic data
sampled from an evolutionary history. In this paper, we develop theoretical
foundations for this approach. First, we present a novel formulation of the
underlying inference problem. Specifically, we introduce and study the novelty
profile, a simple, stable statistic of an evolutionary history which not only
counts recombination events but also quantifies how recombination creates
genetic diversity. We propose that the (hitherto implicit) goal of the
topological approach to recombination is the estimation of novelty profiles. We
then study the problem of obtaining a lower bound on the novelty profile using
barcodes. We focus on a low-recombination regime, where the evolutionary
history can be described by a directed acyclic graph called a galled tree,
which differs from a tree only by isolated topological defects. We show that in
this regime, under a complete sampling assumption, the $1^\mathrm{st}$ barcode
yields a lower bound on the novelty profile, and hence on the number of
recombination events. For $i>1$, the $i^{\mathrm{th}}$ barcode is empty. In
addition, we use a stability principle to strengthen these results to ones
which hold for any subsample of an arbitrary evolutionary history. To establish
these results, we describe the topology of the Vietoris--Rips filtrations
arising from evolutionary histories indexed by galled trees. As a step towards
a probabilistic theory, we also show that for a random history indexed by a
fixed galled tree and satisfying biologically reasonable conditions, the
intervals of the $1^{\mathrm{st}}$ barcode are independent random variables.
Using simulations, we explore the sensitivity of these intervals to
recombination.


  The spreading of epidemics is very much determined by the structure of the
contact network, which may be impacted by the mobility dynamics of the
individuals themselves. In confined scenarios where a small, closed population
spends most of its time in localized environments and has easily identifiable
mobility patterns --such as workplaces, university campuses or schools-- it is
of critical importance to identify the factors controlling the rate of disease
spread. Here we present a discrete-time, metapopulation-based model to describe
the transmission of SIS-like diseases that take place in confined scenarios
where the mobilities of the individuals are not random but, rather, follow
clear recurrent travel patterns. This model allows analytical determination of
the onset of epidemics, as well as the ability to discern which contact
structures are most suited to prevent the infection to spread. It thereby
determines whether common prevention mechanisms, as isolation, are worth
implementing in such a scenario and their expected impact.


  Generative Adversarial Networks (GANs) represent an attractive and novel
approach to generate realistic data, such as genes, proteins, or drugs, in
synthetic biology. Here, we apply GANs to generate synthetic DNA sequences
encoding for proteins of variable length. We propose a novel feedback-loop
architecture, called Feedback GAN (FBGAN), to optimize the synthetic gene
sequences for desired properties using an external function analyzer. The
proposed architecture also has the advantage that the analyzer need not be
differentiable. We apply the feedback-loop mechanism to two examples: 1)
generating synthetic genes coding for antimicrobial peptides, and 2) optimizing
synthetic genes for the secondary structure of their resulting peptides. A
suite of metrics demonstrate that the GAN generated proteins have desirable
biophysical properties. The FBGAN architecture can also be used to optimize
GAN-generated datapoints for useful properties in domains beyond genomics.


  During development of biological organisms, multiple complex structures are
formed. In many instances, these structures need to exhibit a high degree of
order to be functional, although many of their constituents are intrinsically
stochastic. Hence, it has been suggested that biological robustness ultimately
must rely on complex gene regulatory networks and clean-up mechanisms. Here we
explore developmental processes that have evolved inherent robustness against
stochasticity. In the context of the Drosophila eye disc, multiple optical
units, ommatidia, develop into crystal-like patterns. During the larva-to-pupa
stage of metamorphosis, the centers of the ommatidia are specified initially
through the diffusion of morphogens, followed by the specification of R8 cells.
Establishing the R8 cell is crucial in setting up the geometric, and
functional, relationships of cells within an ommatidium and among neighboring
ommatidia. Here we study a mathematical model of these spatio-temporal
processes in the presence of stochasticity, defining and applying measures that
quantify order within the resulting spatial patterns. We observe a universal
sigmoidal response to increasing transcriptional noise. Ordered patterns
persist up to a threshold noise level in the model parameters. As the noise is
further increased past a threshold point of no return, these ordered patterns
rapidly become disordered. Such robustness in development allows for the
accumulation of genetic variation without any observable changes in phenotype.
We argue that the observed sigmoidal dependence introduces robustness allowing
for sizable amounts of genetic variation and transcriptional noise to be
tolerated in natural populations without resulting in phenotype variation.


  In most studies regarding evolutionary game dynamics, the effective payoff, a
quantity that translates the payoff derived from game interactions into
reproductive success, is usually assumed to be a specific function of the
payoff. Meanwhile, the effect of different function forms of effective payoff
on evolutionary dynamics is always left in the basket. With introducing a
generalized mapping that the effective payoff of individuals is a non-negative
function of two variables on selection intensity and payoff, we study how
different effective payoff functions affect evolutionary dynamics in a
symmetrical mutation-selection process. For standard two-strategy two-player
games, we find that under weak selection the condition for one strategy to
dominate the other depends not only on the classical {\sigma}-rule, but also on
an extra constant that is determined by the form of the effective payoff
function. By changing the sign of the constant, we can alter the direction of
strategy selection. Taking the Moran process and pairwise comparison process as
specific models in well-mixed populations, we find that different fitness or
imitation mappings are equivalent under weak selection. Moreover, the sign of
the extra constant determines the direction of one-third law and risk-dominance
for sufficiently large populations. This work thus helps to elucidate how the
effective payoff function as another fundamental ingredient of evolution affect
evolutionary dynamics.


  The Fair Proportion of a species in a phylogenetic tree is a very simple
measure that has been used to assess its value relative to the overall
phylogenetic diversity represented by the tree. It has recently been proved by
Fuchs and Jin to be equal to the Shapley Value of the coallitional game that
sends each subset of species to its rooted Phylogenetic Diversity in the tree.
We prove in this paper that this result extends to the natural translations of
the Fair Proportion and the rooted Phylogenetic Diversity to rooted
phylogenetic networks. We also generalize to rooted phylogenetic networks the
expression for the Shapley Value of the unrooted Phylogenetic Diversity game on
a phylogenetic tree established by Haake, Kashiwada and Su.


  We study weighted particle systems in which new generations are resampled
from current particles with probabilities proportional to their weights. This
covers a broad class of sequential Monte Carlo (SMC) methods, widely-used in
applied statistics and cognate disciplines. We consider the genealogical tree
embedded into such particle systems, and identify conditions, as well as an
appropriate time-scaling, under which they converge to the Kingman n-coalescent
in the infinite system size limit in the sense of finite-dimensional
distributions. Thus, the tractable n-coalescent can be used to predict the
shape and size of SMC genealogies, as we illustrate by characterising the
limiting mean and variance of the tree height. SMC genealogies are known to be
connected to algorithm performance, so that our results are likely to have
applications in the design of new methods as well. Our conditions for
convergence are strong, but we show by simulation that they do not appear to be
necessary.


  MicroRNAs play important roles in many biological processes. Their aberrant
expression can have oncogenic or tumor suppressor function directly
participating to carcinogenesis, malignant transformation, invasiveness and
metastasis. Indeed, miRNA profiles can distinguish not only between normal and
cancerous tissue but they can also successfully classify different subtypes of
a particular cancer. Here, we focus on a particular class of transcripts
encoding polycistronic miRNA genes that yields multiple miRNA components. We
describe clustered MiRNA Master Regulator Analysis (ClustMMRA), a fully
redesigned release of the MMRA computational pipeline (MiRNA Master Regulator
Analysis), developed to search for clustered miRNAs potentially driving cancer
molecular subtyping. Genomically clustered miRNAs are frequently co-expressed
to target different components of pro-tumorigenic signalling pathways. By
applying ClustMMRA to breast cancer patient data, we identified key miRNA
clusters driving the phenotype of different tumor subgroups. The pipeline was
applied to two independent breast cancer datasets, providing statistically
concordant results between the two analysis. We validated in cell lines the
miR-199/miR-214 as a novel cluster of miRNAs promoting the triple negative
subtype phenotype through its control of proliferation and EMT.


  Phylogenetic networks are becoming of increasing interest to evolutionary
biologists due to their ability to capture complex non-treelike evolutionary
processes. From a combinatorial point of view, such networks are certain types
of rooted directed acyclic graphs whose leaves are labelled by, for example,
species. A number of mathematically interesting classes of phylogenetic
networks are known. These include the biologically relevant class of stable
phylogenetic networks whose members are defined via certain "fold-up" and
"un-fold" operations that link them with concepts arising within the theory of,
for example, graph fibrations. Despite this exciting link, the structural
complexity of stable phylogenetic networks is still relatively poorly
understood. Employing the popular tree-based, reticulation-visible, and
tree-child properties which allow one to gauge this complexity in one way or
another, we provide novel characterizations for when a stable phylogenetic
network satisfies either one of these three properties.


  Automatic design with machine learning and molecular simulations has shown a
remarkable ability to generate new and promising drug candidates. Current
models, however, still have problems in simulation concurrency and molecular
diversity. Most methods generate one molecule at a time and do not allow
multiple simulators to run simultaneously. Additionally, better molecular
diversity could boost the success rate in the subsequent drug discovery
process. We propose a new population-based approach using grammatical evolution
named ChemGE. In our method, a large population of molecules are updated
concurrently and evaluated by multiple simulators in parallel. In docking
experiments with thymidine kinase, ChemGE succeeded in generating hundreds of
high-affinity molecules whose diversity is better than that of known inding
molecules in DUD-E.


  The Moran process is a random process that models the spread of genetic
mutations through graphs. If the graph is connected, the process eventually
reaches "fixation", where every vertex is a mutant, or "extinction", where no
vertex is a mutant.
  Our main result is an almost-tight bound on expected absorption time. For all
epsilon > 0, we show that the expected absorption time on an n-vertex graph is
o(n^(3+epsilon)). In fact, we show that it is at most n^3 * exp(O((log log
n)^3)) and that there is a family of graphs where it is Omega(n^3). In the
course of proving our main result, we also establish a phase transition in the
probability of fixation, depending on the fitness parameter r of the mutation.
We show that no similar phase transition occurs for digraphs, where it is
already known that the expected absorption time can also be exponential.
Finally, we give an improved FPRAS for approximating the probability of
fixation. Its running time is independent of the size of the graph when the
maximum degree is bounded and some basic properties of the graph are given.


  Tsallis' non-extensive entropy is extended to incorporate the dependence on
affinities between the microstates of a system. At the core of our construction
of the extended entropy ($\mathcal{H}$) is the concept of the effective number
of dissimilar states, termed the effective diversity ($\mathit{\Delta}$). It is
a unique integrated measure derived from the probability distribution among
states and the affinities between states. The effective diversity is related to
the extended entropy through the Boltzmann's-equation-like relation,
$\mathcal{H}=\ln_{q}\mathit{\Delta}$, in terms of the Tsallis' $q$-logarithm. A
new principle called the Nesting Principle is established, stating that the
effective diversity remains invariant under an arbitrary grouping of the
constituent states. It is shown that this invariance property holds only for
$q=2$; however, the invariance is recovered for general $q$ in the
zero-affinity limit (i.e. the Tsallis and Boltzmann-Gibbs case). Using the
affinity-based extended Tsallis entropy, the microcanonical and the canonical
ensembles are constructed in the presence of general between-state affinities.
It is shown that the classic postulate of equal a priori probabilities no
longer holds but is modified by affinity-dependent terms. As an illustration, a
two-level system is investigated by the extended canonical method, which
manifests that the thermal behaviours of the thermodynamic quantities at
equilibrium are affected by the between-state affinity. Furthermore, some
applications and implications of the affinity-based extended diversity/entropy
for information theory and biodiversity theory are addressed in appendices.


  A central goal of evolutionary biology is to explain the origins and
distribution of diversity across life. Beyond species or genetic diversity, we
also observe diversity in the circuits (genetic or otherwise) underlying
complex functional traits. However, while the theory behind the origins and
maintenance of genetic and species diversity has been studied for decades,
theory concerning the origin of diverse functional circuits is still in its
infancy. It is not known how many different circuit structures can implement
any given function, which evolutionary factors lead to different circuits, and
whether the evolution of a particular circuit was due to adaptive or
non-adaptive processes. Here, we use digital experimental evolution to study
the diversity of neural circuits that encode motion detection in digital
(artificial) brains. We find that evolution leads to an enormous diversity of
potential neural architectures encoding motion detection circuits, even for
circuits encoding the exact same function. Evolved circuits vary in both
redundancy and complexity (as previously found in genetic circuits) suggesting
that similar evolutionary principles underlie circuit formation using any
substrate. We also show that a simple (designed) motion detection circuit that
is optimally-adapted gains in complexity when evolved further, and that
selection for mutational robustness led this gain in complexity.


  Taking into account an evolutionary model of mutations in term of Levy Fights
that was previously constructed, we designed an algorithm to reproduce the
evolutionary dynamics of the Long-Term Evolution Experiment (LTEE) with E. Coli
bacteria. The algorithm enables us to simulate mutations under natural
selection conditions. The results of simulations on competition of clones, mean
fitness, etc., are compared with experimental data. We attained to reproduce
the behavior of the mean fitness of the bacteria cultures, get our own
interpretations and more tuned descriptions of some phenomena taking part
within the experiment, such as fixation and drift processes, clonal
interference and epistasis.


  In this paper, a system architecture is proposed that approximately models
the functionality of metabolic networks. The AND/OR graph model is used to
represent the metabolic network and each processing element in the system
emulates the functionality of a metabolite. The system is implemented on a
graphics processing unit (GPU) as the hardware platform using CUDA environment.
The proposed architecture takes advantage of the inherent parallelism in the
network structure in terms of both pathway and metabolite traversal. The
function of each element is defined such that it can find flux-balanced
pathways. Pathways in both small and large metabolic networks are applied to
the proposed architecture and the results are discussed.


  The data on molecular phylogeny of the Colias complex of South America (SA),
obtained via barcoding a mitochondrial part of genome, are reported. Barcoding
was trialed employing the Barcoding of Life Database platform and comprised 93
specimens of SA Colias sp. plus 2 outgroup specimens (C. philodice guatemalena
from North America). It is established that Colias sp. from SA form a single
monophyletic clade, characterized by notable interspecific mutational
divergences (3.5-4.5%) and moderate intraspecific ones (0.3-0.9%). It also
reveals the occurrence of 3 well-established (notably diverged genetically)
species (C. dimera, C. vauthierii, and C. alticola, the latter re-raised to sp.
level), plus a complex of 3 'emerging' species (C. lesbia, C. euxanthe, and C.
flaveola) which are weakly differentiated genetically. Furthermore, it
disregards the sp. status of taxa weberbaueri, mossi, nigerrima, erika,
blameyi, and mendozina, as certain evidences are being brought to consider them
ssp. of a clinal type of super-sp. C. flaveola s. lat. Importantly, the clade
formed by SA Colias is found to be sister to the one comprising Colias sp.
occurring outside the region (in Eurasia, North America, and Africa). A direct
comparison of the results obtained for SA Colias with those for worldwide
Colias (GenBank data) leads to the conclusion that these are largely divergent
autochtonous subcomplexes, being in fact subgenera s. lat. of genus Colias. The
first subgenus comprising all SA Colias sp. and named Aucolias (Austral Colias,
or of southern hemisphere stem) is sister to the second one (comprising all
Colias sp. from the Old and New Worlds), for which the name Colias (Colias) s.
str. is preserved.


  It has been an old unsolved puzzle to evolutionary theorists on which
mechanisms would increase large-scale cooperation in human societies. Thus, how
such mechanisms operate in a biological network is still not very understood.
This study addresses these questions with empirical evidence from agent-based
models designed to understand these network interactions. Repeated Prisoner's
Dilemma games were designed to study how costly punishment, diversity, and
density of connectivity interact to influence cooperation in a biological
network. There are 1000 rounds in each game made up of 18 players engaged in
pairwise relationship with their neighbors. This study shows three important
interactions. (1) Introducing diversity to costly punishment favors both
cooperation and defection, but not vice versa. Introducing costly punishment to
diversity disfavors defection but favors cooperation. (2) Costly Punishment,
alone, disfavors defection but decreases average payoff. Decreasing the density
of connectivity, Dc, when there is no costly punishment applied, increases
average payoff. (3) A synergy of diversity and decreasing density of
connectivity favors cooperation in a biological network. Furthermore, this
study also suggests a likelihood from empirical findings that spatial
structures may not be favoring cooperation, as is the widely-accepted notion,
but rather disfavoring defection in the global scale.


  In this paper, association results from genome-wide association studies
(GWAS) are combined with a deep learning framework to test the predictive
capacity of statistically significant single nucleotide polymorphism (SNPs)
associated with obesity phenotype. Our approach demonstrates the potential of
deep learning as a powerful framework for GWAS analysis that can capture
information about SNPs and the important interactions between them. Basic
statistical methods and techniques for the analysis of genetic SNP data from
population-based genome-wide studies have been considered. Statistical
association testing between individual SNPs and obesity was conducted under an
additive model using logistic regression. Four subsets of loci after
quality-control (QC) and association analysis were selected: P-values lower
than 1x10-5 (5 SNPs), 1x10-4 (32 SNPs), 1x10-3 (248 SNPs) and 1x10-2 (2465
SNPs). A deep learning classifier is initialised using these sets of SNPs and
fine-tuned to classify obese and non-obese observations. Using a deep learning
classifier model and genetic variants with P-value < 1x10-2 (2465 SNPs) it was
possible to obtain results (SE=0.9604, SP=0.9712, Gini=0.9817, LogLoss=0.1150,
AUC=0.9908 and MSE=0.0300). As the P-value increased, an evident deterioration
in performance was observed. Results demonstrate that single SNP analysis fails
to capture the cumulative effect of less significant variants and their overall
contribution to the outcome in disease prediction, which is captured using a
deep learning framework.


  We consider a particular instance of the lift of controlled systems recently
proposed in the theory of irreversible thermodynamics and show that it leads to
a variational principle for an optimal control in the sense of Pontryagin. Then
we focus on two important applications: in thermodynamics and in evolutionary
biology. In the thermodynamic context, we show that this principle provides a
dynamical implementation of the Second Law, which stabilizes the equilibrium
manifold of a system. In the evolutionary context, we discuss several
interesting features: it provides a robust scheme for the coevolution of the
population and its fitness landscape; it has a clear informational
interpretation; it recovers Price equation naturally; and finally, it extends
standard evolutionary dynamics to include phenomena such as the emergence of
cooperation and the coexistence of qualitatively different phases of evolution,
which we speculate can be associated with Darwinism and punctuated equilibria.


  An emerging approach to data-limited fisheries stock assessment uses
hierarchical multi-stock assessment models to group stocks together, sharing
information from data-rich to data-poor stocks. In this paper, we simulate
data-rich and data-poor fishery and survey data scenarios for a complex of
dover sole stocks. Simulated data for individual stocks were used to compare
estimation performance for single-stock and hierarchical multi-stock versions
of a Schaefer production model. The single-stock and best performing
multi-stock models were then used in stock assessments for the real dover sole
data. Multi-stock models often had lower estimation errors than single-stock
models when assessment data had low statistical power. Relative errors for
productivity and relative biomass parameters were lower for multi-stock
assessment model configurations. In addition, multi-stock models that estimated
hierarchical priors for survey catchability performed the best under data-poor
scenarios. We conclude that hierarchical multi-stock assessment models are
useful for data-limited stocks and could provide a more flexible alternative to
data-pooling and catch only methods; however, these models are subject to
non-linear side-effects of parameter shrinkage. Therefore, we recommend testing
hierarchical multi-stock models in closed-loop simulations before application
to real fishery management systems.


  The most prominent defence of the unicellular parasite Trypanosoma brucei
against the host immune system is a dense coat that comprises a variant surface
glycoprotein (VSG). Despite the importance of the VSG family, no complete
structure of a VSG has been reported. Making use of high-resolution structures
of individual VSG domains, we employed small-angle X-ray scattering to
elucidate the first two complete VSG structures. The resulting models imply
that the linker regions confer great flexibility between domains, which
suggests that VSGs can adopt two main conformations to respond to obstacles and
changes of protein density, while maintaining a protective barrier at all
times. Single-molecule diffusion measurements of VSG in supported lipid
bilayers substantiate this possibility, as two freely diffusing populations
could be detected. This translates into a highly flexible overall topology of
the surface VSG coat, which displays both lateral movement in the plane of the
membrane and variation in the overall thickness of the coat.


  We study the dynamics of secondary infections on networks, in which only the
individuals currently carrying a certain primary infection are susceptible to
the secondary infection. In the limit of large sparse networks, the model is
mapped to a branching process spreading in a random time-sensitive environment,
determined by the dynamics of the underlying primary infection. When both
epidemics follow the Susceptible-Infective-Recovered model, we show that in
order to survive, it is necessary for the secondary infection to evolve on a
timescale that is closely matched to that of the primary infection on which it
depends.


  This book chapter introduces to the problem to which extent search strategies
of foraging biological organisms can be identified by statistical data analysis
and mathematical modeling. A famous paradigm in this field is the Levy Flight
Hypothesis: It states that under certain mathematical conditions Levy flights,
which are a key concept in the theory of anomalous stochastic processes,
provide an optimal search strategy. This hypothesis may be understood
biologically as the claim that Levy flights represent an evolutionary adaptive
optimal search strategy for foraging organisms. Another interpretation,
however, is that Levy flights emerge from the interaction between a forager and
a given (scale-free) distribution of food sources. These hypotheses are
discussed controversially in the current literature. We give examples and
counterexamples of experimental data and their analyses supporting and
challenging them.


  We present a new non-Archimedean model of evolutionary dynamics, in which the
genomes are represented by p-adic numbers. In this model the genomes have a
variable length, not necessarily bounded, in contrast with the classical models
where the length is fixed. The time evolution of the concentration of a given
genome is controlled by a p-adic evolution equation. This equation depends on a
fitness function f and on mutation measure Q. By choosing a mutation measure of
Gibbs type, and by using a p-adic version of the Maynard Smith Ansatz, we show
the existence of threshold function M_{c}(f,Q), such that the long term
survival of a genome requires that its length grows faster than M_{c}(f,Q).
This implies that Eigen's paradox does not occur if the complexity of genomes
grows at the right pace. About twenty years ago, Scheuring and Poole, Jeffares,
Penny proposed a hypothesis to explain Eigen's paradox. Our mathematical model
shows that this biological hypothesis is feasible, but it requires p-adic
analysis instead of real analysis. More exactly, the Darwin-Eigen cycle
proposed by Poole et al. takes place if the length of the genomes exceeds
M_{c}(f,Q).


  Biological organisms are open, adaptve systems that can respond to changes in
environment in specific ways. Adaptation and response can be posed as an
optimization problem, with a tradeoff between the benefit obtained from a
response and the cost of producing environment-specific responses. Using recent
results in stochastic thermodynamics, we formulate the cost as the mutual
information between the environment and the stochastic response. The problem of
designing an optimally performing network now reduces to a problem in rate
distortion theory -- a branch of information theory that deals with lossy data
compression. We find that as the cost of unit information goes down, the system
undergoes a sequence of transitions, corresponding to the recruitment of an
increasing number of responses, thus improving response specificity as well as
the net payoff. We derive formal equations for the transition points and
exactly solve them for special cases. The first transition point, also called
the {\it coding transition}, demarcates the boundary between a passive response
and an active decision-making by the system. We study this transition point in
detail, and derive three classes of asymptotic behavior, corresponding to the
three limiting distributions of the statistics of extreme values. Our work
points to the necessity of a union between information theory and the theory of
adaptive biomolecular networks, in particular metabolic networks.


  The ecological invasion problem in which a weaker exotic species invades an
ecosystem inhabited by two strongly competing native species is modelled by a
three-species competition-diffusion system. It is known that for a certain
range of parameter values competitor-mediated coexistence occurs and complex
spatio-temporal patterns are observed in two spatial dimensions. In this paper
we uncover the mechanism which generates such patterns. Under some assumptions
on the parameters the three-species competition-diffusion system admits two
planarly stable travelling waves. Their interaction in one spatial dimension
may result in either reflection or merging into a single homoclinic wave,
depending on the strength of the invading species. This transition can be
understood by studying the bifurcation structure of the homoclinic wave. In
particular, a time-periodic homoclinic wave (breathing wave) is born from a
Hopf bifurcation and its unstable branch acts as a separator between the
reflection and merging regimes. The same transition occurs in two spatial
dimensions: the stable regular spiral associated to the homoclinic wave
destabilizes, giving rise first to an oscillating breathing spiral and then
breaking up producing a dynamic pattern characterized by many spiral cores. We
find that these complex patterns are generated by the interaction of two
planarly stable travelling waves, in contrast with many other well known cases
of pattern formation where planar instability plays a central role.


  Trait-mediated indirect effects are increasingly acknowledged as important
components in the dynamics of ecological systems. The hamiltonian form of the
LV equations is traditionally modified by adding density dependence to the prey
variable and functional response to the predator variable. Enriching these
non-linear elements with a trait-mediation added to the carrying capacity of
the prey creates the dynamics of critical transitions and hysteretic zones.


  Motivation: Protein-protein interactions (PPIs) are usually modelled as
networks. These networks have extensively been studied using graphlets, small
induced subgraphs capturing the local wiring patterns around nodes in networks.
They revealed that proteins involved in similar functions tend to be similarly
wired. However, such simple models can only represent pairwise relationships
and cannot fully capture the higher-order organization of protein interactions,
including protein complexes. Results: To model the multi-sale organization of
these complex biological systems, we utilize simplicial complexes from
computational geometry. The question is how to mine these new representations
of PPI networks to reveal additional biological information. To address this,
we define simplets, a generalization of graphlets to simplicial complexes. By
using simplets, we define a sensitive measure of similarity between simplicial
complex network representations that allows for clustering them according to
their data types better than clustering them by using other state-of-the-art
measures, e.g., spectral distance, or facet distribution distance. We model
human and baker's yeast PPI networks as simplicial complexes that capture PPIs
and protein complexes as simplices. On these models, we show that our newly
introduced simplet-based methods cluster proteins by function better than the
clustering methods that use the standard PPI networks, uncovering the new
underlying functional organization of the cell. We demonstrate the existence of
the functional geometry in the PPI data and the superiority of our
simplet-based methods to effectively mine for new biological information hidden
in the complexity of the higher order organization of PPI networks.


  T-cell receptor (TCR) repertoire data contain information about infections
that could be used in disease diagnostics and vaccine development, but
extracting that information remains a major challenge. Here we developed a
statistical framework to detect TCR clone proliferation and contraction from
longitudinal repertoire data. We applied this framework to data from three
pairs of identical twins immunized with the yellow fever vaccine. We identified
500-1500 responding TCRs in each donor and validated them using three
independent assays. While the responding TCRs were mostly private, albeit with
higher overlap between twins, they could be well predicted using a classifier
based on sequence similarity. Our method can also be applied to samples
obtained post-infection, making it suitable for systematic discovery of new
infection-specific TCRs in the clinic.


  Biological systems transduce signals from their surroundings through a myriad
of pathways. In this paper, we describe signal transduction as a communication
system: the signal transduction receptor acts as the receiver in this system,
and can be modeled as a finite-state Markov chain with transition rates
governed by the input signal. Using this general model, we give the mutual
information under IID inputs in discrete time, and obtain the mutual
information in the continuous-time limit. We show that the mutual information
has a concise closed-form expression with clear physical significance. We also
give a sufficient condition under which the Shannon capacity is achieved with
IID inputs. We illustrate our results with three examples: the light-gated
Channelrhodopsin-2 (ChR2) receptor; the ligand-gated nicotinic acetylcholine
(ACh) receptor; and the ligand-gated Calmodulin (CaM) receptor. In particular,
we show that the IID capacity of the ChR2 receptor is equal to its Shannon
capacity. We finally discuss how the results change if only certain properties
of each state can be observed, such as whether an ion channel is open or
closed.


  Nonlinear spectroscopy has revealed long-lasting oscillations in the optical
response of a variety of photosynthetic complexes. Different theoretical models
which involve the coherent coupling of electronic (excitonic) or
electronic-vibrational (vibronic) degrees of freedom have been put forward to
explain these observations. The ensuing debate concerning the relevance of
either one or the other mechanism may have obscured their potential synergy. To
illustrate this synergy, we quantify how the excitonic delocalization in the
LH2 unit of Rhodopseudomonas Acidophila purple bacterium, leads to correlations
of excitonic energy fluctuations, relevant coherent vibronic coupling and,
importantly, a decrease in the excitonic dephasing rates. Combining these
effects, we identify a feasible origin for the long-lasting oscillations
observed in fluorescent traces from time-delayed two-pulse single molecule
experiments performed on this photosynthetic complex.


  Experimental determination of protein function is resource-consuming. As an
alternative, computational prediction of protein function has received
attention. In this context, protein structural classification (PSC) can help,
by allowing for determining structural classes of currently unclassified
proteins based on their features, and then relying on the fact that proteins
with similar structures have similar functions. Existing PSC approaches rely on
sequence-based or direct 3-dimensional (3D) structure-based protein features.
In contrast, we first model 3D structures of proteins as protein structure
networks (PSNs). Then, we use network-based features for PSC. We propose the
use of graphlets, state-of-the-art features in many research areas of network
science, in the task of PSC. Moreover, because graphlets can deal only with
unweighted PSNs, and because accounting for edge weights when constructing PSNs
could improve PSC accuracy, we also propose a deep learning framework that
automatically learns network features from weighted PSNs. When evaluated on a
large set of ~9,400 CATH and ~12,800 SCOP protein domains (spanning 36 PSN
sets), our proposed approaches are superior to existing PSC approaches in terms
of accuracy, with comparable running time.


  We predict vaccine efficacy with a measure of antigenic distance between
influenza A(H3N2) and candidate vaccine viruses based on amino acid
substitutions in the dominant epitopes. In 2016-2017, our model predicts 19%
efficacy compared to 20% observed. This tool assists candidate vaccine
selection by predicting human protection against circulating strains.


  It is shown that quantum sustainability is a universal phenomenon which
emerges during environment-assisted electronic excitation energy transfer (EET)
in photobiological complexes (PBCs), such as photosynthetic reaction centers
and centers of melanogenesis. We demonstrate that quantum photobiological
systems must be sustainable for them to simultaneously endure continuous energy
transfer and keep their internal structure from destruction or critical
instability. These quantum effects occur due to the interaction of PBCs with
their environment which can be described by means of the reduced density
operator and effective non-Hermitian Hamiltonian (NH). Sustainable NH models of
EET predict the coherence beats, followed by the decrease of coherence down to
a small, yet non-zero value. This indicates that in sustainable PBCs, quantum
effects survive on a much larger time scale than the energy relaxation of an
exciton. We show that sustainable evolution significantly lowers the entropy of
PBCs and improves the speed and capacity of EET.


  DNA methylation is a well-studied genetic modification crucial to regulate
the functioning of the genome. Its alterations play an important role in
tumorigenesis and tumor-suppression. Thus, studying DNA methylation data may
help biomarker discovery in cancer. Since public data on DNA methylation become
abundant, and considering the high number of methylated sites (features)
present in the genome, it is important to have a method for efficiently
processing such large datasets. Relying on big data technologies, we propose
BIGBIOCL an algorithm that can apply supervised classification methods to
datasets with hundreds of thousands of features. It is designed for the
extraction of alternative and equivalent classification models through
iterative deletion of selected features. We run experiments on DNA methylation
datasets extracted from The Cancer Genome Atlas, focusing on three tumor types:
breast, kidney, and thyroid carcinomas. We perform classifications extracting
several methylated sites and their associated genes with accurate performance.
Results suggest that BIGBIOCL can perform hundreds of classification iterations
on hundreds of thousands of features in few hours. Moreover, we compare the
performance of our method with other state-of-the-art classifiers and with a
wide-spread DNA methylation analysis method based on network analysis. Finally,
we are able to efficiently compute multiple alternative classification models
and extract, from DNA-methylation large datasets, a set of candidate genes to
be further investigated to determine their active role in cancer. BIGBIOCL,
results of experiments, and a guide to carry on new experiments are freely
available on GitHub.


  A probabilistic reconstruction of genealogies in a polyploid population (from
2x to 4x) is investigated, by considering genetic data analyzed as the
probability of allele presence in a given genotype. Based on the likelihood of
all possible crossbreeding patterns, our model enables us to infer and to
quantify the whole potential genealogies in the population. We explain in
particular how to deal with the uncertain allelic multiplicity that may occur
with polyploids. Then we build an \textit{ad hoc} penalized likelihood to
compare genealogies and to decide whether a particular individual brings
sufficient information to be included in the taken genealogy. This decision
criterion enables us in a next part to suggest a greedy algorithm in order to
explore missing links and to rebuild some connections in the genealogies,
retrospectively. As a by-product, we also give a way to infer the individuals
that may have been favored by breeders over the years. In the last part we
highlight the results given by our model and our algorithm, firstly on a
simulated population and then on a real population of rose bushes. Most of the
methodology relies on the maximum likelihood principle and on graph theory.


  The aim of this paper is to tackle part of the program set by Diekmann et al.
in their seminal paper Diekmann et al. (2001). We quote "It remains to
investigate whether, and in what sense, the nonlinear determin-istic model
formulation is the limit of a stochastic model for initial population size
tending to infinity" We set a precise and general framework for a stochastic
individual based model : it is a piecewise deterministic Markov process defined
on the set of finite measures. We then establish a law of large numbers under
conditions easy to verify. Finally we show how this applies to old and new
examples.


  Molecular interactions have widely been modelled as networks. The local
wiring patterns around molecules in molecular networks are linked with their
biological functions. However, networks model only pairwise interactions
between molecules and cannot explicitly and directly capture the higher order
molecular organisation, such as protein complexes and pathways. Hence, we ask
if hypergraphs (hypernetworks), that directly capture entire complexes and
pathways along with protein-protein interactions (PPIs), carry additional
functional information beyond what can be uncovered from networks of pairwise
molecular interactions. The mathematical formalism of a hypergraph has long
been known, but not often used in studying molecular networks due to the lack
of sophisticated algorithms for mining the underlying biological information
hidden in the wiring patterns of molecular systems modelled as hypernetworks.
  We propose a new, multi-scale, protein interaction hypernetwork model that
utilizes hypergraphs to capture different scales of protein organization,
including PPIs, protein complexes and pathways. In analogy to graphlets, we
introduce hypergraphlets, small, connected, non-isomorphic, induced
sub-hypergraphs of a hypergraph, to quantify the local wiring patterns of these
multi-scale molecular hypergraphs and to mine them for new biological
information. We apply them to model the multi-scale protein networks of baker
yeast and human and show that the higher order molecular organisation captured
by these hypergraphs is strongly related to the underlying biology.
Importantly, we demonstrate that our new models and data mining tools reveal
different, but complementary biological information compared to classical PPI
networks. We apply our hypergraphlets to successfully predict biological
functions of uncharacterised proteins.


  Length-based methods are the cornerstone of many population studies and stock
assessments. This study tested two widely used methods: the Powell-Wetherall
(P-W) plot and the Lmax approach, i.e., estimating Linf directly from Lmax. In
most simulations, P-W estimates of the ratio total mortality / growth (Z/K
ratio) were biased beyond acceptable limits (bias > 30%). Bias in Z/K showed a
complex behavior, without possible corrections. Estimates of asymptotic length
(Linf) were less biased than Z/K, but were very sensible to intra-cohort
variability in growth and to changes in the occurrence of large individuals in
the sample. Exclusion of the largest size classes during the regression
procedure or weighing by abundance does not solve these issues. Perfect
linearization of the data and absurdly narrow confidence intervals for Z/K will
lead users to erroneous overconfidence in outputs. Clearly, the P-W method is
not suitable for the assessment of Z/K ratios of natural populations.
Estimation of Linf may be tentatively possible under very specific conditions,
with necessary external verifications. Also, this study demonstrates that there
is no way to estimate Linf directly from Lmax, since there is no particular
relationship to expect a priori between Linf and Lmax. Errors in estimating
Linf directly affect the estimate of the growth constant K and all other
subsequent calculations in population studies, stock assessments and ecosystem
models. New approaches are urgently needed for length-based studies of body
growth (e.g., unconstrained curve fit with subsequent bootstrapping), that
consider the inherent uncertainty regarding the underlying data and processes.


  Navigation on graphs is the problem how an agent walking on the graph can get
from a source to a target with limited information about the graph. The
information and the way to exploit it can vary. In this paper, we study
navigation on temporal networks -- networks where we have explicit information
about the time of the interaction, not only who interacts with whom. We
contrast a type of greedy navigation -- where agents follow paths that would
have worked well in the past -- with two strategies that do not exploit the
additional information. We test these on empirical temporal network data sets.
The greedy navigation is indeed more efficient than the reference strategies,
meaning that there are correlations in the real temporal networks that can be
exploited. We find that the navigability for individual nodes is most strongly
correlated with degree and burstiness, i.e., both topological and temporal
structures affect the navigation efficiency.


  Random polyampholytes (PAs) contain positively and negatively charged
monomers that are distributed randomly along the polymer chain. The interaction
between charges is assumed to be given by the Debye-Huckel potential. We show
that the size of the PA is determined by an interplay between electrostatic
interactions, giving rise to the polyelectrolyte (PE) effect due to net charge
per monomer ($\sigma$), and an effective attractive PA interaction due to
charge fluctuations, $\delta \sigma$. The interplay between these terms gives
rise to non-monotonic dependence of the radius of gyration, $R_g$ on the
inverse Debye length, $\kappa$ when PA effects are important (${\frac{\delta
\sigma}{\sigma}} > 1$). In the opposite limit, $R_g$ decreases monotonically
with increasing $\kappa$. Simulations of PA chains, using a charged bead-spring
model, further corroborates our theoretical predictions. The simulations
unambiguously show that conformational heterogeneity manifests itself among
sequences that have identical PA parameters. A clear implication is that the
phases of PA sequences, and by inference IDPs, cannot be determined using only
the bare PA parameters ($\sigma$ and $\delta \sigma$). The theory is used to
calculate the changes in $R_g$ on $N$, the number of residues for a set of
Intrinsically Disordered Proteins (IDPs). For a certain class of IDPs, with $N$
between 24 to 441, the size grows as $R_g \sim N^{0.6}$, which agrees with data
from Small Angle X-ray Scattering (SAXS) experiments.


  We consider a stochastic model of gene expression in which transcription
depends on a multistate promoter, including the famous two-state model and
refractory promoters as special cases, and focus on deriving the exact
stationary distribution. Building upon several successful approaches, we
present a more unified viewpoint that enables us to simplify and generalize
existing results. In particular, the original jump process is deeply related to
a multivariate piecewise-deterministic Markov process that may also be of
interest beyond the biological field. In a very particular case of promoter
configuration, this underlying process is shown to have a simple Dirichlet
stationary distribution. In the general case, the corresponding marginal
distributions extend the well-known class of Beta products, involving complex
parameters that directly relate to spectral properties of the promoter
transition matrix. Finally, we illustrate these results with biologically
plausible examples.