  Microbial communities display extreme diversity. A variety of strains or
species coexist even when limited by a single resource. It has been argued that
metabolite secretion creates new niches and facilitates such diversity.
Nonetheless, it is still a controversial topic why cells secrete even essential
metabolites so often; in fact, even under isolation conditions, microbial cells
secrete various metabolites, including those essential for their growth. First,
we demonstrate that leaking essential metabolites can be advantageous. If the
intracellular chemical reactions include multibody reactions like catalytic
reactions, this advantageous leakage of essential metabolites is possible and
indeed typical for most metabolic networks via "flux control" and
"growth-dilution" mechanisms; the later is a result of the balance between
synthesis and growth-induced dilution with autocatalytic reactions.
Counterintuitively, the mechanisms can work even when the supplied resource is
scarce. Next, when such cells are crowded, the presence of another cell type,
which consumes the leaked chemicals is beneficial for both cell types, so that
their coexistence enhances the growth of both. The latter part of the paper is
devoted to the analysis of such unusual form of symbiosis: "consumer" cell
types benefit from the uptake of metabolites secreted by "leaker" cell types,
and such consumption reduces the concentration of metabolites accumulated in
the environment; this environmental change enables further secretion from the
leaker cell types. This situation leads to frequency-dependent coexistence of
several cell types, as supported by extensive simulations. A new look at the
diversity in a microbial ecosystem is thus presented.


  Circadian clocks play a pivotal role in orchestrating numerous physiological
and developmental events. Waveform shapes of the oscillations of protein
abundances can be informative about the underlying biochemical processes of
circadian clocks. We derive a mathematical framework where waveforms do reveal
hidden biochemical mechanisms of circadian timekeeping. We find that the cost
of synthesizing proteins with particular waveforms can be substantially reduced
by rhythmic protein half-lives over time, as supported by previous plant and
mammalian data, as well as our own seedling experiment. We also find that
previously-enigmatic, cyclic expression of positive arm components within the
mammalian and insect clocks allows both a broad range of peak time differences
between protein waveforms and the symmetries of the waveforms about the peak
times. Such various peak-time differences may facilitate tissue-specific or
developmental stage-specific multicellular processes. Our waveform-guided
approach can be extended to various biological oscillators, including
cell-cycle and synthetic genetic oscillators.


  This paper proposes that the distinctively human capacity for cumulative,
adaptive, open-ended cultural evolution came about through two
temporally-distinct cognitive transitions. First, the origin of Homo-specific
culture over two MYA was made possible by the onset of a finer-grained
associative memory that allowed episodes to be encoded in greater detail. This
in turn meant more overlap amongst the distributed representations of these
episodes, such that they could more readily evoke one another through
self-triggered recall (STR). STR enabled representational redescription, the
chaining of thoughts and actions, and the capacity for a stream of thought.
Second, fully cognitive modernity following the appearance of anatomical
modernity after 200,000 BP, was made possible by the onset of contextual focus
(CF): the ability to shift between an explicit convergent mode conducive to
logic and refinement of ideas, and an implicit divergent mode conducive to
free-association, viewing situations from radically new perspectives, concept
combination, analogical thinking, and insight. This paved the way for an
integrated, creative internal network of understandings, and behavioral
modernity. We discuss feasible neural mechanisms for this two-stage proposal,
and outline how STR and CF differ from other proposals. We provide
computational evidence for the proposal obtained with an agent-based model of
cultural evolution in which agents invent ideas for actions and imitate the
fittest of their neighbors' actions. Mean fitness and diversity of actions
across the artificial society increased with STR, and even more so with CF, but
CF was only effective if STR was already in place. CF was most effective
following a change in task, which supports its hypothesized role in escaping
mental fixation. The proposal is discussed in the context of transition theory
in the life sciences.


  The potential benefits of applying machine learning methods to -omics data
are becoming increasingly apparent, especially in clinical settings. However,
the unique characteristics of these data are not always well suited to machine
learning techniques. These data are often generated across different
technologies in different labs, and frequently with high dimensionality. In
this paper we present a framework for combining -omics data sets, and for
handling high dimensional data, making -omics research more accessible to
machine learning applications. We demonstrate the success of this framework
through integration and analysis of multi-analyte data for a set of 3,533
breast cancers. We then use this data-set to predict breast cancer patient
survival for individuals at risk of an impending event, with higher accuracy
and lower variance than methods trained on individual data-sets. We hope that
our pipelines for data-set generation and transformation will open up -omics
data to machine learning researchers. We have made these freely available for
noncommercial use at www.ccg.ai.


  Humans are increasingly stressing ecosystems via habitat destruction, climate
change and global population movements leading to the widespread loss of
biodiversity and the disruption of key ecological services. Ecosystems
characterized primarily by mutualistic relationships between species such as
plant-pollinator interactions may be particularly vulnerable to such
perturbations because the loss of biodiversity can cause extinction cascades
that can compromise the entire network. Here, we develop a general restoration
strategy based on network-science for degraded ecosystems. Specifically, we
show that network topology can be used to identify the optimal sequence of
species reintroductions needed to maximize biodiversity gains following partial
and full ecosystem collapse. This restoration strategy generalizes across
topologically-disparate and geographically-distributed ecosystems.
Additionally, we find that although higher connectance and diversity promote
persistence in pristine ecosystems, these attributes reduce the effectiveness
of restoration efforts in degraded networks. Hence, focusing on restoring the
factors that promote persistence in pristine ecosystems may yield suboptimal
recovery strategies for degraded ecosystems. Overall, our results have
important insights for designing effective ecosystem restoration strategies to
preserve biodiversity and ensure the delivery of critical natural services that
fuel economic development, food security and human health around the globe


  Physical models of biological systems can become difficult to interpret when
they have a large number of parameters. But the models themselves actually
depend on (i.e. are sensitive to) only a subset of those parameters. Rigorously
identifying this subset of "stiff" parameters has been made possible by the
development of parameter space compression (PSC). However, PSC has only been
applied to analytically-solvable physical models. We have generalized this
powerful method by developing a numerical approach to PSC that can be applied
to any computational model. We validated our method against
analytically-solvable models of random walk with drift and protein production
and degradation. We then applied our method to an active area of biophysics
research, namely to a simple computational model of microtubule dynamic
instability. Such models have become increasingly complex, perhaps
unnecessarily. By adding two new parameters that account for prominent
structural features of microtubules, we identify one that can be "compressed
away" (the "seam" in the microtubule) and another that is essential to model
performance (the "tapering" of microtubule ends). Furthermore, we show that the
microtubule model has an underlying, low-dimensional structure that explains
the vast majority of our experimental data. We argue that numerical PSC can
identify the low-dimensional structure of any computational model in
biophysics. The low-dimensional structure of a model is easier to interpret and
identifies the mechanisms and experiments that best characterize the system.


  The collective behavior of biological oscillators has been recognized as an
important problem for several decades, but its control has come into limelight
only recently. Much of the focus for control has been on desynchronization of
an oscillator population, motivated by the pathological neural synchrony
present in essential and parkinsonian tremor. Other applications, such as the
beating of the heart and insulin secretion, require synchronization, and
recently there has been interest in forming clusters within an oscillator
population as well. In this article, we use a formulation that allows us to
devise control frameworks to achieve all of these distinct collective behaviors
observed in biological oscillators. This is based on the Fourier decomposition
of the partial differential equation governing the evolution of the phase
distribution of a population of identical, uncoupled oscillators. Our first two
control algorithms are Lyapunov-based, which work by decreasing a positive
definite Lyapunov function towards zero. Our third control is an optimal
control algorithm, which minimizes the control energy consumption while
achieving the desired collective behavior of an oscillator population.
Motivated by pathological neural synchrony, we apply our control to
desynchronize an initially synchronized neural population. Given the proposed
importance of enhancing spike time dependent plasticity to stabilize neural
clusters and counteract pathological neural synchronization, we formulate the
phase difference distribution in terms of the phase distribution, and prove
some of its fundamental properties, and in turn apply our control to transform
the neural phase distribution to form clusters. Finally, motivated by
eliminating cardiac alternans, we apply our control to phase shift a
synchronous cardiac pacemaker cell population.


  Current techniques in sequencing a genome allow a service provider (e.g. a
sequencing company) to have full access to the genome information, and thus the
privacy of individuals regarding their lifetime secret is violated. In this
paper, we introduce the problem of private DNA sequencing, where the goal is to
keep the DNA sequence private to the sequencer. We propose an architecture,
where the task of reading fragments of DNA and the task of DNA assembly are
separated, the former is done at the sequencer(s), and the later is completed
at a local trusted data collector. To satisfy the privacy constraint at the
sequencer and reconstruction condition at the data collector, we create an
information gap between these two relying on two techniques: (i) we use more
than one non-colluding sequencer, all reporting the read fragments to the
single data collector, (ii) adding the fragments of some known DNA molecules,
which are still unknown to the sequencers, to the pool. We prove that these two
techniques provide enough freedom to satisfy both conditions at the same time.


  Machine learning (ML)-guided directed evolution is a new paradigm for
biological design that enables optimization of complex functions. ML methods
use data to predict how sequence maps to function without requiring a detailed
model of the underlying physics or biological pathways. To demonstrate
ML-guided directed evolution, we introduce the steps required to build ML
sequence-function models and use them to guide engineering, making
recommendations at each stage. This review covers basic concepts relevant to
using ML for protein engineering as well as the current literature and
applications of this new engineering paradigm. ML methods accelerate directed
evolution by learning from information contained in all measured variants and
using that information to select sequences that are likely to be improved. We
then provide two case studies that demonstrate the ML-guided directed evolution
process. We also look to future opportunities where ML will enable discovery of
new protein functions and uncover the relationship between protein sequence and
function.


  Bacterial populations that colonize a host play important roles in host
