  The cellular uptake of nanoparticles or viruses requires that the gain of
adhesion energy overcomes the cost of plasma membrane bending. It is well known
that this leads to a minimal particle size for uptake. Using a simple
deterministic theory for this process, we first show that, for the same radius
and volume, cylindrical particles should be taken up faster than spherical
particles, both for normal and parallel orientations. We then address
stochastic effects, which are expected to be relevant due to small system size,
and show that, now, spherical particles can have a faster uptake, because the
mean first passage time profits from the multiplicative noise induced by the
spherical geometry. We conclude that stochastic effects are strongly geometry
dependent and may favor spherical shapes during adhesion-driven particle
uptake.


  Receptor-mediated endocytosis requires that the energy of adhesion overcomes
the deformation energy of the plasma membrane. The resulting driving force is
balanced by dissipative forces, leading to deterministic dynamical equations.
While the shape of the free membrane does not play an important role for tensed
and loose membranes, in the intermediate regime it leads to an important energy
barrier. Here we show that this barrier is similar to but different from an
effective line tension and suggest a simple analytical approximation for it. We
then explore the rich dynamics of uptake for particles of different shapes and
present the corresponding dynamical state diagrams. We also extend our model to
include stochastic fluctuations, which facilitate uptake and lead to
corresponding changes in the phase diagrams.


  Recently it was shown that, if the subtree and chain reduction rules have
been applied exhaustively to two unrooted phylogenetic trees, the reduced trees
will have at most 15k-9 taxa where k is the TBR (Tree Bisection and
Reconnection) distance between the two trees, and that this bound is tight.
Here we propose five new reduction rules and show that these further reduce the
bound to 11k-9. The new rules combine the ``unrooted generator'' approach
introduced in [Kelk and Linz 2018] with a careful analysis of agreement forests
to identify (i) situations when chains of length 3 can be further shortened
without reducing the TBR distance, and (ii) situations when small subtrees can
be identified whose deletion is guaranteed to reduce the TBR distance by 1. To
the best of our knowledge these are the first reduction rules that strictly
enhance the reductive power of the subtree and chain reduction rules.


  Evolving systems, be it an antibody repertoire in the face of mutating
pathogens or a microbial population exposed to varied antibiotics, constantly
search for adaptive solutions in time-varying fitness landscapes. Generalists
correspond to genotypes that remain fit across diverse selective pressures;
cross-reactive antibodies are much wanted but rare, while multi-drug resistant
microbes are undesired yet prevalent. However, little is known about under what
conditions such solutions with a high capacity to adapt would be efficiently
discovered by evolution, as environmental changes alter the relative fitness
and accessibility of neighboring genotypes. In addition, can epistasis --- the
source of landscape ruggedness and path constraints --- play a different role,
if the environments are correlated in time? We present a generative model to
estimate the propensity of evolving generalists in rugged landscapes that are
tunably related and cycling relatively slowly. We find that environment cycling
can substantially facilitate the search for fit generalists by dynamically
enlarging their effective basins of attraction. Importantly, these high
performers are most likely to emerge at an intermediate level of both
ruggedness and environmental relatedness, trading diversity for fitness and
accessibility. Our work provides a conceptual framework to study evolution in
correlated varying complex environments, and offers statistical understanding
that suggests general strategies for speeding up the generation of broadly
neutralizing antibodies or preventing microbes from evolving multi-drug
resistance.


  In this experimental study fluorescent silver clusters on tRNA matrix were
synthesized for the first time. Obtained complexes have two emission regions in
visible part of spectrum. We have also studied our complexes using FTIR
spectroscopy. The possible binding sites of the clusters have been suggested.


  Consumers face numerous risks that can be minimized by incorporating
different life-history strategies. How much and when a consumer adds to its
energetic reserves or invests in reproduction are key behavioral and
physiological adaptations that structure much of how organisms interact. Here
we develop a theoretical framework that explicitly accounts for stochastic
fluctuations of an individual consumer's energetic reserves while foraging and
reproducing on a landscape with resources that range from uniformly distributed
to highly clustered. First, we show that optimal life-history strategies vary
in response to changes in the mean productivity of the resource landscape,
where depleted environments promote reproduction at lower energetic states,
greater investment in each reproduction event, and smaller litter sizes. We
then show that if resource variance scales with body size due to landscape
clustering, consumers that forage for clustered foods are susceptible to strong
Allee effects, increasing extinction risk. Finally, we show that the proposed
relationship between consumer body size, resource clustering, and Allee
effect-induced population instability offers key ecological insights into the
evolution of large-bodied grazing herbivores from small-bodied browsing
ancestors.


  Chaperone-assisted translocation through a nanopore embedded in membrane
holds a prominent role in the transport of biopolymers. Inspired by classical
Brownian ratchet, we develop a theoretical framework characterizing such
translocation process through a master equation approach. In this framework,
the polymer chain, provided with reversible binding of chaperones, undergoes
forward/backward diffusion, which is rectified by chaperones. We drop the
assumption of timescale separation and keep the length of a polymer chain
finite, both of which happen to be the key points in most of the previous
studies. Our framework makes it accessible to derive analytical expressions for
mean translocation velocity and effective diffusion constant in stationary
state, which is the basis of a comprehensive understanding towards the dynamics
of such process. Generally, the translocation of polymer chain across membrane
consists of three subprocesses: initiation, termination, and translocation of
the main body part of a polymer chain, where the translocation of the main body
part depends on the binding/unbinding kinetics of chaperones. That is the main
concern of this study. Our results show that the increase of forward/backward
diffusion rate of a polymer chain and the binding/unbinding ratio of chaperones
both raise the mean translocation velocity of a polymer chain, and roughly
speaking, the dependence of effective diffusion constant on these two factors
achieves similar behavior.


  We provide a machine learning solution that replaces the traditional methods
for deciding the pesticide application time of Sunn Pest. We correlate climate
data with phases of Sunn Pest in its life-cycle and decide whether the fields
should be sprayed. Our solution includes two groups of prediction models. The
first group contains decision trees that predict migration time of Sunn Pest
from winter quarters to wheat fields. The second group contains random forest
models that predict the nymphal stage percentages of Sunn Pest which is a
criterion for pesticide application. We trained our models on four years of
climate data which was collected from Kir\c{s}ehir and Aksaray. The experiments
show that our promised solution make correct predictions with high accuracies.


  With the multiallelic parent-independent mutation-drift model, the
equilibrium proportions of alleles are known to be Dirichlet distributed. A
special case is the biallelic model, in which the proportions are beta
distributed. A sample taken from these models is then Dirichlet-multinomially
or beta-binomially distributed, respectively. Maximum likelihood (ML)
estimators for the mutation parameters of the biallelic parent-independent
mutation model are available via an expectation maximization algorithm.
Assuming small scaled mutation rates, the distribution of a sample of size $M$
can be expanded in a Taylor series of first order. Then the ML estimators for
the two parameters in the biallelic model can be expressed using the site
frequency spectrum. In this article, we go beyond parent-independent mutation
and analyse a strand-symmetric mutation model with six scaled mutation
parameters that deviates from parent independent mutation and, generally, from
detailed balance. We derive ML estimators for these six parameters assuming
mutation-drift equilibrium and small scaled mutation rates. This is the first
time that ML estimators are provided for a mutation model more complex than
parent-independent mutation.


  Cells can utilize chemical communication to exchange information and
coordinate their behavior in the presence of noise. Communication can reduce
noise to shape a collective response, or amplify noise to generate distinct
phenotypic subpopulations. Here we discuss a moment-based approach to study how
cell-cell communication affects noise in biochemical networks that arises from
both intrinsic and extrinsic sources. We derive a system of approximate
differential equations that captures lower-order moments of a population of
cells, which communicate by secreting and sensing a diffusing molecule. Since
the number of obtained equations grows combinatorially with number of
considered cells, we employ a previously proposed model reduction technique,
which exploits symmetries in the underlying moment dynamics. Importantly, the
number of equations obtained in this way is independent of the number of
considered cells such that the method scales to arbitrary population sizes.
Based on this approach, we study how cell-cell communication affects population
variability in several biochemical networks. Moreover, we analyze the accuracy
and computational efficiency of the moment-based approximation by comparing it
with moments obtained from stochastic simulations.


  In the marine environment biological processes are strongly affected by
oceanic currents, particularly by eddies (vortices) formed by the hydrodynamic
flow field. Employing a kinematic flow field coupled to a population dynamical
model for plankton growth, we study the impact of an intermittent upwelling of
nutrients on triggering harmful algal blooms (HABs). Though it is widely
believed that additional nutrients boost the formation of HABs or algal blooms
in general, we show that the response of the plankton to nutrient plumes
depends crucially on the mesoscale hydrodynamic flow structure. In general
nutrients can either be quickly washed out from the observation area, or can be
captured by the vortices in the flow. The occurrence of either scenario depends
on the relation between the time scales of the vortex formation and nutrient
upwelling as well as the time instants at which upwelling pulse occurs and how
long do they last. We show that these two scenarios result in very different
responses in plankton dynamics which makes it very difficult to predict,
whether nutrient upwelling will lead to a HAB or not. This explains, why
observational data are sometimes inconclusive establishing a correlation
between upwelling events and plankton blooms.


  Since a long time the application of small electric potentials on biological
