  It is well-known that the Hessian of deep loss landscape matters to
optimization, generalization, and even robustness of deep learning. Recent
works empirically discovered that the Hessian spectrum in deep learning has a
two-component structure that consists of a small number of large eigenvalues
and a large number of nearly-zero eigenvalues. However, the theoretical
mechanism or the mathematical behind the Hessian spectrum is still largely
under-explored. To the best of our knowledge, we are the first to demonstrate
that the Hessian spectrums of well-trained deep neural networks exhibit simple
power-law structures. Inspired by the statistical physical theories and the
spectral analysis of natural proteins, we provide a maximum-entropy theoretical
interpretation for explaining why the power-law structure exist and suggest a
spectral parallel between protein evolution and training of deep neural
networks. By conducing extensive experiments, we further use the power-law
spectral framework as a useful tool to explore multiple novel behaviors of deep
learning.


  Gene transcription is a stochastic process that involves thousands of
reactions. The first set of these reactions, which happen near a gene promoter,
are considered to be the most important in the context of stochastic noise. The
most common models of transcription are primarily concerned with the effect of
activators/repressors on the overall transcription rate and approximate the
basal transcription processes as a one step event. According to such effective
models, the Fano factor of mRNA copy distributions is always greater than
(super-Poissonian) or equal to 1 (Poissonian), and the only way to go below
this limit (sub-Poissonian) is via a negative feedback. It is partly due to
this limit that the first stage of transcription is held responsible for most
of the stochastic noise in mRNA copy numbers. However, by considering all major
reactions that build and drive the basal transcription machinery, from the
first protein that binds a promoter to the entrance of the transcription
complex (TC) into productive elongation, it is shown that the first two stages
of transcription, namely the pre-initiation complex (PIC) formation and the
promoter proximal pausing (PPP), is a highly punctual process. In other words,
the time between the first and the last step of this process is narrowly
distributed, which gives rise to sub-Poissonian distributions for the number of
TCs that have entered productive elongation. In fact, having simulated the PIC
formation and the PPP via the Gillespie algorithm using 2000 distinct parameter
sets and 4 different reaction network topologies, it is shown that only 4.4%
give rise to a Fano factor that is > 1 with the upper bound of 1.7, while for
31% of cases the Fano factor is below 0.5, with 0.19 as the lower bound. These
results cast doubt on the notion that most of the stochastic noise observed in
mRNA distributions always originates at the promoter.


  A family of results, referred to as inheritance results, tell us which
enlargements of a chemical reaction network (CRN) preserve its capacity for
nontrivial behaviours such as multistationarity and oscillation. In this paper,
the following inheritance result is proved: under mild assumptions, splitting
chemical reactions and inserting complexes involving some new chemical species
preserves the capacity of a mass action CRN for multiple nondegenerate
equilibria and/or periodic orbits. The claim has been proved previously for
equilibria alone; however, the generalisation to include oscillation involves
extensive development of rather different techniques. Several inheritance
results for multistationarity and oscillation in mass action CRNs, including
the main result of this paper, are gathered into a single theorem. Examples are
presented showing how these results can be used together to make claims about
reaction networks based on knowledge of their subnetworks. The examples include
some networks of biological importance.


  We are concerned with the problem of decomposing the parameter space of a
parametric system of polynomial equations, and possibly some polynomial
inequality constraints, with respect to the number of real solutions that the
system attains. Previous studies apply a two step approach to this problem,
where first the discriminant variety of the system is computed via a Groebner
Basis (GB), and then a Cylindrical Algebraic Decomposition (CAD) of this is
produced to give the desired computation. However, even on some reasonably
small applied examples this process is too expensive, with computation of the
discriminant variety alone infeasible. In this paper we develop new approaches
to build the discriminant variety using resultant methods (the Dixon resultant
and a new method using iterated univariate resultants). This reduces the
complexity compared to GB and allows for a previous infeasible example to be
tackled. We demonstrate the benefit by giving a symbolic solution to a problem
from population dynamics -- the analysis of the steady states of three
connected populations which exhibit Allee effects - which previously could only
be tackled numerically.


  A protein performs biological functions by folding to a particular 3D
structure. To accurately model the protein structures, both the overall
geometric topology and local fine-grained relations between amino acids (e.g.
side-chain torsion angles and inter-amino-acid orientations) should be
carefully considered. In this work, we propose the Directed Weight Neural
Network for better capturing geometric relations among different amino acids.
Extending a single weight from a scalar to a 3D directed vector, our new
framework supports a rich set of geometric operations on both classical and
SO(3)--representation features, on top of which we construct a perceptron unit
for processing amino-acid information. In addition, we introduce an equivariant
message passing paradigm on proteins for plugging the directed weight
perceptrons into existing Graph Neural Networks, showing superior versatility
in maintaining SO(3)-equivariance at the global scale. Experiments show that
our network has remarkably better expressiveness in representing geometric
relations in comparison to classical neural networks and the (globally)
equivariant networks. It also achieves state-of-the-art performance on various
computational biology applications related to protein 3D structures.


  The innocuousness property of a controller is that property that makes the
closed-loop system stable regardless the values of the controller parameters.
In other words, the closed-loop system exhibits some structural stability
property with respect to the parameters of the controller. The innocuousness
property was first emphasized in [Briat, Gupta, and Khammash, Cell Systems,
2016] where it was shown that for stochastic unimolecular networks, the
Antithetic Integral Controller (AIC) is innocuous under very mild conditions on
the controlled network; namely the ergodicity and the output-controllability of
the open-loop network, and the admissibility of the set-point value. We show
here that the class of p-type AIC controllers arising in the use of Antithetic
Integral Rein Controllers (AIRC) also exhibit such a property. It is shown in
the unimolecular reaction network case that the closed-loop network is
structurally stable with respect to the controller parameters provided that the
open-loop dynamics is stable and that the set-point is admissible. Those
results are then extended to the case of so-called output unstable linear
systems and to stable nonlinear networks. Analogous results are obtained for
exponential and logistic integral controllers. Some examples are given for
illustration.


  Gene set analysis methods rely on knowledge-based representations of genetic
interactions in the form of both gene set collections and protein-protein
interaction (PPI) networks. Explicit representations of genetic interactions
often fail to capture complex interdependencies among genes, limiting the
analytic power of such methods. Here we propose an extension of gene set
enrichment analysis to a latent feature space reflecting PPI network topology,
called gene set proximity analysis (GSPA). Compared with existing methods, GSPA
provides improved ability to identify disease-associated pathways in
disease-matched gene expression datasets, while improving reproducibility of
enrichment statistics for similar gene sets. GSPA is statistically
straightforward, reducing to classical gene set enrichment through a single
user-defined parameter. We apply our method to identify novel drug associations
with SARS-CoV-2 viral entry. Finally, we validate our drug association
predictions through retrospective clinical analysis of claims data from 8
million patients, supporting a role for gabapentin as a risk factor and
metformin as a protective factor for COVID-19 hospitalization.


  Antipredator behaviour is a self-preservation strategy present in many
biological systems, where individuals join the effort in a collective reaction
to avoid being caught by an approaching predator. We study a nonhierarchical
tritrophic system, whose predator-prey interactions are described by the
rock-paper-scissors game rules. We performed a set of spatial stochastic
simulations where organisms of one out of the species can resist predation in a
collective strategy. The drop in predation capacity is local, which means that
each predator faces a particular opposition depending on the prey group size
surrounding it. Considering that the interference in a predator action depends
on the prey's physical and cognitive ability, we explore the role of a
conditioning factor that indicates the fraction of the species apt to perform
the antipredator strategy. Because of the local unbalancing of the cyclic
predator-prey interactions, departed spatial domains mainly occupied by a
single species emerge. Unlike the rock-paper-scissors model with a weak species
because a nonlocal reason, our findings show that if the predation probability
of one species is reduced because individuals face local antipredator response,
the species does not predominate. Instead, the local unbalancing of the
rock-paper-scissors model results in the prevalence of the weak species' prey.
Finally, the outcomes show that local unevenness may jeopardise biodiversity,
with the coexistence being more threatened for high mobility.


  We derive sufficient conditions for the existence of a periodic traveling
wave solution to an integro-difference equation with a piecewise constant
growth function exhibiting a stable period2 cycle and strong Allee effect. The
mean traveling wave speed is shown to be the asymptotic spreading speed of
solutions with compactly supported initial data under appropriate conditions.
We then conduct case studies for the Laplace kernel and uniform kernel.


  An epidemic model where disease transmission can occur either through global
contacts or through local, nearest neighbor interactions is considered. The
classical SIR--model describing the global interactions is extended by adding
additional equations for the density of local pairs in different epidemic
states. A locality parameter $p\in [0,1]$ characterizes the probability of
global or local infections. The equilibria of the resulting model are analyzed
in dependence of the locality parameter and the transmission rate of the
pathogen. An explicit expression for the reproduction number in terms of the
locality parameter and the disease parameters is obtained. Transient
simulations confirm these findings. Neighboring pairs of one infected and one
susceptible can be considered as active pairs, since local transmission of the
disease can only occur in that situation. Our analysis shows, that the fraction
of active pairs is minimal for intermediate values of the locality parameter.


  Nowadays there is a big spotlight cast on the development of techniques of
explainable machine learning. Here we introduce a new computational paradigm
based on Group Equivariant Non-Expansive Operators, that can be regarded as the
product of a rising mathematical theory of information-processing observers.
This approach, that can be adjusted to different situations, may have many
advantages over other common tools, like Neural Networks, such as: knowledge
injection and information engineering, selection of relevant features, small
number of parameters and higher transparency. We chose to test our method,
called GENEOnet, on a key problem in drug design: detecting pockets on the
surface of proteins that can host ligands. Experimental results confirmed that
our method works well even with a quite small training set, providing thus a
great computational advantage, while the final comparison with other
state-of-the-art methods shows that GENEOnet provides better or comparable
results in terms of accuracy.


  Patient-derived extracellular vesicles (EVs) that contains a complex
