  The assembly of functional biomolecular condensates often involves
liquid-liquid phase separation (LLPS) of proteins with multiple modular
domains, which can be folded or conformationally disordered to various degrees.
To understand the LLPS-driving domain-domain interactions, a fundamental
question is how readily the interactions in the condensed phase can be inferred
from inter-domain interactions in dilute solutions. In particular, are the
interactions leading to LLPS exclusively those underlying the formation of
discrete inter-domain complexes in homogeneous solutions? We address this
question by developing a mean-field LLPS theory of two stoichiometrically
constrained solute species. The theory is applied to the neuronal proteins
SynGAP and PSD-95, whose complex coacervate serves as a rudimentary model for
neuronal postsynaptic densities (PSDs). The predicted phase behaviors are
compared with experiments. Previously, a three-SynGAP, two-PSD-95 ratio was
determined for SynGAP/PSD-95 complexes in dilute solutions. However, when this
3:2 stoichiometry is uniformly imposed in our theory encompassing both dilute
and condensed phases, the tie-line pattern of the predicted SynGAP/PSD-95 phase
diagram differs drastically from that obtained experimentally. In contrast,
theories embodying alternate scenarios postulating auxiliary SynGAP-PSD-95 as
well as SynGAP-SynGAP and PSD-95-PSD-95 interactions in addition to those
responsible for stoichiometric SynGAP/PSD-95 complexes produce tie-line
patterns consistent with experiment. Hence, our combined
theoretical-experimental analysis indicates that weaker interactions or
higher-order complexes beyond the 3:2 stoichiometry, but not yet documented,
are involved in the formation of SynGAP/PSD-95 condensates, imploring future
efforts to ascertain the nature of these auxiliary interactions in PSD-like
LLPS.


  Compartmental epidemic models have been widely used for predicting the course
of epidemics, from estimating the basic reproduction number to guiding
intervention policies. Studies commonly acknowledge these models' assumptions
but less often justify their validity in the specific context in which they are
being used. Our purpose is not to argue for specific alternatives or
modifications to compartmental models, but rather to show how assumptions can
constrain model outcomes to a narrow portion of the wide landscape of potential
epidemic behaviors. This concrete examination of well-known models also serves
to illustrate general principles of modeling that can be applied in other
contexts.


  We analyze the susceptible-infected-susceptible model for epidemic spreading
in which a fraction of the individuals become immune by vaccination. This
process is understood as a dilution by vaccination, which decreases the
fraction of the susceptible individuals. For a nonzero fraction of vaccinated
individuals, the model predicts a new state in which the disease spreads but
eventually becomes extinct. The new state emerges when the fraction of
vaccinated individuals is greater than a critical value. The model predicts
that this critical value increases as one increases the infection rate reaching
an asymptotic value, which is strictly less than the unity. Above this
asymptotic value, the extinction occurs no matter how large the infection rate
is.


  We investigate possible shapes of the electric field, which oscillating
dipoles in a certain region of biological tissue can produce in a neighboring
region, or outside the tissue boundaries. We find that a wide range of shapes,
including the typical morphology of limbs and appendages, can be generated as a
zone of extremely low field amplitudes embedded in a zone of much larger field
amplitudes. Neutral molecules with a resonance close to the frequency of the
oscillating field may be attracted to this zone or be repelled from it, while
the driving effect on molecules with an electric charge is only extremely weak.
The forces would be sufficient for the controlled deposition of molecules
during growth or regeneration. They could also serve as a method of information
transfer.


  Black-box optimization formulations for biological sequence design have drawn
recent attention due to their promising potential impact on the pharmaceutical
industry. In this work, we propose to unify two seemingly distinct worlds:
likelihood-free inference and black-box optimization, under one probabilistic
framework. In tandem, we provide a recipe for constructing various sequence
design methods based on this framework. We show how previous optimization
approaches can be "reinvented" in our framework, and further propose new
probabilistic black-box optimization algorithms. Extensive experiments on
sequence design application illustrate the benefits of the proposed
methodology.


  COPASI is a popular application for simulation and analysis of biochemical
networks and their dynamics. While this software is widely used, it works as a
standalone application and until now it was not possible for users to interact
with its models through the web. We built ShinyCOPASI, a web-based application
that allows COPASI models to be explored through a web browser. ShinyCOPASI was
written in R with the CoRC package, which provides a high-level R API for
COPASI, and the Shiny package to expose it as a web application. The web view
provided by ShinyCOPASI follows a similar interface to the standalone COPASI
and allows users to explore the details of a model, as well as running a subset
of the tasks available in COPASI from within a browser. A generic version
allows users to load model files from their computer, while another one
pre-loads a specific model from the server and may be useful to provide web
access to published models. The application is available at:
http://shiny.copasi.org/; and the source code is at:
https://github.com/copasi/shinyCOPASI.


  Tree-child networks are an important network class which are used in
phylogenetics to model reticulate evolution. In a recent paper, Pons and Batle
(2021) conjectured a relation between tree-child networks and certain words. In
this short note, we prove their conjecture for the (important) class of
one-component tree-child networks.


  Phylogenetic trees are widely used to understand the evolutionary history of
organisms. Tree shapes provide information about macroevolutionary processes.
However, macroevolutionary models are unreliable for inferring the true
processes underlying empirical trees. Here, we propose a flexible and
biologically plausible macroevolutionary model for phylogenetic trees where
times to speciation or extinction events are drawn from a Coxian phase-type
(PH) distribution. First, we show that different choices of parameters in our
model lead to a range of tree balances as measured by Aldous' $\beta$
statistic. In particular, we demonstrate that it is possible to find parameters
that correspond well to empirical tree balance. Next, we provide a natural
extension of the $\beta$ statistic to sets of trees. This extension produces
less biased estimates of $\beta$ compared to using the median $\beta$ values
from individual trees. Furthermore, we derive a likelihood expression for the
probability of observing any tree with branch lengths under a model with
speciation but no extinction. Finally, we illustrate the application of our
model by performing both absolute and relative goodness-of-fit tests for two
large empirical phylogenies (squamates and angiosperms) that compare models
with Coxian PH distributed times to speciation with models that assume
exponential or Weibull distributed waiting times. In our numerical analysis, we
found that, in most cases, models assuming a Coxian PH distribution provided
the best fit.


  Molecular property prediction is one of the fastest-growing applications of
deep learning with critical real-world impacts. Including 3D molecular
structure as input to learned models improves their performance for many
molecular tasks. However, this information is infeasible to compute at the
scale required by several real-world applications. We propose pre-training a
model to reason about the geometry of molecules given only their 2D molecular
graphs. Using methods from self-supervised learning, we maximize the mutual
information between 3D summary vectors and the representations of a Graph
Neural Network (GNN) such that they contain latent 3D information. During
fine-tuning on molecules with unknown geometry, the GNN still generates
implicit 3D information and can use it to improve downstream tasks. We show
that 3D pre-training provides significant improvements for a wide range of
properties, such as a 22% average MAE reduction on eight quantum mechanical
properties. Moreover, the learned representations can be effectively
transferred between datasets in different molecular spaces.


  Several important aspects related to SARS-CoV-2 transmission are not well
known due to a lack of appropriate data. However, mathematical and
computational tools can be used to extract part of this information from the
available data, like some hidden age-related characteristics. In this paper, we
present a method to investigate age-specific differences in transmission
parameters related to susceptibility to and infectiousness upon contracting
SARS-CoV-2 infection. More specifically, we use panel-based social contact data
from diary-based surveys conducted in Belgium combined with the next generation
principle to infer the relative incidence and we compare this to real-life
incidence data. Comparing these two allows for the estimation of age-specific
transmission parameters. Our analysis implies the susceptibility in children to
be around half of the susceptibility in adults, and even lower for very young
children (preschooler). However, the probability of adults and the elderly to
contract the infection is decreasing throughout the vaccination campaign,
thereby modifying the picture over time.


  The dose-response model has been widely used for quantifying the risk of
infection of airborne diseases like COVID-19. The model has been used in the
room-average analysis of infection risk and analysis using passive scalars as a
proxy for aerosol transport. However, it has not been employed for risk
estimation in numerical simulations of droplet dispersion. In this work, we
develop a framework for the evaluation of the probability of infection in
droplet dispersion simulations using the dose-response model. We introduce a
version of the model that can incorporate the higher transmissibility of
variant strains of SARS-CoV2 and the effect of vaccination in evaluating the
probability of infection. Numerical simulations of droplet dispersion during
speech are carried out to investigate the infection risk over space and time
using the model. The advantage of droplet dispersion simulations for risk
evaluation is demonstrated through the analysis of the effect of humidity on
infection risk.


  In this article, we investigate various numerical methods for computing