membranes (iontophoresis) proved enabling control and improvement of
transdermal delivery of substances across this barrier. In spite of a large
experimental data, the detailed molecular mechanism of iontophoresis is absent.
In the present work the interaction among the external continuous electric
field with the outermost layer of the skin (\textit{stratum corneum}) was
studied by coarse-grained molecular dynamics. Our results pointed out the
occurrence of water-rich vesicles formation depending on the field strength.
The corresponding phase diagram indicated that the large set of phenomena
(vesicle formation, reversibility, phase transition, disruption) could be
completely controlled by tuning external continuous electric fields.
Interestingly, electric field shielding effects are in the origin of observed
effects and followed a general Arrhenius-like time dependence. Direct current
(DC) electric field usage would also have booster diffusion effects due to
vesicles creation and reincorporation which would have direct beneficial
absorption effects on water-soluble topical agents (vitamins) or dermal
jet-injection of drugs by fine needles.


  Molecular data systems have the potential to store information at
dramatically higher density than existing electronic media. Some of the first
experimental demonstrations of this idea have used DNA, but nature also uses a
wide diversity of smaller non-polymeric molecules to preserve, process, and
transmit information. In this paper, we present a general framework for
quantifying chemical memory, which is not limited to polymers and extends to
mixtures of molecules of all types. We show that the theoretical limit for
molecular information is two orders of magnitude denser by mass than DNA,
although this comes with different practical constraints on total capacity. We
experimentally demonstrate kilobyte-scale information storage in mixtures of
small synthetic molecules, and we consider some of the new perspectives that
will be necessary to harness the information capacity available from the vast
non-genomic chemical space.


  Spatial studies of transcriptome provide biologists with gene expression maps
of heterogeneous and complex tissues. However, most experimental protocols for
spatial transcriptomics suffer from the need to select beforehand a small
fraction of genes to be quantified over the entire transcriptome. Standard
single-cell RNA sequencing (scRNA-seq) is more prevalent, easier to implement
and can in principle capture any gene but cannot recover the spatial location
of the cells. In this manuscript, we focus on the problem of imputation of
missing genes in spatial transcriptomic data based on (unpaired) standard
scRNA-seq data from the same biological tissue. Building upon domain adaptation
work, we propose gimVI, a deep generative model for the integration of spatial
transcriptomic data and scRNA-seq data that can be used to impute missing
genes. After describing our generative model and an inference procedure for it,
we compare gimVI to alternative methods from computational biology or domain
adaptation on real datasets and outperform Seurat Anchors, Liger and CORAL to
impute held-out genes.


  Gene interaction graphs aim to capture various relationships between genes