biological cargo is a valuable source of liquid biopsy diagnostics to aid in
early detection, cancer screening, and precision nanotherapeutics. In this
study, we predicted that coupling cancer patient blood-derived EVs to
time-resolved spectroscopy and artificial intelligence (AI) could provide a
robust cancer screening and follow-up tools. Methods: Fluorescence correlation
spectroscopy (FCS) measurements were performed on 24 blood samples-derived EVs.
Blood samples were obtained from 15 cancer patients (presenting 5 different
types of cancers), and 9 healthy controls (including patients with benign
lesions). The obtained FCS autocorrelation spectra were processed into power
spectra using the Fast-Fourier Transform algorithm and subjected to various
machine learning algorithms to distinguish cancer spectra from healthy control
spectra. Results and Applications: The performance of AdaBoost Random Forest
(RF) classifier, support vector machine, and multilayer perceptron, were tested
on selected frequencies in the N=118 power spectra. The RF classifier exhibited
a 90% classification accuracy and high sensitivity and specificity in
distinguishing the FCS power spectra of cancer patients from those of healthy
controls. Further, an image convolutional neural network (CNN), ResNet network,
and a quantum CNN were assessed on the power spectral images as additional
validation tools. All image-based CNNs exhibited a nearly equal classification
performance with an accuracy of roughly 82% and reasonably high sensitivity and
specificity scores. Our pilot study demonstrates that AI-algorithms coupled to
time-resolved FCS power spectra can accurately and differentially classify the
complex patient-derived EVs from different cancer samples of distinct tissue
subtypes.


  It is widely accepted that the number of reported cases during the first
stages of the COVID-19 pandemic severely underestimates the number of actual
cases. We leverage delay embedding theorems of Whitney and Takens and use
Gaussian Process regression to estimate the number of cases during the first
2020 wave based on the second wave of the epidemic in several European
countries, South Korea, and Brazil. We assume that the second wave was more
accurately monitored and hence that it can be trusted. We then construct a
manifold diffeomorphic to that of the implied original dynamical system, using
fatalities or hospitalizations only. Finally, we restrict the diffeomorphism to
the reported cases coordinate of the dynamical system. Our main finding is that
in the European countries studied, the actual cases are under-reported by as
much as 50\%. On the other hand, in South Korea -- which had an exemplary and
proactive mitigation approach -- a far smaller discrepancy between the actual
and reported cases is predicted, with an approximately 17\% predicted
under-estimation. We believe that our backcasting framework is applicable to
other epidemic outbreaks where (due to limited or poor quality data) there is
uncertainty around the actual cases.


  We develop a mathematical model for transferring the vaccine BNT162b2 based
on the heat diffusion equation. Then, we apply optimal control theory to the
proposed generalized SEIR model. We introduce vaccination for the susceptible
population to control the spread of the COVID-19 epidemic. For this, we use the
Pontryagin minimum principle to find the necessary optimality conditions for
the optimal control. The optimal control problem and the heat diffusion
equation are solved numerically. Finally, several simulations are done to study
and predict the spread of the COVID-19 epidemic in Italy. In particular, we
compare the model in the presence and absence of vaccination.


  We propose and study a shifted SICA epidemic model, extending the one of
Silva and Torres (2017) to the stochastic setting driven by both Brownian
motion processes and jump L\'evy noise. L\'evy noise perturbations are usually
ignored by existing works of mathematical modelling in epidemiology, but its
incorporation into the SICA epidemic model is worth to consider because of the
presence of strong fluctuations in HIV/AIDS dynamics, often leading to the
emergence of a number of discontinuities in the processes under investigation.
Our work is organised as follows: (i) we begin by presenting our model, by
clearly justifying its used form, namely the component related to the L\'evy
noise; (ii) we prove existence and uniqueness of a global positive solution by
constructing a suitable stopping time; (iii) under some assumptions, we show
extinction of HIV/AIDS; (iv) we obtain sufficient conditions assuring
persistence of HIV/AIDS; (v) we illustrate our mathematical results through
numerical simulations.


  This study presents new estimates for fractional derivatives without singular
kernels defined by some specific functions. Based on obtained inequalities, we
give a useful method to establish the global stability of steady states for
fractional-order systems and generalize some works existing in the literature.
Finally, we apply our results to prove the global stability of a
fractional-order SEIR model with a general incidence rate.


  Covid-19 vaccines are widely available in the United States, yet our Covid-19
vaccination rates have remained far below 100%. Not only that, but CDC data
shows that even in places where vaccine acceptance was proportionally high at
the outset of the Covid-19 vaccination effort, that willingness has not
necessarily translated into high rates of vaccination over the subsequent
months. We model how such a shift could have arisen, using parameters in
agreement with data from the state of Alabama. The simulations suggest that in
Alabama, local interactions would have favored the emergence of tight consensus
around the initial majority view, which was to accept the Covid-19 vaccine. Yet
this is not what happened. We therefore add to our model the impact of
mega-influencers such as mass media, the governor of the state, etc. Our
simulations show that a single vaccine-hesitant mega-influencer, reaching a
large fraction of the population, can indeed cause the consensus to shift
radically, from acceptance to hesitancy. Surprisingly this is true even when
the mega-influencer only reaches individuals who are already somewhat inclined
to agree with them, and under the conservative assumption that individuals give
no more weight to the mega-influencer than they would give to a single one of
their friends or neighbors. Our simulations also suggest that a competing
mega-influencer with the opposite view can shift the mean population opinion
back, but cannot restore the tightness of consensus around that view. Our code
and data are distributed in the ODyN (Opinion Dynamic Networks) library
available at https://github.com/annahaensch/ODyN.


  Calcium and reactive oxygen species (ROS) interact with each other and play
an important role in cell signaling networks. Based on the existing
mathematical models, we develop an age-dependent feedback control model to
simulate the interaction. The model consists of three subsystems: cytosolic
calcium dynamics, ROS generation from the respiratory chain in mitochondria,
and mitochondrial energy metabolism. In the model, we hypothesized that ROS
induces calcium release from the yeast endoplasmic reticulum , Golgi apparatus,
and vacuoles, and that ROS damages calmodulin and calcineurin by oxidizing
them. The dependence of calcium uptake by Vcx1p on ATP is incorporated into the
model. The model can approximately reproduce the log phase calcium dynamics.
The simulated interaction between the cytosolic calcium and mitochondrial ROS
shows that an increase in calcium results in a decrease in ROS initially (in
log phase), but the increase-decrease relation is changed to an
increase-increase relation when the cell is getting old. This could accord with
the experimental observation that calcium diminishes ROS from complexes I and
III of the respiratory chain under normal conditions, but enhances ROS when the
complex formations are inhibited. The model predicts that the subsystem of the
calcium regulators Pmc1p, Pmr1p, and Vex1p is stable, controllable, and
observable. These structural properties of the dynamical system could
mathematically confirm that cells have evolved delicate feedback control
mechanisms to maintain their calcium homeostasis.


  Motivation: Predicting the drug-target interaction is crucial for drug
discovery as well as drug repurposing. Machine learning is commonly used in
drug-target affinity (DTA) problem. However, machine learning model faces the
cold-start problem where the model performance drops when predicting the
interaction of a novel drug or target. Previous works try to solve the cold
start problem by learning the drug or target representation using unsupervised
learning. While the drug or target representation can be learned in an
unsupervised manner, it still lacks the interaction information, which is
critical in drug-target interaction. Results: To incorporate the interaction
information into the drug and protein interaction, we proposed using transfer
learning from chemical-chemical interaction (CCI) and protein-protein
interaction (PPI) task to drug-target interaction task. The representation
learned by CCI and PPI tasks can be transferred smoothly to the DTA task due to
the similar nature of the tasks. The result on the drug-target affinity
datasets shows that our proposed method has advantages compared to other
pretraining methods in the DTA task.


  Given the severe impact of COVID-19 on several societal levels, it is of
crucial importance to model the impact of restriction measures on the pandemic
evolution, so that governments are able to take informed decisions. Even though
there have been countless attempts to propose diverse models since the raise of
the outbreak, the increase in data availability and start of vaccination
campaigns calls for updated models and studies. Furthermore, most of the works
are focused on a very particular place or application and we strive to attain a
more general model, resorting to data from different countries. In particular,
we compare Great Britain and Israel, two highly different scenarios in terms of
vaccination plans and social structure. We build a network-based model, complex
enough to model different scenarios of government-mandated restrictions, but
generic enough to be applied to any population. To ease the computational load
we propose a decomposition strategy for our model.


  We provide a rigorous definition of free-energy transduction and its
efficiency in arbitrary -- linear or nonlinear -- open chemical reaction
networks (CRNs) operating at steady state. Our method is based on the knowledge
of the stoichiometric matrix and of the chemostatted species (i.e. the species
maintained at constant concentration by the environment) to identify the
fundamental currents and forces contributing to the entropy production.
Transduction occurs when the current of a stoichiometrically balanced process
is driven against its spontaneous direction (set by its force) thanks to other
processes flowing along their spontaneous direction. In these regimes, open
CRNs operate as thermodynamic machines. After exemplifying these general ideas
using toy models, we analyze central energy metabolism. We relate the
fundamental currents to metabolic pathways and discuss the efficiency with
which they are able to transduce free energy.


  Despite significant progress of generative models in the natural sciences,
their controllability remains challenging. One fundamentally missing aspect of
molecular or protein generative models is an inductive bias that can reflect
continuous properties of interest. To that end, we propose the Regression
Transformer (RT), a novel method that abstracts regression as a conditional
sequence modeling problem. This introduces a new paradigm of multitask language
models which seamlessly bridge sequence regression and conditional sequence
generation.
  We thoroughly demonstrate that, despite using a nominal-scale training