health, including serving as a reservoir that transmits to other hosts and from
which invasive strains emerge, thus emphasizing the importance of understanding
rates of acquisition and clearance of colonizing populations. Studies of
colonization dynamics have been based on assessment of whether serial samples
represent a single population or distinct colonization events. A common
solution to estimate acquisition and clearance rates is to use a fixed genetic
distance threshold. However, this approach is often inadequate to account for
the diversity of the underlying within-host evolving population, the time
intervals between consecutive measurements, and the uncertainty in the
estimated acquisition and clearance rates. Here, we summarize recently
submitted work \cite{jarvenpaa2018named} and present a Bayesian model that
provides probabilities of whether two strains should be considered the same,
allowing to determine bacterial clearance and acquisition from genomes sampled
over time. We explicitly model the within-host variation using population
genetic simulation, and the inference is done by combining information from
multiple data sources by using a combination of Approximate Bayesian
Computation (ABC) and Markov Chain Monte Carlo (MCMC). We use the method to
analyse a collection of methicillin resistant Staphylococcus aureus (MRSA)
isolates.


  Bayesian Markov chain Monte Carlo explores tree space slowly, in part because
it frequently returns to the same tree topology. An alternative strategy would
be to explore tree space systematically, and never return to the same topology.
In this paper, we present an efficient parallelized method to map out the high
likelihood set of phylogenetic tree topologies via systematic search, which we
show to be a good approximation of the high posterior set of tree topologies.
Here `likelihood' of a topology refers to the tree likelihood for the
corresponding tree with optimized branch lengths. We call this method
`phylogenetic topographer' (PT). The PT strategy is very simple: starting in a
number of local topology maxima (obtained by hill-climbing from random starting
points), explore out using local topology rearrangements, only continuing
through topologies that are better than than some likelihood threshold below
the best observed topology. We show that the normalized topology likelihoods
are a useful proxy for the Bayesian posterior probability of those topologies.
By using a non-blocking hash table keyed on unique representations of tree
topologies, we avoid visiting topologies more than once across all concurrent
threads exploring tree space. We demonstrate that PT can be used directly to
approximate a Bayesian consensus tree topology. When combined with an accurate
means of evaluating per-topology marginal likelihoods, PT gives an alternative
procedure for obtaining Bayesian posterior distributions on phylogenetic tree
topologies.


  After irradiation of cancer cells in the ion beam therapy method the
concentration of hydrogen peroxide in the cell medium grows significantly. But
the role of hydrogen peroxide molecules in cancer treatment has not been
determined yet. We assume that interaction of peroxide molecules with DNA
atomic groups can block the genetic information of the cancer cell and lead to
its neutralization. To understand the possibility of DNA deactivation in the
cell, in the present study the formation of complexes of hydrogen peroxide
molecules with DNA specific recognition sites (nucleic bases) is considered.
Using atom-atom potential functions method and quantum-chemical approach, based
on density functional theory, the spatial configurations and energy minima for
the complexes of peroxide and water molecules with nucleic bases are studied.
The most probable positions of hydrogen peroxide molecules interacting with
nucleic bases are determined, and the possibility of blocking of genetic
information transfer processes is shown. The obtained data allows us to
formulate a new mechanism of the ion irradiation action on living cells, that
can be useful for cancer treatment.


  When modelling HIV epidemics, it is important to incorporate set-point viral
load and its heritability. As set-point viral load distributions can differ
significantly amongst epidemics, it is imperative to account for the observed
local variation. This can be done by using a heritability model and fitting it
to a local set-point viral load distribution. However, as the fitting procedure
needs to take into account the actual transmission dynamics (i.e., social
network, sexual behaviour), a complex model is required. Furthermore, in order
to use the estimates in subsequent modelling analyses to inform prevention
policies, it is important to assess parameter robustness.
  In order to fit set-point viral load models without the need to capture
explicitly the transmission dynamics, we present a new protocol. Firstly, we
approximate the transmission network from a phylogeny that was inferred from
sequences collected in the local epidemic. Secondly, as this transmission
network only comprises a single instance of the transmission network space, and
our aim is to assess parameter robustness, we infer the transmission network
distribution. Thirdly, we fit the parameters of the selected set-point viral
load model on multiple samples from the transmission network distribution using
approximate Bayesian inference.
  Our new protocol enables researchers to fit set-point viral load models in
their local context, and diagnose the model parameter's uncertainty. Such
parameter estimates are essential to enable subsequent modelling analyses, and
thus crucial to improve prevention policies.


  Antibiotic resistance monitoring is of paramount importance in the face of
this on-going global epidemic. Deep learning models trained with traditional
optimization algorithms (e.g. Adam, SGD) provide poor posterior estimates when
tested against out-of-distribution (OoD) antibiotic resistant/non-resistant
genes. In this paper, we introduce a deep learning model trained with
Stochastic Gradient Langevin Dynamics (SGLD) to classify antibiotic resistant
genes. The model provides better uncertainty estimates when tested against OoD
data compared to traditional optimization methods such as Adam.


  In this paper we develop and test algorithmic techniques to estimate
genotypes fitnesses by analysis of observed daily frequency data monitoring the
long-term evolution of bacterial populations. In particular, we develop a
non-linear least squares approach to estimate selective advantages of emerging
new mutant strains in locked-box stochastic models describing bacterial genetic
evolution similar to the celebrated Lenski experiment on Escherichia Coli. Our
algorithm first analyses emergence of new mutant strains for each individual
trajectory. For each trajectory our analysis is progressive in time, and
successively focuses on the first mutation event before analyzing the second
mutation event. The basic principle applied here is to minimize (for each
trajectory) the mean squared errors of prediction w(t) - W(t) where the
observed white cell frequencies w(t) are predicted by W(t), which is computed
as the conditional expectation of w(t) given the available information at time
(t-1). The pooling of all selective advantages estimates across all
trajectories provides histograms on which we perform a precise peak analysis to
compute final estimates of selective advantages. We validate our approach using
ensembles of simulated trajectories.


  The transportation of the cargoes in biological cells is primarily driven by
the motor proteins on filamentous protein tracks. The stochastic nature of the
motion of motor protein often leads to its spontaneous detachment from the
track. Using the available experimental data, we demonstrate a tradeoff between
the speed of the motor and its rate of spontaneous detachment from the track.
Further, it is also shown that this speed-detachment relation follows a power
law where its exponent dictates the nature of the motor protein processivity.
We utilize this information to study the motion of motor protein on track using
a random-walk model. We obtain the average distance travelled in fixed duration
and average time required for covering a given distance by the motor protein.
These analyses reveal non-monotonic dependence of the motor protein speed on
its transport and, therefore, optimal motor speeds can be identified for the
time and distance controlled conditions.


  The first-passage time (FPT), i.e., the moment when a stochastic process
reaches a given threshold value for the first time, is a fundamental
mathematical concept with immediate applications. In particular, it quantifies
the statistics of instances when biomolecules in a biological cell reach their
specific binding sites and trigger cellular regulation. Typically, the
first-passage properties are given in terms of mean first-passage times.
However, modern experiments now monitor single-molecular binding-processes in
living cells and thus provide access to the full statistics of the underlying
first-passage events, in particular, inherent cell-to-cell fluctuations. We
here present a robust explicit approach for obtaining the distribution of FPTs
to a small partially-reactive target in cylindrical-annulus domains, which
represent typical bacterial and neuronal cell shapes. We investigate various
asymptotic behaviours of this FPT distribution and show that it typically is
very broad in many biological situations: thus, the mean FPT can differ from
the most probable FPT by orders of magnitude. The most probable FPT is shown to
strongly depend only on the starting position within the geometry and to be
almost independent of the target size and reactivity. These findings
demonstrate the dramatic relevance of knowing the full distribution of FPTs and
thus open new perspectives for a more reliable description of many
intracellular processes initiated by the arrival of one or few biomolecules to
a small, spatially localised region inside the cell.


  DNA sequencing to identify genetic variants is becoming increasingly valuable
in clinical settings. Assessment of variants in such sequencing data is
commonly implemented through Bayesian heuristic algorithms. Machine learning
has shown great promise in improving on these variant calls, but the input for
these is still a standardized "pile-up" image, which is not always best suited.
In this paper, we present a novel method for generating images from DNA
sequencing data, which interlaces the human reference genome with personalized
sequencing output, to maximize usage of sequencing reads and improve machine
learning algorithm performance. We demonstrate the success of this in improving
standard germline variant calling. We also furthered this approach to include
somatic variant calling across tumor/normal data with Siamese networks. These
approaches can be used in machine learning applications on sequencing data with
the hope of improving clinical outcomes, and are freely available for
noncommercial use at www.ccg.ai.


  The modeling of atomistic biomolecular simulations using kinetic models such
as Markov state models (MSMs) has had many notable algorithmic advances in
recent years. The variational principle has opened the door for a nearly fully
automated toolkit for selecting models that predict the long-time kinetics from
molecular dynamics simulations. However, one yet-unoptimized step of the
pipeline involves choosing the features, or collective variables, from which
the model should be constructed. In order to build intuitive models, these
collective variables are often sought to be interpretable and familiar
features, such as torsional angles or contact distances in a protein structure.
However, previous approaches for evaluating the chosen features rely on
constructing a full MSM, which in turn requires additional hyperparameters to
be chosen, and hence leads to a computationally expensive framework. Here, we
present a method to optimize the feature choice directly, without requiring the
construction of the final kinetic model. We demonstrate our rigorous
preprocessing algorithm on a canonical set of twelve fast-folding protein
simulations, and show that our procedure leads to more efficient model
selection.


  The marginal likelihood of a model is a key quantity for assessing the