scaled or logarithmic sensitivities of the form $\partial \ln y/\partial \ln
x$. The methods tested include One Point, Two Point, Five Point, and the
Richardson Extrapolation. The different methods were applied to a variety of
mathematical functions as well as a reaction network model. The algorithms were
validated by comparing results with known analytical solutions for functions
and using the Reder method for computing the sensitivities in reaction networks
via the Tellurium package. For evaluation, two aspects were looked at, accuracy
and time taken to compute the sensitivities. Of the four methods, Richardson's
extrapolation was by far the most accurate but also the slowest in terms of
performance. For fast, reasonably accurate estimates, we recommend the
two-point method. For most other cases where the derivatives are changing
rapidly, the five-point method is a good choice, although it is three times
slower than the two-point method. For ultimate accuracy which would apply
particularly to very fast changing derivatives the Richardson method is without
doubt the best, but it is seven-times slower than the two point method. We do
not recommend the one-point method in any circumstance. The Python software
that was used in the study with documentation is available at:
\url{https://github.com/evanyfyip/SensitivityAnalysis}.


  Computational modeling is crucial for understanding and analyzing complex
systems. In biology, model creation is a human dependent task that requires
reading hundreds of papers and conducting wet lab experiments, which would take
days or months. To overcome this hurdle, we propose a novel automated method,
that utilizes the knowledge published in literature to suggest model extensions
by selecting most relevant and useful information in few seconds. In
particular, our novel approach organizes the events extracted from the
literature as a collaboration graph with additional metric that relies on the
event occurrence frequency in literature. Additionally, we show that common
graph centrality metrics vary in the assessment of the extracted events. We
have demonstrated the reliability of the proposed method using three different
selected models, namely, T cell differentiation, T cell large granular
lymphocyte, and pancreatic cancer cell. Our proposed method was able to find
high percent of the desired new events with an average recall of 82%.


  Antibodies are versatile proteins that bind to pathogens like viruses and
stimulate the adaptive immune system. The specificity of antibody binding is
determined by complementarity-determining regions (CDRs) at the tips of these
Y-shaped proteins. In this paper, we propose a generative model to
automatically design the CDRs of antibodies with enhanced binding specificity
or neutralization capabilities. Previous generative approaches formulate
protein design as a structure-conditioned sequence generation task, assuming
the desired 3D structure is given a priori. In contrast, we propose to
co-design the sequence and 3D structure of CDRs as graphs. Our model unravels a
sequence autoregressively while iteratively refining its predicted global
structure. The inferred structure in turn guides subsequent residue choices.
For efficiency, we model the conditional dependence between residues inside and
outside of a CDR in a coarse-grained manner. Our method achieves superior
log-likelihood on the test set and outperforms previous baselines in designing
antibodies capable of neutralizing the SARS-CoV-2 virus.


  The importance of understanding SARS-CoV-2 evolution cannot be
overemphasized. Recent studies confirm that natural selection is the dominating
mechanism of SARS-CoV-2 evolution, which favors mutations that strengthen viral
infectivity. We demonstrate that vaccine-breakthrough or antibody-resistant
mutations provide a new mechanism of viral evolution. Specifically,
vaccine-resistant mutation Y449S in the spike (S) protein receptor-bonding
domain (RBD), which occurred in co-mutation [Y449S, N501Y], has reduced
infectivity compared to the original SARS-CoV-2 but can disrupt existing
antibodies that neutralize the virus. By tracing the evolutionary trajectories
of vaccine-resistant mutations in over 1.9 million SARS-CoV-2 genomes, we
reveal that the occurrence and frequency of vaccine-resistant mutations
correlate strongly with the vaccination rates in Europe and America. We
anticipate that as a complementary transmission pathway, vaccine-resistant
mutations will become a dominating mechanism of SARS-CoV-2 evolution when most
of the world's population is vaccinated. Our study sheds light on SARS-CoV-2
evolution and transmission and enables the design of the next-generation
mutation-proof vaccines and antibody drugs.


  Effective contact tracing is crucial to contain epidemic spreading without
disrupting societal activities especially in the present time of coexistence
with a pandemic outbreak. Large gatherings play a key role, potentially
favouring superspreading events. However, the effects of tracing in large
groups have not been fully assessed so far. We show that beside forward
tracing, which reconstructs to whom disease spreads, and backward tracing,
which searches from whom disease spreads, a third "sideward" tracing is always
present, when tracing gatherings. This is an indirect tracing that detects
infected asymptomatic individuals, even if they have neither been directly
infected by, nor they have directly transmitted the infection to the index
case. We analyse this effect in a model of epidemic spreading for SARS-CoV-2,
within the framework of simplicial activity-driven temporal networks. We
determine the contribution of the three tracing mechanisms to the suppression
of epidemic spreading, showing that sideward tracing induces a non-monotonic
behaviour in the tracing efficiency, as a function of the size of the
gatherings. Based on our results, we suggest an optimal choice for the sizes of
the gatherings to be traced and we test the strategy on an empirical dataset of
gatherings in a University Campus.


  Gene expression is a fundamental process in a living system. The small RNAs
(sRNAs) is widely observed as a global regulator in gene expression. The
inherent nonlinearity in this regulatory process together with the bursty
production of messenger RNA (mRNA), sRNA and protein make the exact solution
for this stochastic process intractable. This is particularly the case when
quantifying the protein noise level, which has great impact on multiple
cellular processes. Here we propose an approximate yet reasonably accurate
solution for the gene expression noise with infrequent burst and strong
regulation by sRNAs. This analytical solution allows us to better analyze the
noise and stochastic deviation of protein level. We find that the regulation
amplifies the noise, reduces the protein level. The stochasticity in the
regulation generates more proteins than what if the stochasticity is removed
from the system. The sRNA level is most important to the relationship between
the noise and stochastic deviation. The results provide analytical tools for
more general studies of gene expression and strengthen our quantitative
understandings of post-transcriptional regulation in controlling gene
expression processes.


  Recombination is a powerful evolutionary process that shapes the genetic
diversity observed in the populations of many species. Reconstructing
genealogies in the presence of recombination from sequencing data is a very
challenging problem, as this relies on mutations having occurred on the correct
lineages in order to detect the recombination and resolve the placement of
edges in the local trees. We investigate the probability of recovering the true
topology of ancestral recombination graphs (ARGs)under the coalescent with
recombination and gene conversion. We explore how sample size and mutation rate
affect the inherent uncertainty in reconstructed ARGs; this sheds light on the
theoretical limitations of ARG reconstruction methods. We illustrate our
results using estimates of evolutionary rates for several biological organisms;
in particular, we find that for parameter values that are realistic for
SARS-CoV-2, the probability of reconstructing genealogies that are close to the
truth is low.


  Mathematical modeling is one of the key factors of the effective control of
newly found infectious diseases, such as COVID-19. Our knowledge about the
parameters and the course of the infection is highly limited in the beginning
of the epidemic, hence computer implementation of the models have to be quick
and flexible. The propensity matrix - update graph method we discuss in this
paper serves as a convenient approach to efficiently implement age structured
stochastic epidemic models. The code base we implemented for our forecasting
work is also included in the attached GitHub repository.


  Recent work suggests that cross-feeding -- the secretion and consumption of
metabolic biproducts by microbes -- is essential for understanding microbial
ecology. Yet how cross-feeding and competition combine to give rise to
ecosystem-level properties remains poorly understood. To address this question,
we analytically analyze the Microbial Consumer Resource Model (MiCRM), a
prominent ecological model commonly used to study microbial communities. Our
mean-field solution exploits the fact that unlike replicas, the cavity method
does not require the existence of a Lyapunov function. We use our solution to
derive new species-packing bounds for diverse ecosystems in the presence of
cross-feeding, as well as simple expressions for species richness and the
abundance of secreted resources as a function of cross-feeding (metabolic
leakage) and competition. Our results show how a complex interplay between
competition for resources and cooperation resulting from metabolic exchange
combine to shape the properties of microbial ecosystems.


  Feedback is ubiquitous in both biological and engineered control systems. In
biology, in addition to typical feedback between plant and controller, we
observe feedback pathways within control systems, which we call internal
feedback pathways (IFPs), that are often very complex. IFPs are most familiar
in neural systems, our primary motivation, but they appear everywhere from
bacterial signal transduction to the human immune system. In this paper, we
describe these very different motivating examples and introduce the concepts
necessary to explain their complex IFPs, particularly the severe speed-accuracy
tradeoffs that constrain the hardware in biology. We also sketch some minimal
theory for extremely simplified toy models that nevertheless highlight the
importance of diversity-enabled sweet spots (DESS) in mitigating the impact of
hardware tradeoffs. For more realistic models, standard modern and robust
control theory can give some insights into previously cryptic IFPs, and the new
System Level Synthesis theory expands this substantially. These additional
theories explaining IFPs will be explored in more detail in several companion
papers.


  Birth-and-death processes (BDPs) form a class of continuous-time Markov
chains that are particularly suited to describing the changes in the size of a
population over time. Population-size-dependent BDPs (PSDBDPs) allow the rate
at which a population grows to depend on the current population size. The main
purpose of our new Python package BirDePy is to provide easy-to-use functions
that allow the parameters of discretely-observed PSDBDPs to be estimated. The
package can also be used to estimate parameters of continuously-observed
PSDBDPs, simulate sample paths, approximate transition probabilities, and
generate forecasts. We describe in detail several methods which have been
incorporated into BirDePy to achieve each of these tasks. The usage and
effectiveness of the package is demonstrated through a variety of examples of
PSDBDPs, as well as case studies involving annual population count data of two
endangered bird species.


  In the genome biology research, regulatory genome modeling is an important
topic for many regulatory downstream tasks, such as promoter classification,
transaction factor binding sites prediction. The core problem is to model how
regulatory elements interact with each other and its variability across
different cell types. However, current deep learning methods often focus on
modeling genome sequences of a fixed set of cell types and do not account for
the interaction between multiple regulatory elements, making them only perform
well on the cell types in the training set and lack the generalizability
required in biological applications. In this work, we propose a simple yet
effective approach for pre-training genome data in a multi-modal and
self-supervised manner, which we call GeneBERT. Specifically, we simultaneously
take the 1d sequence of genome data and a 2d matrix of (transcription factors x
regions) as the input, where three pre-training tasks are proposed to improve
the robustness and generalizability of our model. We pre-train our model on the
ATAC-seq dataset with 17 million genome sequences. We evaluate our GeneBERT on
regulatory downstream tasks across different cell types, including promoter
classification, transaction factor binding sites prediction, disease risk
estimation, and splicing sites prediction. Extensive experiments demonstrate
the effectiveness of multi-modal and self-supervised pre-training for
large-scale regulatory genomics data.


  The Eocene--Oligocene transition (EOT) represents a period of global
environmental changes particularly marked in Europe and coincides with a
dramatic biotic turnover. Here, using an exceptional fossil preservation, we
document and analyse the diversity dynamics of a mammal clade, Cainotherioidea
(Artiodactyla), that survived the EOT and radiated rapidly immediately after.
We infer their diversification history from Quercy Konzentrat--Lagerst{\"a}tte
(south-west France) at the species level using Bayesian birth--death models. We
show that cainotherioid diversity fluctuated through time, with extinction
events at the EOT and in the late Oligocene, and a major speciation burst in
the early Oligocene. The latter is in line with our finding that cainotherioids
had a high morphological adaptability following environmental changes
throughout the EOT, which probably played a key role in the survival and
evolutionary success of this clade in the aftermath. Speciation is positively
associated with temperature and continental fragmentation in a time-continuous
way, while extinction seems to synchronize with environmental change in a
punctuated way. Within-clade interactions negatively affected the cainotherioid
diversification, while inter-clade competition might explain their final
decline during the late Oligocene. Our results provide a detailed dynamic
picture of the evolutionary history of a mammal clade in a context of global
change.


  Interacting particle systems undergoing repeated mutation and selection steps
model genetic evolution, and also describe a broad class of sequential Monte
Carlo methods. The genealogical tree embedded into the system is important in
both applications. Under neutrality, when fitnesses of particles are
independent from those of their parents, rescaled genealogies are known to
converge to Kingman's coalescent. Recent work has established convergence under
non-neutrality, but only for finite-dimensional distributions. We prove weak
convergence of non-neutral genealogies on the space of c\`adl\`ag paths under
standard assumptions, enabling analysis of the whole genealogical tree.


  We report a fast-track computationally-driven discovery of new SARS-CoV2 Main
Protease (M$^{pro}$) inhibitors whose potency range from mM for initial
non-covalent ligands to sub-$\mu$M for the final covalent compound (IC50=830
+/- 50 nM). The project extensively relied on high-resolution all-atom
molecular dynamics simulations and absolute binding free energy calculations
performed using the polarizable AMOEBA force field. The study is complemented
by extensive adaptive sampling simulations that are used to rationalize the
different ligands binding poses through the explicit reconstruction of the
ligand-protein conformation spaces. Machine Learning predictions are also
performed to predict selected compound properties. While simulations
extensively use High Performance Computing to strongly reduce time-to-solution,
they were systematically coupled to Nuclear Magnetic Resonance experiments to
drive synthesis and to in vitro characterization of compounds. Such study
highlights the power of in silico strategies that rely on structure-based
approaches for drug design and allows to address the protein conformational
multiplicity problem. The proposed fluorinated tetrahydroquinolines open routes
for further optimization of M$^{pro}$ inhibitors towards low nM affinities.


  In this work, we develop a kinetic model of tumour growth taking into account
the effects of clinical uncertainties characterising the tumours' progression.
The action of therapeutic protocols trying to steer the tumours' volume towards
a target size is then investigated by means of suitable selective-type controls
acting at the level of cellular dynamics. By means of classical tools of
statistical mechanics for many-agent systems, we are able to prove that it is
possible to dampen clinical uncertainties across the scales. To take into
account the scarcity of clinical data and the possible source of error in the
image segmentation of tumours' evolution, we estimated empirical distributions
of relevant parameters that are considered to calibrate the resulting model
obtained from real cases of primary glioblastoma. Suitable numerical methods
for uncertainty quantification of the resulting kinetic equations are discussed
and, in the last part of the paper, we compare the effectiveness of the
introduced control approaches in reducing the variability in tumours' size due
to the presence of uncertain quantities.


  There is a scientific and ethical imperative to embrace a multidimensional,
continuous view of ancestry and move away from continental ancestry categories


  Biomolecular motor proteins that generate forces by consuming chemical energy
obtained from ATP hydrolysis are pivotal for organizing broad cytoskeletal
structures in living cells. The control of such cytoskeletal structures
benefits programmable protein patterning; however, our current knowledge is
limited owing to the underdevelopment of an engineering approach for
controlling pattern formation. Here, we demonstrate the tailoring of assembled
patterns of microtubules (MTs) driven by kinesin motors by designing the
boundary shape in fabricated microwells. We found an MT bundle structure along
the microwell wall and a bridging structure perpendicular to the wall.
Corroborated by the theory of self-propelled rods, we further showed that the
alignment of MTs defined by the boundary shape determined the transition of the
assembled patterns, providing a blueprint to reconstruct bridge structures in
microchannels. Our findings provide a geometric rule to tailor the
self-organization of cytoskeletons and motor proteins for nanotechnological
applications.


  Gaussian processes are widely used for the analysis of spatial data due to
their nonparametric flexibility and ability to quantify uncertainty, and
recently developed scalable approximations have facilitated application to
massive datasets. For multivariate outcomes, linear models of coregionalization
combine dimension reduction with spatial correlation. However, their
real-valued latent factors and loadings are difficult to interpret because,
unlike nonnegative models, they do not recover a parts-based representation. We
present nonnegative spatial factorization (NSF), a spatially-aware
probabilistic dimension reduction model that naturally encourages sparsity. We
compare NSF to real-valued spatial factorizations such as MEFISTO and
nonspatial dimension reduction methods using simulations and high-dimensional
spatial transcriptomics data. NSF identifies generalizable spatial patterns of
gene expression. Since not all patterns of gene expression are spatial, we also
propose a hybrid extension of NSF that combines spatial and nonspatial
components, enabling quantification of spatial importance for both observations
and features. A TensorFlow implementation of NSF is available from
https://github.com/willtownes/nsf-paper .


  Extremophiles have gained prominence by providing an experimental approach to
astrobiology. Extremophiles gain equal value by being part of a framework for
high-level characterisation of the evolutionary mechanisms that must
necessarily restrict or promote their emergence and presence on solar system
bodies. Thus, extremophiles exist in extreme environments, and therein lies the
paradox: extremophiles can only live in extreme environments but yet are not
able to originate in such environments. Therefore, even though the range of
extremophile capabilities in extreme environments is wider than that in
mesophiles, the range of their emergence possibilities is still equally
restricted. Therefore, even if one locates an extreme exoworld where
terrestrial extremophiles could live here-and-now, it can be predicted that no
extremophile analogues are present anyway. Furthermore, it is possible for a
world to be uninhabited, yet be habitable, and therein arises the extreme
environment paradox: an extreme environment can sustain chemical evolution as
well as arriving non-native life, yet native life cannot be built up in that
very environment. Thus, life may exist on an extraterrestrial extreme world (if
imported there), and chemical evolution may be present on that world. However,
it can be predicted that there is no native life anyway. This situation can be
predicted to function as a chemosignature and eventually as a biosignature.
However, the fact that a non-native extremopile in principle can exist in
extreme environments may demonstrate that the intermediate step between
chemical evolution and extremophiles can still occur in the form of a
statistical deviation.


  The first year of the COVID-19 pandemic put considerable strain on the
national healthcare system in England. In order to predict the effect of the
local epidemic on hospital capacity in England, we used a variety of data
streams to inform the construction and parameterisation of a hospital
progression model, which was coupled to a model of the generalised epidemic. We
named this model EpiBeds. Data from a partially complete patient-pathway
line-list was used to provide initial estimates of the mean duration that
individuals spend in the different hospital compartments. We then fitted
EpiBeds using complete data on hospital occupancy and hospital deaths, enabling
estimation of the proportion of individuals that follow different clinical
pathways, and the reproduction number of the generalised epidemic. The
construction of EpiBeds makes it straightforward to adapt to different patient
pathways and settings beyond England. As part of the UK response to the
pandemic, EpiBeds has provided weekly forecasts to the NHS for hospital bed
occupancy and admissions in England, Wales, Scotland, and Northern Ireland.


  Models of transcriptional regulation that assume equilibrium binding of
transcription factors have been very successful at predicting gene expression
from sequence in bacteria. However, analogous equilibrium models do not perform
as well in eukaryotes, most likely due to the largely out-of-equilibrium nature
of eukaryotic regulatory processes. These processes come with unavoidable
energy expenditure at the molecular level to support precise, reliable, or fast
gene expression responses that could correspond to evolutionarily optimized
regulatory strategies. Unfortunately, the space of possible non-equilibrium
mechanisms is vast and predominantly uninteresting. The key question is
therefore how this space can be navigated efficiently, to focus on mechanisms
and models that are biologically relevant. In this review, we advocate for the
normative role of theory - theory that prescribes rather than just describes -
in providing such a focus. Theory should expand its remit beyond inferring
models from data (by fitting), towards identifying non-equilibrium gene
regulatory schemes (by optimizing) that may have been evolutionarily selected
due to their favorable functional characteristics which outperform regulation
at equilibrium. This approach can help us navigate the expanding complexity of
regulatory architectures, and refocus the questions about the observed
biological mechanisms from how they work towards why they have evolved in the
first place. We illustrate our reasoning by toy examples for which we provide
simulation code.


  In Fall 2020, several European countries reported rapid increases in COVID-19
cases along with growing estimates of the effective reproduction rates. Such an
acceleration in epidemic spread is usually attributed to time-dependent
effects, e.g. human travel, seasonal behavioral changes, mutations of the
pathogen etc. In this case however the acceleration occurred when counter
measures such as testing and contact tracing exceeded their capacity limit.
Considering Austria as an example, here we show that this dynamics can be
captured by a time-independent, i.e. autonomous, compartmental model that
incorporates these capacity limits. In this model, the epidemic acceleration
coincides with the exhaustion of mitigation efforts, resulting in an increasing
fraction of undetected cases that drive the effective reproduction rate
progressively higher. We demonstrate that standard models which does not
include this effect necessarily result in a systematic underestimation of the
effective reproduction rate.


  Algorithmic strategies for the spatio-temporal simulation of multi-cellular
systems are crucial to generate synthetic datasets for bioinformatics tools
benchmarking, as well as to investigate experimental hypotheses on real-world
systems in a variety of in-silico scenarios. In particular, efficient
algorithms are needed to overcome the harsh trade-off between scalability and
expressivity, which typically limits our capability to produce realistic
simulations, especially in the context of cancer evolution. We introduce the
Optimized Gillespie algorithm for simulating Stochastic sPAtial models of
Cancer Evolution (OG-SPACE), a computational framework for the simulation of
the spatio-temporal evolution of cancer subpopulations and of the experimental
procedures of both bulk andsingle-cell sequencing. OG-SPACE relies on an
evolution of the Gillespie algorithm optimized to deal with large numbers of
cells and is designed tohandle a variety of birth-death processes and
interaction rules on arbitrary lattices. As output OG-SPACE returns: the visual
snapshots of the spatial configuration of the system over time, the phylogeny
of the (sampled) cells, the mutational tree, the variant allele frequency
spectrum (for bulk experiments) and the cell genotypes (for single-cell
experiments).OG-SPACE is freely available
at:https://github.com/BIMIB-DISCo/OG-SPACE


  In a viral pandemic, a few important tests are required for successful
containment of the virus and reduction in severity of the infection. Among
those tests, a test for the neutralizing ability of an antibody is crucial for
assessment of population immunity gained through vaccination, and to test
therapeutic value of antibodies made to counter the infections. Here, we report
a sensitive technique to detect the relative neutralizing strength of various
antibodies against the SARS-CoV-2 virus. We used bright, photostable,
background-free, fluorescent upconversion nanoparticles conjugated with
SARS-CoV-2 receptor binding domain as a phantom virion. A glass bottom plate
coated with angiotensin-converting enzyme 2 (ACE-2) protein imitates the target
cells. When no neutralizing IgG antibody was present in the sample, the
particles would bind to the ACE-2 with high affinity. In contrast, a
neutralizing antibody can prevent particle attachment to the ACE-2-coated
substrate. A prototype system consisting of a custom-made confocal microscope
was used to quantify particle attachment to the substrate. The sensitivity of
this assay can reach 4.0 ng/ml and the dynamic range is from 1.0 ng/ml to 3.2
{\mu}g/ml. This is to be compared to 19 ng/ml sensitivity of commercially
available kits.


  Edge Constrained Vertex Coloring (ECVC) problems are defined on a finite
multigraph, their solutions are characterized, and a linear time algorithm is
given for solving $N$ ECVCs on the same underlying multigraph. Using ECVC
problems we develop a novel family-based haplotype reconstruction method which
has linear-time complexity in both number of markers and family size and has
many other desirable properties. To do so, we define a multigraph given a
genomic interval on which a family is recombination-free.


  The chromosomal DNA of bacteria is folded into a compact body called the
nucleoid, which is composed essentially of DNA (80%), RNA (10%), and a number
of different proteins (10%). These nucleoid proteins act as regulators of gene
expression and influence the organization of the nucleoid by bridging, bending,
or wrapping the DNA. These so-called architectural properties of nucleoid
proteins are still poorly understood. For example, the reason why certain
proteins compact the DNA coil in certain environments but make instead the DNA
more rigid in other environments is the matter of ongoing debates. Here, we
address the question of the impact of the self-association of nucleoid proteins
on their architectural properties and try to determine whether differences in
self-association are sufficient to induce large changes in the organization of
the DNA coil. More specifically, we developed two coarse-grained models of
proteins, which interact identically with the DNA but self-associate
differently by forming either clusters or filaments in the absence of the DNA.
We showed through Brownian dynamics simulations that self-association of the
proteins increases dramatically their ability to shape the DNA coil. Moreover,
we observed that cluster-forming proteins compact significantly the DNA coil
(similar to the DNA-bridging mode of H-NS proteins), whereas filament-forming
proteins increase instead significantly the stiffness of the DNA chain (similar
to the DNA-stiffening mode of H-NS proteins). This work consequently suggests
that the knowledge of the DNA-binding properties of the proteins is in itself
not sufficient to understand their architectural properties. Rather, their
self-association properties must also be investigated in detail, because they
might actually drive the formation of different DNA/protein complexes.


  Consider a Markovian SIR epidemic model in a homogeneous community. To this
model we add a rate at which individuals are tested, and once an infectious
individual tests positive it is isolated and each of their contacts are traced
and tested independently with some fixed probability. If such a traced
individual tests positive it is isolated, and the contact tracing is iterated.
This model is analysed using large population approximations, both for the
early stage of the epidemic when the "to-be-traced components" of the epidemic
behaves like a branching process, and for the main stage of the epidemic where
the process of to-be-traced components converges to a deterministic process
defined by a system of differential equations. These approximations are used to
quantify the effect of testing and of contact tracing on the effective
reproduction numbers (for the components as well as for the individuals), the
probability of a major outbreak, and the final fraction getting infected. Using
numerical illustrations when rates of infection and natural recovery are fixed,
it is shown that Test-and-Trace strategy is effective in reducing the
reproduction number. Surprisingly, the reproduction number for the branching
process of components is not monotonically decreasing in the tracing
probability, but the individual reproduction number is conjectured to be
monotonic as expected. Further, in the situation where individuals also
self-report for testing, the tracing probability is more influential than the
screening rate (measured by the fraction infected being screened).


  We introduce a system of differential equations to assess the impact of
(self-)quarantine of symptomatic infectious individuals on disease dynamics. To
this end we depart from using the classic bilinear infection process, but
remain still within the framework of the mass-action assumption. From the
mathematical point of view our model is interesting due to the lack of
continuous differentiability at disease free steady states, which implies also
that the basic reproductive number cannot be computed following established
approaches for certain parameter values. However, we parametrise our
mathematical model using published values from the COVID-19 literature, and
analyse the model simulations. We also contrast model simulations against
publicly available COVID-19 test data focusing on the first wave of the
pandemic during March - July 2020 in the UK. Our simulations indicate that
actual peak case numbers might have been as much as 200 times higher than the
reported positive test cases during the first wave in the UK. We find that very
strong adherence to self-quarantine rules yields (only) a reduction of 22$\%$
of peak numbers and delays the onset of the peak by approximately 30-35 days.
However, during the early phase of the outbreak the impact of (self)-quarantine
is much more significant. We also take into account the effect of a national
lockdown in a simplistic way by reducing the effective susceptible population
size. We find that in case of a 90$\%$ reduction of the effective susceptible
population size, strong adherence to self-quarantine still only yields a 25$\%$
reduction of peak infectious numbers when compared to low adherence. This is
due to the significant number of asymptomatic infectious individuals in the
population.


  Messenger RNA-based medicines hold immense potential, as evidenced by their
rapid deployment as COVID-19 vaccines. However, worldwide distribution of mRNA
molecules has been limited by their thermostability, which is fundamentally
limited by the intrinsic instability of RNA molecules to a chemical degradation
reaction called in-line hydrolysis. Predicting the degradation of an RNA
molecule is a key task in designing more stable RNA-based therapeutics. Here,
we describe a crowdsourced machine learning competition ("Stanford
OpenVaccine") on Kaggle, involving single-nucleotide resolution measurements on
6043 102-130-nucleotide diverse RNA constructs that were themselves solicited
through crowdsourcing on the RNA design platform Eterna. The entire experiment
was completed in less than 6 months, and 41% of nucleotide-level predictions
from the winning model were within experimental error of the ground truth
measurement. Furthermore, these models generalized to blindly predicting
orthogonal degradation data on much longer mRNA molecules (504-1588
nucleotides) with improved accuracy compared to previously published models.
Top teams integrated natural language processing architectures and data
augmentation techniques with predictions from previous dynamic programming
models for RNA secondary structure. These results indicate that such models are
capable of representing in-line hydrolysis with excellent accuracy, supporting
their use for designing stabilized messenger RNAs. The integration of two
crowdsourcing platforms, one for data set creation and another for machine
learning, may be fruitful for other urgent problems that demand scientific
discovery on rapid timescales.


  Human and animal tissues consist of heterogeneous cell types that organize
and interact in highly structured manners. Bulk and single-cell sequencing
technologies remove cells from their original microenvironments, resulting in a
loss of spatial information. Spatial transcriptomics is a recent technological
innovation that measures transcriptomic information while preserving spatial
information. Spatial transcriptomic data can be generated in several ways. RNA
molecules are measured by in situ sequencing, in situ hybridization, or spatial
barcoding to recover original spatial coordinates. The inclusion of spatial
information expands the range of possibilities for analysis and visualization,
and spurred the development of numerous novel methods. In this review, we
summarize the core concepts of spatial genomics technology and provide a
comprehensive review of current analysis and visualization methods for spatial
transcriptomics.


  We introduce cellanneal, a python-based software for deconvolving bulk RNA
sequencing data. cellanneal relies on the optimization of Spearman's rank
correlation coefficient between experimental and computational mixture gene
expression vectors using simulated annealing. cellanneal can be used as a
python package or via a command line interface, but importantly also provides a
simple graphical user interface which is distributed as a single executable
file for user convenience. The python package is available at
https://github.com/LiBuchauer/cellanneal , the graphical software can be
downloaded at http://shalevlab.weizmann.ac.il/resources .


  Deep learning has become popular because of its potential to achieve high
accuracy in prediction tasks. However, accuracy is not always the only goal of
statistical modelling, especially for models developed as part of scientific
research. Rather, many scientific models are developed to facilitate scientific
discovery, by which we mean to abstract a human-understandable representation
of the natural world. Unfortunately, the opacity of deep neural networks limit
their role in scientific discovery, creating a new demand for models that are
transparently interpretable. This article is a field guide to transparent model
design. It provides a taxonomy of transparent model design concepts, a
practical workflow for putting design concepts into practice, and a general
template for reporting design choices. We hope this field guide will help
researchers more effectively design transparently interpretable models, and
thus enable them to use deep learning for scientific discovery.


  The concept of individual admixture (IA) assumes that the genome of
individuals is composed of alleles inherited from $K$ ancestral populations.
Each copy of each allele has the same chance $q_k$ to originate from population
$k$, and together with the allele frequencies $p$ in all populations at all $M$
markers, comprises the admixture model. Here, we assume a supervised scheme,
i.e.\ allele frequencies $p$ are given through a reference database of size
$N$, and $q$ is estimated via maximum likelihood for a single sample. We study
laws of large numbers and central limit theorems describing effects of
finiteness of both, $M$ and $N$, on the estimate of $q$. We recall results for
the effect of finite $M$, and provide a central limit theorem for the effect of
finite $N$, introduce a new way to express the uncertainty in estimates in
standard barplots, give simulation results, and discuss applications in
forensic genetics.


  Current metagenome assemblers developed for short sequence reads or noisy
long readswere not optimized for accurate long reads. Here we describe
hifiasm-meta, a new metagenome assembler that exploits the high accuracy of
recent data. Evaluated on seven empirical datasets, hifiasm-meta reconstructed
tens to hundreds of complete circular bacterial genomes per dataset,
consistently outperforming other metagenome assemblers.


  We consider the problem of projecting a vector onto the so-called k-capped
simplex, which is a hyper-cube cut by a hyperplane. For an n-dimensional input
vector with bounded elements, we found that a simple algorithm based on
Newton's method is able to solve the projection problem to high precision with
a complexity roughly about O(n), which has a much lower computational cost
compared with the existing sorting-based methods proposed in the literature. We
provide a theory for partial explanation and justification of the method.
  We demonstrate that the proposed algorithm can produce a solution of the
projection problem with high precision on large scale datasets, and the
algorithm is able to significantly outperform the state-of-the-art methods in
terms of runtime (about 6-8 times faster than a commercial software with
respect to CPU time for input vector with 1 million variables or more).
  We further illustrate the effectiveness of the proposed algorithm on solving
sparse regression in a bioinformatics problem. Empirical results on the GWAS
dataset (with 1,500,000 single-nucleotide polymorphisms) show that, when using
the proposed method to accelerate the Projected Quasi-Newton (PQN) method, the
accelerated PQN algorithm is able to handle huge-scale regression problem and
it is more efficient (about 3-6 times faster) than the current state-of-the-art
methods.


  As a classical problem in the field of complex networks, link prediction has
attracted much attention from researchers, which is of great significance to
help us understand the evolution and dynamic development mechanisms of
networks. Although various network type-specific algorithms have been proposed
to tackle the link prediction problem, most of them suppose that the network
structure is dominated by the Triadic Closure Principle. We still lack an
adaptive and comprehensive understanding of network formation patterns for
predicting potential links. In addition, it is valuable to investigate how
network local information can be better utilized. To this end, we proposed a
novel method named Link prediction using Multiple Order Local Information
(MOLI) that exploits the local information from the neighbors of different
distances, with parameters that can be a prior-driven based on prior knowledge,
or data-driven by solving an optimization problem on observed networks. MOLI
defined a local network diffusion process via random walks on the graph,
resulting in better use of network information. We show that MOLI outperforms
the other 11 widely used link prediction algorithms on 11 different types of
simulated and real-world networks. We also conclude that there are different
patterns of local information utilization for different networks, including
social networks, communication networks, biological networks, etc. In
particular, the classical common neighbor-based algorithm is not as adaptable
to all social networks as it is perceived to be; instead, some of the social
networks obey the Quadrilateral Closure Principle which preferentially connects
paths of length three.


  The infection dynamics of a population under stationary isolation conditions
is modeled. It is underlined that the stationary character of the isolation
measures can be expected to imply that an effective SIR model with constant
parameters should describe the infection process. Then, a derivation of this
property is presented, assuming that the statistical fluctuations in the number
of infection and recovered cases are disregarded. This effective SIR model
shows a reduced population number and a constant $\beta$ parameter. The effects
of also including the retardation between recovery and infection process is
also considered. Next, it is shown that any solution of the effective SIR also
solves the linear problem to which the SIR equations reduce when the total
population is much larger than the number of the infected cases. Then, it is
also argued that this equivalence follows for a specific contagious parameter
$\beta(t)$ which time dependence is analytically derived here. Then, two
equivalent predictive calculational methods for the infection dynamics under
stationary isolation measures are proposed. The results represent a solutions
for the known and challenging problem of defining the time dependence of the
contagion parameter, when the SIR parameter $N$ is assumed to be the whole
population number. Finally, the model is applied to describe the known
infection curves for countries that already had passed the epidemic process
under strict stationary isolation measures. The cases of Iceland, New Zealand,
Korea and Cuba were considered. Although, non subject to stationary isolation
measures the cases of U.S.A. and Mexico are also examined due to their
interest. The results support the argued validity of SIR model including
retardation.


  Monte Carlo estimation in plays a crucial role in stochastic reaction
networks. However, reducing the statistical uncertainty of the corresponding
estimators requires sampling a large number of trajectories. We propose control
variates based on the statistical moments of the process to reduce the
estimators' variances. We develop an algorithm that selects an efficient subset
of infinitely many control variates. To this end, the algorithm uses resampling
and a redundancy-aware greedy selection. We demonstrate the efficiency of our
approach in several case studies.


  We solve MacArthur's resource-competition model with random species-resource
couplings in the `thermodynamic' limit of infinitely many species and resources
using dynamical path-integrals a la De Domincis. We analyze how the steady
state picture changes upon modifying several parameters, including the degree
of heterogeneity of metabolic strategies (encoding the preferences of species)
and of maximal resource levels (carrying capacities), and discuss its
stability. Ultimately, the scenario obtained by other approaches is recovered
by analyzing an effective one-species-one-resource ecosystem that is fully
equivalent to the original multi-species one. The technique used here can be
applied for the analysis of other model ecosystems related to the version of
MacArthur's model considered here.


  Single-cell sequencing has a significant role to explore biological processes
such as embryonic development, cancer evolution, and cell differentiation.
These biological properties can be presented by a two-dimensional scatter plot.
However, single-cell sequencing data generally has very high dimensionality.
Therefore, dimensionality reduction should be used to process the high
dimensional sequencing data for 2D visualization and subsequent biological
analysis. The traditional dimensionality reduction methods, which do not
consider the structure characteristics of single-cell sequencing data, are
difficult to reveal the data structure in the 2D representation. In this paper,
we develop a 2D feature representation method based on graph convolutional
networks (GCN) for the visualization of single-cell data, termed single-cell
sequencing graph embedding networks (SGEN). This method constructs the graph by
the similarity relationship between cells and adopts GCN to analyze the
neighbor embedding information of samples, which makes the similar cell closer
to each other on the 2D scatter plot. The results show SGEN achieves obvious 2D
distribution and preserves the high-dimensional relationship of different
cells. Meanwhile, similar cell clusters have spatial continuity rather than
relying heavily on random initialization, which can reflect the trajectory of
cell development in this scatter plot.


  In recent studies, it has been shown that a cooperative interaction in a
co-infection spread can lead to a discontinuous transition at a decreased
threshold. Here, we investigate effects of immunization with a rate
proportional to the extent of the infection on phase transitions of a
cooperative co-infection. We use the mean-field approximation to illustrate how
measures that remove a portion of the susceptible compartment, like
vaccination, with high enough rates can change discontinuous transitions in two
coupled susceptible-infectious-recovered dynamics into continuous ones while
increasing the threshold of transitions. First, we introduce vaccination with a
fixed rate into a symmetric spread of two diseases and investigate the
numerical results. Second, we set the rate of measures proportional to the size
of the infectious compartment and scrutinize the dynamics. We solve the
equations numerically and analytically and probe the transitions for a wide
range of parameters. We also determine transition points from the analytical
solutions. Third, we adopt a heterogeneous mean-field approach to include
heterogeneity and asymmetry in the dynamics and see if the results
corresponding to homogeneous symmetric case stand.


  In order to model an epidemic, different approaches can be adopted. Mainly,
the deterministic approach and the stochastic one. Recently, a large amount of
literature has been published using the two approaches. The aim of this paper
is to illustrate the usual framework used for commonly adopted compartmental
models in epidemiology and introduce variant analytic and numerical tools that
interfere on each one of those models, as well as the general related types of
existing, ongoing and future possible contributions.


  Enzyme reactions are highly dependent on reaction conditions. To ensure
reproducibility of enzyme reaction parameters, experiments need to be carefully
designed and kinetic modelling meticulously executed. Furthermore, to enable
the judgement of the quality of enzyme reaction parameters, the experimental
conditions, the modelling process as well as the raw data need to be reported
comprehensively. By taking these steps, enzyme reaction parameters can be open
and FAIR (findable, accessible, interoperable, re-usable) as well as
repeatable, replicable and reproducible. This review discusses these issues and
provides a practical guide to designing initial rate experiments for the
determination of enzyme reaction parameters and gives an open, FAIR and
re-editable example of the kinetic modelling of an enzyme reaction. Both the
guide and example are scripted with Python in Jupyter Notebooks and are
publicly available (https://fairdomhub.org/investigations/483). Finally, the
prerequisites of automated data analysis and machine learning algorithms are
briefly discussed to provide further motivation for the comprehensive, open and
FAIR reporting of enzyme reaction parameters.


  Protein dynamics has been investigated on a wide range of time scales. Nano-
and picosecond dynamics have been assigned to local fluctuations, while slower
dynamics have been attributed to larger conformational changes. However, it is
largely unknown how local fluctuations can lead to global allosteric changes.
Here we show that molecule-spanning dynamics on the 100 ns time scale precede
larger allosteric changes. We assign global real-space movements to dynamic
modes on the 100 ns time scales, which became possible by a combination of
single-molecule fluorescence, quasi-elastic neutron scattering and all-atom MD
simulations. Additionally, we demonstrate the effect of Sba1, a co-chaperone of
Hsp90, on these molecule-spanning dynamics, which implies functional importance
of such dynamics. Our integrative approach provides comprehensive insights into
molecule-spanning dynamics on the nanosecond time scale for a multi-domain
protein and indicates that such dynamics are the molecular basis for allostery
and large conformational changes in proteins.


  A recent article by Weidner et al. [2021] presents a method to extract graph
properties that are predictive of the dynamical behavior of multivariate,
discrete models of biochemical regulation. In other words, a method that uses
only features from the structure of network interactions to predict which nodes
are most involved in automata network dynamics. However, the authors claim that
dynamical analysis of large automata network models is "not even feasible." To
make sure that others are not discouraged from working on this problem, it is
important to clarify that effective dynamical analysis of automata network
models, to the contrary, is feasible. Unlike what is suggested in the article,
graph-based analysis of static features is not the only analytical avenue for
large systems biology models of regulation and signaling dynamics because there
are dynamical methods that are, indeed, scalable.


  The automated assembly and extension of dynamic network models using
information extracted from literature are challenging due to the amount and
inconsistency in published literature. Recently, efforts have been made to
automatically and efficiently assemble the information extracted from
literature into models. In this review, we summarize the basic concept,
performance, advantages, and limitations of five automated extension methods.
Each method was tested for its ability to reconstruct a model of T-cell
differentiation as compared against a number of predefined system properties.


  Consider a cooperation game on a spatial network of habitat patches, where
players can relocate between patches if they judge the local conditions to be
unfavorable. In time, the relocation events may lead to a homogeneous state
where all patches harbor the same relative densities of cooperators and
defectors or they may lead to self-organized patterns, where some patches
become safe havens that maintain an elevated cooperator density. Here we
analyze the transition between these states mathematically. We show that safe
havens form once a certain threshold in connectivity is crossed. This threshold
can be analytically linked to the structure of the patch network and
specifically to certain network motifs. Surprisingly, a forgiving defector
avoidance strategy may be most favorable for cooperators. Our results
demonstrate that the analysis of cooperation games in ecological metacommunity
models is mathematically tractable and has the potential to link topics such as
macroecological patterns, behavioral evolution, and network topology.


  Assembling optimal microbial communities is key for various applications in
biofuel production, agriculture, and human health. Finding the optimal
community is challenging because the number of possible communities grows
exponentially with the number of species, and so an exhaustive search cannot be
performed even for a dozen species. A heuristic search that improves community
function by adding or removing one species at a time is more practical, but it
is unknown whether this strategy can discover an optimal or nearly optimal
community. Using consumer-resource models with and without cross-feeding, we
investigate how the efficacy of search depends on the distribution of
resources, niche overlap, cross-feeding, and other aspects of community
ecology. We show that search efficacy is determined by the ruggedness of the
appropriately-defined ecological landscape. We identify specific ruggedness
measures that are both predictive of search performance and robust to noise and
low sampling density. The feasibility of our approach is demonstrated using
experimental data from a soil microbial community. Overall, our results
establish conditions necessary for the success of the heuristic search and
provide concrete design principles for building high-performing microbial
consortia.


  Attractant-based trap networks targeting insects are ubiquitous worldwide.
These networks have diverse targets, goals, and efficiencies, but all are
constrained by practical considerations like cost and available lures. An
important way to balance goals and constrains is through quantitative
mathematical modeling. Here we describe an extension of a computer model of
trapping networks known as "TrapGrid" to include an alternative mode of
calculating the probability of capture over time in a trapping network: Strict
detection ("capture of one or more") compared with the average probability of
capture as implemented in the original version. We suggest that this new
calculation may be useful in situations of low prevalence where trap network
operators wish to interpret the meaning of zero captures at a small scale. The
original remains preferred for comparing the sensitivity and suitability of
alternate trap networks (i.e. density of traps, their placement, lure
attractiveness, etc.).


  Preserving the privacy and security of big data in the context of cloud
computing, while maintaining a certain level of efficiency of its processing
remains to be a subject, open for improvement. One of the most popular
applications epitomizing said concerns is found to be useful in genome
analysis. This work proposes a secure multi-label tumor classification method
using homomorphic encryption, whereby two different machine learning
algorithms, SVM and XGBoost, are used to classify the encrypted genome data of
different tumor types.


  Phylogenetic networks extend phylogenetic trees to model non-vertical
inheritance, by which a lineage inherits material from multiple parents. The
computational complexity of estimating phylogenetic networks from genome-wide
data with likelihood-based methods limits the size of networks that can be
handled. Methods based on pairwise distances could offer faster alternatives.
We study here the information that average pairwise distances contain on the
underlying phylogenetic network, by characterizing local and global features
that can or cannot be identified. For general networks, we clarify that the
root and edge lengths adjacent to reticulations are not identifiable, and then
focus on the class of zipped-up semidirected networks. We provide a criterion
to swap subgraphs locally, such as 3-cycles, resulting in indistinguishable
networks. We propose the "distance split tree", which can be constructed from
pairwise distances, and prove that it is a refinement of the network's tree of
blobs, capturing the tree-like features of the network. For level-1 networks,
this distance split tree is equal to the tree of blobs refined to separate
polytomies from blobs, and we prove that the mixed representation of the
network is identifiable. The information loss is localized around 4-cycles, for
which the placement of the reticulation is unidentifiable. The mixed
representation combines split edges for 4-cycles, regular tree and hybrid edges
from the semidirected network, and edge parameters that encode all information
identifiable from average pairwise distances.


  Understanding the conditions of feasibility and stability in ecological
systems is a major challenge in theoretical ecology. The seminal work of May in
1972 and recent developments based on the theory of random matrices have shown
the existence of emergent universal patterns of both stability and feasibility
in ecological dynamics. However, only a few studies have investigated the role
of delay coupled with population dynamics in the emergence of feasible and
stable states. In this work, we study the effects of delay on Generalized
Loka-Volterra population dynamics of several interacting species in closed
ecological environments. First, we investigate the relation between feasibility
and stability of the modeled ecological community in the absence of delay and
find a simple analytical relation when intra-species interactions are dominant.
We then show how, by increasing the time delay, there is a transition in the
stability phases of the population dynamics: from an equilibrium state to a
stable non-point attractor phase. We calculate analytically the critical delay
of that transition and show that it is in excellent agreement with numerical
simulations. Finally, we introduce a measure of stability that holds for out of
equilibrium dynamics and we show that in the oscillatory regime induced by the
delay stability increases for increasing ecosystem diversity.


  The native structures of proteins, except for notable exceptions of
intrinsically disordered proteins, in general take their most stable
conformation in the physiological condition to maintain their structural
framework so that their biological function can be properly carried out.
Experimentally, the stability of a protein can be measured by several means,
among which the pulling experiment using the atomic force microscope (AFM)
stands as a unique method. AFM directly measures the resistance from unfolding,
which can be quantified from the observed force-extension profile. It has been
shown that key features observed in an AFM pulling experiment can be well
reproduced by computational molecular dynamics simulations. Here, we applied
computational pulling for estimating the accuracy of computational protein
structure models under the hypothesis that the structural stability would
positively correlated with the accuracy, i.e. the closeness to the native, of a
model. We used in total 4,929 structure models for 24 target proteins from the
Critical Assessment of Techniques of Structure Prediction (CASP) and
investigated if the magnitude of the break force, i.e., the force required to
rearrange the model structure, from the force profile was sufficient
information for selecting near-native models. We found that near-native models
can be successfully selected by examining their break forces suggesting that
high break force indeed indicates high stability of models. On the other hand,
there were also near-native models that had relatively low peak forces. The
mechanisms of the stability exhibited by the break forces were explored and
discussed.


  While analyzing the persistence length of DNA atomistic simulations Wu et al.
[J. Chem. Phys. 142, 125103 (2015)] introduced an empirical formula to account
for the observed length-dependence. In particular they found that the
persistence length increases with the distance. Here, we derive the formula by
Wu et al. using a non-local twistable wormlike chain which introduces couplings
between distal sites. Finally, we show that the same formula can account for
the length-scale dependence of the torsional persistence length and is, in
fact, applicable to any kind of polymer model with non-local couplings.


  Simulations of protein folding and protein association happen on timescales
that are orders of magnitude larger than what can typically be covered in
all-atom molecular dynamics simulations. Use of low-resolution models
alleviates this problem but may reduce the accuracy of the simulations. We
introduce a replica-exchange-based multiscale sampling technique that combines
the faster sampling in coarse-grained simulations with the potentially higher
accuracy of all-atom simulations. After testing the efficiency of our
Resolution Exchange with Tunneling (ResET) in simulations of the Trp-cage
protein, an often used model to evaluate sampling techniques in protein
simulations, we use our approach to compare the landscape of wild type and A2T
mutant Abeta 1-42 peptides. Our results suggest a mechanism by that the
mutation of a small hydrophobic Alanine (A) into a bulky polar Threonine (T)
may interfere with the self-assembly of Abeta - fibrils.


  The increasing prevalence of SARS-CoV-2 variants with spike mutations has
raised concerns owing to higher transmission rates, disease severity, and
escape from neutralizing antibodies. Rapid and accurate detection of SARS-CoV-2
variants provides crucial information concerning the outbreaks of SARS-CoV-2
variants and possible lines of transmission. This information is vital for
infection prevention and control. We used a Cas12a-based RT-PCR combined with
CRISPR on-site rapid detection system (RT-CORDS) platform to detect the key
mutations in SARS-COV-2 variants, such as 69/70 deletion, N501Y, and D614G. We
used type-specific CRISPR RNAs (crRNAs) to identify wild-type (crRNA-W) and
mutant (crRNA-M) sequences of SARS-CoV-2. We successfully differentiated mutant
variants from wild-type SARS-CoV-2 with a sensitivity of $10^{-17}$ M
(approximately 6 copies/$\mu$L). The assay took just 10 min with the
Cas12a/crRNA reaction after a simple RT-PCR using a fluorescence reporting
system. In addition, a sensitivity of $10^{-16}$ M could be achieved when
lateral flow strips were used as readouts. The accuracy of RT-CORDS for
SARS-CoV-2 variant detection was 100% consistent with the sequencing data. In
conclusion, using the RT-CORDS platform, we accurately, sensitively,
specifically, and rapidly detected SARS-CoV-2 variants. This method may be used
in clinical diagnosis.


  Mitochondria in plant cells form strikingly dynamic populations of largely
individual organelles. Each mitochondrion contains on average less than a full
copy of the mitochondrial DNA (mtDNA) genome. Here, we asked whether
mitochondrial dynamics may allow individual mitochondria to `collect' a full
copy of the mtDNA genome over time, by facilitating exchange between
individuals. Akin to trade on a social network, exchange of mtDNA fragments
across organelles may lead to the emergence of full `effective' genomes in
individuals over time. We characterise the collective dynamics of mitochondria
in \emph{Arabidopsis thaliana} hypocotyl cells using a recent approach
combining single-cell timelapse microscopy, video analysis, and network
science. We then use a quantitative model to predict the capacity for the
sharing and accumulation of genetic information through the networks of
encounters between mitochondria. We find that biological encounter networks are
strikingly well predisposed to support the collection of full genomes over
time, outperforming a range of other networks generated from theory and
simulation. Using results from the coupon collector's problem, we show that the
upper tail of the degree distribution is a key determinant of an encounter
network's performance at this task and discuss how features of mitochondrial
dynamics observed in biology facilitate the emergence of full effective
genomes.


  Applications of single-cell RNA sequencing in various biomedical research
areas have been blooming. This new technology provides unprecedented
opportunities to study disease heterogeneity at the cellular level. However,
unique characteristics of scRNA-seq data, including large dimensionality, high
dropout rates, and possibly batch effects, bring great difficulty into the
analysis of such data. Not appropriately addressing these issues obstructs true
scientific discovery. Herein, we propose a unified Regularized Zero-inflated
Mixture Model framework designed for scRNA-seq data (RZiMM-scRNA) to
simultaneously detect cell subgroups and identify gene differential expression
based on a developed importance score, accounting for both dropouts and batch
effects. We conduct extensive simulation studies in which we evaluate the
performance of RZiMM-scRNA and compare it with several popular methods,
including Seurat, SC3, K-Means, and Hierarchical Clustering. Simulation results
show that RZiMM-scRNA demonstrates superior clustering performance and enhanced
biomarker detection accuracy compared to alternative methods, especially when
cell subgroups are less distinct, verifying the robustness of our method. Our
empirical investigations focus on two brain tumor studies dealing with
astrocytoma of various grades, including the most malignant of all brain
tumors, glioblastoma multiforme (GBM). Our goal is to delineate cell
heterogeneity and identify driving biomarkers associated with these tumors.
Notably, RZiMM-scNRA successfully identifies a small group of oligodendrocyte
cells which has drawn much attention in biomedical literature on brain cancers.


  Divergence time estimation is crucial to provide temporal signals for dating
biologically important events, from species divergence to viral transmissions
in space and time. With the advent of high-throughput sequencing, recent
Bayesian phylogenetic studies have analyzed hundreds to thousands of sequences.
Such large-scale analyses challenge divergence time reconstruction by requiring
inference on highly-correlated internal node heights that often become
computationally infeasible. To overcome this limitation, we explore a ratio
transformation that maps the original N - 1 internal node heights into a space
of one height parameter and N - 2 ratio parameters. To make analyses scalable,
we develop a collection of linear-time algorithms to compute the gradient and
Jacobian-associated terms of the log-likelihood with respect to these ratios.
We then apply Hamiltonian Monte Carlo sampling with the ratio transform in a
Bayesian framework to learn the divergence times in four pathogenic virus
phylogenies: West Nile virus, rabies virus, Lassa virus and Ebola virus. Our
method both resolves a mixing issue in the West Nile virus example and improves
inference efficiency by at least 5-fold for the Lassa and rabies virus
examples. Our method also makes it now computationally feasible to incorporate
mixed-effects molecular clock models for the Ebola virus example, confirms the
findings from the original study and reveals clearer multimodal distributions
of the divergence times of some clades of interest.


  Both conventional wisdom and empirical evidence suggests that arranging a
prior commitment or agreement before an interaction enhances the chance of
reaching mutual cooperation. Yet it is not clear what mechanisms can promote
the participation in and compliance with such a commitment, especially when the
former is costly and deviating from the latter is profitable. Prior work either
considers regimented commitments where compensation is assumed enforceable from
dishonest committers, or assume implicit commitments from every individual (so
they are all in and thus being treated as such). Here we develop a theory of
participation and compliance with respect to an explicit prior commitment under
institutional incentives where individuals, at first, decide whether or not to
join a cooperative agreement to play a one-shot social dilemma game. Using a
mathematical model, we determine when participating in a costly commitment and
complying with it, is an evolutionary stable strategy (ESS) when playing
against all other possible strategies, and results in high levels of
cooperation in the population. We show that, given a sufficient budget for
providing incentives, reward of commitment compliant behaviours better promotes
cooperation than punishment of non-compliant ones. Moreover, by sparing part of
this budget for rewarding those who are willing to participate in a commitment,
the overall frequency of cooperation can be significantly enhanced, for both
reward and punishment. Finally, we find that, surprisingly, the presence of
errors in a participation decision favours evolutionary stability of commitment
compliant strategies and higher levels of cooperation.


  Interactions in natural communities can be highly heterogeneous, with any
given species interacting appreciably with only some of the others, a situation
commonly represented by sparse interaction networks. We study the consequences
of sparse competitive interactions, in a theoretical model of a community
assembled from a species pool. We find that communities can be in a number of
different regimes, depending on the interaction strength. When interactions are
strong, the network of coexisting species breaks up into small subgraphs, while
for weaker interactions these graphs are larger and more complex, eventually
encompassing all species. This process is driven by emergence of new allowed
subgraphs as interaction strength decreases, leading to sharp changes in
diversity and other community properties, and at weaker interactions to two
distinct collective transitions: a percolation transition, and a transition
between having a unique equilibrium and having multiple alternative equilibria.
Understanding community structure is thus made up of two parts: first, finding
which subgraphs are allowed at a given interaction strength, and secondly, a
discrete problem of matching these structures over the entire community. In a
shift from the focus of many previous theories, these different regimes can be
traversed by modifying the interaction strength alone, without need for
heterogeneity in either interaction strengths or the number of competitors per
species.


  National responses to the Covid-19 pandemic varied markedly across countries,
from business-as-usual to complete shutdowns. Policies aimed at disrupting the
viral transmission cycle and preventing the healthcare system from being
overwhelmed, simultaneously exact an economic toll. We developed a intervention
policy model that comprised the relative human, economic and healthcare costs
of non-pharmaceutical epidemic intervention and arrived at the optimal strategy
using the neuroevolution algorithm. The proposed model finds the minimum
required reduction in contact rates to maintain the burden on the healthcare
system below the maximum capacity. We find that such a policy renders a sharp
increase in the control strength at the early stages of the epidemic, followed
by a steady increase in the subsequent ten weeks as the epidemic approaches its
peak, and finally control strength is gradually decreased as the population
moves towards herd immunity. We have also shown how such a model can provide an
efficient adaptive intervention policy at different stages of the epidemic
without having access to the entire history of its progression in the
population. This work emphasizes the importance of imposing intervention
measures early and provides insights into adaptive intervention policies to
minimize the economic impacts of the epidemic without putting an extra burden
on the healthcare system.


  We present a simple model for the spread of an infection that incorporates
spatial variability in population density. Starting from first principle
considerations, we explore how a novel PDE with state-dependent diffusion can
be obtained. This model exhibits higher infection rates in the areas of higher
population density, a feature that we argue to be consistent with
epidemiological observations. The model also exhibits an infection wave whose
speed varies with population density. In addition, we demonstrate the
possibility that an infection can ``jump'' (i.e., tunnel) across areas of low
population density towards the areas of high population density. We briefly
touch upon the data reported for coronavirus spread in the Canadian province of
Nova Scotia as a case example with a number of qualitatively similar features
as our model. Lastly, we propose a number of generalizations of the model
towards future studies.


  We consider a natural class of reaction networks which consist of reactions
where either two species can inactivate each other (i.e., sequestration), or
some species can be transformed into another (i.e., transmutation), in a way
that gives rise to a feedback cycle. We completely characterize the capacity of
multistationarity of these networks. This is especially interesting because
such networks provide simple examples of "atoms of multistationarity", i.e.,
minimal networks that can give rise to multiple positive steady states


  Models of stochastic processes are widely used in almost all fields of
science. Theory validation, parameter estimation, and prediction all require
model calibration and statistical inference using data. However, data are
almost always incomplete observations of reality. This leads to a great
challenge for statistical inference because the likelihood function will be
intractable for almost all partially observed stochastic processes. This
renders many statistical methods, especially within a Bayesian framework,
impossible to implement. Therefore, computationally expensive likelihood-free
approaches are applied that replace likelihood evaluations with realisations of
the model and observation process. For accurate inference, however,
likelihood-free techniques may require millions of expensive stochastic
simulations. To address this challenge, we develop a new method based on recent
advances in multilevel and multifidelity. Our approach combines the multilevel
Monte Carlo telescoping summation, applied to a sequence of approximate
Bayesian posterior targets, with a multifidelity rejection sampler to minimise
the number of computationally expensive exact simulations required for accurate
inference. We present the derivation of our new algorithm for likelihood-free
Bayesian inference, discuss practical implementation details, and demonstrate
substantial performance improvements. Using examples from systems biology, we
demonstrate improvements of more than two orders of magnitude over standard
rejection sampling techniques. Our approach is generally applicable to
accelerate other sampling schemes, such as sequential Monte Carlo, to enable
feasible Bayesian analysis for realistic practical applications in physics,
chemistry, biology, epidemiology, ecology and economics.


  Cox analysis is a common clinical data analysis technique to link valuable
variables to clinical outcomes including dead and relapse. In the omics era,
Cox model batch processing is a basic strategy for screening clinically
relevant variables, biomarker discovery and gene signature identification.
However, all such analyses have been implemented with homebrew code in research
community, thus lack of transparency and reproducibility. Here, we present
ezcox, the first R/CRAN package for Cox model batch processing and
visualization. ezcox is an open source R package under GPL-3 license and it is
free available at https://github.com/ShixiangWang/ezcox and
https://cran.r-project.org/package=ezcox.


  The vaccination game is a social dilemma that refers to the conundrum
individuals face (to get immunized or not) when the population is exposed to an
infectious disease. The model has recently gained much traction due to the
COVID-19 pandemic since the public perception of vaccines plays a significant
role in disease dynamics. This paper studies the vaccination game in the
thermodynamic limit with an analytical method derived from the 1D Ising model
called Nash equilibrium mapping. The individual dilemma regarding Vaccination
comes from an internal conflict wherein one tries to balance the perceived
advantages of immunizing with the apparent risks associated with Vaccination
which they hear through different news media. We compare the results of Nash
equilibrium(NE) mapping to other 1D Ising-based models, namely Darwinian
evolution and agent-based simulation. This study aims to analyze the behaviour
of an infinite population regarding what fraction of people choose to vaccinate
or not vaccinate. While Nash equilibrium mapping and agent-based simulation
agree mostly, Darwinian evolution strays far from the two models. It fails to
predict the equilibrium behaviour of players in the population reasonably. We
apply the results of our study to analyze the Astra-Zeneca(AZ) COVID-19 vaccine
risk versus disease deaths debate, both via NE mapping and agent-based method.
Both predict near 100% AZ vaccine coverage for people above 40, notwithstanding
the risk. At the same time, younger people show a slight reluctance. We predict
that while government intervention via vaccination mandates and advertisement
campaigns is unnecessary for the older population, for the younger population
(ages: 20-39), some encouragement from the government via media campaigns and
vaccine mandates may be necessary.


  It is essential to understand the dynamics of epidemics in the presence of
coexisting pathogens. There are various phenomenon that can effect the
dynamics. In this paper, we formulate a mathematical model using different
assumptions to capture the effect of these additional phenomena such as partial
cross immunity, density dependence in each class and a role of recovered
population in the dynamics. We found the basic reproduction number for each
model which is the threshold that describes the invasion of disease in
population. The basic reproduction number in each model shows that the
persistence of disease or strains depends on the carrying capacity. In the
model of this paper, we present the local stability analysis of the boundary
equilibrium points and observed that the recovered population is not uniformly
bounded with respect to the carrying capacity.


  Uncertainty in data is certainly one of the main problems in epidemiology, as
shown by the recent COVID-19 pandemic. The need for efficient methods capable
of quantifying uncertainty in the mathematical model is essential in order to
produce realistic scenarios of the spread of infection. In this paper, we
introduce a bi-fidelity approach to quantify uncertainty in spatially dependent
epidemic models. The approach is based on evaluating a high-fidelity model on a
small number of samples properly selected from a large number of evaluations of
a low-fidelity model. In particular, we will consider the class of multiscale
transport models recently introduced in Bertaglia, Boscheri, Dimarco &
Pareschi, Math. Biosci. Eng. (2021) and Boscheri, Dimarco & Pareschi, Math.
Mod. Meth. App. Scie. (2021) as the high-fidelity reference and use simple
two-velocity discrete models for low-fidelity evaluations. Both models share
the same diffusive behavior and are solved with ad-hoc asymptotic-preserving
numerical discretizations. A series of numerical experiments confirm the
validity of the approach.


  We apply the theory of learning to physically renormalizable systems in an
attempt to develop a theory of biological evolution, including the origin of
life, as multilevel learning. We formulate seven fundamental principles of
evolution that appear to be necessary and sufficient to render a universe
observable and show that they entail the major features of biological
evolution, including replication and natural selection. These principles also
follow naturally from the theory of learning. We formulate the theory of
evolution using the mathematical framework of neural networks, which provides
for detailed analysis of evolutionary phenomena. To demonstrate the potential
of the proposed theoretical framework, we derive a generalized version of the
Central Dogma of molecular biology by analyzing the flow of information during
learning (back-propagation) and predicting (forward-propagation) the
environment by evolving organisms. The more complex evolutionary phenomena,
such as major transitions in evolution, in particular, the origin of life, have
to be analyzed in the thermodynamic limit, which is described in detail in the
accompanying paper.


  Various machine-learning models, including deep neural network models, have
already been developed to predict deleteriousness of missense (non-synonymous)
mutations. Potential improvements to the current state of the art, however, may
still benefit from a fresh look at the biological problem using more
sophisticated self-adaptive machine-learning approaches. Recent advances in the
natural language processing field show transformer models-a type of deep neural
network-to be particularly powerful at modeling sequence information with
context dependence. In this study, we introduce MutFormer, a transformer-based
model for the prediction of deleterious missense mutations, which uses
reference and mutated protein sequences from the human genome as the primary
features. MutFormer takes advantage of a combination of self-attention layers
and convolutional layers to learn both long-range and short-range dependencies
between amino acid mutations in a protein sequence. In this study, we first
pre-trained MutFormer on reference protein sequences and mutated protein
sequences resulting from common genetic variants observed in human populations.
We next examined different fine-tuning methods to successfully apply the model
to deleteriousness prediction of missense mutations. Finally, we evaluated
MutFormer's performance on multiple testing data sets. We found that MutFormer
showed similar or improved performance over a variety of existing tools,
including those that used conventional machine-learning approaches. We conclude
that MutFormer successfully considers sequence features that are not explored
in previous studies and could potentially complement existing computational
predictions or empirically generated functional scores to improve our
understanding of disease variants.


  We introduce an extended SEIR infectious disease model with data assimilation
for the study of the spread of COVID-19. In this framework, undetected
asymptomatic and pre-symptomatic cases are taken into account, and the impact
of their uncertain proportion is fully investigated. The standard SEIR model
does not consider these populations, while their role in the propagation of the
disease is acknowledged. An ensemble Kalman filter is implemented to assimilate
reliable observations of three compartments in the model. The system tracks the
evolution of the effective reproduction number and estimates the unobservable
subpopulations. The analysis is carried out for three main prefectures of Japan
and for the entire population of Japan. For these four populations, our
estimated effective reproduction numbers are more stable than the corresponding
ones estimated by a different method (Toyokeizai). We also perform sensitivity
tests for different values of some uncertain medical parameters, like the
relative infectivity of symptomatic / asymptomatic cases. The regional analysis
results suggest the decreasing efficiency of the states of emergency.


  Keeping a memory of evolving stimuli is ubiquitous in biology, an example of
which is immune memory for evolving pathogens. However, learning and memory
storage for dynamic patterns still pose challenges in machine learning. Here,
we introduce an analytical energy-based framework to address this problem. By
accounting for the tradeoff between utility in keeping a high-affinity memory
and the risk in forgetting some of the diverse stimuli, we show that a moderate
tolerance for risk enables a repertoire to robustly classify evolving patterns,
without much fine-tuning. Our approach offers a general guideline for learning
and memory storage in systems interacting with diverse and evolving stimuli.


  We outline a phenomenological theory of evolution and origin of life by
combining the formalism of classical thermodynamics with a statistical
description of learning. The maximum entropy principle constrained by the
requirement for minimization of the loss function is employed to derive a
canonical ensemble of organisms (population), the corresponding partition
function (macroscopic counterpart of fitness) and free energy (macroscopic
counterpart of additive fitness). We further define the biological counterparts
of temperature (biological temperature) as the measure of stochasticity of the
evolutionary process and of chemical potential (evolutionary potential) as the
amount of evolutionary work required to add a new trainable variable (such as
an additional gene) to the evolving system. We then develop a phenomenological
approach to the description of evolution, which involves modeling the grand
potential as a function of the biological temperature and evolutionary
potential. We demonstrate how this phenomenological approach can be used to
study the "ideal mutation" model of evolution and its generalizations. Finally,
we show that, within this thermodynamics framework, major transitions in
evolution, such as the transition from an ensemble of molecules to an ensemble
of organisms, that is, the origin of life, can be modeled as a special case of
bona fide physical phase transitions that are associated with the emergence of
a new type of grand canonical ensemble and the corresponding new level of
description


  We present quantitative measurements on a classic experiment proposed for the
first time in 1947 to illustrate the phenomenon of a chain reaction in nuclear
fission. The experiment involves a number of mousetraps loaded with solid
balls. Once one trap is made to snap, it releases two balls that cause other
traps to snap more traps and the result is a chain reaction. We report for the
first time that the mousetrap experiment can be fitted by a simple dynamic
model. We also discuss the significance of this experiment beyond nuclear chain
reactions, as providing insight in a variety of fields involving complex,
adaptive systems.