objective, the RT matches or surpasses the performance of conventional
regression models in property prediction tasks of small molecules, proteins and
chemical reactions. Critically, priming the same model with continuous
properties yields a highly competitive conditional generative model that
outperforms specialized approaches in a substructure-constrained,
property-driven molecule generation benchmark. Our dichotomous approach is
facilitated by a novel, alternating training scheme that enables the model to
decorate seed sequences by desired properties, e.g., to optimize reaction
yield.
  In sum, the RT is the first report of a multitask model that concurrently
excels at predictive and generative tasks in biochemistry. This finds
particular application in property-driven, local exploration of the chemical or
protein space and could pave the road toward foundation models in material
design.
  The code to reproduce all experiments of the paper is available at:
https://github.com/IBM/regression-transformer


  To control African swine fever (ASF) efficiently, easily interpretable
metrics of the outbreak dynamics are needed to plan and adapt the required
measures. We found that the spread pattern of African Swine Fever cases in wild
boar follows the mechanics of a diffusion process, at least in the early phase,
for the cases that occurred in Germany. Following incursion into a previously
unaffected area, infection disseminates locally within a naive and abundant
wild boar population. Using real case data for Germany, we derive statistics
about the time differences and distances between consecutive case reports. With
the use of these statistics, we generate an ensemble of random walkers
(continuous time random walks, CTRW) that resemble the properties of the
observed outbreak pattern as one possible realization of all possible disease
dissemination patterns. The trained random walker ensemble yields the diffusion
constant, the affected area, and the outbreak velocity of early ASF spread in
wild boar. These quantities are easy to interpret, robust, and may be
generalized or adapted to different regions. Therefore, diffusion metrics can
be useful descriptors of early disease dynamics and help facilitate efficient
control of African Swine Fever.


  Early literature on genome rearrangement modelling views the problem of
computing evolutionary distances as an inherently combinatorial one. In
particular, attention was given to estimating distances using the minimum
number of events required to transform one genome into another. In hindsight,
this approach is analogous to early methods for inferring phylogenetic trees
from DNA sequences such as maximum parsimony -- both are motivated by the
principle that the true distance minimises evolutionary change, and both are
effective if this principle is a true reflection of reality. Recent literature
considers genome rearrangement under statistical models, continuing this
parallel with DNA-based methods; the goal here is to use model-based methods
(for example maximum likelihood techniques) to compute distance estimates that
incorporate the large number of rearrangement paths that can transform one
genome into another. Crucially, this approach requires one to decide upon a set
of feasible rearrangement events and, in this paper, we focus on characterising
well-motivated models for signed, uni-chromosomal circular genomes, where the
number of regions remains fixed. Since rearrangements are often mathematically
described using permutations, we isolate the sets of permutations representing
rearrangements that are biologically reasonable in this context, for example
inversions and translocations. We provide precise mathematical expressions for
these rearrangements, and then describe them in terms of the set of cuts made
in the genome when they are applied. We directly compare cuts to breakpoints,
and use this concept to count the distinct rearrangement actions which apply a
given number of cuts. Finally, we provide some examples of rearrangement
models, and include a discussion of some questions that arise when defining
plausible models.


  Ecological and evolutionary processes show various population dynamics
depending on internal interactions and environmental changes. While crucial in
predicting biological processes, discovering general relations for such
nonlinear dynamics has remained a challenge. Here, we derive a universal
information-theoretical constraint on a broad class of nonlinear dynamical
systems represented as population dynamics. The constraint is interpreted as a
generalization of Fisher's fundamental theorem of natural selection.
Furthermore, the constraint indicates nontrivial bounds for the speed of
critical relaxation around bifurcation points, which we argue are universally
determined only by the type of bifurcation. Our theory is verified for an
evolutionary model and an epidemiological model, which exhibit the
transcritical bifurcation, as well as for an ecological model, which undergoes
limit-cycle oscillation. This work paves a way to predict biological dynamics
in light of information theory, by providing fundamental relations in
nonequilibrium statistical mechanics of nonlinear systems.


  We study the global dynamics of a reaction-diffusion SEIR infection model
with distributed delay and nonlinear incidence rate. The well-posedness of the
proposed model is proved. By means of Lyapunov functionals, we show that the
disease free equilibrium state is globally asymptotically stable when the basic
reproduction number is less or equal than one, and that the disease endemic
equilibrium is globally asymptotically stable when the basic reproduction
number is greater than one. Numerical simulations are provided to illustrate
the obtained theoretical results.


  In this study, we explore interactions between cancer cells by using the
hawk-dove game. We analyze the heterogeneity of tumors by considering games
with populations composed of 2 or 3 types of cells. We determine what
strategies are evolutionarily stable in the 2-type and 3-type population games
and the corresponding expected payoffs. Our results show that the payoff of the
best-off cell in the 2-type population game is higher than that of the best-off
cell in the 3-type population game. Translating these mathematical findings to
the field of oncology suggests that a branching-type tumor pursues a less
aggressive course than a punctuated-type one. Some histological and genomic
data of clear cell renal cell carcinomas are consistent with these results.


  Understanding the evolutionary stability of cooperation is a central problem
in biology, sociology, and economics. There exist only a few known mechanisms
that guarantee the existence of cooperation and its robustness to cheating.
Here, we introduce a new mechanism for the emergence of cooperation in the
presence of fluctuations. We consider agents whose wealth change stochastically
in a multiplicative fashion. Each agent can share part of her wealth as public
good, which is equally distributed among all the agents. We show that, when
agents operate with long time-horizons, cooperation produce an advantage at the
individual level, as it effectively screens agents from the deleterious effect
of environmental fluctuations.


  Pre-trained protein models (PTPMs) represent a protein with one fixed
embedding and thus are not capable for diverse tasks. For example, protein
structures can shift, namely protein folding, between several conformations in
various biological processes. To enable PTPMs to produce task-aware
representations, we propose to learn interpretable, pluggable and extensible
protein prompts as a way of injecting task-related knowledge into PTPMs. In
this regard, prior PTPM optimization with the masked language modeling task can
be interpreted as learning a sequence prompt (Seq prompt) that enables PTPMs to
capture the sequential dependency between amino acids. To incorporate
conformational knowledge to PTPMs, we propose an interaction-conformation
prompt (IC prompt) that is learned through back-propagation with the
protein-protein interaction task. As an instantiation, we present a
conformation-aware pre-trained protein model that learns both sequence and
interaction-conformation prompts in a multi-task setting. We conduct
comprehensive experiments on nine protein datasets. Results confirm our
expectation that using the sequence prompt does not hurt PTPMs' performance on
sequence-related tasks while incorporating the interaction-conformation prompt
significantly improves PTPMs' performance on tasks where conformational
knowledge counts. We also show the learned prompts can be combined and extended
to deal with new complex tasks.


  Computational analysis methods including machine learning have a significant
impact in the fields of genomics and medicine. High-throughput gene expression
analysis methods such as microarray technology and RNA sequencing produce
enormous amounts of data. Traditionally, statistical methods are used for
comparative analysis of gene expression data. However, more complex analysis
for classification of sample observations, or discovery of feature genes
requires sophisticated computational approaches. In this review, we compile
various statistical and computational tools used in analysis of expression
microarray data. Even though the methods are discussed in the context of
expression microarrays, they can also be applied for the analysis of RNA
sequencing and quantitative proteomics datasets. We discuss the types of
missing values, and the methods and approaches usually employed in their
imputation. We also discuss methods of data normalization, feature selection,
and feature extraction. Lastly, methods of classification and class discovery
along with their evaluation parameters are described in detail. We believe that
this detailed review will help the users to select appropriate methods for
preprocessing and analysis of their data based on the expected outcome.


  Surface assays such as ELISA are pervasive in clinics and research and
predominantly standardized in microtiter plates (MTP). MTPs provide many
advantages but are often detrimental to surface assay efficiency due to
inherent mass transport limitations. Microscale flows can overcome these and
largely improve assay kinetics. However, the disruptive nature of microfluidics
with existing labware and protocols has narrowed its transformative potential.
We present WellProbe, a novel microfluidic concept compatible with MTPs. With
it, we show that immunoassays become more sensitive at low concentrations (up
to 9-fold signal improvement in 12-fold less time), richer in information with
3-4 different kinetic conditions, and can be used to estimate kinetic
parameters, minimize washing steps and non-specific binding, and identify
compromised results. We further multiplex single-well assays combining
WellProbes kinetic regions with tailored microarrays. Finally, we demonstrate
our system in a context of immunoglobulin subclass evaluation, increasingly
regarded as clinically relevant.


  Seagrasses are vital organisms in coastal waters, and the drastic demise of
their population in the last decades has worrying implications for marine
ecosystems. Spatial models for seagrass meadows provide a mathematical
framework to study their dynamical processes and emergent collective behavior.
These models are crucial to predict the response of seagrasses to different
global warming scenarios, analyze the resilience of existing seagrass
distributions, and optimize restoration strategies. In this article, we propose
a model that includes interactions among different species based on the clonal
growth of seagrasses. We present a theoretical analysis of the model
considering the specific case of the seagrass-seaweed interaction between {\it
Cymodocea nodosa} and {\it Caulerpa prolifera}. Our simulations successfully
reproduce field observations of shoot densities in mixed meadows in the Ebro
River Delta in the Mediterranean Sea. Besides, the proposed model allows us to
investigate the possible underlying mechanisms that mediate the interaction
among the two macrophytes.


  COVID-19 has been prevalent worldwide for about 2 years now and has brought