evidence provided by the data in support of a model. The marginal likelihood is
the normalizing constant for the posterior density, obtained by integrating the
product of the likelihood and the prior with respect to model parameters. Thus,
the computational burden of computing the marginal likelihood scales with the
dimension of the parameter space. In phylogenetics, where we work with tree
topologies that are high-dimensional models, standard approaches to computing
marginal likelihoods are very slow. Here we study methods to quickly compute
the marginal likelihood of a single fixed tree topology. We benchmark the speed
and accuracy of 19 different methods to compute the marginal likelihood of
phylogenetic topologies on a suite of real datasets. These methods include
several new ones that we develop explicitly to solve this problem, as well as
existing algorithms that we apply to phylogenetic models for the first time.
Altogether, our results show that the accuracy of these methods varies widely,
and that accuracy does not necessarily correlate with computational burden. Our
newly developed methods are orders of magnitude faster than standard
approaches, and in some cases, their accuracy rivals the best established
estimators.


  The process of RNA base fraying (i.e. the transient opening of the termini of
a helix) is involved in many aspects of RNA dynamics. We here use molecular
dynamics simulations and Markov state models to characterize the kinetics of
RNA fraying and its sequence and direction dependence. In particular, we first
introduce a method for determining biomolecular dynamics employing core-set
Markov state models constructed using an advanced clustering technique. The
method is validated on previously reported simulations. We then use the method
to analyze extensive trajectories for four different RNA model duplexes.
Results obtained using D. E. Shaw research and AMBER force fields are compared
and discussed in detail, and show a non-trivial interplay between the stability
of intermediate states and the overall fraying kinetics.


  Protein activation and deactivation is central to a variety of biological
mechanisms, including cellular signaling and transport. Unimolecular
fluorescent resonance energy transfer (FRET) probes are a class of fusion
protein sensors that allow biologists to visualize using an optical microscope
whether specific proteins are activated due to the presence nearby of small
drug-like signaling molecules, ligands or analytes. Often such probes comprise
a donor fluorescent protein attached to a ligand binding domain, a sensor or
reporter domain attached to the acceptor fluorescent protein, with these ligand
binding and sensor domains connected by a protein linker. Various choices of
linker type are possible ranging from highly flexible proteins to hinge-like
proteins. It is also possible to select donor and acceptor pairs according to
their corresponding F\"oster radius, or even to mutate binding and sensor
domains so as to change their binding energy in the activated or inactivated
states. The focus of the present work is the exploration through simulation of
the impact of such choices on sensor performance.


  The inapplicability of amino acid covariation methods to small protein
families has limited their use for structural annotation of whole genomes.
Recently, deep learning has shown promise in allowing accurate residue-residue
contact prediction even for shallow sequence alignments. Here we introduce
DMPfold, which uses deep learning to predict inter-atomic distance bounds, the
main chain hydrogen bond network, and torsion angles, which it uses to build
models in an iterative fashion. DMPfold produces more accurate models than two
popular methods for a test set of CASP12 domains, and works just as well for
transmembrane proteins. Applied to all Pfam domains without known structures,
confident models for 25% of these so-called dark families were produced in
under a week on a small 200 core cluster. DMPfold provides models for 16% of
human proteome UniProt entries without structures, generates accurate models
with fewer than 100 sequences in some cases, and is freely available.


  The effect of ethanol-water media on the synthesis of Ag/TiO2 nanocomposite
was investigated with 0.05, 0.1 and 0.5 (wt.%) of Ag content. Ethanol was used
as hole-scavenger enhancing the photodecomposition of Ag+ ions under
illumination of near-UV light. The nanocomposites were further calcined to 33C
and 400C under controlled atmosphere. The synthesized nanocomposites were
tested for photocatalytic efficiency decomposing formic acid as an organic
pollutant under irradiation of a 4W near UV lamp with lmax at 365nm and the
results were compared with that of nan-modified commercially available Degussa
TiO2 (P25). The synthesized nanocomposites were characterized using XRD and SEM
micrographs. The photocatalytic studies have revealed a constant overall
performance for all synthesized nanocomposites. The initial rate of
decomposition was observed to increase with decreasing Ag content, in the case
of 0.05 wt.% having the best initial rate. Calcination of the nanocomposites
was found to have activating properties on the synthesized nanocomposites where
nanocomposite with 0.05 (wt.%) Ag content calcined at 300{\deg}C was
distinctively advantageous over other calcined nanocomposites.


  Intratumor heterogeneity (ITH), referring to coexistence of different cell
subpopulations in a single tumor, has been a major puzzle in cancer research
for almost half a century. The lack of understanding of the underlying
mechanism of ITH hinders progress in developing effective therapies for
cancers. Based on the findings in a recent quantitative experiment on
pancreatic cancer, we developed a general evolutionary model for one type of
cancer, accounting for interactions between different cell populations through
paracrine or juxtacrine factors. We show that the emergence of a stable
heterogeneous state in a tumor requires an unequal allocation of paracrine
growth factors ("public goods") between cells that produce them and those that
merely consume them. Our model provides a quantitative explanation of recent
{\it in vitro} experimental studies in pancreatic cancer in which insulin
growth factor (IGF-II) plays the role of public goods. The calculated phase
diagrams as a function of exogenous resources and fraction of growth factor
producing cells show ITH persists only in a narrow range of concentration of
exogenous IGF-II. Remarkably, maintenance of ITH requires cooperation among
tumor cell subpopulations in harsh conditions, specified by lack of exogenous
IGF-II, whereas surplus exogenous IGF-II elicits competition. Our theory also
quantitatively accounts for measured {\it in vivo} tumor growth in glioblastoma
multiforme (GBM). The predictions for GBM tumor growth as a function of the
fraction of tumor cells are amenable to experimental tests. The mechanism for
ITH also provides hints for devising efficacious therapies.


  A popular approach to analyze the dynamics of high-dimensional many-body
systems, such as macromolecules, is to project the trajectories onto a space of
slowly-varying collective variables, where subsequent analyses are made, such
as clustering or estimation of free energy profiles or Markov state models
(MSMs). However, existing "dynamical" dimension reduction methods, such as the
time-lagged independent component analysis (TICA) are only valid if the
dynamics obeys detailed balance (microscopic reversibility) and typically
require long, equilibrated simulation trajectories. Here we develop a dimension
reduction method for non-equilibrium dynamics based on the recently developed
Variational Approach for Markov Processes (VAMP) by Wu and No\'e. VAMP is
illustrated by obtaining a low-dimensional description of a single file ion
diffusion model and by identifying long-lived states from molecular dynamics
simulations of the KcsA channel protein in an external electrochemical
potential. This analysis provides detailed insights into the coupling of
conformational dynamics, the configuration of the selectivity filter, and the
conductance of the channel. We recommend VAMP as a replacement for the less
general TICA method.


  Biomolecular networks have already found great utility in characterizing
complex biological systems arising from pair-wise interactions amongst
biomolecules. Here, we review how graph theoretical approaches can be applied
not only for a better understanding of various proximate (mechanistic)
relations, but also, ultimate (evolutionary) structures encoded in such
networks. A central question deals with the evolutionary dynamics by which
different topologies of biomolecular networks might have evolved, as well as
the biological principles that can be hypothesized from a deeper understanding
of the induced network dynamics. We emphasize the role of gene duplication in
terms of signaling game theory, whereby sender and receiver gene players accrue
benefit from gene duplication, leading to a preferential attachment mode of
network growth. Information asymmetry between sender and receiver genes is
hypothesized as a key driver of the resulting network topology. The study of
the resulting dynamics suggests many mathematical/computational problems, the
majority of which are intractable but yield to efficient approximation
algorithms, when studied through an algebraic graph theoretic lens.


  Nowadays, evidence is mounting that the race of living organisms for
adaptation to the chemicals synthesized by their neighbours may drive community
structures. Particularly, some bacterial infections and plant invasions
disruptive of the native community rely on the release of allelochemicals that
inhibit or kill sensitive strains or individuals from their own or other
species. In this report, an eco-evolutionary model for community assembly
through resource competition, allelophatic interactions, and evolutionary
branching is presented and studied by numerical analysis. Our major findings
are that stable communities with increasing biodiversity can emerge at weak
allelopathic suppression, but stronger allelophaty is negatively correlated
with community diversity. In the former regime the allelopathic interaction
networks exhibit Gaussian degree distributions, while in the later one the
network degrees are Weibull distributed.


  We investigate the effect of filament-obstacle interactions on the
force-velocity relation of growing biopolymers, via calculations explicitly
treating obstacle diffusion and stochastic addition and subtraction of
subunits. We first show that the subunit on- and off-rates satisfy a rigorous
thermodynamic relationship determined by the filament-obstacle interaction
potential. Both the on- and off-rates depend not only on the average force on
the obstacle, but also on the shape of the potential on the nanometer length
scale. Basing obstacle-induced reduction of the on-rate entirely on the force
overestimates the stall force when there are fluctuations in the force exerted
on a filament tip. We then perform simulations and analytic calculations using
the thermodynamic relationship. We find, consistent with expectations from
general thermodynamic relations, that the "Brownian-Ratchet" model is an upper
bound to the growth velocity and that for purely repulsive potentials the
growth velocity is essentially that predicted by the Brownian-Ratchet model.
For shallow potential wells of depth ~ 5k_BT, which might correspond to
transient filament-membrane attachments, the zero-force velocity is a
substantial fraction of the free-filament velocity. In this case, the growth
velocity can depend strongly on the obstacle diffusion coefficient even when
the dimensionless diffusion coefficient is large. The velocity also drops more
rapidly than predicted by the Brownian-ratchet model, in some cases by as much
as a factor of 50 at an opposing force of 1 pN. For deep potential wells, as
might result from strong filament-membrane links, both the on- and off-rates
are reduced significantly, slowing polymerization. Such potentials can sustain
pulling forces while polymerizing, but only if the attractive well has a
"shelf" comparable to or greater than the monomer size.


  The Extended Evolutionary Synthesis (EES) is beginning to fulfill the whole