and can represent decades of biology research. When trying to make predictions
from genomic data, those graphs could be used to overcome the curse of
dimensionality by making machine learning models sparser and more consistent
with biological common knowledge. In this work, we focus on assessing how well
those graphs capture dependencies seen in gene expression data to evaluate the
adequacy of the prior knowledge provided by those graphs. We propose a
condition graphs should satisfy to provide good prior knowledge and test it
using `Single Gene Inference' tasks. We also compare with randomly generated
graphs, aiming to measure the true benefit of using biologically relevant
graphs in this context, and validate our findings with five clinical tasks. We
find some graphs capture relevant dependencies for most genes while being very
sparse. Our analysis with random graphs finds that dependencies can be captured
almost as well at random which suggests that, in terms of gene expression
levels, the relevant information about the state of the cell is spread across
many genes.


  The behavior of ecological systems mainly relies on the interactions between
the species it involves. We consider the problem of inferring the species
interaction network from abundance data. To be relevant, any network inference
methodology needs to handle count data and to account for possible
environmental effects. It also needs to distinguish between direct interactions
and indirect associations and graphical models provide a convenient framework
for this purpose. We introduce a generic statistical model for network
inference based on abundance data. The model includes fixed effects to account
for environmental covariates and sampling efforts, and correlated random
effects to encode species interactions. The inferred network is obtained by
averaging over all possible tree-shaped (and therefore sparse) networks, in a
computationally efficient manner. An output of the procedure is the probability
for each edge to be part of the underlying network. A simulation study shows
that the proposed methodology compares well with state-of-the-art approaches,
even when the underlying graph strongly differs from a tree. The analysis of
two datasets highlights the influence of covariates on the inferred network.
Accounting for covariates is critical to avoid spurious edges. The proposed
approach could be extended to perform network comparison or to look for missing
species.


  Shiga toxin-producing Escherichia coli O157: H7 is a food-borne pathogen and
the major cause of hemorrhagic colitis. Pseudomonas is the genus most frequent
psychrotrophic spoilage microorganisms present in milk. Two-species bacterial
systems with Escherichia coli O157: H7, non-pathogenic Escherichia coli, and
Pseudomonas fluorescens in skimmed milk at 7, 13, 19, or 25 C was studied.
Bacterial interactions were modelled after applying a Bayesian approach. No
direct correlation between Pseudomonas fluorescent's growth rate and its effect
on the maximum population densities of Escherichia coli species was found. The
results show the complexity of the interactions between two species into a food
model. The use of natural microbiota members to control foodborne pathogens
could be useful to improve food safety during the processing and storage of
refrigerated foods.


  In population dynamics, mathematical models often contain too many parameters
to be easily testable. A way to reliably estimate parameters for a broad range
of systems would help us obtain clearer predictions from theory. In this paper,
we examine how the allometric scaling of a number of biological quantities with
animal mass may be useful to parameterise population dynamical models. Using
this allometric scaling, we make predictions about the ratio of prey to
predators in real ecosystems, and we attempt to estimate the length of animal
population cycles as a function of mass. Our analytical and numerical results
turn out to compare reasonably to data from a number of ecosystems. This paves
the way for a wider usage of allometric scaling to simplify mathematical models
in population dynamics and make testable predictions.


  The aims of this work is to analyse of the global stability of the extended
model of hepatitis C virus(HCV) infection with cellular proliferation,
spontaneous cure and hepatocyte homeostasis. We first give general information
about hepatitis C. Secondly, We prove the existence of local, maximal and
global solutions of the model and establish some properties of this solution as
positivity and asymptotic behaviour. Thirdly we show, by the construction of an
appropriate Lyapunov function, that the uninfected equilibrium and the unique
infected equilibrium of the model of HCV are globally asymptotically stable
respectively when the threshold number $\mathcal{R}_{0}<1-\frac{q}{d_{I}+q}$
and when $\mathcal{R}_{0}>1$. Finally, some numerical simulations are carried
out using Maple software confirm these theoretical results.


  In this paper we analyze the equilibrium properties of a large class of
stochastic processes describing the fundamental biological process within
bacterial cells, {\em the production process of proteins}. Stochastic models
classically used in this context to describe the time evolution of the numbers
of mRNAs and proteins are presented and discussed. An extension of these
models, which includes elongation phases of mRNAs and proteins, is introduced.
A convergence result to equilibrium for the process associated to the number of
proteins and mRNAs is proved and a representation of this equilibrium as a
functional of a Poisson process in an extended state space is obtained.
Explicit expressions for the first two moments of the number of mRNAs and
proteins at equilibrium are derived, generalizing some classical formulas.
Approximations used in the biological literature for the equilibrium
distribution of the number of proteins are discussed and investigated in the
light of these results. Several convergence results for the distribution of the
number of proteins at equilibrium are in particular obtained under different
scaling assumptions.


  Background: Chromatin immunoprecipitation coupled to next generation
sequencing (ChIP-Seq) is a widely used technique to investigate the function of
chromatin-related proteins in a genome-wide manner. ChIP-Seq generates large
quantities of data which can be difficult to process and analyse, particularly
for organisms with contig based genomes. Contig-based genomes often have poor
annotations for cis-elements, for example enhancers, that are important for
gene expression. Poorly annotated genomes make a comprehensive analysis of
ChIP-Seq data difficult and as such standardized analysis pipelines are
lacking. Methods: We report a computational pipeline that utilizes traditional
High-Performance Computing techniques and open source tools for processing and
analysing data obtained from ChIP-Seq. We applied our computational pipeline
"Rapid Analysis of ChIP-Seq data" (RACS) to ChIP-Seq data that was generated in
the model organism Tetrahymena thermophila, an example of an organism with a
genome that is available in contigs. Results: To test the performance and
efficiency of RACs, we performed control ChIP-Seq experiments allowing us to
rapidly eliminate false positives when analyzing our previously published data
set. Our pipeline segregates the found read accumulations between genic and
intergenic regions and is highly efficient for rapid downstream analyses.
Conclusions: Altogether, the computational pipeline presented in this report is
an efficient and highly reliable tool to analyze genome-wide ChIP-Seq data
generated in model organisms with contig-based genomes.
  RACS is an open source computational pipeline available to download from:
https://bitbucket.org/mjponce/racs --or--
https://gitrepos.scinet.utoronto.ca/public/?a=summary&p=RACS


  A growing number of experimental evidence shows that it is general for a
ligand binding protein to have a potential for allosteric regulation and for
further evolution. In addition, such proteins generically change their
conformation upon binding. O. Rivoire has recently proposed an evolutionary
scenario that explains these properties as a generic byproduct of selection for
exquisite discrimination between very similar ligands. The initial claim was
supported by two classes of basic examples: continuous protein models with
small numbers of degrees of freedom, on which the development of a
conformational switch was established, and a 2-dimensional spin glass model
supporting the rest of the statement. This work aimed to clarify the
implication of the exquisite discrimination for smooth models with large number
of degrees of freedom, the situation closer to real biological systems. With
the help of differential geometry, jet-space analysis, and transversality
theorems, it is shown that the claim holds true for any generic flexible system
that can be described in terms of smooth manifolds. The result suggests that,
indeed, evolutionary solutions to the exquisite discrimination problem, if
exist, are located near a codimension-1 subspace of the appropriate genotypical
space. This constraint, in turn, gives rise to a potential for the allosteric
regulation of the discrimination via generic conformational changes upon
binding.


  Genomic data visualization is essential for interpretation and hypothesis
generation as well as a valuable aid in communicating discoveries. Visual tools
bridge the gap between algorithmic approaches and the cognitive skills of
investigators. Addressing this need has become crucial in genomics, as
biomedical research is increasingly data-driven and many studies lack
well-defined hypotheses. A key challenge in data-driven research is to discover
unexpected patterns and to formulate hypotheses in an unbiased manner in vast
amounts of genomic and other associated data. Over the past two decades, this
has driven the development of numerous data visualization techniques and tools
for visualizing genomic data. Based on a comprehensive literature survey, we
propose taxonomies for data, visualization, and tasks involved in genomic data
visualization. Furthermore, we provide a comprehensive review of published
genomic visualization tools in the context of the proposed taxonomies.


  Extensive post reproductive lifespan (PRLS) is observed only in a few
species, such as humans or resident killer whales, and its origin is under
debate. Hypotheses like mother-care and grandmother-care invoke strategies of
investment--provision to one's descendants to enhance one's overall
reproductive success--to explain PRLS. The contribution of an investment
strategy varies with the age of the caregiver, as the number of care-receiving
descendant changes with age. Here we simulated an agent based model, which is
sensitive to age-specific selection, to examine how the investment strategies
in different hypotheses affect survival and reproduction across different
stages of life. We found that extensive PRLS emerges if we combine multiple
investment strategies, including grandmother-care but not mother-care, which
allow an individual to have an increasing contribution as it ages. We also
found that, if mother-care is further introduced to the PRLS-enabling
strategies, it will let contribution at mid-life to substitute contribution at
late life, which consequently terminates extensive PRLS.


  Free energy landscapes provide insights into conformational ensembles of
biomolecules. In order to analyze these landscapes and elucidate mechanisms
underlying conformational changes, there is a need to extract metastable states
with limited noise. This has remained a formidable task, despite a plethora of
existing clustering methods. We present InfleCS, a novel method for extracting
well-defined core states from free energy landscapes. The method is based on a
Gaussian mixture free energy estimator and exploits the shape of the estimated
density landscape. The core states that naturally arise from the clustering
allow for detailed characterization of the conformational ensemble. The
clustering quality is evaluated on three toy models with different properties,
where the method is shown to consistently outperform other conventional and
state-of-the-art clustering methods. Finally, the method is applied to a
temperature enhanced molecular dynamics simulation of Ca2+ bound Calmodulin.
Through the free energy landscape, we discover a pathway between a canonical
and a compact state, revealing conformational changes driven by electrostatic
interactions.


  Immunotherapy has recently shown important clinical successes in a
substantial number of oncology indications. Additionally, the tumor somatic
mutation load has been shown to associate with response to these therapeutic
agents, and specific mutational signatures are hypothesized to improve this
association, including signatures related to pathogen insults. We sought to
study in silico the validity of these observations and how they relate to each
other. We first addressed whether somatic mutations typically involved in
cancer may increase, in a statistically meaningful manner, the similarity
between common pathogens and the human exome. Our study shows that common
mutagenic processes increase, in the upper range of biologically plausible
frequencies, the similarity between cancer exomes and pathogen DNA at a scale
of 12-16 nucleotide sequences and established that this increased similarity is
due to the specific mutation distribution of the considered mutagenic
processes. Next, we studied the impact of mutation rate and showed that
increasing mutation rate generally results in an increased similarity between
the cancer exome and pathogen DNA, at a scale of 4-5 amino acids. Finally, we
investigated whether the considered mutational processes result in amino-acid
changes with functional relevance that are more likely to be immunogenic. We
showed that functional tolerance to mutagenic processes across species
generally suggests more resilience to mutagenic processes that are due to
exposure to elements of nature than to mutagenic processes that are due to
exposure to cancer-causing artificial substances. These results support the
idea that recognition of pathogen sequences as well as differential functional
tolerance to mutagenic processes may play an important role in the immune
recognition process involved in tumor infiltration by lymphocytes.


  We consider the dynamics of spatial stochastic May-Leonard models with mutual
predation interactions of equal strength between any two individuals of
different species. Using two-dimensional simulations, with two and three
pecies, we investigate the dynamical impact of the death of individuals after a
given threshold number of successive unsuccessful predation attempts. We find
that the death of these individuals can have a strong impact on the dynamics of
population networks and provide a crucial contribution to the preservation of
coexistence.


  The fracture and severing of polymer chains plays a critical role in the
failure of fibrous materials and the regulated turnover of intracellular
filaments. Using continuum wormlike chain models, we investigate the fracture
of semiflexible polymers via thermal bending fluctuations, focusing on the role
of filament flexibility and dynamics. Our results highlight a previously
unappreciated consequence of mechanical heterogeneity in the filament, which
enhances the rate of thermal fragmentation particularly in cases where
constraints hinder the movement of the chain ends. Although generally
applicable to semiflexible chains with regions of different bending stiffness,
the model is motivated by a specific biophysical system: the enhanced severing
of actin filaments at the boundary between stiff bare regions and mechanically
softened regions that are coated with cofilin regulatory proteins. The results
presented here point to a potential mechanism for disassembly of filament
networks in general and cytoskeletal actin networks in particular by the
introduction of locally softened chain regions, as occurs with cofilin binding.


  Thyroid hormone (TH) signalling, an evolutionary conserved pathway, is
crucial for brain function and cognition throughout life, from early
development to ageing. In humans, TH deficiency during pregnancy alters
offspring brain development, increasing the risk of cognitive disorders. How TH
regulates neurogenesis and subsequent behaviour and cognitive functions remains
a major research challenge. Cellular and molecular mechanisms underlying TH
signalling on proliferation, survival, determination, migration,
differentiation and maturation have been studied in mammalian animal models for
over a century. However, recent data show that THs also influence embryonic and
adult neurogenesis throughout vertebrates (from mammals to teleosts). These
latest observations raise the question of how TH availability is controlled
during neurogenesis and particularly in specific neural stem cell populations.
This review deals with the role of TH in regulating neurogenesis in the
developing and the adult brain across different vertebrate species. Such
evo-devo approaches can shed new light on (i) the evolution of the nervous
system and (ii) the evolutionary control of neurogenesis by TH across animal
phyla. We also discuss the role of thyroid disruptors on brain development in
an evolutionary context.


  In social-ecological systems (SESs), social and biophysical dynamics interact
within and between structural levels separated by spatial and temporal scales.
Cross-scale interactions (CSIs) are interdependences between processes at
different scales, generating behaviour unpredictable at single scales.
Understanding CSIs is important for improving SES governance but they remain
understudied. Theoretical models are needed, which capture essential features
while being simple enough to yield insights into mechanisms. In a stylised
model, we study CSIs in a two-level system of weakly interacting communities
harvesting a common-pool resource. Community members adaptively conform to, or
defect from, a norm of socially optimal harvesting, enforced through social
sanctioning both within and between communities. Each subsystem's dynamics
depend sensitively on the other despite interactions being much weaker between
subsystems than within them. When interaction is purely biophysical, stably
high cooperation in one community can cause cooperation in the other to
collapse. However, even weak social interaction can prevent collapse of
cooperation and instead cause collapse of defection. We identify conditions
under which subsystem-level cooperation produces desirable system-level
outcomes. Our findings expand evidence that collaboration is important for
sustainably managing shared resources, showing its importance even when
resource sharing and social relationships are weak.


  Rock-scissors-paper game, as the simplest model of intransitive relation
between competing agents, is a frequently quoted model to explain the stable
diversity of competitors in the race of surviving. When increasing the number
of competitors we may face a novel situation because beside the mentioned
unidirectional predator-prey-like dominance a balanced or peer relation can
emerge between some competitors. By utilizing this possibility in the present
work we generalize a four-state predator-prey type model where we establish two
groups of species labeled by even and odd numbers. In particular, we introduce
different invasion probabilities between and within these groups, which results
in a tunable intensity of bidirectional invasion among peer species. Our study
reveals an exceptional richness of pattern formations where five quantitatively
different phases are observed by varying solely the strength of the mentioned
inner invasion. The related transition points can be identified with the help
of appropriate order parameters based on the spatial autocorrelation decay, on
the fraction of empty sites, and on the variance of the species density.
Furthermore, the application of diverse, alliance-specific inner invasion rates
for different groups may result in the extinction of the pair of species where
this inner invasion is moderate. These observations highlight that beyond the
well-known and intensively studied cyclic dominance there is an additional
source of complexity of pattern formation that has not been explored earlier.


  The recent realization that entire communities fuse and separate (community
coalescence) has led to a reappraisal of the forces determining species
diversity and dynamics, especially in microbial communities where coalescence
is likely widespread. To understand if connectedness by coalescence results in
different outcomes from connectedness by individual dispersal, we investigated
chance processes leading to loss of species diversity using a model of a
neutral two-species metacommunity. Two scenarios were investigated: pairwise
coalescence where the communities coalesce in pairs, intermix and then
separate; and diffuse coalescence where several communities mix as a pool and
are re-distributed to their original patches. When standardized for the same
net movement, both types of coalescence led to a longer time to single species
dominance than dispersal. Coalescence therefore may be an important process
contributing to the surprisingly high microbial species diversity in nature.


  Genetic systems with multiple loci can have complex dynamics. For example,
mean fitness need not always increase and stable cycling is possible. Here, we
study the dynamics of a genetic system inspired by the molecular biology of
recognition-dependent double strand breaks and repair as it happens in
recombination hotspots. The model shows slow-fast dynamics in which the system
converges to the quasi-linkage equilibrium (QLE) manifold. On this manifold,
sustained cycling is possible as the dynamics approach a heteroclinic cycle, in
which allele frequencies alternate between near extinction and near fixation.
We find a closed-form approximation for the QLE manifold and use it to simplify
the model. For the simplified model, we can analytically calculate the
stability of the heteroclinic cycle. In the discrete-time model the cycle is
always stable; in a continuous-time approximation, the cycle is always
unstable. This demonstrates that complex dynamics are possible under
quasi-linkage equilibrium.


  RNA-Seq analysis has revolutionized researchers' understanding of the
transcriptome in biological research. Assessing the differences in
transcriptomic profiles between tissue samples or patient groups enables
researchers to explore the underlying biological impact of transcription.
RNA-Seq analysis requires multiple processing steps and huge computational
capabilities. There are many well-developed R packages for individual steps;
however, there are few R/Bioconductor packages that integrate existing software
tools into a comprehensive RNA-Seq analysis and provide fundamental end-to-end
results in pure R environment so that researchers can quickly and easily get
fundamental information in big sequencing data. To address this need, we have
developed the open source R/Bioconductor package, RNASeqR. It allows users to
run an automated RNA-Seq analysis with only six steps, producing essential
tabular and graphical results for further biological interpretation. The
features of RNASeqR include: six-step analysis, comprehensive visualization,
background execution version, and the integration of both R and command-line
software. RNASeqR provides fast, light-weight, and easy-to-run RNA-Seq analysis
pipeline in pure R environment. It allows users to efficiently utilize popular
software tools, including both R/Bioconductor and command-line tools, without
predefining the resources or environments. RNASeqR is freely available for
Linux and macOS operating systems from Bioconductor
(https://bioconductor.org/packages/release/bioc/html/RNASeqR.html).


  Tumor microenvironment has complex effects on tumorigenesis and metastasis.
However, there is still a lack of comprehensive understanding of the
relationship among molecular and cellular characteristics in tumor
microenvironment, clinical prognosis and immunotherpy response. In this study,
the immune and stromal (non-immune) signatures of tumor microenvironment were
integrated to identify novel subgroups of lung adenocarcinoma by
eigendecomposition and extraction algorithms of bioinformatics and machine
learning, such as non-negative matrix factorization and multitask learning.
Tumors were classified into 4 groups according to the activation of immunity
and stroma by novel signatures. The 4 groups had different mutation landscape,
molecular, cellular characteristics and prognosis, which have been validation
in 6 independent data sets containing 1551 patients. High-immune and
low-stromal activation group links to high immunocyte infiltration, high
immunocompetence, low fibroblasts, endothelial cells, collagen, laminin, tumor
mutation burden, and better overall survival. We developed a novel model based
on tumor microenvironment by integrating immune and stromal activation, namely
PMBT (prognostic model based on tumor microenvironment). The PMBT showed the
value to predict overall survival and immunotherapy responses.


  The Moran discrete process and the Wright-Fisher modelare the most popular
models in population genetics. It is common tounderstand the dynamics of these
models to use an approximating diffusionprocess, called Wright-Fisher
diffusion. Here, we give a quantitativelarge population limit of the error
committed by using the approximationdiffusion in the presence of weak selection
and weak immigrationin one dimension. The approach is robust enough to consider
the casewhere selection and immigration are Markovian processes, with
limitsjump or diffusion processes.


  Gender differences in human behaviour have attracted generations of social
scientists, who have explored whether males and females act differently in
domains involving competition, risk taking, cooperation, altruism, honesty, as
well as many others. Yet, little is known about gender differences in the
equity-efficiency trade-off. It has been suggested that females are more
equitable than males, but the empirical evidence is weak and inconclusive. This
gap is particularly important, because people in power of redistributing
resources often face a conflict between equity and efficiency, to the point
that this trade-off has been named as "the central problem in distributive
justice". The recently introduced Trade-Off Game (TOG) - in which a
decision-maker has to unilaterally choose between being equitable or being
efficient - offers a unique opportunity to fill this gap. To this end, I
analyse gender differences on a large dataset including all N=5,056 TOG
decisions collected by my research group since we introduced this game. The
results show that females prefer equity over efficiency to a greater extent
than males do. These findings suggest that males and females have different
preferences for resource distribution, and point to new avenues for future
research.


  Cilia are ubiquitous organelles involves in eukaryotic motility. They are
long, slender, and motile protrusions from the cell body. They undergo active
regular oscillatory beating patterns that can propel cells, such as the algae
Chlamydomonas, through fluids. When many cilia beat in synchrony they can also
propel fluid along the surfaces of cells, as is the case of nodal cilia.
  The main structural elements inside the cilium are microtubules. There are
also molecular motors of the dynein family that actively power the motion of
the cilium. These motors transform chemical energy in the form of ATP into
mechanical forces that produce sliding displacement between the microtubules.
This sliding is converted to bending by constraints at the base and/or along
the length of the cilium. Forces and displacements within the cilium can
regulate dyneins and provide a feedback mechanism: the dyneins generate forces,
deforming the cilium; the deformations, in turn, regulate the dyneins. This
feedback is believed to be the origin of the coordination of dyneins in space
and time which underlies the regularity of the beat pattern.
  Goals and approach. While the mechanism by which dyneins bend the cilium is
understood, the feedback mechanism is much less clear. The two key questions
are: which forces and displacements are the most relevant in regulating the
beat? and how exactly does this regulation occur?
  In this thesis we develop a framework to describe the spatio-temporal
patterns of a cilium with different mechanisms of motor regulation.
Characterizing and comparing the predicted shapes and beat patterns of these
different mechanisms to those observed in experiments provides us with further
understanding on how dyneins are regulated. This comparison is done both, with
a linear model that can be analytically solved, as with a non-linear model that
we solve numerically.


  We present a new whole-genome based approach to infer large-scale phylogenies
that is alignment- and reference-free. In contrast to other methods, it does
not rely on pairwise comparisons to determine distances to infer edges in a
tree. Instead, a colored de-Bruijn graph is constructed, and information on
common subsequences is extracted to infer phylogenetic splits. Application to
different datasets confirms robustness of the approach. A comparison to other
state-of-the-art whole-genome based methods indicates comparable or higher
accuracy and efficiency.


  The coding space of protein sequences is shaped by evolutionary constraints
set by requirements of function and stability. We show that the coding space of
a given protein family--the total number of sequences in that family--can be
estimated using models of maximum entropy trained on multiple sequence
alignments of naturally occuring amino acid sequences. We analyzed and
calculated the size of three abundant repeat proteins families, whose members
are large proteins made of many repetitions of conserved portions of ~ 30 amino
acids. While amino acid conservation at each position of the alignment explains
most of the reduction of diversity relative to completely random sequences, we
found that correlations between amino acid usage at different positions
significantly impact that diversity. We quantified the impact of different
types of correlations, functional and evolutionary, on sequence diversity.
Analysis of the detailed structure of the coding space of the families revealed
a rugged landscape, with many local energy minima of varying sizes with a
hierarchical structure, reminiscent of fustrated energy landscapes of spin
glass in physics. This clustered structure indicates a multiplicity of subtypes
within each family, and suggests new strategies for protein design.


  Given a set of species whose evolution is represented by a species tree, a
gene family is a group of genes having evolved from a single ancestral gene. A
gene family evolves along the branches of a species tree through various
mechanisms, including - but not limited to - speciation, gene duplication, gene
loss, horizontal gene transfer. The reconstruction of a gene tree representing
the evolution of a gene family constrained by a species tree is an important
problem in phylogenomics. However, unlike in the multispecies coalescent
evolutionary model, very little is known about the search space for gene family
histories accounting for gene duplication, gene loss and horizontal gene
transfer (the DLT-model). We introduce the notion of evolutionary histories
defined as a binary ordered rooted tree describing the evolution of a gene
family, constrained by a species tree in the DLT-model. We provide formal
grammars describing the set of all evolutionary histories that are compatible
with a given species tree, whether it is ranked or unranked. These grammars
allow us, using either analytic combinatorics or dynamic programming, to
efficiently compute the number of histories of a given size, and also to
generate random histories of a given size under the uniform distribution. We
apply these tools to obtain exact asymptotics for the number of gene family
histories for two species trees, the rooted caterpillar and the complete binary
tree, as well as estimates of the range of the exponential growth factor of the
number of histories for random species trees of size up to 25. Our results show
that including horizontal gene transfer induce a dramatic increase of the
number of evolutionary histories. We also show that, within ranked species
trees, the number of evolutionary histories in the DLT-model is almost
independent of the species tree topology.


  The interaction between aphids, ants and ladybirds has been investigated from
an ecological point of view since many decades, while there are no attempts to
describe it from a mathematical point of view. This paper introduces a new
mathematical model to describe the within-season population dynamics in an
ecological patch of a system composed by aphids, ants and ladybirds, through a
set of four differential equations. The proposed model is based on the
Kindlmann and Dixon set of differential equations, focused on the prediction of
the aphids-ladybirds population densities, that share a prey-predator
relationship. The population of ants, in mutualistic relationship with aphids
and in interspecific competition with ladybirds, is described according to the
Holland and De Angelis mathematical model, in which the authors faced the
problem of mutualistic interactions in general terms. The set of differential
equations proposed here is discretized by means the Nonstandard Finite
Difference scheme, successfully applied by Gabbriellini to the mutualistic
model. The constructed finite-difference scheme is positivity-preserving and
characterized by four nonhyperbolic steady-states, as highlighted by the
phase-space and time-series analyses. Particular attention is dedicated to the
steady-state most interesting from an ecological point of view, whose
asymptotic stability is demonstrated via the Centre Manifold Theory. The model
allows to numerically confirm that mutualistic relationship effectively
influences the population dynamic, by increasing the peaks of the aphids and
ants population densities. Nonetheless, it is showed that the asymptotical
populations of aphids and ladybirds collapse for any initial condition, unlike
that of ants that, after the peak, settle on a constant asymptotic value.


  We investigate diffusive search on planar networks, motivated by tubular
organelle networks in cell biology that contain molecules searching for
reaction partners and binding sites. Exact calculation of the diffusive mean
first-passage time on a spatial network is used to characterize the typical
search time as a function of network connectivity. We find that global
structural properties --- the total edge length and number of loops --- are
sufficient to largely determine network exploration times for a variety of both
synthetic planar networks and organelle morphologies extracted from living
cells. For synthetic networks on a lattice, we predict the search time
dependence on these global structural parameters by connecting with percolation
theory, providing a bridge from irregular real-world networks to a simpler
physical model. The dependence of search time on global network structural
properties suggests that network architecture can be designed for efficient
search without controlling the precise arrangement of connections.
Specifically, increasing the number of loops substantially decreases search
times, pointing to a potential physical mechanism for regulating reaction rates
within organelle network structures.


  We study a variant of the cyclic Lotka-Volterra model with three-agent
interactions. Inspired by a multiplayer variation of the Rock-Paper-Scissors
game, the model describes an ideal ecosystem in which cyclic competition among
three species develops through cooperative predation. Its rate equations in a
well-mixed environment display a degenerate Hopf bifurcation, occurring as
reactions involving two predators plus one prey have the same rate as reactions
involving two preys plus one predator. We estimate the magnitude of the
stochastic noise at the bifurcation point, where finite size effects turn
neutrally stable orbits into erratically diverging trajectories. In particular,
we compare analytic predictions for the extinction probability, derived in the
Fokker-Planck approximation, with numerical simulations based on the Gillespie
stochastic algorithm. We then extend the analysis of the phase portrait to
heterogeneous rates. In a well-mixed environment, we observe a continuum of
degenerate Hopf bifurcations, generalizing the above one. Neutral stability
ensues from a complex equilibrium between different reactions. Remarkably, on a
two-dimensional lattice, all bifurcations disappear as a consequence of the
spatial locality of the interactions. In the second part of the paper, we
investigate the effects of mobility in a lattice metapopulation model with
patches hosting several agents. We find that strategies propagate along the
arms of rotating spirals, as they usually do in models of cyclic dominance. We
observe propagation instabilities in the regime of large wavelengths. We also
examine three-agent interactions inducing nonlinear diffusion.


  In this study, considering the importance of how to exploit renewable natural
resources, we analyze a fishing model with nonlinear harvesting function in
which the players at the equilibrium point do a static game with complete
information that, according to the calculations, will cause a waste of energy
for both players and so the selection of cooperative strategies along with the
agreement between the players is the result of this research.


  Oscillatory gene circuits are ubiquitous to biology and are involved in
fundamental processes of cell cycle, circadian rhythms and developmental
systems. The synthesis of small, non-natural oscillatory genetic circuits have
been increasingly used to test fundamental principles of genetic network
dynamics. A recently developed fast, tunable genetic oscillator by Stricker et
al.[23] has demonstrated robustness and tunability of oscillatory behavior by
combining positive and negative feedback loops. This oscillator combining lacI
(negative) and araC (positive) feedback loops, was however modeled using
multiple layers of differential equations to capture the molecular complexity
of regulation, in order to explain the experimentally measured oscillations. We
have developed a reduced model based on delay differential equations (DDEs) of
this dual feedback loop oscillator, that reproduces the tunability of
oscillator period and amplitude based on the concentration of the two inducers
isopropyl b-D-1-thiogalactopyranoside (IPTG) and arabinose. Previous work had
predicted a need for an asymmetry in copy numbers of activator (araC) and
repressor (lacI) genes encoded on plasmids. We use our reduced model to
redesign the network by comparing the effect of asymmetry in gene expression at
the level of (a) DNA copy numbers and the rates of (b) mRNA translation and (c)
degradation. We find the minimal period of the oscillator is sensitive to DNA
copy number asymmetry, but translation rate asymmetry has an identical effect
as plasmid copy numbers, while modulating the asymmetry in mRNA degradation can
improve the tunability of period of the oscillator, together with increased
robustness to replication 'noise' and influence of the host cell cycle. Thus,
our model predicts experimentally testable principles to redesign a potentially
more robust oscillatory genetic network.


  Mutation is a critical mechanism by which evolution explores the functional
landscape of proteins. Despite our ability to experimentally inflict mutations
at will, it remains difficult to link sequence-level perturbations to
systems-level responses. Here, we present a framework centered on measuring
changes in the free energy of the system to link individual mutations in an
allosteric transcriptional repressor to the parameters which govern its
response. We find the energetic effects of the mutations can be categorized
into several classes which have characteristic curves as a function of the
inducer concentration. We experimentally test these diagnostic predictions
using the well-characterized LacI repressor of Escherichia coli, probing
several mutations in the DNA binding and inducer binding domains. We find that
the change in gene expression due to a point mutation can be captured by
modifying only a subset of the model parameters that describe the respective
domain of the wild-type protein. These parameters appear to be insulated, with
mutations in the DNA binding domain altering only the DNA affinity and those in
the inducer binding domain altering only the allosteric parameters. Changing
these subsets of parameters tunes the free energy of the system in a way that
is concordant with theoretical expectations. Finally, we show that the
induction profiles and resulting free energies associated with pairwise double
mutants can be predicted with quantitative accuracy given knowledge of the
single mutants, providing an avenue for identifying and quantifying epistatic
interactions.


  As human spaceflight seeks to expand beyond low-Earth orbit, NASA and its
international partners face numerous challenges related to ensuring the safety
of their astronauts, including the need to provide a safe and effective
pharmacy for long-duration spaceflight. Historical missions have relied upon
frequent resupply of onboard pharmaceuticals; as a result, there has been
little study into the effects of long-term exposure of pharmaceuticals to the
space environment. Of particular concern are the long-term effects of space
radiation on drug stability, especially as missions venture away from the
protective proximity of the Earth. Here we highlight the risk of space
radiation to pharmaceuticals during exploration spaceflight, identifying the
limitations of current understanding. We further seek to identify ways in which
these limitations could be addressed through dedicated research efforts aimed
towards the rapid development of an effective pharmacy for future spaceflight
endeavors.


  Comparative graph and network analysis play an important role in both systems
biology and pattern recognition, but existing surveys on the topic have
historically ignored or underserved one or the other of these fields. We
present an integrative introduction to the key objectives and methods of graph
and network comparison in each field, with the intent of remaining accessible
to relative novices in order to mitigate the barrier to interdisciplinary idea
crossover.
  To guide our investigation, and to quantitatively justify our assertions
about what the key objectives and methods of each field are, we have
constructed a citation network containing 5,793 vertices from the full
reference lists of over two hundred relevant papers, which we collected by
searching Google Scholar for ten different network comparison-related search
terms. We investigate its basic statistics and community structure, and frame
our presentation around the papers found to have high importance according to
five different standard centrality measures.


  Non-coding RNA (ncRNA) are RNA sequences which don't code for a gene but
instead carry important biological functions. The task of ncRNA classification
consists in classifying a given ncRNA sequence into its family. While it has
been shown that the graph structure of an ncRNA sequence folding is of great
importance for the prediction of its family, current methods make use of
machine learning classifiers on hand-crafted graph features. We improve on the
state-of-the-art for this task with a graph convolutional network model which
achieves an accuracy of 85.73% and an F1-score of 85.61% over 13 classes.
Moreover, our model learns in an end-to-end fashion from the raw RNA graphs and
removes the need for expensive feature extraction. To the best of our
knowledge, this also represents the first successful application of graph
convolutional networks to RNA folding data.


  Slime mould plasmodia can adjust their behaviour in response to chemical
trails left by themselves and other Physarum plasmodia. This simple feedback
process increases their foraging efficiency. We still do not know whether other
factors influence plasmodium behaviour in realistic competition settings. Here
we designed a competition experiment where two plasmodia had to find one food
source in a common environment. As previously shown, the time it took plasmodia
to find food depended on their hunger motivation. However, the time it took a
plasmodium to start looking for food depended on its motivation and the
motivation of its competitor. Plasmodia always initiated foraging quicker if
they were in the presence of a competitor and the quickest if they were hungry
and in the presence of a satiated competitor. The time it took to arrive to the
food was not influenced by whether they were alone or with a competitor.
Ultimately, this complex competition response benefited the hungry plasmodia as
they had a 4:1 chance of finding the food first. The sensory ecology of
Physarum polycephalum is more complex than previously thought and yields
complex behaviour in a simple organism.


  Body size or mass is one of the main factors underlying food webs structure.
A large number of evolutionary models have shown that indeed, the adaptive
evolution of body size (or mass) can give rise to hierarchically organised
trophic levels with complex between and within trophic interactions. However,
these models generally make strong arbitrary assumptions on how traits evolve,
casting doubts on their robustness. In particular, biomass conversion
efficiency is always considered independent of the predator and prey size,
which contradicts with the literature. In this paper, we propose a general
model encompassing most previous models which allows to show that relaxing
arbitrary assumptions gives rise to unrealistic food webs. We then show that
considering biomass conversion efficiency dependent on species size is
certainly key for food webs adaptive evolution because realistic food webs can
evolve, making obsolete the need of arbitrary constraints on traits' evolution.
We finally conclude that, on the one hand, ecologists should pay attention to
how biomass flows into food webs in models. On the other hand, we question more
generally the robustness of evolutionary models for the study of food webs.


  Species networks generalize the notion of species trees to allow for
hybridization or other lateral gene transfer. Under the Network Multispecies
Coalescent Model, individual gene trees arising from a network can have any
topology, but arise with frequencies dependent on the network structure and
numerical parameters. We propose a new algorithm for statistical inference of a
level-1 species network under this model, from data consisting of gene tree
topologies, and provide the theoretical justification for it. The algorithm is
based in an analysis of quartets displayed on gene trees, combining several
statistical hypothesis tests with combinatorial ideas such as a quartet-based
intertaxon distance appropriate to networks, the NeighborNet algorithm for
circular split systems, and the Circular Network algorithm for constructing a
splits graph.


  Differential diffusion is a source of instability in population dynamics
systems when species diffuse with different rates. Predator-prey systems show
this instability only under certain specific conditions, usually requiring
Holling-type functionals involved. Here we study the effects of intraspecific
cooperation and competition on diffusion-driven instability in a predator-prey
system with a different structure. We conduct the analysis on a generalized
population dynamics that bounds intraspecific and interspecific interactions
with Verhulst-type saturation terms instead of Holling-type functionals. We
find that instability occurs due to the intraspecific saturation or
intraspecific interactions, both cooperative and competitive. We present
numerical simulations and show spatial patterns due to diffusion.


  Diverse physical systems are characterized by their response to small
perturbations. Near thermodynamic equilibrium, the fluctuation-dissipation
theorem provides a powerful theoretical and experimental tool to determine the
nature of response by observing spontaneous equilibrium fluctuations. In this
spirit, we derive here a collection of equalities and inequalities valid
arbitrarily far from equilibrium that constrain the response of nonequilibrium
steady states in terms of the strength of nonequilibrium driving. Our work
opens new avenues for characterizing nonequilibrium response. As illustrations,
we show how our results rationalize the energetic requirements of two common
biochemical motifs.


  We report the effect of conformational dynamics on the fluctuations of
electric fields in the active site of the enzyme Ketosteroid Isomerase (KSI).
While KSI is considered rigid with little conformational variation to support
different stages of the catalytic cycle, we show that KSI utilizes cooperative
side chain motions of the entire protein scaffold outside the active site,
which contribute negligibly to the electric fields on the substrate, by
progressively stabilizing electric field contributions by particular active
site residues at different timescales. The design of synthetic enzymes could
benefit from strategies that can take advantage of the dynamics by using
electric fields fluctuations as a guide.


  The observed architecture of ecological and socio-economic networks differs
significantly from that of random networks. From a network science standpoint,
non-random structural patterns observed in real networks call for an
explanation of their emergence and an understanding of their potential systemic
consequences. This article focuses on one of these patterns: nestedness. Given
a network of interacting nodes, nestedness can be described as the tendency for
nodes to interact with subsets of the interaction partners of better-connected
nodes. Known since more than $80$ years in biogeography, nestedness has been
found in systems as diverse as ecological mutualistic organizations, world
trade, inter-organizational relations, among many others. This review article
focuses on three main pillars: the existing methodologies to observe nestedness
in networks; the main theoretical mechanisms conceived to explain the emergence
of nestedness in ecological and socio-economic networks; the implications of a
nested topology of interactions for the stability and feasibility of a given
interacting system. We survey results from variegated disciplines, including
statistical physics, graph theory, ecology, and theoretical economics.
Nestedness was found to emerge both in bipartite networks and, more recently,
in unipartite ones; this review is the first comprehensive attempt to unify
both streams of studies, usually disconnected from each other. We believe that
the truly interdisciplinary endeavour -- while rooted in a complex systems
perspective -- may inspire new models and algorithms whose realm of application
will undoubtedly transcend disciplinary boundaries.


  For much of biology, the manner in which genotype maps to phenotype remains a
fundamental mystery. The few maps that are known tend to show modular
pleiotropy: sets of phenotypes are determined by distinct sets of genes. One
key map that has evaded discovery is that of the human brain's network
architecture. Here, we determine the form of this map for gene coexpression and
single nucleotide polymorphisms. We discover that mostly non-overlapping sets
of genes encode the connectivity of brain network modules (or so-called
communities), suggesting that brain network communities demarcate genetic
transitions. We find that these clean boundaries break down at connector hubs,
whose integrative connectivity is encoded by pleiotropic genes from mostly
non-overlapping sets. Broadly, this study opens fundamentally new directions in
the study of genetic encoding of brain development, evolution, and disease.


  RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and
their three-dimensional (3D) structures and stability play essential roles in a
variety of biological functions. Therefore, to predict 3D structures and
stability of RNA pseudoknots is essential for understanding their functions. In
the work, we employed our previously developed coarse-grained model with
implicit salt to make extensive predictions and comprehensive analyses on the
3D structures and stability for RNA pseudoknots in monovalent/divalent ion
solutions. The comparisons with available experimental data show that our model
can successfully predict the 3D structures of RNA pseudoknots from their
sequences, and can also make reliable predictions for the stability of RNA
pseudoknots with different lengths and sequences over a wide range of
monovalent/divalent ion concentrations. Furthermore, we made comprehensive
analyses on the unfolding pathway for various RNA pseudoknots in ion solutions.
Our analyses for extensive pseudokonts and the wide range of
monovalent/divalent ion concentrations verify that the unfolding pathway of RNA
pseudoknots is mainly dependent on the relative stability of unfolded
intermediate states, and show that the unfolding pathway of RNA pseudoknots can
be significantly modulated by their sequences and solution ion conditions.


  The multi-population replicator dynamics (RD) can be considered a dynamic
approach to the study of multi-player games, where it was shown to be related
to Cross' learning, as well as of systems of coevolving populations. However,
not all of its equilibria are Nash equilibria (NE) of the underlying game, and
neither convergence to an NE nor convergence in general are guaranteed.
Although interior equilibria are guaranteed to be NE, no interior equilibrium
can be asymptotically stable in the multi-population RD, resulting, e.g., in
cyclic orbits around a single interior NE. We introduce a new notion of
equilibria of RD, called mutation limits, which is based on the inclusion of a
naturally arising, simple form of mutation, but is invariant under the specific
choice of mutation parameters. We prove the existence of such mutation limits
for a large range of games, and consider a subclass of particular interest,
that of attracting mutation limits. Attracting mutation limits are approximated
by asymptotically stable equilibria of the (mutation-)perturbed RD, and hence,
offer an approximate dynamic solution of the underlying game, especially if the
original dynamic has no asymptotically stable equilibria. In this sense,
mutation stabilises the system in certain cases and makes attracting mutation
limits near-attainable. Furthermore, the relevance of attracting mutation
limits as a game theoretic equilibrium concept is emphasised by a similarity of
(mutation-)perturbed RD to the Q-learning algorithm in the context of
multi-agent reinforcement learning. In contrast to the guaranteed existence of
mutation limits, attracting mutation limits do not exist in all games, raising
the question of their characterization.


  Abeta is a disordered peptide central to Alzheimer's Disease. Aggregation of
Abeta has been widely explored, but its molecular crowding less so. The
synaptic cleft where Abeta locates only holds 60-70 water molecules along its
width. We subjected Abeta40 to 100 different simulations with variable water
cell size. We show that even for this disordered aggregation-prone peptide,
many properties are not cell-size dependent, i.e. a small cell is easily
justified. The radius of gyration, intra-peptide, and peptide-water hydrogen
bonds are well-sampled by short (50 ns) time scales at any cell size. Abeta is
mainly disordered with 0-30% alpha helix but undergoes consistent alpha-beta
transitions up to 14% strand in 5-10% of the simulations regardless of cell
size. The similar prevalence in long and short simulations indicate small
diffusion barriers for structural transitions in contrast to folded globular
proteins, which we suggest is a defining hallmark of intrinsically disordered
proteins. Importantly, the hydrophobic surface increases significantly in small
cells (confidence level 95%, two-tailed t-test), as does the variation in
exposure and backbone conformations (>40% and >27% increased standard
deviations). Whereas hydrophilic exposure dominates hydrophobic exposure in
large cells, this tendency breaks down at low water concentration. We interpret
these findings as a concentration-dependent hydrophobic effect, with the small
water layer unable to keep the protein unexposed, an effect mainly caused by
the layered water-water interactions, not by the peptide dynamics. The exposure
correlates with radius of gyration (R2 0.35-0.50) and could be important in
crowded environments, e.g. the synaptic cleft.


  Protein phosphorylation cycles are important mechanisms of the post
translational modification of a protein and as such an integral part of
intracellular signaling and control. We consider the sequential phosphorylation
and dephosphorylation of a protein at two binding sites. While it is known that
proteins where phosphorylation is processive and dephosphorylation is
distributive admit oscillations (for some value of the rate constants and total
concentrations) it is not known whether or not this is the case if both
phosphorylation and dephosphorylation are distributive. We study four
simplified mass action models of sequential and distributive phosphorylation
and show that for each of those there do not exist rate constants and total
concentrations where a Hopf bifurcation occurs. To arrive at this result we use
convex parameters to parameterize the steady state and Hurwitz matrices.


  We derive a reduction formula for singularly perturbed ordinary differential
equations (in the sense of Tikhonov and Fenichel) with a known parameterization
of the critical manifold. No a priori assumptions concerning separation of slow
and fast variables are made, or necessary.We apply the theoretical results to
chemical reaction networks with mass action kinetics admitting slow and fast
reactions. For some relevant classes of such systems there exist canonical
parameterizations of the variety of stationary points, hence the theory is
applicable in a natural manner. In particular we obtain a closed form
expression for the reduced system when the fast subsystem admits complex
balanced steady states.


  Motivation: Due to the central role of protein structure in molecular
recognition, great computational efforts are devoted to modeling protein
structures and motions that mediate structural rearrangements. The size,
dimensionality, and non-linearity of the protein structure space present
outstanding challenges. Such challenges also arise in robot motion planning,
and robotics-inspired treatments of protein structure and motion are
increasingly showing high exploration capability. Encouraged by such findings,
we debut here ROMEO, which stands for Robotics prOtein Motion ExplOration
framework. ROMEO is an open-source, object-oriented platform that allows
researchers access to and reproducibility of published robotics-inspired
algorithms for modeling protein structures and motions, as well as facilitates
novel algorithmic design via its plug-and-play architecture.
  Availability and implementation: ROMEO is written in C++ and is available in
GitLab (https://github.com/). This software is freely available under the
Creative Commons license (Attribution and Non-Commercial).
  Contact: amarda@gmu.edu


  Networks of filamentous actin (F-actin) are important for the mechanics of
most animal cells. These cytoskeletal networks are highly dynamic, with a
variety of actin-associated proteins that control cross-linking, polymerization
and force generation in the cytoskeleton. Inspired by recent rheological
experiments on reconstituted solutions of dynamic actin filaments, we report a
theoretical model that describes stress relaxation behavior of these solutions
in the presence of severing proteins. We show that depending on the kinetic
rates of assembly, disassembly, and severing, one can observe both
length-dependent and length-independent relaxation behavior.


  Comprehensive and unambiguous identification of small molecules in complex
samples will revolutionize our understanding of the role of metabolites in
biological systems. Existing and emerging technologies have enabled measurement
of chemical properties of molecules in complex mixtures and, in concert, are
sensitive enough to resolve even stereoisomers. Despite these experimental
advances, small molecule identification is inhibited by (i) chemical reference
libraries representing <1% of known molecules, limiting the number of possible
identifications, and (ii) the lack of a method to generate candidate matches
directly from experimental features (i.e. without a library). To this end, we
developed a variational autoencoder (VAE) to learn a continuous numerical, or
latent, representation of molecular structure to expand reference libraries for
small molecule identification. We extended the VAE to include a chemical
property decoder, trained as a multitask network, in order to shape the latent
representation such that it assembles according to desired chemical properties.
The approach is unique in its application to small molecule identification,
with its focus on m/z and CCS, paired with its training paradigm, which
involved a cascade of transfer learning iterations. This allows the network to
learn as much as possible at each stage, enabling success with progressively
smaller datasets without overfitting. Once trained, the network can rapidly
predict chemical properties directly from structure, as well as generate
candidate structures with desired chemical properties. Additionally, the
ability to generate novel molecules along manifolds, defined by chemical
property analogues, positions DarkChem as highly useful in a number of
application areas, including metabolomics and small molecule identification,
drug discovery and design, chemical forensics, and beyond.


  The selectivity of L-lysine acylation by lauric acid catalysed by Candida
antarctica lipase B (CALB) was investigated combining experimental and
theoretical methodologies. Experiments showed the near-exclusive acylation of
lysine $\epsilon$-amino group; only traces of product resulting from the
acylation of lysine $\alpha$-amino group were observed fleetingly. Molecular
modelling simulations were performed aiming to understand the molecular rules
for selectivity. Flexible docking simulations combined with structural
investigations into lysine/CALB binding modes also suggested the preferential
acylation of lysine $\epsilon$-amino group without, however, excluding the
acylation of the lysine $\alpha$-amino group. Electrostatic interaction energy
between lysine and the residues covering the catalytic cavity was calculated in
order to understand the discrimination between the two lysine amino groups. The
results suggests that the proximity of the carboxylate group hinders the
binding of the substrate in configurations enabling the N$\alpha$-acylation.
Key interactions with the polar region covering the catalytic triad were
identified and a plausible explanation for selectivity was proposed.


  Human contact networks are constituted by a multitude of individuals and
pairwise contacts among them. However, the dynamic nature, which generates the
evolution of human contact networks, of contact patterns is not known yet.
Here, we analyse three empirical datasets and identify two crucial mechanisms
of the evolution of temporal human contact networks, i.e. the activity state
transition laws for an individual to be socially active, and the contact
establishment mechanism that active individuals adopt. We consider both of the
two mechanisms to propose a temporal network model, named the memory driven
(MD) model, of human contact networks. Then we study the susceptible-infected
(SI) spreading processes on empirical human contact networks and four
corresponding temporal network models, and compare the full prevalence time of
SI processes with various infection rates on the networks. The full prevalence
time of SI processes in the MD model is the same as that in real-world human
contact networks. Moreover, we find that the individual activity state
transition promotes the spreading process, while, the contact establishment of
active individuals suppress the prevalence. Apart from this, we observe that
even a small percentage of individuals to explore new social ties is able to
induce an explosive spreading on networks. The proposed temporal network
framework could help the further study of dynamic processes in temporal human
contact networks, and offer new insights to predict and control the diffusion
processes on networks.


  Estimates of population abundance are fundamental to wildlife management and
conservation, but are difficult to obtain across large geographic scales and
for cryptic species. We used a state-space model with age-at-harvest data in a
Bayesian framework to model American black bear (Ursus americanus) abundance
and demographic parameters in four management zones in Wisconsin from
2011-2017. We had limited demographic data available from the population, and
relied upon a) the model, b) age-at-harvest data, and c) informative prior
distributions from a literature review. The estimated posterior means and
distributions for abundance and demographic parameters from our models were
reasonable for each management zone, and indicated a decreasing trend in zones
A and B, and a generally stable trend in zones C and D. The age-at-harvest data
updated the posterior distribution and means for initial population size,
harvest season survival, and non-harvest season survival, with a notable
increase in precision for the survival values. A strength of the model for
managers is the formalized process for providing biologically supported
information as prior distributions to transparently accommodate expert opinion
with measures of confidence when estimating wildlife populations, which can
then be updated by the age-at-harvest data and model structure. The integration
of prior information and age-at-harvest state-space models in a Bayesian
framework efficiently leverage all available information for making
zone-specific abundance estimates for the management of harvested species. This
may create more informative data for decision-makers when setting harvest
quotas, and could lead to more effective monitoring, conservation, and
management of cryptic carnivore species.


  Bacteriorhodopsin (bR) is a light-driven proton pump. We use time-resolved
crystallography at an X-ray free-electron laser to follow the structural
changes in multiphoton-excited bR from 250 femtoseconds to 10 picoseconds.
Quantum chemistry and ultrafast spectroscopy allow identifying a sequential
two-photon absorption process, leading to excitation of a tryptophan residue
flanking the retinal chromophore, as a first manifestation of multi-photon
effects. We resolve distinct stages in the structural dynamics of the all-trans
retinal in photoexcited bR to a highly twisted 13-cis conformation. Other
active site sub-picosecond rearrangements include correlated vibrational
motions of the electronically excited retinal chromophore, the surrounding
amino acids and water molecules as well as their hydrogen bonding network.
These results show that this extended photo-active network forms an
electronically and vibrationally coupled system in bR, and most likely in all
retinal proteins.


  Optimising discrete data for a desired characteristic using gradient-based
methods involves projecting the data into a continuous latent space and
carrying out optimisation in this space. Carrying out global optimisation is
difficult as optimisers are likely to follow gradients into regions of the
latent space that the model has not been exposed to during training; samples
generated from these regions are likely to be too dissimilar to the training
data to be useful. We propose Constrained Optimisation with Latent
Distributions (COLD), a constrained global optimisation procedure to find
samples with high values of a desired property that are similar to yet distinct
from the training data. We find that on MNIST, our procedure yields optima for
each of three different objectives, and that enforcing tighter constraints
improves the quality and increases the diversity of the generated images. On
the ChEMBL molecular dataset, our method generates a diverse set of new
molecules with drug-likeness scores similar to those of the highest-scoring
molecules in the training data. We also demonstrate a computationally efficient
way to approximate the constraint when evaluating it exactly is computationally
expensive.


  A guiding principle in self-assembly is that, for high production yield,
nucleation of structures must be significantly slower than their growth.
However, details of the mechanism that impedes nucleation are broadly
considered irrelevant. Here, we analyze self-assembly into finite-sized target
structures employing mathematical modeling. We investigate two key scenarios to
delay nucleation: (i) by introducing a slow activation step for the assembling
constituents and, (ii) by decreasing the dimerization rate. These scenarios
have widely different characteristics. While the dimerization scenario exhibits
robust behavior, the activation scenario is highly sensitive to demographic
fluctuations. These demographic fluctuations ultimately disfavor growth
compared to nucleation and can suppress yield completely. The occurrence of
this stochastic yield catastrophe does not depend on model details but is
generic as soon as number fluctuations between constituents are taken into
account. On a broader perspective, our results reveal that stochasticity is an
important limiting factor for self-assembly and that the specific
implementation of the nucleation process plays a significant role in
determining the yield.


  Analysis of large continuous-time stochastic systems is a computationally
intensive task. In this work we focus on population models arising from
chemical reaction networks (CRNs), which play a fundamental role in analysis
and design of biochemical systems. Many relevant CRNs are particularly
challenging for existing techniques due to complex dynamics including
stochasticity, stiffness or multimodal population distributions. We propose a
novel approach allowing not only to predict, but also to explain both the
transient and steady-state behaviour. It focuses on qualitative description of
the behaviour and aims at quantitative precision only in orders of magnitude.
Firstly, we abstract the CRN into a compact model preserving rough timing
information, distinguishing only signifcinatly different populations, but
capturing relevant sequences of behaviour. Secondly, we approximately analyse
the most probable temporal behaviours of the model through most probable
transitions. As demonstrated on complex CRNs from literature, our approach
reproduces the known results, but in contrast to the state-of-the-art methods,
it runs with virtually no computational cost and thus offers
unprecedented~scalability.


  The environment has a strong influence on a population's evolutionary
dynamics. Driven by both intrinsic and external factors, the environment is
subject to continual change in nature. To capture an ever-changing environment,
we consider a model of evolutionary dynamics with game transitions, where
individuals' behaviors together with the games they play in one time step
influence the games to be played next time step. Within this model, we study
the evolution of cooperation in structured populations and find a simple rule:
weak selection favors cooperation over defection if the ratio of the benefit
provided by an altruistic behavior, $b$, to the corresponding cost, $c$,
exceeds $k-k'$, where $k$ is the average number of neighbors of an individual
and $k'$ captures the effects of the game transitions. Even if cooperation
cannot be favored in each individual game, allowing for a transition to a
relatively valuable game after mutual cooperation and to a less valuable game
after defection can result in a favorable outcome for cooperation. In
particular, small variations in different games being played can promote
cooperation markedly. Our results suggest that simple game transitions can
serve as a mechanism for supporting prosocial behaviors in highly-connected
populations.


  We model evolution of plants in a world, made up of different locations, with
multiple environments (mutually exclusive and collectively exhaustive subsets
of locations). Each environment (landmass) has temperature, rainfall, and other
attributes that directly affect plant growth and reproduction. Each plant has
preferences for environment attributes. Depending on how suitable the
environment is to the plants, seeds are released or death occurs. With every
reproductive cycle, genetic mutations occur. To model competition, plants in
compete for survival, and success is stochastically dependent on environmental
fitness. Our model determines whether and how evolution occurs, and how the
attributes of plants change and possibly converge over time in relation to the
attributes of the environment.


  It is well established that gene expression can be modeled as a Markovian
stochastic process and hence proper observables might be subjected to large
fluctuations and rare events. Since dynamics is often more than statics, one
can work with ensembles of trajectories for long but fixed times, instead of
states or configurations, to study dynamics of these Markovian stochastic
processes and glean more information. In this paper we aim to show that the
concept of ensemble of trajectories can be applied to a variety of stochastic
models of gene expression ranging from a simple birth-death process to a more
sophisticate model containing burst and switch. By considering the protein
numbers as a relevant dynamical observable, apart from asymptotic behavior of
remote tails of probability distribution, generating function for the cumulants
of this observable can also be obtained. We discuss the unconditional
stochastic Markov processes which generate the statistics of rare events in
these models.


  In $\textit{Escherichia coli}$ bacterium, the molecular compounds involved in
protein synthesis, messenger RNAs (mRNAs) and ribosomes, show marked
intracellular localization patterns. Yet a quantitative understanding of the
physical principles which would allow one to control protein synthesis by
designing, bioengineering, and optimizing these localization patterns is still
lacking. In this study, we consider a scenario where a synthetic modification
of mRNA reaction-diffusion properties allows for controlling the localization
and stoichiometry of mRNAs and polysomes$\mathrm{-}$complexes of multiple
ribosomes bound to mRNAs. Our analysis demonstrates that protein synthesis can
be controlled, e.g., optimally enhanced or inhibited, by leveraging mRNA
spatial localization and stoichiometry only, without resorting to alterations
of mRNA expression levels. We identify the physical mechanisms that control the
protein-synthesis rate, highlighting the importance of colocalization between
mRNAs and freely diffusing ribosomes, and the interplay between polysome
stoichiometry and excluded-volume effects due to the DNA nucleoid. The
genome-wide, quantitative predictions of our work may allow for a direct
verification and implementation in cell-biology experiments, where localization
patterns and protein-synthesis rates may be monitored by fluorescence
microscopy in single cells and populations.


  Atomic Force Microscopy was utilized to study the morphology of Gag,
{\Psi}RNA, and their binding complexes with lipids in a solution environment
with 0.1{\AA} vertical and 1nm lateral resolution. TARpolyA RNA was used as a
RNA control. The lipid used was phospha-tidylinositol-(4,5)-bisphosphate
(PI(4,5)P2). The morphology of specific complexes Gag-{\Psi}RNA, Gag-TARpolyA
RNA, Gag-PI(4,5)P2 and PI(4,5)P2-{\Psi}RNA-Gag were studied. They were imaged
on either positively or negatively charged mica substrates depending on the net
charges carried. Gag and its complexes consist of monomers, dimers and
tetramers, which was confirmed by gel electrophoresis. The addition of specific
{\Psi}RNA to Gag is found to increase Gag multimerization. Non-specific
TARpolyA RNA was found not to lead to an increase in Gag multimerization. The
addition PI(4,5)P2 to Gag increases Gag multimerization, but to a lesser extent
than {\Psi}RNA. When both {\Psi}RNA and PI(4,5)P2 are present Gag undergoes
comformational changes and an even higher degree of multimerization.


  While cooperativity in ligand-induced receptor dimerization has been linked
with receptor-receptor couplings via minimal representations of physical
observables, effects arising from higher-order oligomer (e.g., trimer and
tetramer) formations of unobserved receptors have received less attention.
Here, we propose a dimerization model of ligand-induced receptors in
multivalent form representing physical observables under basis vectors of
various aggregated receptor-states. Our simulations of multivalent models not
only reject Wofsy-Goldstein parameter conditions for cooperativity, but show
higher-order oligomer formations can shift cooperativity from positive to
negative.


  We introduce here a multi-type bootstrap percolation model, which we call
T-Bootstrap Percolation (T-BP), and apply it to study information propagation
in social networks. In this model, a social network is represented by a graph G
whose vertices have different labels corresponding to the type of role the
person plays in the network (e.g. a student, an educator, etc.). Once an
initial set of vertices of G is randomly selected to be carrying a gossip (e.g.
to be infected), the gossip propagates to a new vertex provided it is
transmitted by a minimum threshold of vertices with different labels. By
considering random graphs, which have been shown to closely represent social
networks, we study different properties of the T-BP model through numerical
simulations, and describe its implications when applied to rumour spread, fake
news, and marketing strategies.


  This paper describes an investigation into part of the mechanical mechanisms
underlying the formation of mitotic spindle, the cellular machinery responsible
for chromosomal separation during cell division. In normal eukaryotic cells,
spindles are composed of microtubule filaments that radiate outward from two
centrosomes. In many transformed cells, however, centrosome number is
misregulated resulting in cells with more than two centrosomes. Addressing the
question of how these cells accommodate these additional structures by
coalescing supernumerary centrosomes to form normal spindles will provide a
powerful insight toward understanding the proliferation of cancer cells and
developing new therapeutics. The process of centrosome coalescence is thought
to involve motor proteins that function to slide microtubules relative to one
another. Here we use in vitro motility assays combined with fluorescence
microscopy to visualize, characterize and quantify microtubule-microtubule
interactions. After segmenting the microtubules, their speed and direction of
movement are the extracted features to cluster their interaction type. In order
to evaluate the potential of our processing algorithm, we created a simulated
dataset similar to the time-lapse series. Once our procedure has been optimized
using the simulated data, we will apply it to the real data. Results of our
analyses will provide a quantitative description of interaction among
microtubules. This is a potentially important step toward more thorough
understanding of cancer.


  Recent work in graph neural networks (GNNs) has led to improvements in
molecular activity and property prediction tasks. Unfortunately, GNNs often
fail to capture the relative importance of interactions between molecular
substructures, in part due to the absence of efficient intermediate pooling
steps. To address these issues, we propose LaPool (Laplacian Pooling), a novel,
data-driven, and interpretable hierarchical graph pooling method that takes
into account both node features and graph structure to improve molecular
representation. We benchmark LaPool on molecular graph prediction and
understanding tasks and show that it outperforms recent GNNs. Interestingly,
LaPool also remains competitive on non-molecular tasks. Both quantitative and
qualitative assessments are done to demonstrate LaPool's improved
interpretability and highlight its potential benefits in drug design. Finally,
we demonstrate LaPool's utility for the generation of valid and novel molecules
by incorporating it into an adversarial autoencoder.


  Native contacts between residues could be predicted from the amino acid
sequence of proteins, and the predicted contact information could assist the de
novo protein structure prediction. Here, we present a novel pipeline of a
residue contact predictor AmoebaContact and a contact-assisted folder GDFold
for rapid protein structure prediction. Unlike mainstream contact predictors
that utilize human-designed neural networks, AmoebaContact adopts a set of
network architectures that are found as optimal for contact prediction through
automatic searching and predicts the residue contacts at a series of cutoffs.
Different from conventional contact-assisted folders that only use top-scored
contact pairs, GDFold considers all residue pairs from the prediction results
of AmoebaContact in a differentiable loss function and optimizes the atom
coordinates using the gradient descent algorithm. Combination of AmoebaContact
and GDFold allows quick reconstruction of the protein structure, with
comparable model quality to the state-of-the-art protein structure prediction
methods.


  Living species, ranging from bacteria to animals, exist in environmental
conditions that exhibit spatial and temporal heterogeneity which requires them
to adapt. Risk-spreading through spontaneous phenotypic variations is a known
concept in ecology, which is used to explain how species may survive when faced
with the evolutionary risks associated with temporally varying environments. In
order to support a deeper understanding of the adaptive role of spontaneous
phenotypic variations in fluctuating environments, we consider a system of
non-local partial differential equations modelling the evolutionary dynamics of
two competing phenotype-structured populations in the presence of periodically
oscillating nutrient levels. The two populations undergo spontaneous phenotypic
variations at different rates. The phenotypic state of each individual is
represented by a continuous variable, and the phenotypic landscape of the
populations evolves in time due to variations in the nutrient level. Exploiting
the analytical tractability of our model, we study the long-time behaviour of
the solutions to obtain a detailed mathematical depiction of evolutionary
dynamics. The results suggest that when nutrient levels undergo small and slow
oscillations, it is evolutionarily more convenient to rarely undergo
spontaneous phenotypic variations. Conversely, under relatively large and fast
periodic oscillations in the nutrient levels, which bring about alternating
cycles of starvation and nutrient abundance, higher rates of spontaneous
phenotypic variations confer a competitive advantage. We discuss the
implications of our results in the context of cancer metabolism.