unprecedented challenges to our society. Before vaccines were available, the
main disease intervention strategies were non-pharmaceutical. Starting December
2020, in Ontario, Canada, vaccines were approved for administering to
vulnerable individuals and gradually expanded to all individuals above the age
of 12. As the vaccine coverage reached a satisfactory level among the eligible
population, normal social activities resumed and schools reopened starting
September 2021. However, when schools reopen for in-person learning, children
under the age of 12 are unvaccinated and are at higher risks of contracting the
virus. We propose an age-stratified model based on the age and vaccine
eligibility of the individuals. We fit our model to the data in Ontario, Canada
and obtain a good fitting result. The results show that a relaxed between-group
contact rate may trigger future epidemic waves more easily than an increased
within-group contact rate. An increasing mixed contact rate of the older group
quickly amplifies the daily incidence numbers for both groups whereas an
increasing mixed contact rate of the younger group mainly leads to future waves
in the younger group alone. The results indicate the importance of accelerating
vaccine rollout for younger individuals in mitigating disease spread.


  Phylogenetic placement refers to a family of tools and methods to analyze,
visualize, and interpret the tsunami of metagenomic sequencing data generated
by high-throughput sequencing. Compared to alternative (e. g.,
similarity-based) methods, it puts metabarcoding sequences into a phylogenetic
context using a set of known reference sequences and taking evolutionary
history into account. Thereby, one can increase the accuracy of metagenomic
surveys and eliminate the requirement for having exact or close matches with
existing sequence databases. Phylogenetic placement constitutes a valuable
analysis tool per se, but also entails a plethora of downstream tools to
interpret its results. A common use case is to analyze species communities
obtained from metagenomic sequencing, for example via taxonomic assignment,
diversity quantification, sample comparison, and identification of correlations
with environmental variables. In this review, we provide an overview over the
methods developed during the first ten years. In particular, the goals of this
review are (i) to motivate the usage of phylogenetic placement and illustrate
some of its use cases, (ii) to outline the full workflow, from raw sequences to
publishable figures, including best practices, (iii) to introduce the most
common tools and methods and their capabilities, (iv) to point out common
placement pitfalls and misconceptions,(v) to showcase typical placement-based
analyses, and how they can help to analyze, visualize, and interpret
phylogenetic placement data.


  The potential of virtual reality (VR) to contribute to drug design and
development has been recognised for many years. Hardware and software
developments now mean that this potential is beginning to be realised, and VR
methods are being actively used in this sphere. A recent advance is to use VR
not only to visualise and interact with molecular structures, but also to
interact with molecular dynamics simulations of 'on the fly' (interactive
molecular dynamics in VR, IMD-VR), which is useful not only for flexible
docking but also to examine binding processes and conformational changes.
iMD-VR has been shown to be useful for creating complexes of ligands bound to
target proteins, e.g., recently applied to peptide inhibitors of the SARS-CoV-2
main protease. In this review, we use the term 'interactive VR' to refer to
software where interactivity is an inherent part of the user VR experience
e.g., in making structural modifications or interacting with a physically
rigorous molecular dynamics (MD) simulation, as opposed to simply using VR
controllers to rotate and translate the molecule for enhanced visualisation.
Here, we describe these methods and their application to problems relevant to
drug discovery, highlighting the possibilities that they offer in this arena.
We suggest that the ease of viewing and manipulating molecular structures and
dynamics, and the ability to modify structures on the fly (e.g., adding or
deleting atoms) makes modern interactive VR a valuable tool to add to the
armoury of drug development methods.


  We analyze a within-host model of virus infection with antibody and CD8+
cytotoxic T lymphocyte (CTL) responses proposed by Schwartz et al. (2013). The
goal of this work is to gain an overview of the stability of the
biologically-relevant equilibria as a function of the model's immune response
parameters. We show that the equilibria undergo at most two forward
transcritical bifurcations. The model is also explored numerically and results
are applied to equine infectious anemia virus infection. In order to arrive at
stability of the biologically-relevant endemic equilibrium characterized by
coexistence of antibody and CTL responses, the parameters promoting CTL
responses need to be boosted over parameters promoting antibody production.
This result may seem counter-intuitive (in that a weaker antibody response is
better) but can be understood in terms of a balance between CTL and antibody
responses that is needed to permit existence of CTLs. In conclusion, an
intervention such as a vaccine that is intended to control a persistent viral
infection with both immune responses should moderate the antibody response to
allow for stimulation of the CTL response.


  Protein engineers conventionally use tools such as Directed Evolution to find
new proteins with better functionalities and traits. More recently,
computational techniques and especially machine learning approaches have been
recruited to assist Directed Evolution, showing promising results. In this
paper, we propose POET, a computational Genetic Programming tool based on
evolutionary computation methods to enhance screening and mutagenesis in
Directed Evolution and help protein engineers to find proteins that have better
functionality. As a proof-of-concept we use peptides that generate MRI contrast
detected by the Chemical Exchange Saturation Transfer contrast mechanism. The
evolutionary methods used in POET are described, and the performance of POET in
different epochs of our experiments with Chemical Exchange Saturation Transfer
contrast are studied. Our results indicate that a computational modelling tool
like POET can help to find peptides with 400% better functionality than used
before.


  We demonstrate from first principles a core fallacy employed by a coterie of
authors who claim that data from the Vaccine Adverse Reporting System (VAERS)
show that hundreds of thousands of U.S. deaths are attributable to COVID
vaccination.


  Resource are often not uniformly distributed within a population. Spatial