promise of Darwinian insight through its extension of evolutionary
understanding from the biological domain to include cultural information
evolution. Several decades of important foundation-laying work took a social
Darwinist approach and exhibited and ecologically-deterministic elements. This
is not the case with more recent developments to the evolutionary study of
culture, which emphasize non-Darwinian processes such as self-organization,
potentiality, and epigenetic change.


  Diagnosing an inherited disease often requires identifying the pattern of
inheritance in a patient's family. We represent family trees with genetic
patterns of inheritance using hypergraphs and latent state space models to
provide explainable inheritance pattern predictions. Our approach allows for
exact causal inference over a patient's possible genotypes given their
relatives' phenotypes. By design, inference can be examined at a low level to
provide explainable predictions. Furthermore, we make use of human intuition by
providing a method to assign hypothetical evidence to any inherited gene
alleles. Our analysis supports the application of latent state space models to
improve patient care in cases of rare inherited diseases where access to
genetic specialists is limited.


  For precision medicine and personalized treatment, we need to identify
predictive markers of disease. We focus on Alzheimer's disease (AD), where
magnetic resonance imaging scans provide information about the disease status.
By combining imaging with genome sequencing, we aim at identifying rare genetic
markers associated with quantitative traits predicted from convolutional neural
networks (CNNs), which traditionally have been derived manually by experts.
Kernel-based tests are a powerful tool for associating sets of genetic
variants, but how to optimally model rare genetic variants is still an open
research question. We propose a generalized set of kernels that incorporate
prior information from various annotations and multi-omics data. In the
analysis of data from the Alzheimer's Disease Neuroimaging Initiative (ADNI),
we evaluate whether (i) CNNs yield precise and reliable brain traits, and (ii)
the novel kernel-based tests can help to identify loci associated with AD. The
results indicate that CNNs provide a fast, scalable and precise tool to derive
quantitative AD traits and that new kernels integrating domain knowledge can
yield higher power in association tests of very rare variants.


  Modeling disease progression in healthcare administrative databases is
complicated by the fact that patients are observed only at irregular intervals
when they seek healthcare services. In a longitudinal cohort of 76,888 patients
with chronic obstructive pulmonary disease (COPD), we used a continuous-time
hidden Markov model with a generalized linear model to model healthcare
utilization events. We found that the fitted model provides interpretable
results suitable for summarization and hypothesis generation.


  There are many models of the Calvin cycle of photosynthesis in the
literature. When investigating the dynamics of these models one strategy is to
look at the simplest possible models in order to get the most detailed
insights. We investigate a minimal model of the Calvin cycle introduced by Hahn
while he was pursuing this strategy. In a variant of the model not including
photorespiration it is shown that there exists exactly one positive steady
state and that this steady state is unstable. For generic initial data either
all concentrations tend to infinity at lates times or all concentrations tend
to zero at late times. In a variant including photorespiration it is shown that
for suitable values of the parameters of the model there exist two positive
steady states, one stable and one unstable. For generic initial data either the
solution tends to the stable steady state at late times or all concentrations
tend to zero at late times. Thus we obtain rigorous proofs of mathematical
statements which together confirm the intuitive idea proposed by Hahn that
photorespiration can stabilize the operation of the Calvin cycle. In the case
that the concentrations tend to infinity we derive formulae for the leading
order asymptotics using the Poincar\'e compactification.


  The covalent functionalisation of type I atelocollagen with either
4-vinylbenzyl or methacrylamide residues is presented as a simple synthetic
strategy to achieve customisable, cell-friendly UV-cured hydrogel networks with
widespread clinical applicability. Molecular parameters, i.e. the type of
monomer, degree of atelocollagen functionalisation and UV curing solution, have
been systematically varied and their effect on gelation kinetics, swelling
behaviour, elastic properties and enzymatic degradability investigated.
UV-cured hydrogel networks deriving from atelocollagen precursors
functionalised with equivalent molar content of 4-vinylbenzyl and
methacrylamide adducts proved to display remarkably-different swelling ratio,
storage modulus and collagenase resistance, similarly to the case of UV-cured
hydrogel networks obtained with the same type of methacrylamide adduct, but
varied degree of functionalisation. UV-induced network formation of 4VBC
functionalised atelocollagen molecules yielded hydrogels with increased
stiffness and enzymatic stability, attributed to the molecular rigidity of
resulting aromatised crosslinking segment, whilst no toxic response was
observed with osteosarcoma G292 cells. Although to a lesser extent, the pH of
the UV-curing solution also proved to affect macroscopic hydrogel properties,
likely due to the altered organisation of atelocollagen molecules during
network formation. By leveraging the knowledge gained with classic synthetic
networks, this study highlights how the type of monomer can be conveniently
exploited to realise customisable atelocollagen hydrogels with controlled
structure-property relationships to meet the requirements of unmet clinical
applications.


  Collective behaviour is known to be the result of diverse dynamics and is
sometimes likened to a living system. Although many studies have revealed the
dynamics of various collective behaviours, their main focus was on the
information process inside the collective, not on the whole system itself. For
example, the qualitative difference between two elements and three elements as
a system has rarely been investigated. Tononi et al. have proposed Integrated
Information Theory (IIT) to measure the degree of consciousness $\Phi$. IIT
postulates that the amount of information loss caused by certain partitions is
equivalent to the degree of information integration in the system. This measure
is not only useful for estimating the degree of consciousness but can also be
applied to more general network systems. Here we applied IIT (in particular,
IIT 3.0 using PyPhi) to analyse real fish schools ({\it Plecoglossus
altivelis}). Our hypothesis in this study is a very simple one: a living system
evolves to raise its $\Phi$ value. If we accept this hypothesis, IIT reveals
the existence of continuous and discontinuous properties as group size varies.
For example, leadership in the fish school emerged for a school size of four or
above; but not below three. Furthermore, this transition was not observed by
measuring mutual information or in a simple Boids model. This result suggests
that integrated information $\Phi$ can reveal some inherent properties which
cannot be observed using other measures. We also discuss how the fish
recognition of the figure-ground relation, that is, what determines the
relevant ON and OFF states, may reveal various optimal paths for obtaining the
functional evolution of collective behaviour.


  We present two globally convergent Levenberg-Marquardt methods for finding
zeros of H\"{o}lder metrically subregular mappings that may have non-isolated
zeros. The first method unifies the Levenberg- Marquardt direction and an
Armijo-type line search, while the second incorporates this direction with a
nonmonotone trust-region technique. For both methods, we prove the global
convergence to a first-order stationary point of the associated merit function.
Furthermore, the worst-case global complexity of these methods are provided,
indicating that an approximate stationary point can be computed in at most
$\mathcal{O}(\varepsilon^{-2})$ function and gradient evaluations, for an
accuracy parameter $\varepsilon>0$. We also study the conditions for the
proposed methods to converge to a zero of the associated mappings. Computing a
moiety conserved steady state for biochemical reaction networks can be cast as
the problem of finding a zero of a H\"{o}lder metrically subregular mapping. We
report encouraging numerical results for finding a zero of such mappings
derived from real-world biological data, which supports our theoretical
foundations.


  We consider the mutation--selection differential equation with pairwise
interaction (or, equivalently, the diploid mutation--selection equation) and
establish the corresponding ancestral process, which is a random tree and a
variant of the ancestral selection graph. The formal relation to the forward
model is given via duality. To make the tree tractable, we prune branches upon
mutations, thus reducing it to its informative parts. The hierarchies inherent
in the tree are encoded systematically via tripod trees with weighted leaves;
this leads to the stratified ancestral selection graph. The latter also
satisfies a duality relation with the mutation--selection equation. Each of the
dualities provides a stochastic representation of the solution of the
differential equation. This allows us to connect the equilibria and their
bifurcations to the long-term behaviour of the ancestral process. Furthermore,
with the help of the stratified ancestral selection graph, we obtain explicit
results about the ancestral type distribution in the case of unidirectional
mutation.


  De novo protein structure prediction from amino acid sequence is one of the
most challenging problems in computational biology. As one of the extensively
explored mathematical models for protein folding, Hydrophobic-Polar (HP) model
enables thorough investigation of protein structure formation and evolution.
Although HP model discretizes the conformational space and simplifies the
folding energy function, it has been proven to be an NP-complete problem. In
this paper, we propose a novel protein folding framework FoldingZero,
self-folding a de novo protein 2D HP structure from scratch based on deep
reinforcement learning. FoldingZero features the coupled approach of a two-head
(policy and value heads) deep convolutional neural network (HPNet) and a
regularized Upper Confidence Bounds for Trees (R-UCT). It is trained solely by
a reinforcement learning algorithm, which improves HPNet and R-UCT iteratively
through iterative policy optimization. Without any supervision and domain
knowledge, FoldingZero not only achieves comparable results, but also learns
the latent folding knowledge to stabilize the structure. Without exponential
computation, FoldingZero shows promising potential to be adopted for real-world
protein properties prediction.


  We introduce a framework for analyzing ordinary differential equation (ODE)
models of biological networks using statistical model checking (SMC). A key
aspect of our work is the modeling of single-cell variability by assigning a
probability distribution to intervals of initial concentration values and
kinetic rate constants. We propagate this distribution through the system
dynamics to obtain a distribution over the set of trajectories of the ODEs.
This in turn opens the door for performing statistical analysis of the ODE
system's behavior. To illustrate this we first encode quantitative data and
qualitative trends as bounded linear time temporal logic (BLTL) formulas. Based
on this we construct a parameter estimation method using an SMC-driven
evaluation procedure applied to the stochastic version of the behavior of the
ODE system. We then describe how this SMC framework can be generalized to
hybrid automata by exploiting the given distribution over the initial states
and the--much more sophisticated--system dynamics to associate a Markov chain
with the hybrid automaton. We then establish a strong relationship between the
behaviors of the hybrid automaton and its associated Markov chain.
Consequently, we sample trajectories from the hybrid automaton in a way that
mimics the sampling of the trajectories of the Markov chain. This enables us to
verify approximately that the Markov chain meets a BLTL specification with high
probability. We have applied these methods to ODE based models of Toll-like
receptor signaling and the crosstalk between autophagy and apoptosis, as well
as to systems exhibiting hybrid dynamics including the circadian clock pathway
and cardiac cell physiology. We present an overview of these applications and
summarize the main empirical results. These case studies demonstrate that the
our methods can be applied in a variety of practical settings.


  Predicting the structure of a protein from its sequence is a cornerstone task
of molecular biology. Established methods in the field, such as homology
modeling and fragment assembly, appeared to have reached their limit. However,
this year saw the emergence of promising new approaches: end-to-end protein
structure and dynamics models, as well as reinforcement learning applied to
protein folding. For these approaches to be investigated on a larger scale, an
efficient implementation of their key computational primitives is required. In
this paper we present a library of differentiable mappings from two standard
dihedral-angle representations of protein structure (full-atom representation
"$\phi,\psi,\omega,\chi$" and backbone-only representation
"$\phi,\psi,\omega$") to atomic Cartesian coordinates. The source code and
documentation can be found at https://github.com/lupoglaz/TorchProteinLibrary.


  Genetic regulatory networks (GRNs) have been widely studied, yet there is a
lack of understanding with regards to the final size and properties of these
networks, mainly due to no network currently being complete. In this study, we
analyzed the distribution of GRN structural properties across a large set of
distinct prokaryotic organisms and found a set of constrained characteristics
such as network density and number of regulators. Our results allowed us to
estimate the number of interactions that complete networks would have, a
valuable insight that could aid in the daunting task of network curation,
prediction, and validation. Using state-of-the-art statistical approaches, we
also provided new evidence to settle a previously stated controversy that
raised the possibility of complete biological networks being random and
therefore attributing the observed scale-free properties to an artifact
emerging from the sampling process during network discovery. Furthermore, we
identified a set of properties that enabled us to assess the consistency of the
connectivity distribution for various GRNs against different alternative
statistical distributions. Our results favor the hypothesis that highly
connected nodes (hubs) are not a consequence of network incompleteness.
Finally, an interaction coverage computed for the GRNs as a proxy for
completeness revealed that high-throughput based reconstructions of GRNs could
yield biased networks with a low average clustering coefficient, showing that
classical targeted discovery of interactions is still needed.


  In this paper, we propose a sex-structured entomological model that serves as
a basis for design of control strategies relying on releases of sterile male
mosquitoes (Aedes spp) and aiming at elimination of the wild vector population
in some target locality. We consider different types of releases (constant and
periodic impulsive), providing necessary conditions to reach elimination.
However, the main part of the paper is focused on the study of the periodic
impulsive control in different situations. When the size of wild mosquito
population cannot be assessed in real time, we propose the so-called open-loop
control strategy that relies on periodic impulsive releases of sterile males
with constant release size. Under this control mode, global convergence towards
the mosquito-free equilibrium is proved on the grounds of sufficient condition
that relates the size and frequency of releases. If periodic assessments
(either synchronized with releases or more sparse) of the wild population size
are available in real time, we propose the so-called closed-loop control
strategy, which is adjustable in accordance with reliable estimations of the
wild population sizes. Under this control mode, global convergence to the
mosquito-free equilibrium is proved on the grounds of another sufficient
condition that relates not only the size and frequency of periodic releases but
also the frequency of sparse measurements taken on wild populations. Finally,
we propose a mixed control strategy that combines open-loop and closed-loop
strategies. This control mode renders the best result, in terms of overall time
needed to reach elimination and the number of releases to be effectively
carried out during the whole release campaign, while requiring for a reasonable
amount of released sterile insects.


  We present a new Markov chain Monte Carlo algorithm, implemented in software
Arbores, for inferring the history of a sample of DNA sequences. Our principal
innovation is a bridging procedure, previously applied only for simple
stochastic processes, in which the local computations within a bridge can
proceed independently of the rest of the DNA sequence, facilitating large-scale
parallelisation.


  In this article, we show how the recent statistical techniques developed in
Topological Data Analysis for the Mapper algorithm can be extended and
leveraged to formally define and statistically quantify the presence of
topological structures coming from biological phenomena in datasets of CCC
contact maps.


  One of the major issues in the theoretical modeling of epidemic spreading is
the development of methods to control the transmission of an infectious agent.
Human behavior plays a fundamental role in the spreading dynamics and can be
used to stop a disease from spreading or to reduce its burden, as individuals
aware of the presence of a disease can take measures to reduce their exposure
to contagion. In this paper, we propose a mathematical model for the spread of
diseases with awareness in complex networks. Unlike previous models, the
information is propagated following a generalized Maki-Thompson rumor model.
Flexibility on the timescale between information and disease spreading is also
included. We verify that the velocity characterizing the diffusion of
information awareness greatly influences the disease prevalence. We also show
that a reduction in the fraction of unaware individuals does not always imply a
decrease of the prevalence, as the relative timescale between disease and
awareness spreading plays a crucial role in the systems' dynamics. This result
is shown to be independent of the network topology. We finally calculate the
epidemic threshold of our model, and show that it does not depend on the
relative timescale. Our results provide a new view on how information influence
disease spreading and can be used for the development of more efficient methods
for disease control.


  Proteins display generic properties that are challenging to explain by direct
selection, notably allostery, the capacity to be regulated through long-range
effects, and evolvability, the capacity to adapt to new selective pressures. An
evolutionary scenario is proposed where proteins acquire these two features
indirectly as a by-product of their selection for a more fundamental property,
exquisite discrimination, the capacity to bind discriminatively very similar
ligands. Achieving this task is shown to typically require proteins to undergo
a conformational change. We argue that physical and evolutionary constraints
impel this change to be controlled by a group of sites extending from the
binding site. Proteins can thus acquire a latent potential for allosteric
regulation and evolutionary adaptation because of long-range effects that
initially arise as evolutionary spandrels. This scenario accounts for the
groups of conserved and coevolving residues observed in multiple sequence
alignments. However, we propose that most pairs of coevolving and contacting
residues inferred from such alignments have a different origin, related to
thermal stability. A physical model is presented that illustrates this
evolutionary scenario and its implications. The scenario can be implemented in
experiments of protein evolution to directly test its predictions.


  A recent study in the 2017 Living to 100 Monograph published by the Society
of Actuaries suggests, in contrast to previous research, that the risk of death
after 110 increases with age. By fitting a Gompertz model to estimated central
death rates for the oldest old, the authors challenge existing theory and
empirical research indicating a deceleration of mortality at older ages and the
emergence of a plateau. We argue that their results are inconclusive for three
reasons: (1) the data selection was arbitrary; (2) the statistical analysis was
inappropriate; and (3) the presentation of the results is misleading and
inadequate. We therefore claim that the hypothesis that the human force of
mortality increases after age 110 has not been proved.


  In this work we prove occurrence of a super-critical Hopf bifurcation in a
model of white blood cell formation structured by three maturation stages. We
provide an explicit analytical expression for the bifurcation point depending
on model parameters. The Hopf bifurcation is a unique feature of the
multi-compartment structure as it does not exist in the corresponding
two-compartment model. It appears for a parameter set different from the
parameters identified for healthy hematopoiesis and requires changes in at
least two cell properties. Model analysis allows identifying a range of
biologically plausible parameter sets that can explain persistent oscillations
of white blood cell counts observed in some hematopoietic diseases. Relating
the identified parameter sets to recent experimental and clinical findings
provides insights into the pathological mechanisms leading to oscillating blood
cell counts.


  Simulating the population dynamics of a stage-structured population requires
the knowledge of the biodemographic functions characterizing the species,
namely development, mortality and fecundity. In general, development and
fecundity can satisfactorily be estimated starting from literature data.
Unfortunately, this is often not the case of the mortality function, because of
the lack of experimental data. To overcome this problem we estimate the
mortality rate function from field data on the abundance of the species. The
mortality is expressed as a linear combination of cubic splines and the
estimation method allows to determine its coefficients taking into account the
observations measurement error. Moreover, the variability in the estimate is
quantified by means of the confidence bands for both mortality and dynamics.
The presented method allows to obtain a more flexible shape for the mortality
rate functions compared with previous methods applied to the same pest. The
method is applied to the case of Lobesia botrana, the main pest in the European
vineyards, with abundance data collected in a location in the North of Italy
for five consecutive years. Data collected over three years are used to
estimate the mortality and to analyze the variability in the estimate and its
effects on the population dynamics. Other two datasets are used to validate the
model simulating the dynamics using the estimated mortality.


  Pest phenological models describe the cumulative flux of the individuals into
each stage of the life cycle of a stage-structured population. Phenological
models are widely used tools in pest control decision making. Despite the fact
that these models do not provide information on population abundance, they
share some advantages with respect to the more sophisticated and complex
demographic models. The main advantage is that they do not require data
collection to define the initial conditions of model simulation, reducing the
effort for field sampling and the high uncertainty affecting sample estimates.
Phenological models are often built considering the developmental rate function
only. To the aim of adding more realism to phenological models, in this paper
we explore the consequences of improving these models taking into consideration
three additional elements: the age distribution of individuals which exit from
the overwintering phase, the age- and temperature-dependent profile of the
fecundity rate function and the consideration of a temperature-dependent
mortality rate function. Numerical simulations are performed to investigate the
effect of these elements with respect to phenological models considering
development rate functions only. To further test the implications of different
models formulation, we compare results obtained from different phenological
models to the case study of the codling moth (Cydia pomonella) a primary pest
of the apple orchard. The results obtained from model comparison are discussed
in view of their potential application in pest control decision support.


  Cooperation is prevalent in nature, not only in the context of social