variations of concentration of a resource, change the fitness of competing
strategies locally. The notion of fitness varying with respect to both genotype
and environment is important in modeling cancer initiation, microbial evolution
and evolution of drug resistance. Environmental interactions can be asymmetric,
that is, they affect the fitness of one type more than the other. The question
is how local environmental variations in network population structures change
the selection dynamics in a finite population setting. We consider
one-dimensional lattice population structures with spatial fitness
distributions with a periodic pattern. Heterogeneity is determined by standard
deviation of fitnesses and period. The model covers biologically relevant
limits of two-habitat subdivided populations and randomly-distributed resources
in high- and low-periods. We numerically calculate fixation probability and
fixation times for a constant population birth-death process as fitness
heterogeneity and period vary. We identify levels of heterogeneity for which a
previously deleterious mutant, in a uniform environment, becomes beneficial. In
other regimes of the problem we observe unexpected behavior where the fixation
probability of both types are larger than their neutral value at the same time.
This coincides with an exponential increase in time to fixation as a function
of population size, which points to significant slow-down in selection process
and the potential for coexistence between types in realistic time scales. We
also discuss `fitness shift' model where the fitness function of one type is
identical to the other up to a constant spatial shift. This leads to
significant increase (or decrease) in the fixation probability of the mutant
depending the value of the shift.


  Let $W_1,\ldots,W_N$ be a sample of $\mathrm{Pareto}(\alpha)$ random
variables normalized by their sum, such that $\sum_i W_i=1$. The $W_i$ may
represent the weights of valleys in a spin glass (if $0<\alpha<1$), or the
frequency of different lineages (families) in a genealogy. This paper considers
a population in which there are $N$ individuals reproducing with
$\mathrm{Pareto}(\alpha)$ offspring-number distribution ($1<\alpha<2$). The
probability of two randomly-chosen individuals being siblings, $Y_2=\sum_i
W_i^2$, gives the sample mean of the normalized size of families, and its
reciprocal gives the effective number of families (or reproducing lineages) in
the population, $N_{\mathrm{e}}=1/Y_2$. The typical sample mean is very
different from the average over all possible samples, i.e. $Y_2$ is not a
self-averaging quantity. The typical $Y_2$ and its reciprocal do not vary with
$N$ in opposite ways. Non-self-averaging effects are crucial in understanding
genetic diversity in mass spawning species such as marine fishes.


  The global pandemic due to the outbreak of COVID-19 ravages the whole world
for more than two years in which all the countries are suffering a lot since
December 2019. In order to control this ongoing waves of epidemiological
infections, attempts have been made to understand the dynamics of this pandemic
in deterministic approach with the help of several mathematical models. In this
article characteristics of a multi-wave SIR model have been studied which
successfully explains the features of this pandemic waves in India. Stability
of this model has been studied by identifying the equilibrium points as well as
by finding the eigen values of the corresponding Jacobian matrices. Complex
eigen values are found which ultimately give rise to the oscillatory solutions
for the three categories of populations, say, susceptible, infected and
removed. In this model, a finite probability of the recovered people for
becoming susceptible again is introduced which eventually lead to the
oscillatory solution in other words. The set of differential equations has been
solved numerically in order to obtain the variation for numbers of susceptible,
infected and removed people with time. In this phenomenological study, finally
an additional modification is made in order to explain the aperiodic
oscillation which is found necessary to capture the feature of epidemiological
waves particularly in India.


  We study a prey-predator model based on the classical Lotka-Volterra system
with Leslie-Gower and Holling IV schemes and a constant-effort harvesting. Our
goal is twofold: to present the model proposed by Cheng and Zhang in 2021,
pointing out some inconsistencies; to analyse the number and type of
equilibrium points of the model. We end by proving the stability of the
meaningful equilibrium point, according to the distribution of the eigenvalues.


  The Omicron variant of severe acute respiratory syndrome coronavirus 2
(SARS-CoV-2) has rapidly replaced the Delta variant as a dominating SARS-CoV-2
variant because of natural selection, which favors the variant with higher
infectivity and stronger vaccine breakthrough ability. Omicron has three
lineages or subvariants, BA.1 (B.1.1.529.1), BA.2 (B.1.1.529.2), and BA.3
(B.1.1.529.3). Among them, BA.1 is the currently prevailing subvariant. BA.2
shares 32 mutations with BA.1 but has 28 distinct ones. BA.3 shares most of its
mutations with BA.1 and BA.2 except for one. BA.2 is found to be able to
alarmingly reinfect patients originally infected by Omicron BA.1. An important
question is whether BA.2 or BA.3 will become a new dominating "variant of
concern". Currently, no experimental data has been reported about BA.2 and
BA.3. We construct a novel algebraic topology-based deep learning model trained
with tens of thousands of mutational and deep mutational data to systematically
evaluate BA.2's and BA.3's infectivity, vaccine breakthrough capability, and
antibody resistance. Our comparative analysis of all main variants namely,
Alpha, Beta, Gamma, Delta, Lambda, Mu, BA.1, BA.2, and BA.3, unveils that BA.2
is about 1.5 and 4.2 times as contagious as BA.1 and Delta, respectively. It is
also 30% and 17-fold more capable than BA.1 and Delta, respectively, to escape
current vaccines. Therefore, we project that Omicron BA.2 is on its path to
becoming the next dominating variant. We forecast that like Omicron BA.1, BA.2
will also seriously compromise most existing mAbs, except for sotrovimab
developed by GlaxoSmithKline.


  The Amazon floodplains represent important surfaces of highly valuable
ecosystems, yet they remain neglected from protected areas. While the
efficiency of the protected area network of the Amazon basin may be jeopardised
by climate change, floodplains are exposed to important consequences of climate
change but are omitted from species distribution models and protection gap
analyses. We modelled the current and future (2070) distribution of the giant
bony-tongue fish Arapaima sp. accounting for climate and habitat requirements,
with consideration of dam presence (already existing and planned constructions)
and hydroperiod (high- and low-water stages). We further quantified the amount
of suitable environment which falls inside and outside the current network of
protected areas to identify spatial conservation gaps. We predict climate
change to cause the decline of environmental suitability by 16.6% during the
high-water stage, and by 19.4% during the low-water stage. We found that about
70% of the suitable environments of Arapaima sp. remain currently unprotected,
which is likely to increase by 5% with future climate change effects. Both
current and projected dam constructions may hamper population flows between the
central and the Bolivian and Peruvian parts of the basin. We highlight
protection gaps mostly in the southwestern part of the basin and recommend the
extension of the current network of protected areas in the floodplains of the
upper Ucayali, Juru\`a and Purus Rivers and their tributaries. This study
showed the importance of taking into account hydroperiods and dispersal
barriers in forecasting the distribution of freshwater fish species, and
stresses the urgent need to integrate floodplains to the protected area
networks.


  Drug development is time-consuming and expensive. Repurposing existing drugs
for new therapies is an attractive solution that accelerates drug development
at reduced experimental costs, specifically for Coronavirus Disease 2019
(COVID-19), an infectious disease caused by severe acute respiratory syndrome
coronavirus 2 (SARS-CoV-2). However, comprehensively obtaining and productively
integrating available knowledge and big biomedical data to effectively advance
deep learning models is still challenging for drug repurposing in other complex
diseases. In this review, we introduce guidelines on how to utilize deep
learning methodologies and tools for drug repurposing. We first summarized the
commonly used bioinformatics and pharmacogenomics databases for drug
repurposing. Next, we discuss recently developed sequence-based and graph-based
representation approaches as well as state-of-the-art deep learning-based
methods. Finally, we present applications of drug repurposing to fight the
COVID-19 pandemic, and outline its future challenges.


  Predicting how a drug-like molecule binds to a specific protein target is a
core problem in drug discovery. An extremely fast computational binding method
would enable key applications such as fast virtual screening or drug
engineering. Existing methods are computationally expensive as they rely on
heavy candidate sampling coupled with scoring, ranking, and fine-tuning steps.
We challenge this paradigm with EquiBind, an SE(3)-equivariant geometric deep
learning model performing direct-shot prediction of both i) the receptor
binding location (blind docking) and ii) the ligand's bound pose and
orientation. EquiBind achieves significant speed-ups and better quality
compared to traditional and recent baselines. Further, we show extra
improvements when coupling it with existing fine-tuning techniques at the cost
of increased running time. Finally, we propose a novel and fast fine-tuning
model that adjusts torsion angles of a ligand's rotatable bonds based on
closed-form global minima of the von Mises angular distance to a given input
atomic point cloud, avoiding previous expensive differential evolution
strategies for energy minimization.


  During the early stage of the COVID-19 pandemic, many countries implemented
non-pharmaceutical interventions (NPIs) to control the transmission of
SARS-CoV-2, the causative pathogen of COVID-19. Among those NPIs, quarantine
measures were widely adopted and enforced through stay-at-home and
shelter-in-place orders. Understanding the effectiveness of quarantine measures
can inform decision-making and control planning during the ongoing COVID-19
pandemic and for future disease outbreaks. In this study, we use mathematical
models to evaluate the impact of quarantine measures on COVID-19 spread in four
cities that experienced large-scale outbreaks in the spring of 2020: Wuhan, New
York, Milan, and London. We develop a susceptible-exposed-infected-removed
(SEIR)-type model with a component of quarantine and couple this disease
transmission model with a data assimilation method. By calibrating the model to
case data, we estimate key epidemiological parameters before lockdown in each
city. We further examine the impact of quarantine rates on COVID-19 spread
after lockdown using model simulations. Results indicate that quarantine of
susceptible and exposed individuals and undetected infections is necessary to
contain the outbreak; however, the quarantine rates for these populations can
be reduced through faster isolation of confirmed cases. We generate
counterfactual simulations to estimate effectiveness of quarantine measures.
Without quarantine measures, the cumulative confirmed cases could be 73, 22, 43
and 93 times higher than reported numbers within 40 days after lockdown in
Wuhan, New York, Milan, and London. Our findings underscore the essential role
of quarantine during the early phase of the pandemic.


  Apparent competition is an indirect interaction between species that share
natural resources without any mutual aggression but negatively affect each
other if there is a common enemy. The negative results of the apparent
competition are reflected in the species spatial segregation, which impacts the
dynamics of their populations. Performing a series of stochastic simulations,
we study a model where organisms of two prey species do not compete for space
but share a common predator. Our outcomes elucidate the central role played by
the predator in the pattern formation and coarsening dynamics in apparent
competition models. Investigating the effects of predator mortality on the
persistence of the species, we find a crossover between a curvature driven
scaling regime and a coexistence scenario. For low predator mortality, spatial
domains mainly inhabited by one type of prey arise, surrounded by interfaces
that mostly contain predators. We demonstrate that the dynamics of the
interface network are curvature driven whose coarsening follows a scaling law
common to other nonlinear systems. The effects of the apparent competition
decrease for high predator mortality, allowing organisms of two prey species to
share a more significant fraction of lattice. Finally, our results reveal that
predation capacity in single-prey domains influences the scaling power law that
characterises the coarsening dynamics. Our findings may be helpful to
biologists to understand the pattern formation and dynamics of biodiversity in
systems with apparent competition.


  We discuss the criticality in the stochastic SIR model for infectious
diseases. We adopt the path-integral formalism for the propagation of
infections among susceptible, infectious, and removed individuals, and perform
the perturbative and nonperturbative analyses to evaluate the critical value of
the basic reproduction number ${\cal R}$. In the perturbation theory, we
calculate the mean values and the variances of the number of infectious
individuals near the initial time, and find that the critical value ${\cal
R}_{\text{c}}=1/3$-$2/3$ should be adopted in order to suppress the stochastic
spread of infections sufficiently. In the nonperturbative approach, we derive
the effective potential by integrating out the stochastic fluctuations, and
obtain the effective Euler-Lagrange equations for the time-evolution of the
numbers of susceptible, infectious, and removed individuals. From the
asymptotic behaviors for a long time, we find that the critical value ${\cal
R}_{\text{c}}=2/3$ should be adopted for the sufficient convergence of
infections. We also find that the endemic state can be generated dynamically by
the stochastic fluctuation which is absent in the conventional SIR model. Those
analyses show that the critical value of the basic reproduction number should
be less than one, against the usually known critical value ${\cal
R}_{\text{c}}=1$, when the stochastic fluctuations are taken into account in
the SIR model.


  Monitoring the evolution of the Covid19 pandemic constitutes a critical step
in sanitary policy design. Yet, the assessment of the pandemic intensity within
the pandemic period remains a challenging task because of the limited quality
of data made available by public health authorities (missing data, outliers and
pseudoseasonalities, notably), that calls for cumbersome and ad-hoc
preprocessing (denoising) prior to estimation. Recently, the estimation of the
reproduction number, a measure of the pandemic intensity, was formulated as an
inverse problem, combining data-model fidelity and space-time regularity
constraints, solved by nonsmooth convex proximal minimizations. Though
promising, that formulation lacks robustness against the limited quality of the
Covid19 data and confidence assessment. The present work aims to address both
limitations: First, it discusses solutions to produce a robust assessment of
the pandemic intensity by accounting for the low quality of the data directly
within the inverse problem formulation. Second, exploiting a Bayesian
interpretation of the inverse problem formulation, it devises a Monte Carlo
sampling strategy, tailored to a nonsmooth log-concave a posteriori
distribution, to produce relevant credibility intervalbased estimates for the
Covid19 reproduction number. Clinical relevance Applied to daily counts of new
infections made publicly available by the Health Authorities for around 200
countries, the proposed procedures permit robust assessments of the time
evolution of the Covid19 pandemic intensity, updated automatically and on a
daily basis.


  Since the last century, deterministic compartmental models have emerged as
powerful tools to predict and control epidemic outbreaks, in many cases helping
to mitigate their impacts. A key quantity for these models is the so-called
Basic Reproduction Number, that measures the number of secondary infections
produced by an initial infected individual in a fully susceptible population.
Some methods have been developed to allow the direct computation of this
quantity provided that some conditions are fulfilled, such that the model has a
pre-pandemic disease-free equilibrium state. This condition is only fulfilled
when the populations are stationary. In the case of vector-borne diseases, this
implies that the vector birth and death rates need to be balanced, what is not
fulfilled in many realistic cases in which the vector population grow or
decrease. Here we develop a vector-borne epidemic model with growing and
decaying vector populations and study the conditions under which the standard
methods to compute $R_0$ work and discuss an alternative when they fail. We
also show that growing vector populations produce a delay in the epidemic
dynamics when compared to the case of the stationary vector population.
Finally, we discuss the conditions under which the model can be reduced to the
SIR model with fewer compartments and parameters, which helps in solving the
problem of parameter unidentifiability of many vector-borne epidemic models.


  It is well known that Drug Design is often a costly process both in terms of
time and economic effort. While good Quantitative Structure-Activity
Relationship models (QSAR) can help predicting molecular properties without the
need to synthesize them, it is still required to come up with new molecules to
be tested. This is mostly done in lack of tools to determine which
modifications are more promising or which aspects of a molecule are more
influential for the final activity/property. Here we present an automatic
process which involves Graph Convolutional Network models and input-attribution
methods to generate new molecules. We also explore the problems of
over-optimization and applicability, recognizing them as two important aspects
in the practical use of such automatic tools.


  Since the introduction of artificial intelligence in medicinal chemistry, the
necessity has emerged to analyse how molecular property variation is modulated
by either single atoms or chemical groups. In this paper, we propose to train
graph-to-graph neural network using semi-supervised learning for attributing
structure-property relationships. As initial case studies we apply the method
to solubility and molecular acidity while checking its consistency in
comparison with known experimental chemical data. As final goal, our approach
could represent a valuable tool to deal with problems such as activity cliffs,
lead optimization and de-novo drug design.


  Animals typically forage in groups. Social foraging can help animals avoid
predation and decrease their uncertainty about the richness of food resources.
Despite this, theoretical mechanistic models of patch foraging have
overwhelmingly focused on the behavior of single foragers. In this study, we
develop a mechanistic model describing the behavior of individuals foraging
together and departing food patches following an evidence accumulation process.
Each individual's belief about patch quality is represented by a stochastically
accumulating variable coupled to others' belief, representing the transfer of
information. We consider a cohesive group, and model information sharing as
either intermittent pulsatile coupling (communicate decision to leave) or
continuous diffusive coupling (communicate online belief). Foraging efficiency
under pulsatile coupling has a stronger dependence on the coupling strength
parameter compared to diffusive. Despite employing minimal information
transfer, pulsatile coupling can still provide similar or higher foraging
efficiency compared to diffusive coupling. Conversely, diffusive coupling is
more robust to parameter detuning and performs better when individuals have
heterogeneous departure criteria and social information weighting. Efficiency
is measured by a reward rate function that balances the amount of energy
accumulated against the time spent in a patch, computed by solving an ordered
first passage time problem for the patch departures of each individual. Using
synthetic data we show that we can distinguish between the two modes of
communication and identify the model parameters. Our model establishes a social
patch foraging framework to parse and identify deliberative decision
strategies, to distinguish different forms of social communication, and to
allow model fitting to real world animal behavior data.


  In this paper, we discuss the embedding problem for centrosymmetric matrices,
which are higher order generalizations of the matrices occurring in Strand
Symmetric Models. These models capture the substitution symmetries arising from
the double helix structure of the DNA. Deciding whether a transition matrix is
embeddable or not enables us to know if the observed substitution probabilities
are consistent with a homogeneous continuous time substitution model, such as
the Kimura models, the Jukes-Cantor model or the general time-reversible model.
On the other hand, the generalization to higher order matrices is motivated by
the setting of synthetic biology, which works with different sizes of genetic
alphabets.


  In this article, we develop a mathematical model for the rotary bacterial
flagellar motor (BFM) based on the recently discovered structure of the stator
complex (MotA$_5$MotB$_2$). The structure suggested that the stator also
rotates. The BFM is modeled as two rotating nano-rings that interact with each
other. Specifically, translocation of protons through the stator complex drives
rotation of the MotA pentamer ring, which in turn drives rotation of the FliG
ring in the rotor via interactions between the MotA ring of the stator and the
FliG ring of the rotor. Preliminary results from the structure-informed model
are consistent with the observed torque-speed relation. More importantly, the
model predicts distinctive rotor and stator dynamics and their load dependence,
which may be tested by future experiments. Possible approaches to verify and
improve the model to further understanding of the molecular mechanism for
torque generation in BFM are also discussed.


  Behavioral experiments on the Ultimatum Game have shown that we human beings
have remarkable preference in fair play, contradicting the predictions by the
game theory. Most of the existing models seeking for explanations, however,
strictly follow the assumption of \emph{Homo economicus} in orthodox Economics
that people are self-interested and fully rational to maximize their earnings.
Here we relax this assumption by allowing that people probabilistically choose
to be "good Samaritans", acting as fair players from time to time. For
well-mixed and homogeneously structured populations, we numerically show that
as this probability increases the level of fairness undergoes from the low
scenario abruptly to the full fairness state, where occasional fair behaviors
($\sim5\%$) are sufficient to drive the whole population to behave in the
half-half split manner. We also develop a mean-field theory, which correctly
reproduces the first-order phase transition and points out that the bistability
is an intrinsic property of this game and small fair acts lead to dramatical
change due to its bifurcation structure. Heterogeneously structured
populations, however, display continuous fairness transition; surprisingly,
very few hub nodes acting as fair players are able to entrain the whole
population to the full fairness state. Our results thus reveal the unexpected
strength of "good Samaritans", which may constitute a new explanation for the
emergence of fairness in our society.


  Seagrass meadows play an essential role in blue carbon and aquatic ecosystem
services. However, methods for the flourishing of seagrass are still being
explored. Here, data from 49 public coastal surveys on the distribution of
seagrass and seaweed around the onshore aquaculture facilities are revalidated,
and an exceptional area where the seagrass Zostera marina thrives was found.
The bacterial population there showed an apparent decrease in the pathogen
candidates belonging to the order Flavobacteriales. Moreover, structure
equation modeling and a linear non-Gaussian acyclic model based on the machine
learning data estimated an optimal symbiotic bacterial group candidate for
seagrass growth as follows: the Lachnospiraceae and Ruminococcaceae families as
gut-inhabitant bacteria, Rhodobacteraceae as photosynthetic bacteria, and
Desulfobulbaceae as cable bacteria modulating oxygen or nitrate reduction and
oxidation of sulfide. These observations confer a novel perspective on the
seagrass symbiotic bacterial structures critical for blue carbon and
low-pathogenic marine ecosystems in aquaculture.


  This paper studies the scaling properties of recruitment fluctuations in
randomly varying environments, for abundant marine species with extreme
reproductive behavior. Fisheries stock-recruitment data from the North Atlantic
display fluctuation scaling, a proportionality between the standard deviation
and the average recruitment among stocks. The proportionality covers over five
orders of magnitude in the range studied. A linear-scaling behavior can be a
sign of a universal distribution of the normalized data across stocks. In light
of this conjecture, it is demonstrated that the L\'evy-stable model offers a
better effective description of the recruitment distribution than the
log-normal model. Care is devoted to the problem of random sums of random
variables. Recruitment is calculated by summing random offspring numbers with
infinite variance, where the number of summands (i.e. spawning population size)
is also a random process with infinite variance.


  Almost every country in the world is battling to limit the spread of
COVID-19. As the world strives to get an effective medication to control the
disease, appropriate intervention measures, for now, remains one of the
effective methods to reduce the spread of the disease. Optimal control
strategies have proven to be an effective method in curtailing the spread of
infectious diseases. In this paper, a model has been formulated to study
transmission dynamics of the disease. Basic properties of the model such as the
basic reproduction number, equilibrium points and stability of the equilibrium
points have been determined. Sensitivity analysis was carried on to determine
the impact of the model parameters on the basic reproduction number. We also
introduced a compartment for the deceased and examined the behaviour of
COVID-19 related deaths. The numerical simulation prediction is consistent with
real data from Ghana for the period March 2020 to March 2021. The simulation
revealed the disease had less impact on the population during the first seven
months of the outbreak. To help contain the spread of the disease, time
dependent optimal controls were incorporated into the model and Pontryagin
maximum principle was used to characterize vital conditions of the optimal
control model. Numerical simulations of the optimal control model showed that,
combination of optimal preventive strategies such as nose mask and vaccination
are effective to significantly decrease the number of COVID-19 cases in
different compartments of the model. Vaccination decreases the susceptibility
to the disease whereas mask usage preserved the susceptible population from
extinction.


  Existing studies comparing individual-based models of growing cell
populations and their continuum counterparts have mainly focused on homogeneous
populations, in which all cells have the same phenotypic characteristics.
However, significant intercellular phenotypic variability is commonly observed
in cellular systems. Therefore, we develop here an individual-based model for
the growth of phenotypically heterogeneous cell populations. In this model, the
phenotypic state of each cell is described by a structuring variable that
captures intercellular variability in cell proliferation and migration rates.
The model tracks the spatial evolutionary dynamics of single cells, which
undergo pressure-dependent proliferation, heritable phenotypic changes and
directional movement in response to pressure differentials. We formally show
that the continuum limit of this model comprises a non-local partial
differential equation for the cell population density, which generalises
earlier models of growing cell populations. Results of the individual-based
model illustrate how proliferation-migration tradeoffs shaping the evolution of
single cells can lead to the formation of travelling waves at the population
level where highly-mobile cells locally dominate at the invasive front, while
more-proliferative cells are found at the rear. We demonstrate that there is an
excellent quantitative agreement between these results and the results of
numerical simulations and formal travelling-wave analysis of the continuum
model, when sufficiently large cell numbers are considered. We provide
numerical evidence of scenarios in which the predictions of the two models may
differ due to demographic stochasticity, which cannot be captured by the
continuum model. This indicates the importance of integrating individual-based
and continuum approaches when modelling the growth of phenotypically
heterogeneous cell populations.


  We present our blind predictions for the Statistical Assessment of the
Modeling of Proteins and Ligands (SAMPL), 9th challenge, focusing on binding of
WP6 (carboxy-pillar[6]arene) with ammonium/diammonium cationic guests.
Host-guest binding free energies have been calculated using the recently
developed virtual double system single box approach, based on the enhanced
sampling of the bound and unbound end-states followed by fast switching
nonequilibrium alchemical simulations [M Macchiagodena, M Pagliai, M
Karrenbrock, G Guarnieri, F
  Iannone, Procacci, J Chem Theory Comput, 16, 7260, 2020]. As far as Pearson
and Kendall coefficients are concerned, performances were acceptable and, in
general, better than those we submitted for calixarenes, cucurbituril-like open
cavitand, and beta-cyclodextrines in previous SAMPL host-guest challenges,
confirming the reliability of nonequilibrium approaches for absolute binding
free energy calculations. In comparison with previous submissions, we found,
nonetheless, a rather large mean signed error that we attribute to the way the
finite charge correction was addressed through the assumption of a neutralizing
background plasma.


  It seems paradoxical to have observed the absence of reduced effective
population sizes $N_{\mathrm{e}}$ under marine hatchery practices. This paper
studies the Ryman-Laikre, or two-demographic-component, model of the hatchery
impact related to inbreeding in a population with power-law family-size
distribution, where hatchery inputs are represented by a Dirac delta function.
By examining the asymptotic (i.e. large-population limit) behavior of the
normalized sizes (or weights) of families of the mixture population, I derive
the distribution properties of the average weight of families (i.e. the sum of
the squared weights, $Y$) over the population existing at any given time. The
reciprocal of the average weight $Y$ gives the effective number of families (or
reproducing lineages) in the population, $N_{\mathrm{e}}=1/Y$. When the
specific production in the hatchery (i.e. the number of offspring per
broodstock) is low, the most probable value of $N_{\mathrm{e}}$ is close to the
lower bound of the $N_{\mathrm{e}}$-distribution. When the specific
hatchery-production is increased to a critical value with fixed mixing
proportion of hatchery fish, a discontinuous transition takes place, so that
the most probable $N_{\mathrm{e}}$ jumps to the upper extreme of the
distribution. This hatchery-induced transition is attributed to the breaking of
reciprocal symmetry, i.e. the fact that the typical value of $Y$ and its
reciprocal (the typical $N_{\mathrm{e}}$) do not vary with the population size
in opposite ways. At a high specific hatchery-production, the symmetry breaking
disappears.


  CRISPR-Cas systems are an adaptive immunity that protects prokaryotes against
foreign genetic elements. Genetic templates acquired during past infection
events enable DNA-interacting enzymes to recognize foreign DNA for destruction.
Due to the programmability and specificity of these genetic templates,
CRISPR-Cas systems are potential alternative antibiotics that can be engineered
to self-target antimicrobial resistance genes on the chromosome or plasmid.
However, several fundamental questions remain to repurpose these tools against
drug-resistant bacteria. For endogenous CRISPR-Cas self-targeting,
antimicrobial resistance genes and functional CRISPR-Cas systems have to
co-occur in the target cell. Furthermore, these tools have to outplay DNA
repair pathways that respond to the nuclease activities of Cas proteins, even
for exogenous CRISPR-Cas delivery. Here, we conduct a comprehensive survey of
CRISPR-Cas genomes. First, we address the co-occurrence of CRISPR-Cas systems
and antimicrobial resistance genes in the CRISPR-Cas genomes. We show that the
average number of these genes varies greatly by the CRISPR-Cas type, and some
CRISPR-Cas types (IE and IIIA) have over 20 genes per genome. Next, we
investigate the DNA repair pathways of these CRISPR-Cas genomes, revealing that
the diversity and frequency of these pathways differ by the CRISPR-Cas type.
The interplay between CRISPR-Cas systems and DNA repair pathways is essential
for the acquisition of new spacers in CRISPR arrays. We conduct simulation
studies to demonstrate that the efficiency of these DNA repair pathways may be
inferred from the time-series patterns in the RNA structure of CRISPR repeats.
This bioinformatic survey of CRISPR-Cas genomes elucidates the necessity to
consider multifaceted interactions between different genes and systems to
design effective CRISPR-based antimicrobials.


  We study how expanding the scope for selecting game and learning (adaptation)
partners affects the evolution of cooperation in public goods games on
scale-free networks. We show the following three results. (i) Expanding the
scope for selecting game partners suppresses cooperation. (ii) Expanding the
scope for selecting learning partners promotes cooperation when cooperation
evolution is difficult. (iii) When cooperation is more likely to evolve,
slightly expanding the scope for selecting learning partners causes a
significant drop in the cooperation rate, but expanding the scope further
causes the cooperation rate to recover. Although (i) is explained by the
hub-centered mechanism, the well-known dynamic that promotes cooperation on
scale-free networks, (ii) and (iii) are caused by a completely different
mechanism that has heretofore been rarely mentioned.


  The COVID-19 pandemic has had high mortality rates in the elderly and frail
worldwide, particularly in care homes. This is driven by the difficulty of
isolating care homes from the wider community, the large population sizes
within care facilities (relative to typical households), and the age/frailty of
the residents. To quantify the mortality risk posed by disease, the case
fatality risk (CFR) is an important tool. This quantifies the proportion of
cases that result in death. Throughout the pandemic, CFR amongst care home
residents in England has been monitored closely. To estimate CFR, we apply both
novel and existing methods to data on deaths in care homes, collected by Public
Health England and the Care Quality Commission. We compare these different
methods, evaluating their relative strengths and weaknesses. Using these
methods, we estimate temporal trends in the instantaneous CFR (at both daily
and weekly resolutions) and the overall CFR across the whole of England, and
dis-aggregated at regional level. We also investigate how the CFR varies based
on age and on the type of care required, dis-aggregating by whether care homes
include nursing staff and by age of residents. This work has contributed to the
summary of measures used for monitoring the UK epidemic.


  We model hypothetical bio-dispersal within a single Galactic region using the
stochastic infection dynamics process, which is inspired by these local
properties of life dispersal on Earth. We split the population of stellar
systems into different categories regarding habitability and evolved them
through time using probabilistic cellular automata rules analogous to the
model. As a dynamic effect, we include the existence of natural dispersal
vectors (e.g., dust, asteroids) in a way that avoids assumptions about their
agency. By assuming that dispersal vectors have a finite velocity and range,
the model includes the parameter of 'optical depth of life spreading'. The
effect of the oscillatory infection rate on the long-term behavior of the
dispersal flux, which adds a diffusive component to its progression, is also
taken into account. We found that phase space is separated into subregions of
long-lasting transmission, rapidly terminated transmission, and a transition
region between the two. We observed that depending on the amplitude of the
oscillatory life spreading rate, life-transmission in the Galactic patch might
take on different geometrical shapes. Even if some host systems are
uninhabited, life transmission has a certain threshold, allowing a patch to be
saturated with viable material over a long period. Although stochastic
fluctuations in the local density of habitable systems allow for clusters that
can continuously infect one another, the spatial pattern disappears when life
transmission is below the observed threshold, so that transmission process is
not permanent in time. Both findings suggest that a habitable planet in a
densely populated region may remain uninfected.


  The sequential probability ratio test (SPRT) from statistics is known to have
the least mean decision time compared to other sequential or fixed-time tests
for given error rates. In some circumstances, cells need to make decisions
accurately and quickly, therefore it has been suggested the SPRT may be used to
understand the speed-accuracy tradeoff in cellular decision making. It is
generally thought that in order for cells to make use of the SPRT, it is
necessary to find biochemical circuits that can compute the log-likelihood
ratio needed for the SPRT. However, this paper takes a different approach. We
recognise that the high-level behaviour of the SPRT is defined by its positive
detection or hit rate, and the computation of the log-likelihood ratio is just
one way to realise this behaviour. In this paper, we will present a method
which uses a transcription-based detector to emulate the hit rate of the SPRT
without computing the exact log-likelihood ratio. We consider the problem of
using a promoter with multiple binding sites to accurately and quickly detect
whether the concentration of a transcription factor is above a target level. We
show that it is possible to find binding and unbinding rates of the
transcription factor to the promoter's binding sites so that the probability
that the amount of mRNA produced will be higher than a threshold is
approximately equal to the hit rate of the SPRT detector. Moreover, we show
that the average time that this transcription-based detector needs to make a
positive detection is less than or equal to that of the SPRT for a wide range
of concentrations. We remark that the last statement does not contradict Wald's
optimality result because our transcription-based detector uses an open-ended
test.


  The effectiveness of the first and second dose vaccinations are different for
COVID-19; therefore, a
susceptible-infected-recovered-vaccination1-vaccination2-death (SIRVVD) model
that can represent the states of the first and second vaccination doses has
been proposed. By the previous study, we can carry out simulating the spread of
infectious disease considering the effects of the first and second doses of the
vaccination based on the SIRVVD model. However, theoretical analysis of the
SIRVVD Model is insufficient. Therefore, we obtained an analytical expression
of the infectious number, by treating the numbers of susceptible persons and
vaccinated persons as parameters. We used the solution to determine the target
rate of the vaccination for decreasing the infection numbers of the COVID-19
Delta variant (B.1.617) in Japan. Further, we investigated the target
vaccination rates for cases with strong or weak variants by comparison with the
COVID-19 Delta variant (B.1.617). This study contributes to the mathematical
development of the SIRVVD model and provides insight into the target rate of
the vaccination to decrease the number of infections.


  The GPRO suite is an in-progress bioinformatic project for -omic data
analyses. As part of the continued growth of this project, we introduce a
client side & server side solution for comparative transcriptomics and analysis
of variants. The client side consists of two Java applications called "RNASeq"
and "VariantSeq" to manage workflows for RNA-seq and Variant-seq analysis,
respectively, based on the most common command line interface tools for each
topic. Both applications are coupled with a Linux server infrastructure (named
GPRO Server Side) that hosts all dependencies of each application (scripts,
databases, and command line interface tools). Implementation of the server side
requires a Linux operating system, PHP, SQL, Python, bash scripting, and
third-party software. The GPRO Server Side can be deployed via a Docker
container that can be installed in the user's PC using any operating system or
on remote servers as a cloud solution. The two applications are available as
desktop and cloud applications and provide two execution modes: a Step-by-Step
mode enables each step of a workflow to be executed independently and a
Pipeline mode allows all steps to be run sequentially. The two applications
also feature an experimental support system called GENIE that consists of a
virtual chatbot/assistant and a pipeline jobs panel coupled with an expert
system. The chatbot can troubleshoot issues with the usage of each tool, the
pipeline job panel provides information about the status of each task executed
in the GPRO Server Side, and the expert provides the user with a potential
recommendation to identify or fix failed analyses. The two applications and the
GPRO Server Side combine the user-friendliness and security of client software
with the efficiency of front-end & back-end solutions to manage command line
interface software for RNA-seq and variant-seq analysis via interface
environments.


  I derive uncertainty relations that bound the rate of evolutionary processes
driven by natural selection, mutations, or by genetic drift. These rate limits
imply that the variability -- or statistical uncertainty -- in a population
allows for faster evolutionary rates. In particular, the variability of a given
quantitative trait and of the fitness function of a population are singled out
as necessary for fast changes in the trait of interest. At the same time, lower
variability makes a trait less susceptible to changes due to stochastic genetic
drift. These results generalize Fisher's fundamental theorem of natural
selection to dynamics that allow for mutations and genetic drift in terms of
uncertainty relations that constrain the evolutionary rates of arbitrary
traits.


  We have studied the effect of high hydrostatic pressure and temperature on
the steady state fluorescence anisotropy of Green Fluorescent Protein (GFP). We
find that the fluorescence anisotropy of GFP at a constant temperature
decreases with increasing pressure. At atmospheric pressure, anisotropy
decreases with increasing temperature but exhibits a maximum with temperature
for pressure larger than 20 MPa. The temperature corresponding to the maximum
of anisotropy increases with increasing pressure. By taking into account of the
rotational correlation time changes of GFP with the pressure-temperature
dependent viscosity of the solvent, we argue that viscosity increase with
pressure is not a major contributing factor to the decrease in anisotropy with
pressure. Our results suggest that the decrease of fluorescence anisotropy with
pressure may result from changes in H-bonding environment around the
chromophore.


  Recruitment is calculated by summing random offspring-numbers entering the
population, where the number of summands (i.e. spawning population size) is
also a random process. A priori, it is not clear that individual reproductive
variability would have a significant impact on aggregate measures for
monitoring populations. Usually these variations are averaged out in a large
population, and the aggregate output is merely influenced by population-wide
environmental disturbances such as climate and fisheries. However, such
arguments break down if the distribution of the individual offspring numbers is
heavy-tailed. In a world with power-law offspring-number distribution with
exponent $1<\alpha<2$, the recruitment distribution has a putative power-law
regime in the tail with the same $\alpha$. The question is to what extent
individual reproductive variability can have a noticeable impact on the
recruitment under environmentally driven population fluctuations. This question
is answered by considering the L\'evy-stable fluctuations as embedded in a
randomly varying environment. I report fluctuation scaling and asymmetric
fluctuations in recruitment of commercially exploited fish stocks throughout
the North Atlantic. The linear scaling of recruitment standard deviation with
recruitment level implies that the individual reproductive variability is
dominated by population fluctuations. The totally asymmetric (skewed to the
right) character is a sign of idiosyncratic variation in reproductive success.


  Pretrained Language Models (PLMs) such as BERT have revolutionized the
landscape of Natural Language Processing (NLP). Inspired by their
proliferation, tremendous efforts have been devoted to Pretrained Graph Models
(PGMs). Owing to the powerful model architectures of PGMs, abundant knowledge
from massive labeled and unlabeled graph data can be captured. The knowledge
implicitly encoded in model parameters can benefit various downstream tasks and
help to alleviate several fundamental issues of learning on graphs. In this
paper, we provide the first comprehensive survey for PGMs. We firstly present
the limitations of graph representation learning and thus introduce the
motivation for graph pre-training. Then, we systematically categorize existing
PGMs based on a taxonomy from four different perspectives. Next, we present the
applications of PGMs in social recommendation and drug discovery. Finally, we
outline several promising research directions that can serve as a guideline for
future research.


  Star-shaped branching patterns of genealogies are common in marine species.
High-fecundity marine populations are characterized by low ratios of effective
to actual population size, which reflect high variance in reproductive success
among parents in mass spawns. When extreme reproduction events occur, offspring
from very few parents dominate the population (whereby multiple mergers, or
subsets of lineages with star-like trees, arise) and thus, the loss of genetic
diversity is significant. Under high reproductive-variance conditions (assuming
that reproduction occurs by sampling from the Pareto distribution), this paper
explores the distribution of heterozygosity across generations. The result
shows that zero heterozygosity is not achieved, implying that the populations
may decline without evident loss of genetic variation. It is also found that
there are singularities in the heterozygosity distributions. However, in the
case of high reproductive variance, the locations of the singular points subtly
deviate from those of the case where reproduction occurs by Wright-Fisher
sampling.


  Since steep declines in a population's size also typically alter its
composition, population bottlenecks are considered highly important for
evolution. However, despite such significance, the mechanisms governing the
impact of a given population bottleneck remain poorly understood. In this
context, we show that long-term post-bottleneck outcomes can depend crucially
on the rate at which a diminished population grows whilst recovering. That is,
two otherwise identical populations, each having undergone the same bottleneck,
can fixate on dramatically different demographics due to different rates of
post-bottleneck growth. This behaviour is moreover shown to change
non-trivially with different levels of mutation. The underlying mechanism can
be traced to intrinsic fluctuations whose standard deviation decreases in time
with the inverse square root of the growing population size. Crucially, such
decreasing fluctuations couple to the underlying dynamics and result in
distinct regimes of (non-equilibrium) demographic behaviour, delimited by
abrupt transitions at critical population sizes. Seen through this lens, growth
and mutation alter evolutionary outcomes by changing the duration and character
of dynamical phases; a feature that we speculate may be of generic importance
across many systems.


  Understanding the structure of supramolecular complexes provides insight into
their functional capabilities and how they can be modulated in the context of
disease. Super-resolution microscopy (SRM) excels in performing this task by
resolving ultrastructural details at the nanoscale with molecular specificity.
However, technical limitations, such as underlabelling, preclude its ability to
provide complete structures. Single-particle analysis (SPA) overcomes this
limitation by combining information from multiple images of identical
structures and producing an averaged model, effectively enhancing the
resolution and coverage of image reconstructions. This review highlights
important studies using SRM-SPA, demonstrating how it broadens our knowledge by
elucidating features of key biological structures with unprecedented detail.


  The origin and early evolution of eukaryotes are one of the major transitions
in the evolution of life on earth. One of its most interesting aspects is the
emergence of cellular organelles, their dynamics, their functions, and their
divergence. Cell compartmentalization and architecture in prokaryotes is a less
understood complex property. In eukaryotes it is related to cell size, specific
genomic architecture, evolution of cell cycles, biogenesis of membranes and
endosymbiotic processes. Explaining cell evolution through form and function
demands an interdisciplinary approach focused on microbial diversity,
phylogenetic and functional cell biology. Two centuries of views on eukaryotic
origin have completed the disciplinary tools necessarily to answer these
questions. We have moved from Haeckel SCALA NATURAE to the un-rooted tree of
life. However, the major relations among cell domains are still elusive and
keep the nature of eukaryotic ancestor enigmatic. Here we present a review on
state of art views of eukaryogenesis; the background and perspectives of
different disciplines involved in this topic


  The diversity and quality of natural systems have been a puzzle and
inspiration for communities studying artificial life. It is now widely admitted
that the adaptation mechanisms enabling these properties are largely influenced
by the environments they inhabit. Organisms facing environmental variability
have two alternative adaptation mechanisms operating at different timescales:
\textit{plasticity}, the ability of a phenotype to survive in diverse
environments and \textit{evolvability}, the ability to adapt through mutations.
Although vital under environmental variability, both mechanisms are associated
with fitness costs hypothesized to render them unnecessary in stable
environments. In this work, we study the interplay between environmental
dynamics and adaptation in a minimal model of the evolution of plasticity and
evolvability. We experiment with different types of environments characterized
by the presence of niches and a climate function that determines the fitness
landscape. We empirically show that environmental dynamics affect plasticity
and evolvability differently and that the presence of diverse ecological niches
favors adaptability even in stable environments. We perform ablation studies of
the selection mechanisms to separate the role of fitness-based selection and
niche-limited competition. Results obtained from our minimal model allow us to
propose promising research directions in the study of open-endedness in
biological and artificial systems.


  We argue that information from countries who had earlier COVID-19 surges can
be used to inform another country's current model, then generating what we call
back-to-the-future (BTF) projections. We show that these projections can be
used to accurately predict future COVID-19 surges prior to an inflection point
of the daily infection curve. We show, across 12 different countries from all
populated continents around the world, that our method can often predict future
surges in scenarios where the traditional approaches would always predict no
future surges. However, as expected, BTF projections cannot accurately predict
a surge due to the emergence of a new variant. To generate BTF projections, we
make use of a matching scheme for asynchronous time series combined with a
response coaching SIR model.