interactions within the animal kingdom, but also on the cellular level. In
cancer for example, tumour cells can cooperate by producing growth factors. The
evolution of cooperation has traditionally been studied for well-mixed
populations under the framework of evolutionary game theory, and more recently
for structured populations using evolutionary graph theory. The population
structures arising due to cellular arrangement in tissues however are dynamic
and thus cannot be accurately represented by either of these frameworks. In
this work we compare the conditions for cooperative success in an epithelium
modelled using evolutionary graph theory, to those in a mechanical model of an
epithelium =- the Voronoi tessellation model. Crucially, in this latter model
cells are able to move, and birth and death are not spatially coupled. We
calculate fixation probabilities in the Voronoi tessellation model through
simulation and an approximate analytic technique and show that this leads to
stronger promotion of cooperation in comparison with the evolutionary graph
theory model.


  Classical understanding of the outcome of the struggle for existence results
in the Darwinian survival of the fittest. Here we show that the situation may
be different, more complex and arguably more interesting. Specifically, we show
that different versions of non-homogeneous logistic-like models with a
distributed Malthusian parameter imply non-Darwinian survival of everybody. In
contrast, the non-homogeneous logistic equation with distributed carrying
capacity shows Darwinian survival of the fittest. We also consider an
non-homogeneous birth-and-death equation and give a simple proof that this
equation results in the survival of the fittest. In addition to this known
result, we find an exact limit distribution of the parameters of this equation.
We also consider frequency-dependent non-homogeneous models and show that
although some of these models show Darwinian survival of the fittest, there is
not enough time for selection of the fittest species. We discuss the well-known
Gauze Competitive exclusion principle that states that Complete competitors
cannot coexist. While this principle is often considered as a direct
consequence of the Darwinian survival of the fittest, we show that from the
point of view of developed mathematical theory complete competitors can in fact
coexist indefinitely.


  Numerous studies have demonstrated the important role of noise in the
dynamical behaviour of a complex system. The most probable trajectories of
nonlinear systems under the influence of Gaussian noise have recently been
studied already. However, there has been only a few works that examine how most
probable trajectories in the two-dimensional system (MeKS network) are
influenced under non-Gaussian stable L\'evy noise. Therefore, we discuss the
most probable trajectories of a two-dimensional model depicting the competence
behaviour in B. subtilis under the influence of stable L\'evy noise. On the
basis of the Fokker-Planck equation, we describe the noise-induced most
probable trajectories of the MeKS network from the low ComK protein
concentration (vegetative state) to the high ComK protein concentration
(competence state) under stable L\'evy noise. We demonstrate choices of the
non-Gaussianity index $\alpha$ and the noise intensity $\epsilon$ which
generate the ComK protein escape from the low concentration to the high
concentration. We also reveal the optimal combination of both parameters
$\alpha$ and $\epsilon$ making the tipping time shortest. Moreover, we find
that different initial concentrations around the low ComK protein concentration
evolve to a metastable state, and provide the optimal $\alpha$ and $\epsilon$
such that the distance between the deterministic competence state and the
metastable state is smallest.


  Cyclic dominant systems, like rock-paper-scissors game, are frequently used
to explain biodiversity in nature, where mobility, reproduction and
intransitive competition are on stage to provide the coexistence of
competitors. A significantly new situation emerges if we introduce an apex
predator who can superior all members of the mentioned three-species system. In
the latter case the evolution may terminate into three qualitatively different
destinations depending on the apex predator decaying ratio $q$. In particular,
the whole population goes extinct or all four species survive or only the
original three-species system remains alive as we vary the control parameter.
These solutions are separated by a discontinuous and a continuous phase
transitions at critical $q$ values. Our results highlight that cyclic dominant
competition can offer a stable way to survive even in a predator-prey-like
system that can be maintained for large interval of critical parameter values.


  Mutualistic interactions are vital constituents of ecological and
socio-economic systems. Empirical studies have found that the patterns of
reciprocal relations among the participants often shows the salient features of
being simultaneously nested and modular. Whether and how these two structural
properties of mutualistic networks can emerge out of a common mechanism however
remains unclear. We propose a unified dynamic model based on the adaptation of
niche relations that gives rise to both structural features. We apply
Hutchinson's concept of niche interaction to networked cooperative species.
Their niche relation evolves under the assumption of fitness maximization.
Modularity and nestedness emerge concurrently through the accumulated local
advantages in the structural and demographic distribution. A rich ensemble of
key dynamical behaviors are unveiled in the dynamical framework. We demonstrate
that mutualism can exhibit either a stabilizing or destabilizing effect on the
evolved network, which undergoes a drastic transition with the overall
competition level. Most strikingly, the adaptive network may exhibit a profound
nature of history-dependency in response to environmental changes, allowing it
to be found in alternative stable structures. The adaptive nature of niche
interactions, as captured in our framework, can underlie a broad class of
ecological relations and also socio-economic networks that engage in bipartite
cooperation.


  By combining analytical results and simulations of various coarse-grained
models we investigate the minimal energy shape of DNA minicircles which are
torsionally constrained by an imposed over or undertwist. We show that
twist-bend coupling, a cross interaction term discussed in the recent DNA
literature, induces minimal energy shapes with a periodic alternance of parts
with high and low curvature resembling rounded polygons. We briefly discuss the
possible experimental relevance of these findings. We finally show that the
twist and bending energies of minicircles are governed by renormalized
stiffness constants, not the bare ones. This has important consequences for the
analysis of experiments involving circular DNA meant to determine DNA elastic
constants.


  Phylogeny is the field of modelling the temporal discrete dynamics of
speciation. Complex models can nowadays be studied using the Approximate
Bayesian Computation approach which avoids likelihood calculations. The field's
progression is hampered by the lack of robust software to estimate the numerous
parameters of the speciation process. In this work we present an R package,
pcmabc, based on Approximate Bayesian Computations, that implements three novel
phylogenetic algorithms for trait-dependent speciation modelling. Our
phylogenetic comparative methodology takes into account both the simulated
traits and phylogeny, attempting to estimate the parameters of the processes
generating the phenotype and the trait. The user is not restricted to a
predefined set of models and can specify a variety of evolutionary and
branching models. We illustrate the software with a simulation-reestimation
study focused around the branching Ornstein-Uhlenbeck process, where the
branching rate depends non-linearly on the value of the driving
Ornstein-Uhlenbeck process. Included in this work is a tutorial on how to use
the software.


  We introduce a new Wright-Fisher type model for seed banks incorporating
"simultaneous switching", which is motivated by recent work on microbial
dormancy. We show that the simultaneous switching mechanism leads to a new
jump-diffusion limit for the scaled frequency processes, extending the
classical Wright-Fisher and seed bank diffusion limits. We further establish a
new dual coalescent structure with multiple activation and deactivation events
of lineages. While this seems reminiscent of multiple merger events in general
exchangeable coalescents, it actually leads to an entirely new class of
coalescent processes with unique qualitative and quantitative behaviour. To
illustrate this, we provide a novel kind of condition for coming down from
infinity for these coalescents using recent results of Griffiths.


  Life-history tradeoffs among species are major drivers of community assembly.
Most studies investigate how tradeoffs promote deterministic coexistence of
species. It remains unclear how tradeoffs may instead promote historically
contingent exclusion of species, where species dominance is affected by initial
abundances, causing alternative community states. Focusing on the
establishment-longevity tradeoff, we study the transient dynamics and
equilibrium outcomes of competitive interactions in a simulation model of plant
community assembly. We show that, in this model, the establishment-longevity
tradeoff is a necessary but not sufficient condition for alternative stable
equilibria that require also low fecundity for both species. An analytical
approximation of our simulation model demonstrates that alternative stable
equilibria are driven by demographic stochasticity in the number of seeds
arriving at each establishment site. This site-scale stochasticity is only
affected by fecundity and therefore occurs even in infinitely large
communities. In many cases where the establishment-longevity tradeoff does not
cause alternative stable equilibria, it still decreases the rate of convergence
toward the single equilibrium, resulting in decades of transient dynamics that
can appear indistinguishable from alternative stable equilibria in empirical
studies.


  Consider the standard stochastic reaction network model where the dynamics is
given by a continuous-time Markov chain over a discrete lattice. For such
models, estimation of parameter sensitivities is an important problem, but the
existing computational approaches to solve this problem usually require
time-consuming Monte Carlo simulations of the reaction dynamics. Therefore
these simulation-based approaches can only be expected to work over finite
time-intervals, while it is often of interest in applications to examine the
sensitivity values at the steady-state after the Markov chain has relaxed to
its stationary distribution. The aim of this paper is to present a
computational method for the estimation of steady-state parameter
sensitivities, which instead of using simulations, relies on the recently
developed stationary Finite State Projection (sFSP) algorithm [J. Chem. Phys.
147, 154101 (2017)] that provides an accurate estimate of the stationary
distribution at a fixed set of parameters. We show that sensitivity values at
these parameters can be estimated from the solution of a Poisson equation
associated with the infinitesimal generator of the Markov chain. We develop an
approach to numerically solve the Poisson equation and this yields an efficient
estimator for steady-state parameter sensitivities. We illustrate this method
using several examples.


  Bacterial growth presents many beautiful phenomena that pose new theoretical
challenges to statistical physicists, and are also amenable to laboratory
experimentation. This review provides some of the essential biological
background, discusses recent applications of statistical physics in this field,
and highlights the potential for future research.


  Natural evolution gives the impression of leading to an open-ended process of
increasing diversity and complexity. If our goal is to produce such
open-endedness artificially, this suggests an approach driven by evolutionary
metaphor. On the other hand, techniques from machine learning and artificial
intelligence are often considered too narrow to provide the sort of exploratory
dynamics associated with evolution. In this paper, we hope to bridge that gap
by reviewing common barriers to open-endedness in the evolution-inspired
approach and how they are dealt with in the evolutionary case - collapse of
diversity, saturation of complexity, and failure to form new kinds of
individuality. We then show how these problems map onto similar issues in the
machine learning approach, and discuss how the same insights and solutions
which alleviated those barriers in evolutionary approaches can be ported over.
At the same time, the form these issues take in the machine learning
formulation suggests new ways to analyze and resolve barriers to
open-endedness. Ultimately, we hope to inspire researchers to be able to
interchangeably use evolutionary and gradient-descent-based machine learning
methods to approach the design and creation of open-ended systems.


  Allosteric signaling in biological molecules, which may be viewed as specific
action at a distance due to localized perturbation upon binding of ligands or
changes in environmental cues, is pervasive in biology. Phenomenological MWC
and KNF models galvanized research for over five decades, and these models
continue to be the basis for describing the allostery in an array of systems.
However, understanding allosteric signaling and the associated dynamics between
distinct allosteric states is challenging. In this review, we first describe
symmetry and rigidity as essential requirements for allosteric proteins. The
general features, with MWC and KNF as two extreme scenarios, emerge when
allosteric signaling is viewed from an energy landscape perspective. To go
beyond the general theories, we describe computational tools to predict the
network of residues that carry allosteric signals. Methods to obtain molecular
insights into the dynamics of allosteric transitions are briefly mentioned. The
utility of the methods is illustrated by applications to systems ranging from
monomeric proteins to multisubunit proteins. Finally, the role allostery plays
in the functions of ATP-consuming molecular machines, bacterial chaperonin
GroEL and molecular motors, is described. Although universal molecular
principles governing allosteric signaling do not exist, we can draw the
following general conclusions. (1) Multiple pathways connecting allosteric
states are highly heterogeneous. (2) Allosteric signaling is exquisitely
sensitive to the architecture of the system, which implies that the capacity
for allostery is encoded in the structure itself. (3) The mechanical modes that
connect distinct allosteric states are robust to sequence variations. (4)
Extensive investigations of allostery in Hemoglobin and more recently GroEL,
show that a network of salt-bridge rearrangements serves as allosteric
switches.


  The transition distribution of a sample taken from a Wright-Fisher diffusion
with general small mutation rates is found using a coalescent approach. The
approximation is equivalent to having at most one mutation in the coalescent
tree of the sample up to the most recent common ancestor with additional
mutations occurring on the lineage from the most recent common ancestor to the
time origin if complete coalescence occurs before the origin. The sampling
distribution leads to an approximation for the transition density in the
diffusion with small mutation rates. This new solution has interest because the
transition density in a Wright-Fisher diffusion with general mutation rates is
not known.


  Cell growth is determined by substrate availability and the cell's metabolic
capacity to assimilate substrates into building blocks. Metabolic genes that
determine growth rate may interact synergistically or antagonistically, and can
accelerate or slow growth, depending on the genetic background and
environmental conditions. We evolved a diverse set of Escherichia coli
single-gene deletion mutants with a spectrum of growth rates and identified
mutations that generally increase growth rate. Despite the metabolic
differences between parent strains, mutations that enhanced growth largely
mapped to the core transcription machinery, including the $\beta$ and $\beta'$
subunits of RNA polymerase (RNAP) and the transcription elongation factor,
NusA. The structural segments of RNAP that determine enhanced growth have been
previously implicated in antibiotic resistance and in the control of
transcription elongation and pausing. We further developed a computational
framework to characterize how the transcriptional changes that occur upon
acquisition of these mutations affect growth rate across strains. Our
experimental and computational results provide evidence for cases in which RNAP
mutations shift the competitive balance between active transcription and gene
silencing. This study demonstrates that mutations in specific regions of RNAP
are a convergent adaptive solution that can enhance the growth rate of cells
from distinct metabolic states.


  Evolutionary graph theory models the effects of natural selection and random
drift on structured populations of competing mutant and non-mutant individuals.
Recent studies have found that fixation times in such systems often have
right-skewed distributions. Little is known, however, about how these
distributions and their skew depend on mutant fitness. Here we calculate the
fitness dependence of the fixation-time distribution for the Moran Birth-death
process in populations modeled by two extreme networks: the complete graph and
the one-dimensional ring lattice, obtaining exact solutions in the limit of
large network size. We find that with non-neutral fitness, the Moran process on
the ring has normally distributed fixation times, independent of the relative
fitness of mutants and non-mutants. In contrast, on the complete graph, the
fixation-time distribution is a fitness-weighted convolution of two Gumbel
distributions. When fitness is neutral the fixation-time distribution jumps
discontinuously and becomes highly skewed on both the complete graph and the
ring. Even on these simple networks, the fixation-time distribution exhibits
rich fitness dependence, with discontinuities and regions of universality.
Extensions of our results to two-fitness Moran models, times to partial
fixation, and evolution on random networks are discussed.


  Anthropogenic changes of the biota and human hyper-dominance are modulating
the evolution of life on our planet. Humankind has spread worldwide supported
by cultural and technological knowledge, and has already modified uncountable
biological interactions. While numerous species have been extinguished by human
actions, others are directly favored, such as alien species, hybrids, and
genetically modified organisms. These biodiversity shifts have generated new
interactions among all living organisms in anthropized or anthropogenic
ecosystems, with the consequent establishment of novel evolutionary pathways.
Thus, humans have created a strong evolutionary bias on Earth, leading to
unexpected and irreversible outcomes. Anthropogenic changes and novelty
organisms are shifting the evolutionary paths of all organisms towards the
Biological Anthropocene, a new concept of our imprint on biodiversity and
evolution.


  Stochasticity is a key characteristic of intracellular processes such as gene
regulation and chemical signalling. Therefore, characterising stochastic
effects in biochemical systems is essential to understand the complex dynamics
of living things. Mathematical idealisations of biochemically reacting systems
must be able to capture stochastic phenomena. While robust theory exists to
describe such stochastic models, the computational challenges in exploring
these models can be a significant burden in practice since realistic models are
analytically intractable. Determining the expected behaviour and variability of
a stochastic biochemical reaction network requires many probabilistic
simulations of its evolution. Using a biochemical reaction network model to
assist in the interpretation of time course data from a biological experiment
is an even greater challenge due to the intractability of the likelihood
function for determining observation probabilities. These computational
challenges have been subjects of active research for over four decades. In this
review, we present an accessible discussion of the major historical
developments and state-of-the-art computational techniques relevant to
simulation and inference problems for stochastic biochemical reaction network
models. Detailed algorithms for particularly important methods are described
and complemented with MATLAB implementations. As a result, this review provides
a practical and accessible introduction to computational methods for stochastic
models within the life sciences community.


  The plane of bacterial cell division must be precisely positioned. In the
bacterium Myxococcus xanthus, the proteins PomX and PomY form a large cluster,
which is tethered to the nucleoid by the ATPase PomZ and moves in a stochastic,
but biased manner towards midcell, where it initiates cell division.
Previously, a positioning mechanism based on the fluxes of PomZ on the nucleoid
was proposed. However, the cluster dynamics was analyzed in a reduced,
one-dimensional geometry. Here we introduce a mathematical model that accounts
for the three-dimensional shape of the nucleoid, such that nucleoid-bound PomZ
dimers can diffuse past the cluster without interacting with it. Using
stochastic simulations, we find that the cluster still moves to and localizes
at midcell. Redistribution of PomZ by diffusion in the cytosol is essential for
this cluster dynamics. Our mechanism also positions two clusters equidistantly
on the nucleoid. We conclude that a flux-based mechanism allows for cluster
positioning in a biologically realistic three-dimensional cell geometry.


  We present a detailed heuristic method to quantify the degree of local
energetic frustration manifested by protein molecules. Current applications are
realized in computational experiments where a protein structure is visualized
highlighting the energetic conflicts or the concordance of the local
interactions in that structure. Minimally frustrated linkages highlight the
stable folding core of the molecule. Sites of high local frustration, in
contrast, often indicate functionally relevant regions such as binding, active
or allosteric sites.


  Thermogenesis by uncoupling protein (UCP) has traditionally been explained as
the dissipation of proton gradient across the inner mitochondrial membrane into
heat. Herein, we propose that UCPs, aided by the large pore and positively
charged amino acids of suspended loops, enable protonation and transport of
DROS. Thus, UCP facilitates DROS-reactions amongst themselves, forming water
and liberating heat around the inner mitochondrial membrane. Thereby, the
simple murburn scheme for biothermogenesis integrates structural information of
UCP with its attributed physiological function.


  Analytical expressions for the basic reproduction number, R0, have been
obtained in the past for a wide variety of mathematical models for infectious
disease spread, along with expressions for the expected final size of an
outbreak. However, what has so far not been studied is the average number of
infections that descend down the chains of infection begun by each of the
individuals infected in an outbreak (we refer to this quantity as the "average
number of descendant infections" per infectious individual, or ANDI). ANDI
includes not only the number of people that an individual directly contacts and
infects, but also the number of people that those go on to infect, and so on
until that particular chain of infection dies out. Quantification of ANDI has
relevance to the vaccination debate, since with ANDI one can calculate the
probability that one or more people are hospitalised (or die) from a disease
down an average chain of infection descending from an infected un-vaccinated
individual. Here we obtain estimates of ANDI using both deterministic and
stochastic modelling formalisms. With both formalisms we find that even for
relatively small community sizes and under most scenarios for R0 and initial
fraction vaccinated, ANDI can be surprisingly large when the effective
reproduction number is >1, leading to high probabilities of adverse outcomes
for one or more people down an average chain of infection in outbreaks of
diseases like measles.


  Ribonucleic Acid (RNA) can fold into shapes that perform functions in the
cell. These foldings are governed by Watson-Crick base pairing, namely Adenine
to Uracil and Cytosine to Guanine (A-U and G-C). The properties of the H-P
(hydrophobic-hydrophilic) model of protein folding has been well studied in the
two dimensional orthogonal case, and we attempt to achieve similar results. We
prove that (1) there is an infinite family of even-length sequences with unique
optimal foldings, (2) there are two ideal foldings for an even length sequence
and given a sequence is is quickly verifiable if both, one, or neither are
optimal, (3) finding an optimal foldings under this model is NP-hard, and (4)
we give a constant-factor approximation algorithm for optimally folding RNA
sequences.


  While G-protein coupled receptors (GPCRs) constitute the largest class of
membrane proteins, structures and endogenous ligands of a large portion of
GPCRs remain unknown. Due to the involvement of GPCRs in various signaling
pathways and physiological roles, the identification of endogenous ligands as
well as designing novel drugs is of high interest to the research and medical
communities. Along with highlighting the recent advances in structure-based
ligand discovery, including docking and molecular dynamics, this article
focuses on the latest advances for automating the discovery of bioactive
ligands using machine learning. Machine learning is centered around the
development and applications of algorithms that can learn from data
automatically. Such an approach offers immense opportunities for bioactivity
prediction as well as quantitative structure-activity relationship studies.
This review describes the most recent and successful applications of machine
learning for bioactive ligand discovery, concluding with an outlook on deep
learning methods that are capable of automatically extracting salient
information from structural data as a promising future direction for rapid and
efficient bioactive ligand discovery.


  In lieu of an abstract here is the first paragraph: No other species remotely
approaches the human capacity for the cultural evolution of novelty that is
accumulative, adaptive, and open-ended (i.e., with no a priori limit on the
size or scope of possibilities). By culture we mean extrasomatic
adaptations--including behavior and technology--that are socially rather than
sexually transmitted. This chapter synthesizes research from anthropology,
psychology, archaeology, and agent-based modeling into a speculative yet
coherent account of two fundamental cognitive transitions underlying human
cultural evolution that is consistent with contemporary psychology. While the
chapter overlaps with a more technical paper on this topic (Gabora & Smith
2018), it incorporates new research and elaborates a genetic component to our
overall argument. The ideas in this chapter grew out of a non-Darwinian
framework for cultural evolution, referred to as the Self-other Reorganization
(SOR) theory of cultural evolution (Gabora, 2013, in press; Smith, 2013), which
was inspired by research on the origin and earliest stage in the evolution of
life (Cornish-Bowden & C\'ardenas 2017; Goldenfeld, Biancalani, & Jafarpour,
2017, Vetsigian, Woese, & Goldenfeld 2006; Woese, 2002). SOR bridges
psychological research on fundamental aspects of our human nature such as
creativity and our proclivity to reflect on ideas from different perspectives,
with the literature on evolutionary approaches to cultural evolution that
aspire to synthesize the behavioral sciences much as has been done for the
biological scientists. The current chapter is complementary to this effort, but
less abstract; it attempts to ground the theory of cultural evolution in terms
of cognitive transitions as suggested by archaeological evidence.


  Neuroinflammation in utero may result in lifelong neurological disabilities.
Astrocytes play a pivotal role, but the mechanisms are poorly understood. No
early postnatal treatment strategies exist to enhance neuroprotective potential
of astrocytes. We hypothesized that agonism on alpha7 nicotinic acetylcholine
receptor (alpha7nAChR) in fetal astrocytes will augment their neuroprotective
transcriptome profile, while the antagonistic stimulation of alpha7nAChR will
achieve the opposite. Using an in vivo - in vitro model of developmental
programming of neuroinflammation induced by lipopolysaccharide (LPS), we
validated this hypothesis in primary fetal sheep astrocytes cultures re-exposed
to LPS in presence of a selective alpha7nAChR agonist or antagonist. Our RNAseq
findings show that a pro-inflammatory astrocyte transcriptome phenotype
acquired in vitro by LPS stimulation is reversed with alpha7nAChR agonistic
stimulation. Conversely, antagonistic alpha7nAChR stimulation potentiates the
pro-inflammatory astrocytic transcriptome phenotype. Furthermore, we conduct a
secondary transcriptome analysis against the identical alpha7nAChR experiments
in fetal sheep primary microglia cultures and discuss the implications for
neurodevelopment.


  Antifragility is a property that enhances the capability of a system in
response to external perturbations. Although the concept has been applied in
many areas, a practical measure of antifragility has not been developed yet.
Here we propose a simply calculable measure of antifragility, based on the
change of "satisfaction" before and after adding perturbations, and apply it to
random Boolean networks (RBNs). Using the measure, we found that ordered RBNs
are the most antifragile. Also, we demonstrated that seven biological systems
are antifragile. Our measure and results can be used in various applications of
Boolean networks (BNs) including creating antifragile engineering systems,
identifying the genetic mechanism of antifragile biological systems, and
developing new treatment strategies for various diseases.


  We used a hybrid fluorescence spectroscopic toolkit to monitor T4 Lysozyme
(T4L) in action. By unraveling the kinetic and dynamic interplay of the
conformational states, we sought to elucidate the dynamic structural biology of
T4L. In particular, by combining single-molecule and ensemble multiparameter
fluorescence detection, EPR spectroscopy, mutagenesis, and FRET-positioning and
screening, we characterized three short-lived conformational states within the
conformational landscape of the T4L over the ns-ms timescale. The use of 33
FRET-derived distance sets, to screen known T4L structures, revealed that T4L
in solution mainly adopts the known open and closed states in exchange at 4
$\mu$s. A newly found minor state, undisclosed by at present more than 500
crystal structures of T4L and sampled at 230 $\mu$s, may be actively involved
in the product release step in catalysis. The presented fluorescence
spectroscopic toolkit is anticipated to accelerate the development of dynamic
structural biology.


  We present an exhaustive description of a stochastic, event-driven,
hierarchical agent-based model designed to reproduce the infectious state of
the cattle disease called Bovine Viral Diarrhea, for which the livestock-trade
network is the main route of spreading. For the farm-node dynamics, it takes
into account a vast number of breeding, infectious and animal movement
mechanisms via a susceptible-infected-recovered type of dynamics with an
additional permanently infectious class. The interaction between the farms is
described by a supply-demand farm manager mechanism governing the network
structure and dynamics. The model includes realistic disease and breeding
dynamics and allows to study numerous mitigation strategies of present and past
government regulations, including different testing and vaccination scenarios.


  We study the global dynamics of a SICA infection model with general incidence
rate. The proposed model is calibrated with cumulative cases of infection by
HIV-AIDS in Morocco from 1986 to 2015. We first prove that our model is
biologically and mathematically well-posed. Stability analysis of different
steady states is performed and threshold parameters are identified where the
model exhibits clearance of infection or maintenance of a chronic infection.
Furthermore, we examine the robustness of the model to some parameter values by
examining the sensitivity of the basic reproduction number. Finally, using
numerical simulations with real data from Morocco, we show that the model
predicts well such reality.


  p27$^{Kip1}$ (p27) is an intrinsically disordered protein (IDP) that folds
upon binding to cyclin-dependent kinase (Cdk)$/$cyclin complexes (e.g.,
Cdk2$/$cyclin A), inhibiting their catalytic activity and causing cell cycle
arrest. However, cell division progresses when stably Cdk2$/$cyclin A-bound p27
is phosphorylated on one or two structurally occluded tyrosine residues
$[$tyrosines 88 (Y88) and 74 (Y74)$]$ and a distal threonine residue
$[$threonine 187 (T187)$]$. These events trigger ubiquitination and degradation
of p27, fully activating Cdk2$/$cyclin A to drive cell division. Using an
integrated approach comprising structural, biochemical, biophysical and
single-molecule fluorescence methods, we show that Cdk2$/$cyclin A-bound p27
samples lowly-populated conformations that dynamically anticipate the
sequential steps of this signaling cascade. "Dynamic anticipation" provides
access to the non-receptor tyrosine kinases, BCR-ABL and Src, which
sequentially phosphorylate Y88 and Y74 and promote intra-assembly
phosphorylation (of p27) on distal T187. Tyrosine phosphorylation also
allosterically relieves p27-dependent inhibition of substrate binding to
Cdk2$/$cyclin A, a phenomenon we term "cross-complex allostery". Even when
tightly bound to Cdk2$/$cyclin A, intrinsic flexibility enables p27 to
integrate and process signaling inputs, and generate outputs including altered
Cdk2 activity, p27 stability, and, ultimately, cell cycle progression.
Intrinsic dynamics within multi-component assemblies may be a general mechanism
of signaling by regulatory IDPs, which can be subverted in human disease, as
exemplified by hyper-active BCR-ABL and Src in certain cancers.


  Here we explore a new algorithmic framework for multi-agent reinforcement
learning, called Malthusian reinforcement learning, which extends self-play to
include fitness-linked population size dynamics that drive ongoing innovation.
In Malthusian RL, increases in a subpopulation's average return drive
subsequent increases in its size, just as Thomas Malthus argued in 1798 was the
relationship between preindustrial income levels and population growth.
Malthusian reinforcement learning harnesses the competitive pressures arising
from growing and shrinking population size to drive agents to explore regions
of state and policy spaces that they could not otherwise reach. Furthermore, in
environments where there are potential gains from specialization and division
of labor, we show that Malthusian reinforcement learning is better positioned
to take advantage of such synergies than algorithms based on self-play.