  The existing concept of the "fitness value of information" provides a
theoretical upper bound on the fitness advantage of using information
concerning a fluctuating environment. Using concepts from rate-distortion
theory, we develop a theoretical framework to answer a different pair of
questions: What is the minimal amount of information needed for a population to
achieve a certain growth rate? What is the minimal amount of information gain
needed for one sub-population to achieve a certain average selection
coefficient over another? We introduce a correspondence between fitness and
distortion and solve for the rate-distortion functions of several systems using
analytical and numerical methods. Because accurate information processing is
energetically costly, our approach provides a theoretical basis for
understanding evolutionary "design principles" underlying information-cost
trade-offs.


  Nucleic acids physical properties have been investigated by theoretical
methods based both on fully atomistic representations and on coarse grained
models, e.g. the worm-like-chain, taken from polymer physics. In this article,
I present an intermediate (mesoscopic) approach and show how to build a three
dimensional Hamiltonian model which accounts for the main interactions
responsible for the stability of the helical molecules. While the 3D mesoscopic
model yields a sufficiently detailed description of the helix at the level of
the base pair, it also allows one to predict the thermodynamical and structural
properties of molecules in solution. Relying on the idea that the base pair
fluctuations can be conceived as trajectories, I have built a computational
method based on the time dependent path integral formalism to derive the
partition function. While the main features of the method are presented, I
focus here in particular on a newly developed statistical method to set the
maximum amplitude of the base pair fluctuations, a key parameter of the theory.
Some applications to the calculation of DNA flexibility properties are
discussed together with the available experimental data.


  A major goal shared by neuroscience and collective behavior is to understand
how dynamic interactions between individual elements give rise to behaviors in
populations of neurons and animals, respectively. This goal has recently become
within reach thanks to techniques providing access to the connectivity and
activity of neuronal ensembles as well as to behaviors among animal
collectives. The next challenge using these datasets is to unravel network
mechanisms generating population behaviors. This is aided by network theory, a
field that studies structure-function relationships in interconnected systems.
Here we review studies that have taken a network view on modern datasets to
provide unique insights into individual and collective animal behaviors.
Specifically, we focus on how analyzing signal propagation, controllability,
symmetry, and geometry of networks can tame the complexity of collective system
dynamics. These studies illustrate the potential of network theory to
accelerate our understanding of behavior across ethological scales.


  We consider the $SEIRS$ epidemiology model with such features of the COVID-19
outbreak as: abundance of unidentified infected individuals, limited time of
immunity and a possibility of vaccination. Within a compartmental realization
of this model, we found the disease-free and the endemic stationary states.
They exist in their respective restricted regions of the le via linear
stability analysis. The expression for the basic reproductive number is
obtained as well. The positions and heights of a first peak for the fractions
of infected individuals are obtained numerically and are fitted to simple
algebraic forms, that depend on model rates. Computer simulations of a
lattice-based realization for this model was performed by means of the cellular
automaton algorithm. These allowed to study the effect of the quarantine
measures explicitly, via changing the neighbourhood size. The attempt is made
to match both quarantine and vaccination measures aimed on balanced solution
for effective suppression of the pandemic.


  Bacteria are ubiquitous in our daily lives, either as motile planktonic cells
or as immobilized surface-attached biofilms. These different phenotypic states
play key roles in agriculture, environment, industry, and medicine; hence, it
is critically important to be able to predict the conditions under which
bacteria transition from one state to the other. Unfortunately, these
transitions depend on a dizzyingly complex array of factors that are determined
by the intrinsic properties of the individual cells as well as those of their
surrounding environments, and are thus challenging to describe. To address this
issue, here, we develop a generally-applicable biophysical model of the
interplay between motility-mediated dispersal and biofilm formation under
positive quorum sensing control. Using this model, we establish a universal
rule predicting how the onset and extent of biofilm formation depend
collectively on cell concentration and motility, nutrient diffusion and
consumption, chemotactic sensing, and autoinducer production. Our work thus
provides a key step toward quantitatively predicting and controlling biofilm
formation in diverse and complex settings.


  Understanding the thermodynamics of Darwinian evolution has important
implications for biophysics, evolutionary biology, and the study of the origin
of life. We show that for autocatalytic replicators in a nonequilibrium steady
state, the critical selection coefficient (minimal fitness difference visible
to selection) is lower bounded by the Gibbs free energy dissipated per
replication event. This bound presents a fundamental thermodynamic threshold
for Darwinian evolution, which is complementary to other thresholds that may
arise from finite population sizes, large mutation rates, etc. Our results
apply to a large class of molecular replicators, including many types of
autocatalytic sets, polymer-based replicators, and multistep autocatalytic
mechanisms. We illustrate our approach on a model of simple replicators in a
chemostat.


  Understanding the emergence of prosocial behaviours (e.g., cooperation and
trust) among self-interested agents is an important problem in many
disciplines. Network structure and institutional incentives (e.g., punishing
antisocial agents) are known to promote prosocial behaviours, when acting in
isolation, one mechanism being present at a time. Here we study the interplay
between these two mechanisms to see whether they are independent, interfering
or synergetic. Using evolutionary game theory, we show that punishing
antisocial agents and a regular networked structure not only promote prosocial
behaviours among agents playing the trust game, but they also interplay with
each other, leading to interference or synergy, depending on the game
parameters. Synergy emerges on a wider range of parameters than interference
does. In this domain, the combination of incentives and networked structure
improves the efficiency of incentives, yielding prosocial behaviours at a lower
cost than the incentive does alone. This has a significant implication in the
promotion of prosocial behaviours in multi-agent systems.


  Infectious disease epidemiologists routinely fit stochastic epidemic models
to time series data to elucidate infectious disease dynamics, evaluate
interventions, and forecast epidemic trajectories. To improve computational
tractability, many approximate stochastic models have been proposed. In this
paper, we focus on one class of such approximations -- time series
Susceptible-Infectious-Removed (TSIR) models. Infectious disease modeling often
starts with a homogeneous mixing assumption, postulating that the rate of
disease transmission is proportional to a product of the numbers of susceptible
and infectious individuals. One popular way to relax this assumption proceeds
by raising the number of susceptible and/or infectious individuals to some
positive powers. We show that when this technique is used within the TSIR
models they cannot be interpreted as approximate SIR models, which has
important implications for TSIR-based statistical inference. Our simulation
study shows that TSIR-based estimates of infection and mixing rates are
systematically biased in the presence of non-homogeneous mixing, suggesting
that caution is needed when interpreting TSIR model parameter estimates when
this class of models relaxes the homogeneous mixing assumption.


  Understanding kinetics and thermodynamics profile of biomolecules is
necessary to understand their functional roles which has a major impact in
mechanism driven drug discovery. Molecular dynamics simulation has been
routinely used to understand conformational dynamics and molecular recognition
in biomolecules. Statistical analysis of high-dimensional spatiotemporal data
generated from molecular dynamics simulation requires identification of few
low-dimensional variables which can describe essential dynamics of a system
without significant loss of informations. In physical chemistry, these
low-dimensional variables often called collective variables. Collective
variables are used to generated reduced representation of free energy surface
and calculate transition probabilities between different metastable basins.
However the choice of collective variables is not trivial for complex systems.
Collective variables ranges from geometric criteria's such as distances,
dihedral angles to abstract ones such as weighted linear combinations of
multiple geometric variables. Advent of machine learning algorithms led to
increasing use of abstract collective variables to represent biomolecular
dynamics. In this review, I will highlight several nuances of commonly used
collective variables ranging from geometric to abstract ones. Further, I will
put forward some cases where machine learning based collective variables were
used to describe simple systems which in principle could have been described by
geometric ones. Finally, I will put forward my thoughts on artificial general
intelligence and how it can be used to discover and predict collective
variables from spatiotemporal data generated by molecular dynamics simulations.


  Muilti-modality data are ubiquitous in biology, especially that we have
entered the multi-omics era, when we can measure the same biological object
(cell) from different aspects (omics) to provide a more comprehensive insight
into the cellular system. When dealing with such multi-omics data, the first
step is to determine the correspondence among different modalities. In other
words, we should match data from different spaces corresponding to the same
object. This problem is particularly challenging in the single-cell multi-omics
scenario because such data are very sparse with extremely high dimensions.
Secondly, matched single-cell multi-omics data are rare and hard to collect.
Furthermore, due to the limitations of the experimental environment, the data
are usually highly noisy. To promote the single-cell multi-omics research, we
overcome the above challenges, proposing a novel framework to align and
integrate single-cell RNA-seq data and single-cell ATAC-seq data. Our approach
can efficiently map the above data with high sparsity and noise from different
spaces to a low-dimensional manifold in a unified space, making the downstream
alignment and integration straightforward. Compared with the other
state-of-the-art methods, our method performs better in both simulated and real
single-cell data. The proposed method is helpful for the single-cell
multi-omics research. The improvement for integration on the simulated data is
significant.


  In this work we study the stability properties of the equilibrium points of
deterministic epidemic models with nonconstant population size. Models with
nonconstant population have been studied in the past only in particular cases,
two of which we review and combine. Our main result shows that for simple
"matrix epidemic models" introduced in [1], an explicit general formula for the
reproduction number and the corresponding "weak stability alternative" still
holds, under small modifications, for models with nonconstant population size,
and even when the model allows for vaccination and loss of immunity. The
importance of this result is clear once we note that the models of [1] include
a large number of viral and bacterial models of epidemic propagation, including
for example the totality of homogeneous COVID-19 models. To better understand
the nature of the result, we emphasize that the models proposed in [1] and
considered here are extensions of the SIR-PH model, which is essentially
characterized by a phase-type distribution that models transitions between the
"disease/infectious compartments". In these cases, the reproduction number and
a certain Lyapunov function for the disease free equilibrium are explicitly
expressible. Not surprisingly, accounting for varying demography, loss of
immunity, and vaccinations lead to several challenges. One of the most
important is that a varying population size leads to multiple endemic
equilibrium points: this is in contrast with "classic models" which in general
admit unique disease-free and endemic equilibria. As a special case of our
analysis, we consider a "first approximation" (FA) of our model, which
coincides with the constant-demography model often studied in the literature,
and for which more explicit results are available. Furthermore, we propose a
second heuristic approximation named "intermediate approximation" (IA).


  We present DEIMoS: Data Extraction for Integrated Multidimensional
Spectrometry, a Python application programming interface (API) and command-line
tool for high-dimensional mass spectrometry data analysis workflows that offers
ease of development and access to efficient algorithmic implementations.
Functionality includes feature detection, feature alignment, collision cross
section (CCS) calibration, isotope detection, and MS/MS spectral deconvolution,
with the output comprising detected features aligned across study samples and
characterized by mass, CCS, tandem mass spectra, and isotopic signature.
Notably, DEIMoS operates on N-dimensional data, largely agnostic to acquisition
instrumentation; algorithm implementations simultaneously utilize all
dimensions to (i) offer greater separation between features, thus improving
detection sensitivity, (ii) increase alignment/feature matching confidence
among datasets, and (iii) mitigate convolution artifacts in tandem mass
spectra. We demonstrate DEIMoS with LC-IMS-MS/MS data to illustrate the
advantages of a multidimensional approach in each data processing step.


  A Caputo-type fractional-order mathematical model for "metapopulation cholera
transmission" was recently proposed in [Chaos Solitons Fractals 117 (2018),
37--49]. A sensitivity analysis of that model is done here to show the accuracy
relevance of parameter estimation. Then, a fractional optimal control (FOC)
problem is formulated and numerically solved. A cost-effectiveness analysis is
performed to assess the relevance of studied control measures. Moreover, such
analysis allows us to assess the cost and effectiveness of the control measures
during intervention. We conclude that the FOC system is more effective only in
part of the time interval. For this reason, we propose a system where the
derivative order varies along the time interval, being fractional or classical
when more advantageous. Such variable-order fractional model, that we call a
'FractInt' system, shows to be the most effective in the control of the
disease.


  In this article, we propose a theory to explain the reduction in vaccine
effectiveness (VE) against the Delta SARS-CoV-2 variant and decreasing VE over
time reported in recent studies. Using a model illustration, we show that
in-host viral replication dynamics and delays in immune response could play a
key role in VE. Given this, current laboratory approaches solely measuring
reductions in neutralizing ability cannot fully represent the potential impact
of new SARS-CoV-2 variants. We instead propose an alternative approach that
incorporates viral replication dynamics into evaluations of SARS-CoV-2 variant
impact on immunity and VE. This more robust assessment may better inform public
health response to new variants like the newly detected Omicron variant.


  Accurate prediction of the number of daily or weekly confirmed cases of
COVID-19 is critical to the control of the pandemic. Existing mechanistic
models nicely capture the disease dynamics. However, to forecast the future,
they require the transmission rate to be known, limiting their prediction
power. Typically, a hypothesis is made on the form of the transmission rate
with respect to time. Yet the real form is too complex to be mechanistically
modeled due to the unknown dynamics of many influential factors. We tackle this
problem by using a hypothesis-free machine-learning algorithm to estimate the
transmission rate from data on non-pharmaceutical policies, and in turn
forecast the confirmed cases using a mechanistic disease model. More
specifically, we build a hybrid model consisting of a mechanistic ordinary
differential equation (ODE) model and a generalized boosting model (GBM). To
calibrate the parameters, we develop an "inverse method" that obtains the
transmission rate inversely in time from the other variables in the ODE model
and then feed it into the GBM to connect with the policy data. The resulting
model forecasted the number of daily confirmed cases up to 35 days in the
future in the United States with an averaged mean absolute percentage error of
27%. Being partly data-driven, the method is more accurate than typical
mechanistic models and meanwhile more intuitive, and possibly reliable, than
purely data-based machine learning models. Moreover, it can identify the most
informative predictive variables, which can be helpful in designing improved
forecasters as well as informing policymakers.


  Protein-protein interactions (PPIs) are essentials for many biological
processes where two or more proteins physically bind together to achieve their
functions. Modeling PPIs is useful for many biomedical applications, such as
vaccine design, antibody therapeutics, and peptide drug discovery. Pre-training
a protein model to learn effective representation is critical for PPIs. Most
pre-training models for PPIs are sequence-based, which naively adopt the
language models used in natural language processing to amino acid sequences.
More advanced works utilize the structure-aware pre-training technique, taking
advantage of the contact maps of known protein structures. However, neither
sequences nor contact maps can fully characterize structures and functions of
the proteins, which are closely related to the PPI problem. Inspired by this
insight, we propose a multimodal protein pre-training model with three
modalities: sequence, structure, and function (S2F). Notably, instead of using
contact maps to learn the amino acid-level rigid structures, we encode the
structure feature with the topology complex of point clouds of heavy atoms. It
allows our model to learn structural information about not only the backbones
but also the side chains. Moreover, our model incorporates the knowledge from
the functional description of proteins extracted from literature or manual
annotations. Our experiments show that the S2F learns protein embeddings that
achieve good performances on a variety of PPIs tasks, including cross-species
PPI, antibody-antigen affinity prediction, antibody neutralization prediction
for SARS-CoV-2, and mutation-driven binding affinity change prediction.


  By generating transient encounters between individuals beyond their immediate
social environment, transport can have a profound impact on the spreading of an
epidemic. In this work, we consider epidemic dynamics in the presence of the
transport process that gives rise to a multiplex network model. In addition to
a static layer, the (multiplex) epidemic network consists of a second dynamic
layer in which any two individuals are connected for the time they occupy the
same site during a random walk they perform on a separate transport network. We
develop a mean-field description of the stochastic network model and study the
influence the transport process has on the epidemic threshold. We show that any
transport process generally lowers the epidemic threshold because of the
additional connections it generates. In contrast, considering also random walks
of fractional order that in some sense are a more realistic model of human
mobility, we find that these non-local transport dynamics raise the epidemic
threshold in comparison to a classical local random walk. We also test our
model on a realistic transport network (the Munich U-Bahn network), and
carefully compare mean-field solutions with stochastic trajectories in a range
of scenarios.


  Synthetic genetic circuits are implemented in living cells for their
operation. During gene expression, proteins are produced from the respective
genes, by formation of complexes through the process of transcription and
translation. In transcription the circuit uses RNAP, etc. as resource from the
host cell and in translation, ribosome, tRNA and other cellular resources are
supplied to the operating circuit. As the cell contains these resources in
limited number, the circuit can suffer from unprecedented resource competition
which might destroy the circuit functionality, or introduce some emergent
responses. In this paper, we have studied a three-gene motif under resource
competition where interesting behaviour, similar to regulatory responses occur
due to limited supply of necessary resources. The system of interest exhibits
prominent changes in behaviour which can be observed experimentally. We focus
on two specific aspects, namely, dynamic range and half-induction point, which
inherently describe the circuit functionalities, and can be affected by
corresponding resource affinity and availability.


  Explaining biodiversity is a fundamental issue in ecology. A long-standing
puzzle lies in the paradox of the plankton: many species of plankton feeding on
a limited type of resources coexist, apparently flouting the competitive
exclusion principle (CEP), which holds that the number of predator (consumer)
species cannot exceed that of the resources at steady state. Here, we present a
mechanistic model and show that the intraspecific interference among the
consumers enables a plethora of consumer species to coexist at constant
population densities with only one or a handful of resource species. The
facilitated biodiversity is resistant to stochasticity, either with the
stochastic simulation algorithm or individual-based modeling. Our model
naturally explains the classical experiments that invalidate CEP,
quantitatively illustrates the universal S-shaped pattern of the rank-abundance
curves across a wide range of ecological communities, and can be broadly used
to resolve the mystery of biodiversity in many natural ecosystems.


  Recent advances in genomic sequencing technology have resulted in an
abundance of genome sequence data. Despite the progress in interpreting those
data, there remains a broad scope for their translation into clinical and
societal benefits. Loss-of-function variations in the human genome can be
causal in disease development. Precise identification of such variations and
pathogenicity prediction may lead to better drug targeting, among other
benefits. Machine learning comes across as a promising method for its proven
predictive ability. We have curated a novel dataset for the classification of
LOF variants using high-quality databases of genetic variation. We trained and
validated seven different classification algorithms using the new dataset to
classify the variants as Benign, Pathogenic and Likely pathogenic. We recorded
the best overall performance using the XG-Boost algorithm with an F1-score of
0.88 on the test set. We observed fair performance on Pathogenic samples with
high recall and moderate precision and subpar performance on Likely pathogenic
class, albeit with moderate precision. Overall, the encouraging results make
our final model a promising candidate for further real-world tests.


  The paper presents classification and analysis of the mathematical models of
COVID-19 spread in different groups of populations such as the family, school,
office (3-100 people), neighborhood (100-5000 people), city, region (0.5-15
million people), country, continent and the world. The classification covers
the main types of models including time-series, differential, imitation ones,
and their combinations. The time-series models are built from analysis of the
time series derived using filtration, regression and network methods (Section
2). The differential models include those derived from systems of ordinary and
stochastic differential equations as well as partial-derivative equations
(Section 3). The imitation models include cellular automata and agent-based
models (Section 4). The fourth group in the classification is combinations of
nonlinear Markov chains and optimal control theory, derived within the
framework of the mean-field game theory. Due to the novelty of the disease and
the difficulties it causes, the parameters of most models are, as a rule,
unknown, which necessitates one to solve the inverse problem, so the paper also
analyses the main algorithms to solve the inverse problem such as stochastic
optimization; nature-like algorithms (genetic; differential evolution; particle
swarm, etc.); the understanding method; big-data analysis, and machine
learning.


  The purpose of this paper is to present a model of a phenomenon of plant stem
morphogenesis observed by Cesar Gomez-Campo in 1970. We consider a simplified
model of auxin dynamics in plant stems, showing that, after creation of the
original primordium, it can represent random, distichous and spiral phyllotaxis
(leaf arrangement) just by varying one parameter, the rate of diffusion. The
same analysis extends to the $n$-jugate case where $n$ primordia are initiated
at each plastochrone. Having validated the model, we consider how it can give
rise to the Gomez-Campo phenomenon, showing how a stem with spiral phyllotaxis
can produce branches of the same or opposite chirality. And finally, how the
relationship can change from discordant to concordant over the course of a
growing season.


  The basic reproduction number, $R_0$ is an important and widely used concept
in the study of infectious diseases. We briefly review the recent trend of
calculating the average of various $R_0$ estimates in systematic reviews aimed
at estimating the basic reproduction number of SARS-CoV-2, and discuss the
drawbacks and implications of using such averaging methods. Additionally, we
argue that even a theoretically grounded approach such as next generation
matrix could have practical impediments in its use. More generally, the
practice of associating an infectious disease with a single value of $R_0$ is
problematic, when the disease can, in fact have different reproduction numbers
in various populations.


  Intrinsically disordered proteins (IDPs) and multidomain proteins with
flexible linkers show a high level of structural heterogeneity and are best
described by ensembles consisting of multiple conformations with associated
thermodynamic weights. Determining conformational ensembles usually involves
integration of biophysical experiments and computational models. In this
review, we discuss current approaches to determining conformational ensembles
of IDPs and multidomain proteins, including the choice of biophysical
experiments, computational models used to sample protein conformations, models
to calculate experimental observables from protein structure, and methods to
refine ensembles against experimental data. We also provide examples of recent
applications of integrative conformational ensemble determination to study IDPs
and multidomain proteins and suggest future directions for research in the
field.


  The SARS-CoV-2 pandemic has so far unfolded diversely across the fifty United
States of America, reflected both in different time progressions of infection
"waves" and in magnitudes of local infection rates. Despite a marked diversity
of presentations, most U.S. states experienced their single greatest surge in
daily new cases during the transition from Fall 2020 to Winter 2021. Popular
media also cite additional similarities between states -- often despite
disparities in governmental policies, reported mask-wearing compliance rates,
and vaccination percentages. Here, we identify a set of robust, low-dimensional
clusters that 1) summarize the timings and relative heights of four historical
COVID-19 "wave opportunities" accessible to all 50 U.S. states, 2) correlate
with geographical and intervention patterns associated with those groups of
states they encompass, and 3) predict aspects of the "fifth wave" of new
infections in the late Summer of 2021. In particular, we argue that clustering
elucidates a negative relationship between vaccination rates and subsequent
case-load variabilities within state groups. We advance the hypothesis that
vaccination acts as a ``seat belt," in effect constraining the likely range of
new-case upticks, even in the context of the Summer 2021, variant-driven surge.


  The analysis of single particle trajectories plays an important role in
elucidating dynamics within complex environments such as those found in living
cells. However, the characterization of intracellular particle motion is often
confounded by confinement of the particles within non-trivial subcellular
geometries. Here, we focus specifically on the case of particles undergoing
Brownian motion within a tubular network, as found in some cellular organelles.
An unraveling algorithm is developed to uncouple particle motion from the
confining network structure, allowing for an accurate extraction of the
diffusion coefficient, as well as differentiating between Brownian and
fractional Brownian dynamics. We validate the algorithm with simulated
trajectories and then highlight its application to an example system: analyzing
the motion of membrane proteins confined in the tubules of the peripheral
endoplasmic reticulum in mammalian cells. We show that these proteins undergo
diffusive motion with a well-characterized diffusivity. Our algorithm provides
a generally applicable approach for disentangling geometric morphology and
particle dynamics in networked architectures.


  A manuscript identified bat sarbecoviruses with high sequence homology to
SARS-CoV-2 found in caves in Laos that can directly infect human cells via the
human ACE2 receptor (Coronaviruses with a SARS-CoV-2-like receptor binding
domain allowing ACE2-mediated entry into human cells isolated from bats of
Indochinese peninsula, Temmam S., et al.). Here, I examine the genomic sequence
of one of these viruses, BANAL-236, and show it has 5-UTR and 3-UTR secondary
structures that are non-canonical and, in fact, have never been seen in an
infective coronavirus. Specifically, the 5-UTR has a 177 nt copy-back extension
which forms an extended, highly stable duplex RNA structure. Because of this
copy-back, the four obligate Stem Loops (SL) -1, -2, -3, and -4 cis-acting
elements found in all currently known replicating coronaviruses are buried in
the extended duplex. The 3-UTR has a similar fold-back duplex of 144 nts and is
missing the obligate poly-A tail. Taken together, these findings demonstrate
BANAL-236 is missing eight obligate UTR cis-acting elements; each one of which
has previously been lethal to replication when modified individually. Neither
duplex copyback has ever been observed in an infective sarbecovirus, although
some of the features have been seen in defective interfering particles, which
can be found in co-infections with non-defective, replicating viruses. They are
also a common error seen during synthetic genome assembly in a laboratory.
BANAL-236 must have evolved an entirely unique mechanism for replication, RNA
translation, and RNA packaging never seen in a coronavirus and because it is a
bat sarbecovirus closely related to SARS-CoV-2, it is imperative that we
understand its unique mode of infectivity by a collaborative, international
research effort.


  Clonal microbes can switch between different phenotypes and recent
theoretical work has shown that stochastic switching between these
subpopulations can stabilize microbial communities. This phenotypic switching
need not be stochastic, however, but could also be in response to environmental
factors, both biotic and abiotic. Here, motivated by the bacterial persistence
phenotype, we explore the ecological effects of such responsive switching by
analyzing phenotypic switching in response to competing species. We show that
the stability of microbial communities with responsive switching differs
generically from that of communities with stochastic switching only. To
understand the mechanisms by which responsive switching stabilizes coexistence,
we go on to analyze simple two-species models. Combining exact results and
numerical simulations, we extend the classical stability results for the
competition of two species without phenotypic variation to the case in which
one species switches, stochastically and responsively, between two phenotypes.
In particular, we show that responsive switching can stabilize coexistence even
when stochastic switching on its own does not affect the stability of the
community.


  Reliable short term forecasting can provide potentially lifesaving insights
into logistical planning, and in particular, into the optimal allocation of
resources such as hospital staff and equipment. By reinterpreting COVID-19
daily cases in terms of candlesticks, we are able to apply some of the most
popular stock market technical indicators to obtain predictive power over the
course of the pandemics. By providing a quantitative assessment of MACD, RSI,
and candlestick analyses, we show their statistical significance in making
predictions for both stock market data and WHO COVID-19 data. In particular, we
show the utility of this novel approach by considering the identification of
the beginnings of subsequent waves of the pandemic. Finally, our new methods
are used to assess whether current health policies are impacting the growth in
new COVID-19 cases.


  Adoption of recently developed methods from machine learning has given rise
to creation of drug-discovery knowledge graphs (KG) that utilize the
interconnected nature of the domain. Graph-based modelling of the data,
combined with KG embedding (KGE) methods, are promising as they provide a more
intuitive representation and are suitable for inference tasks such as
predicting missing links. One common application is to produce ranked lists of
genes for a given disease, where the rank is based on the perceived likelihood
of association between the gene and the disease. It is thus critical that these
predictions are not only pertinent but also biologically meaningful. However,
KGs can be biased either directly due to the underlying data sources that are
integrated or due to modeling choices in the construction of the graph, one
consequence of which is that certain entities can get topologically
overrepresented. We demonstrate the effect of these inherent structural
imbalances, resulting in densely-connected entities being highly ranked no
matter the context. We provide support for this observation across different
datasets, models as well as predictive tasks. Further, we present various graph
perturbation experiments which yield more support to the observation that KGE
models can be more influenced by the frequency of entities rather than any
biological information encoded within the relations. Our results highlight the
importance of data modeling choices, and emphasizes the need for practitioners
to be mindful of these issues when interpreting model outputs and during KG
composition.


  Understanding the structure of the human T-cell receptor repertoire is a
crucial precondition to understand the ability of the immune system to
recognize and respond to antigens. T-cells are often compared via the
complementarity determining region 3 (CDR3) of their respective T-cell receptor
beta chains. Nevertheless, previous studies often simply compared if CDR3beta
sequences were equal, while network theory studies were usually limited to
several thousand sequences due to the high computational effort of constructing
the network. To overcome that hurdle, we introduce the GPU-based algorithm
TCR-NET to construct large-scale CDR3beta similarity networks using
model-generated and empirical sequence data with up to 800,000 CDR3beta
sequences on a normal computer for the first time. Using network analysis
methods we study the structural properties of these networks and conclude that
(i) the fraction of public TCRs depends on the size of the TCR repertoire,
along with the exact (not unified) definition of "public" sequences, (ii) the
TCR network is assortative with the average neighbor degree being proportional
to the squareroot of the degree of a node and (iii) the repertoire is robust
against losses of TCRs. Moreover, we analyze the networks of antigen-specific
TCRs for different antigen families and find differing clustering coefficients
and assortativities. TCR-NET offers better access to assess large-scale TCR
repertoire networks, opening the possibility to quantify their structure and
quantitatively distinguish their ability to react to antigens, which we
anticipate to become a useful tool in a time of increasingly large amounts of
repertoire sequencing data becoming available.


  Species-rich communities, such as the microbiota or microbial ecosystems,
provide key functions for human health and climatic resilience. Increasing
effort is being dedicated to design experimental protocols for selecting
community-level functions of interest. These experiments typically involve
selection acting on populations of communities, each of which is composed of
multiple species. If numerical simulations started to explore the evolutionary
dynamics of this complex, multi-scale system, a comprehensive theoretical
understanding of the process of artificial selection of communities is still
lacking. Here, we propose a general model for the evolutionary dynamics of
communities composed of a large number of interacting species, described by
disordered generalised Lotka-Volterra equations. Our analytical and numerical
results reveal that selection for scalar community functions leads to the
emergence, along an evolutionary trajectory, of a low-dimensional structure in
an initially featureless interaction matrix. Such structure reflects the
combination of the properties of the ancestral community and of the selective
pressure. Our analysis determines how the speed of adaptation scales with the
system parameters and the abundance distribution of the evolved communities.
Artificial selection for larger total abundance is thus shown to drive
increased levels of mutualism and interaction diversity. Inference of the
interaction matrix is proposed as a method to assess the emergence of
structured interactions from experimentally accessible measures.


  Microtubules are dynamic intracellular fibers that have been observed
experimentally to undergo spontaneous self-alignment. We formulate a 3D
mean-field theory model to analyze the nematic phase transition of microtubules
growing and interacting within a 3D space then make a comparison with
computational simulations. We identify a control parameter $G_\text{eff}$ and
predict a unique critical value $G_\text{eff}=1.56$ for which a phase
transition can occur. Furthermore, we show both analytically and using
simulations that this predicted critical value does not depend on the presence
of zippering. The mean-field theory developed here provides an analytical
estimate of microtubule patterning characteristics without running
time-consuming simulations and is a step towards bridging scales from
microtubule behavior to multicellular simulations.


  The Susceptible-Infectious-Recovered (SIR) equations and their extensions
comprise a commonly utilized set of models for understanding and predicting the
course of an epidemic. In practice, it is of substantial interest to estimate
the model parameters based on noisy observations early in the outbreak, well
before the epidemic reaches its peak. This allows prediction of the subsequent
course of the epidemic and design of appropriate interventions. However,
accurately inferring SIR model parameters in such scenarios is problematic.
This article provides novel, theoretical insight on this issue of practical
identifiability of the SIR model. Our theory provides new understanding of the
inferential limits of routinely used epidemic models and provides a valuable
addition to current simulate-and-check methods. We illustrate some practical
implications through application to a real-world epidemic data set.


  This study demonstrates how to construct the solutions of a more general form
of population dynamics models via a blend of variational iterative method with
Sumudu transform. In this paper, population growth models are formulated in the
form of delay differential equations of pantograph type which is a general form
for the existing models. Innovative ways are presented for obtaining the
solutions of population growth models where other analytic methods fail.
Stimulating procedures for finding patterns and regularities in seemingly
chaotic processes have been elucidated in this paper. How, when and why the
changes in population sizes occur can be deduced through this study.


  The state of the art of epidemic modelling in terrestrial ecosystems is the
compartmental SIR model and its extensions from the now classical work of
Kermack-Mackendrick. In contrast, epidemic modelling of marine ecosystems is a
bit behind, and compartmental models have been introduced only recently. One of
the reasons is that many epidemic processes in terrestrial ecosystems can be
described through a contact process, while modelling marine epidemics is more
subtle in many cases. Here we present a model describing disease outbreaks
caused by parasites in bivalve populations. The SIRP model is a
multicompartmental model with four compartments, three of which describe the
different states of the host, susceptible (i.e. healthy), S, infected, I, and
removed (dead), R, and one compartment for the parasite in the marine medium,
P, written as a 4-dimensional dynamical system. Even if this is the simplest
model one can write to describe this system, it is still too complicated for
both direct analytical manipulation and direct comparison with experimental
observations, as it depends on four parameters to be fitted. We show that it is
possible to simplify the model, including a reduction to the standard SIR model
if the parameters fulfil certain conditions. The model is validated with
available data for the recent Mass Mortality Event of the noble pen shell Pinna
nobilis, a disease caused by the parasite Haplosporidium pinnae, showing that
the reduced SIR model is able to fit the data. So, we show that a model in
which the species that suffers the epidemics (host) cannot move, and contagion
occurs through parasites, can be reduced to the standard SIR model that
represents epidemic transmission between mobile hosts. The fit indicates that
the assumptions made to simplify the model are reasonable in practice, although
it leads to an indeterminacy in three of the original parameters.


  PEACH Tree is an easy-to-use, online tool for displaying multiple sequence
alignments and phylogenetic trees side-by-side. PEACH Tree is powerful for
rapidly tracing evolutionary and transmission histories by filtering invariant
sites out of the display, and allowing samples to readily be filtered out of
the display. These features, coupled with the ability to display
epidemiological metadata, make the tool suitable for infectious disease
epidemiology. PEACH Tree further enables much needed communication between the
fields of genomics and infectious disease epidemiology, as exemplified by the
COVID-19 pandemic.


  Large scale self-supervised pre-training of Transformer language models has
advanced the field of Natural Language Processing and shown promise in
cross-application to the biological `languages' of proteins and DNA. Learning
effective representations of DNA sequences using large genomic sequence
corpuses may accelerate the development of models of gene regulation and
function through transfer learning. However, to accurately model cell
type-specific gene regulation and function, it is necessary to consider not
only the information contained in DNA nucleotide sequences, which is mostly
invariant between cell types, but also how the local chemical and structural
`epigenetic state' of chromosomes varies between cell types. Here, we introduce
a Bidirectional Encoder Representations from Transformers (BERT) model that
learns representations based on both DNA sequence and paired epigenetic state
inputs, which we call Epigenomic BERT (or EBERT). We pre-train EBERT with a
masked language model objective across the entire human genome and across 127
cell types. Training this complex model with a previously prohibitively large
dataset was made possible for the first time by a partnership with Cerebras
Systems, whose CS-1 system powered all pre-training experiments. We show
EBERT's transfer learning potential by demonstrating strong performance on a
cell type-specific transcription factor binding prediction task. Our fine-tuned
model exceeds state of the art performance on 4 of 13 evaluation datasets from
ENCODE-DREAM benchmarks and earns an overall rank of 3rd on the challenge
leaderboard. We explore how the inclusion of epigenetic data and task specific
feature augmentation impact transfer learning performance.


  Comparative studies of cancer-related genes allow us to gain novel
information about the evolution and function of these genes, but also to
understand cancer as a driving force in biological systems and species life
histories. So far, comparative studies of cancer genes have focused on mammals.
Here, we provide the first comparative study of cancer-related gene copy number
variation in fish. As fish are evolutionarily older and genetically more
diverse than mammals, their tumour suppression mechanisms should not only
include most of the mammalian mechanisms, but also reveal novel (but
potentially phylogenetically older) previously undetected mechanisms. We have
matched the sequenced genomes of 65 fish species from the Ensemble database
with the cancer gene information from the COSMIC database. By calculating the
number of gene copies across species using the Ensembl CAFE data (providing
species trees for gene copy number counts), we were able to develop a novel,
less resource demanding method for ortholog identification. Our analysis
demonstrates a masked relationship with cancer-related gene copy number
variation (CNV) and maximum lifespan in fish species, suggesting that higher
tumour suppressor gene CNV lengthens and oncogene CNV shortens lifespan, when
both traits are added to the model. Based on the correlation between tumour
suppressor and oncogene CNV, we were able to show which species have more
tumour suppressors in relation to oncogenes. It could therefore be suggested
that these species have stronger genetic defences against oncogenic processes.
Fish studies could yet be a largely unexplored treasure trove for understanding
the evolution and ecology of cancer, by providing novel insights into the study
of cancer and tumour suppression, in addition to the study of fish evolution,
life-history trade-offs, and ecology.


  The control of biopolymer length is mediated by proteins that localize to
polymer ends and regulate polymerization dynamics. Several mechanisms have been
proposed to achieve end localization. Here, we propose a novel mechanism by
which a protein that binds to a shrinking polymer and slows its shrinkage will
be spontaneously enriched at the shrinking end through a "herding" effect. We
formalize this process using both a lattice-gas model and a continuum
description, and we present experimental evidence that the microtubule
regulator spastin employs this mechanism. Our findings extend to more general
problems involving diffusion within shrinking domains.


  The increasing spread of COVID-19, caused by the virus SARS-CoV-2, raises
concerns about the extent to which mutations have occurred across the viral
genome. We present a partial replication of an earlier 2021 study by Wang, R.
et al. that determined the presence of four substrains and eleven top mutations
in the United States. We analyze a portion of the authors' data set in order to
recreate Figure S1 from the paper, recapitulating the same features observed in
the original figure. We further generate a summary of mutation characteristics
for each of the 26 named proteins and confirm the significance of the spike
protein at roughly 24% of all recorded mutations. Our analysis suggests that
additional factors may affect per-protein mutation rate besides protein length.


  Living things benefit from exquisite molecular sensitivity in many of their
key processes, including DNA replication, transcription and translation,
chemical sensing, and morphogenesis. At thermodynamic equilibrium, the basic
biophysical mechanism for sensitivity is cooperative binding, for which it can
be shown that the Hill coefficient, a sensitivity measure, cannot exceed the
number of binding sites. Generalizing this fact, we find that for any kinetic
scheme, at or away from thermodynamic equilibrium, a very simple structural
quantity, the size of the support of a perturbation, always limits the
effective Hill coefficient. This support bound sheds light on and unifies
diverse sensitivity mechanisms, ranging from kinetic proofreading to a
nonequilibrium Monod-Wyman-Changeux (MWC) model proposed for the E. coli
flagellar motor switch, and represents a simple, precise bridge between
experimental observations and the models we write down. In pursuit of
mechanisms that saturate the support bound, we find a nonequilibrium binding
mechanism, nested hysteresis, with sensitivity exponential in the number of
binding sites, with implications for our understanding of models of gene
regulation.


  In response to pathogens, the adaptive immune system generates specific
antibodies that bind and neutralize foreign antigens. Understanding the
composition of an individual's immune repertoire can provide insights into this
process and reveal potential therapeutic antibodies. In this work, we explore
the application of antibody-specific language models to aid understanding of
immune repertoires. We introduce AntiBERTy, a language model trained on 558M
natural antibody sequences. We find that within repertoires, our model clusters
antibodies into trajectories resembling affinity maturation. Importantly, we
show that models trained to predict highly redundant sequences under a multiple
instance learning framework identify key binding residues in the process. With
further development, the methods presented here will provide new insights into
antigen binding from repertoire sequences alone.


  In conservation biology, phylogenetic diversity (PD) provides a way to
quantify the impact of the current rapid extinction of species on the
evolutionary `Tree of Life'. This approach recognises that extinction not only
removes species but also the branches of the tree on which unique features
shared by the extinct species arose. In this paper, we investigate three
questions that are relevant to PD. The first asks how many sets of species of
given size $k$ preserve the maximum possible amount of PD in a given tree. The
number of such maximum PD sets can be very large, even for moderate-sized
phylogenies. We provide a combinatorial characterisation of maximum PD sets,
focusing on the setting where the branch lengths are ultrametric (e.g.
proportional to time). This leads to a polynomial-time algorithm for
calculating the number of maximum PD sets of size $k$ by applying a generating
function; we also investigate the types of tree shapes that harbour the most
(or fewest) maximum PD sets of size $k$. Our second question concerns
optimising a linear function on the species (regarded as leaves of the
phylogenetic tree) across all the maximum PD sets of a given size. Using the
characterisation result from the first question, we show how this optimisation
problem can be solved in polynomial time, even though the number of maximum PD
sets can grow exponentially. Our third question considers a dual problem: If
$k$ species were to become extinct, then what is the largest possible {\em
loss} of PD in the resulting tree? For this question, we describe a
polynomial-time solution based on dynamical programming.


  Predicting drug-drug interactions (DDI) is the problem of predicting side
effects (unwanted outcomes) of a pair of drugs using drug information and known
side effects of many pairs. This problem can be formulated as predicting labels
(i.e. side effects) for each pair of nodes in a DDI graph, of which nodes are
drugs and edges are interacting drugs with known labels. State-of-the-art
methods for this problem are graph neural networks (GNNs), which leverage
neighborhood information in the graph to learn node representations. For DDI,
however, there are many labels with complicated relationships due to the nature
of side effects. Usual GNNs often fix labels as one-hot vectors that do not
reflect label relationships and potentially do not obtain the highest
performance in the difficult cases of infrequent labels. In this paper, we
formulate DDI as a hypergraph where each hyperedge is a triple: two nodes for
drugs and one node for a label. We then present CentSmoothie, a hypergraph
neural network that learns representations of nodes and labels altogether with
a novel central-smoothing formulation. We empirically demonstrate the
performance advantages of CentSmoothie in simulations as well as real datasets.


  Affinity maturation is crucial for improving the binding affinity of
antibodies to antigens. This process is mainly driven by point substitutions
caused by somatic hypermutations of the immunoglobulin gene. It also includes
deletions and insertions of genomic material known as indels. While the
landscape of point substitutions has been extensively studied, a detailed
statistical description of indels is still lacking. Here we present a
probabilistic inference tool to learn the statistics of indels from repertoire
sequencing data, which overcomes the pitfalls and biases of standard annotation
methods. The model includes antibody-specific maturation ages to account for
variable mutational loads in the repertoire. After validation on synthetic
data, we applied our tool to a large dataset of human immunoglobulin heavy
chains. The inferred model allows us to identify universal statistical features
of indels in heavy chains. We report distinct insertion and deletion hotspots,
and show that the distribution of lengths of indels follows a geometric
distribution, which puts constraints on future mechanistic models of the
hypermutation process.


  This theory seeks to define species and to explore evolutionary forces and
genetic elements in speciation and species maintenance. The theory explains how
speciation and species maintenance are caused by natural selection acting on
non-Mendelian and Mendelian variation, respectively. The emergence and
maintenance of species as groups of populations are balanced by evolutionary
forces including complementary mechanisms of gene flow within and between
populations at population-specific rates: sexual and asexual reproduction,
recombining and nonrecombining genome regions, vertical and horizontal DNA
transfer, and transposon proliferation and control. While recombining genome
regions carry conserved genes and are subjected to meiotic recombination,
nonrecombining genome regions carry accessory genes and are not subjected to
such structural restrain. Sexual reproduction, vertical DNA transfer,
recombining genome regions and transposon control keep species in existence by
maintaining recombining chromosome number and structure, while asexual
reproduction, horizontal DNA transfer, nonrecombining genome regions and
transposon proliferation help species emerge by promoting reproductive
isolation and changes in chromosome number and structure. The theory is based
on the analysis of the genome sequences of isolates in the Fusarium oxysporum
complex. The rate of horizontal supernumerary chromosome transfer in this
complex was estimated to be 0.1 per genome per year.


  The intelligibility of speech relies on the ability of interlocutors to
dynamically align their expectations about the rates at which informative
changes in signals occur. Exactly how this is achieved remains an open
question. We propose that speaker alignment is supported by the statistical
structure of spoken signals and show how pauses offer a time-invariant template
for structuring speech sequences. Consistent with this, we show that pause
distributions in conversational English and Korean provide a memoryless
information source. We describe how this can facilitate both the initial
structuring and maintenance of predictability in spoken signals over time, and
show how the properties of this signal change predictably with speaker
experience. These results indicate that pauses provide a structuring signal
that interacts with the morphological and rhythmical structure of languages,
allowing speakers at all stages of lifespan development to distinguish signal
from noise and maintain mutual predictability in time.


  Accurate drug response prediction (DRP) is a crucial yet challenging task in
precision medicine. This paper presents a novel Attention-Guided Multi-omics
Integration (AGMI) approach for DRP, which first constructs a Multi-edge Graph
(MeG) for each cell line, and then aggregates multi-omics features to predict
drug response using a novel structure, called Graph edge-aware Network (GeNet).
For the first time, our AGMI approach explores gene constraint based
multi-omics integration for DRP with the whole-genome using GNNs. Empirical
experiments on the CCLE and GDSC datasets show that our AGMI largely
outperforms state-of-the-art DRP methods by 8.3%--34.2% on four metrics. Our
data and code are available at https://github.com/yivan-WYYGDSG/AGMI.


  The size and shape of the region affected by an outbreak is relevant to
understand the dynamics of a disease and help to organize future actions to
mitigate similar events. A simple extension of the SIR model is considered,
where agents diffuse on a regular lattice and the disease may be transmitted
when an infected and a susceptible agents are nearest neighbors. We study the
geometric properties of both the connected cluster of sites visited by infected
agents (outbreak cluster) and the set of clusters with sites that have not been
visited. By changing the density of agents, our results show that there is a
mixed-order (hybrid) transition where the region affected by the disease is
finite in one phase but percolates through the system beyond the threshold.
Moreover, the outbreak cluster seems to have the same exponents of the backbone
of the critical cluster of the ordinary percolation while the clusters with
unvisited sites have a size distribution with a Fisher exponent $\tau<2$.


  Motivation: Identifying drug-target interactions (DTIs) is a key step in drug
repositioning. In recent years, the accumulation of a large number of genomics
and pharmacology data has formed mass drug and target related heterogeneous
networks (HNs), which provides new opportunities of developing HN-based
computational models to accurately predict DTIs. The HN implies lots of useful
information about DTIs but also contains irrelevant data, and how to make the
best of heterogeneous networks remains a challenge. Results: In this paper, we
propose a heterogeneous graph automatic meta-path learning based DTI prediction
method (HampDTI). HampDTI automatically learns the important meta-paths between
drugs and targets from the HN, and generates meta-path graphs. For each
meta-path graph, the features learned from drug molecule graphs and target
protein sequences serve as the node attributes, and then a node-type specific
graph convolutional network (NSGCN) which efficiently considers node type
information (drugs or targets) is designed to learn embeddings of drugs and
targets. Finally, the embeddings from multiple meta-path graphs are combined to
predict novel DTIs. The experiments on benchmark datasets show that our
proposed HampDTI achieves superior performance compared with state-of-the-art
DTI prediction methods. More importantly, HampDTI identifies the important
meta-paths for DTI prediction, which could explain how drugs connect with
targets in HNs.


  Generating the hash values of short subsequences, called seeds, enables
quickly identifying similarities between genomic sequences by matching seeds
with a single lookup of their hash values. However, these hash values can be
used only for finding exact-matching seeds as the conventional hashing methods
assign distinct hash values for different seeds, including highly similar
seeds. Finding only exact-matching seeds causes either 1) increasing the use of
the costly sequence alignment or 2) limited sensitivity.
  We introduce BLEND, the first efficient and accurate mechanism that can
identify both exact-matching and highly similar seeds with a single lookup of
their hash values, called fuzzy seed matches. BLEND 1) utilizes a technique
called SimHash, that can generate the same hash value for similar sets, and 2)
provides the proper mechanisms for using seeds as sets with the SimHash
technique to find fuzzy seed matches efficiently.
  We show the benefits of BLEND when used in read overlapping and read mapping.
For read overlapping, BLEND is faster by 2.4x - 83.9x (on average 19.3x), has a
lower memory footprint by 0.9x - 14.1x (on average 3.8x), and finds higher
quality overlaps leading to accurate de novo assemblies than the
state-of-the-art tool, minimap2. For read mapping, BLEND is faster by 0.8x -
4.1x (on average 1.7x) than minimap2. Source code is available at
https://github.com/CMU-SAFARI/BLEND.


  COVID-19 pandemic is one of the major disasters that humanity has ever faced.
In this paper, we try to model the effect of vaccination in controlling the
pandemic, particularly in context to the third wave which is predicted to hit
globally. Here we have modified the SEIRD model by introducing a vaccination
term. One of our main assumptions is that the infection rate (\b{eta}(t)) is
oscillatory. This oscillatory nature has been discussed earlier in literature
with reference to the seasonality of epidemics. However, in our case we invoke
this nature of the infection rate (\b{eta}(t)) to model the cyclical behavior
of the COVID-19 pandemic within a short period. This study focuses on a
minimalistic approach where we have logically deduced that the infection rate
(\b{eta}(t)) and the vaccination rate ({\lambda}) are the most important
parameters while the other parameters can be assumed to be constants throughout
the simulation. Finally, we have studied the rich interplay between the
infection rate (\b{eta}(t)) and the vaccination rate ({\lambda}) on the
infectious cases of COVID-19 and made some robust conclusions regarding the
global behavior of this pandemic in near future.


  Protein folding is an indispensable process for the majority of proteins
after their synthesis from ribosomes in the cell. Most in vitro protein folding
studies have focused on single-domain proteins. Hence, it is important to
understand the folding process of multi-domain proteins, especially when
domains are discontinuous. We choose the Maltose binding protein (MBP) as a
model system. In particular, we studied a mutant of MBP that folds slowly.
Here, using two- and three-color single-molecule Foerster resonance energy
transfer (smFRET) experiments, we study the refolding of both the domains and
the interaction between the domains of DM-MBP. Initial two-color smFRET
measurements of the N-terminal domain (NTD) reveal the presence of a folding
intermediate. The same folding intermediate is observed in measurements
monitoring the C- terminal domain (CTD) and the NTD-CTD (N-C) interface. The
refolding intermediate is dynamic on the sub-millisecond timescale.
Quantitative analysis on underlying dynamic interconversions revealed a delay
in NTD folding imposed by the entropic barrier being the primary cause for slow
DM-MBP folding. Moreover, CTD folds after NTD completes the folding. Using
three-color smFRET, we could show the NTD folds first and CTD later in a same
protein. Molecular dynamic simulations for temperature-induced unfolding on WT-
and DM-MBP identify a folding nucleus in the NTD, which is rich in hydrophobic
residues, and explains why the two mutations slow down the folding kinetics. In
the presence of the bacterial Hsp60 chaperonin system GroEL/ES, we observe that
DM-MBP is dynamic within the chaperonin cavity but the chaperonin limits the
conformational space of substrate. Hence, confinement aids DM-MBP in overcoming
the entropic barrier. The study reports on the subtle tuning and co-dependency
for protein folding between two-domains with a discontinuous arrangement.


  The "genetic purity" of French wild boar populations has been monitored since
the 1980s based on a cytogenetic difference between wild boars and domestic
pigs (36 and 38 chromosomes, respectively). This difference makes it possible
to identify any boar with 37 or 38 chromosomes as "hybrid", without however
being able to determine the origin (recent or ancient) of the hybridization,
nor guarantee the "purity" of an animal with 36 chromosomes. Analysis of
results of more than 4,600 tests performed over the last 12 years reveals an
average "hybrid" rate of 15.8%, with high variability between populations. To
analyse hybridization in greater detail and overcome inherent limitations of
the cytogenetic approach, 362 wild boars recently collected in different
regions of France were genotyped on a 70K SNP (GeneSeek GGP Porcine HD) chip.
This study showed that for 96.4% of the wild boars analysed, includingmost of
those with 37 or 38 chromosomes, the percentage of the genome of "domestic pig"
origin varied from 0 to 18%. This suggests that hybridization is a fairly
common phenomenon but of moderate intensity, and often ancient. Nevertheless,
higher rates of hybridization have been observed in some regions such as
Ard{\`e}che, and several cases of recent hybridization with domestic pigs were
found in 3.6% of the wild boars analysed, most of them with pigs of Asian
origin.


  The growth and evolution of microbial populations is often subjected to
advection by fluid flows in spatially extended environments, with immediate
consequences for questions of spatial population genetics in marine ecology,
planktonic diversity and origin of life scenarios. Here, we review recent
progress made in understanding this rich problem in the simplified setting of
two competing genetic microbial strains subjected to fluid flows. As a
pedagogical example we focus on antagonsim, i.e., two killer microorganism
strains, each secreting toxins that impede the growth of their competitors
(competitive exclusion), in the presence of stationary fluid flows. By solving
two coupled reaction-diffusion equations that include advection by simple
steady cellular flows composed of characteristic flow motifs in two dimensions
(2d), we show how local flow shear and compressibility effects can interact
with selective advantage to have a dramatic influence on genetic competition
and fixation in spatially distributed populations. We analyze several 1d and 2d
flow geometries including sources, sinks, vortices and saddles, and show how
simple analytical models of the dynamics of the genetic interface can be used
to shed light on the nucleation, coexistence and flow-driven instabilities of
genetic drops. By exploiting an analogy with phase separation with nonconserved
order parameters, we uncover how these genetic drops harness fluid flows for
novel evolutionary strategies, even in the presence of number fluctuations, as
confirmed by agent-based simulations as well.


  The vast majority of multi-cellular organisms are anisogamous, meaning that
male and female sex cells differ in size. It remains an open question how this
asymmetric state evolved, presumably from the symmetric isogamous state where
all gametes are roughly the same size (drawn from the same distribution). Here,
we use tools from the study of nonlinear dynamical systems to develop a simple
mathematical model for this phenomenon. Using theoretical analysis and
numerical simulation, we demonstrate that competition between individuals that
is linked to the mean gamete size will almost inevitably result in a stable
anisogamous equilibrium, and thus isogamy may naturally lead to anisogamy.


  Inference of the marginal likelihood of sample allele configurations using
backward algorithms yields identical results with the Kingman coalescent, the
Moran model, and the diffusion model (up to a scaling of time). For inference
of probabilities of ancestral population allele frequencies at any given point
in the past - either of discrete ancestral allele configurations as in the
coalescent, or of ancestral allele proportions as in the backward diffusion -
backward approaches need to be combined with corresponding forward ones. This
is done in so-called forward-backward algorithms. In this article, we utilize
orthogonal polynomials in forward-backward algorithms. They enable efficient
calculation of past allele configurations of an extant sample and probabilities
of ancestral population allele frequencies in equilibrium and in
non-equilibrium. We show that the genealogy of a sample is fully described by
the backward polynomial expansion of the marginal likelihood of its allele
configuration.


  In this study we investigate computationally tumour-oncolytic virus(OV)
interactions that take place within a heterogeneous ExtraCellular Matrix (ECM).
The ECM is viewed as a mixture of two constitutive phases, namely a fibre phase
and a non-fibre phase. The multiscale mathematical model presented here focuses
on the nonlocal cell-cell and cell-ECM interactions, and how these interactions
might be impacted by the infection of cancer cells with the OV. At macroscale
we track the kinetics of cancer cells, virus particles and the ECM. At
microscale we track (i) the degradation of ECM by matrix degrading enzymes
(MDEs) produced by cancer cells, which further influences the movement of
tumour boundary; (ii) the re-arrangement of the microfibres that influences the
re-arrangement of macrofibres (i.e., fibres at macroscale). With the help of
this new multiscale model, we investigate two questions: (i) whether the
infected cancer cell fluxes are the result of local or non-local advection in
response to ECM density; and (ii) what is the effect of ECM fibres on the the
spatial spread of oncolytic viruses and the outcome of oncolytic virotherapy.


  Detections of mutations of the SARS-CoV-2 virus gave rise to new packages of
interventions. Among them, international travel restrictions have been one of
the fastest and most visible responses to limit the spread of the variants.
While inducing large economic losses, the epidemiological consequences of such
travel restrictions are highly uncertain. They may be poorly effective when the
new highly transmissible strain of the virus already circulates in many
regions. Assessing the effectiveness of travel bans is difficult given the
paucity of data on daily cross-border mobility and on existing variant
circulation. The question is topical and timely as the new omicron variant --
classified as a variant of concern by WHO -- has been detected in Southern
Africa, and perceived as (potentially) more contagious than previous strains.
In this study, we develop a multi-country compartmental model of the SIR type.
We use it to simulate the spread of a new variant across European countries,
and to assess the effectiveness of unilateral and multilateral travel bans.


  Sequencing technologies are prone to errors, making error correction (EC)
necessary for downstream applications. EC tools need to be manually configured
for optimal performance. We find that the optimal parameters (e.g., k-mer size)
are both tool- and dataset-dependent. Moreover, evaluating the performance
(i.e., Alignment-rate or Gain) of a given tool usually relies on a reference
genome, but quality reference genomes are not always available. We introduce
Lerna for the automated configuration of k-mer-based EC tools. Lerna first
creates a language model (LM) of the uncorrected genomic reads; then,
calculates the perplexity metric to evaluate the corrected reads for different
parameter choices. Next, it finds the one that produces the highest alignment
rate without using a reference genome. The fundamental intuition of our
approach is that the perplexity metric is inversely correlated with the quality
of the assembly after error correction. Results: First, we show that the best
k-mer value can vary for different datasets, even for the same EC tool. Second,
we show the gains of our LM using its component attention-based transformers.
We show the model's estimation of the perplexity metric before and after error
correction. The lower the perplexity after correction, the better the k-mer
size. We also show that the alignment rate and assembly quality computed for
the corrected reads are strongly negatively correlated with the perplexity,
enabling the automated selection of k-mer values for better error correction,
and hence, improved assembly quality. Additionally, we show that our
attention-based models have significant runtime improvement for the entire
pipeline -- 18X faster than previous works, due to parallelizing the attention
mechanism and the use of JIT compilation for GPU inferencing.


  The emergence of new variants of SARS-CoV-2 is a major concern given their
potential impact on the transmissibility and pathogenicity of the virus as well
as the efficacy of therapeutic interventions. Here, we predict the mutability
of all positions in SARS-CoV-2 protein domains to forecast the appearance of
unseen variants. Using sequence data from other coronaviruses, pre-existing to
SARS-CoV-2, we build statistical models that do not only capture amino-acid
conservation but more complex patterns resulting from epistasis. We show that
these models are notably superior to conservation profiles in estimating the
already observable SARS-CoV-2 variability. In the receptor binding domain of
the spike protein, we observe that the predicted mutability correlates well
with experimental measures of protein stability and that both are reliable
mutability predictors (ROC AUC ~0.8). Most interestingly, we observe an
increasing agreement between our model and the observed variability as more
data become available over time, proving the anticipatory capacity of our
model. When combined with data concerning the immune response, our approach
identifies positions where current variants of concern are highly
overrepresented. These results could assist studies on viral evolution, future
viral outbreaks and, in particular, guide the exploration and anticipation of
potentially harmful future SARS-CoV-2 variants.


  Vaccination hesitancy is a major obstacle to achieving and maintaining herd
immunity. It is therefore of prime importance for public health authorities to
understand the dynamics of an anti-vaccine opinion in the population. We
introduce a novel mathematical model of opinion dynamics with spatial
reinforcement, which can generate echo chambers, i.e. opinion bubbles in which
information that is incompatible with one's entrenched worldview, is likely
disregarded. In a first mathematical part, we scale the model both to a
deterministic limit and to a weak-effects limit, and obtain bifurcations, phase
transitions, and the invariant measure. In a second part, we fit our model to
measles and meningococci vaccination coverage across 413 districts in Germany.
We reveal that strong echo chambers explain the occurrence and persistence of
the anti-vaccination opinion. We predict and compare the effectiveness of
different policies aimed at influencing opinion dynamics in order to increase
vaccination uptake.


  The design and implementation of regulation motifs ensuring robust perfect
adaptation are challenging problems in synthetic biology. Indeed, the design of
high-yield robust metabolic pathways producing, for instance, drug precursors
and biofuels, could be easily imagined to rely on such a control strategy in
order to optimize production levels and reduce production costs, despite the
presence of environmental disturbance and model uncertainty. We propose here a
motif that ensures tracking and robust perfect adaptation for the controlled
reaction network through integral feedback. Its metabolic load on the host is
fully tunable and can be made arbitrarily close to the constitutive limit, the
universal minimal metabolic load of all possible controllers. A DNA
implementation of the controller network is finally provided. Computer
simulations using realistic parameters demonstrate the good agreement between
the DNA implementation and the ideal controller dynamics.


  Classical (deterministic) and single-enzyme (stochastic) description of the
Michaelis-Menten (MM) kinetics, assume fast equilibration between enzyme and
complex, and identify detailed balance as a sufficient condition for the
hyperbolic substrate dependence of the MM equation (MME). Stochastic MM
kinetics based on the chemical master equation (CME), with no a priori
assumption of fast equilibration, however, unravels an observably long
non-classical transient regime at mesoscopic enzyme concentrations. The
enzymatic velocity in the transient regime is non-hyperbolic for product
turnovers below a critical time, but asymptotically recovers the hyperbolic MME
at long times. Here, we use this description to introduce a new kinetic
measure, the turnover number dependent fractional enzyme velocity. This measure
quantifies the degree of non-hyperbolicity in the non-classical transient
regime with respect to the hyperbolic MME. From this, we obtain a generalized
rate parameter condition for detailed balance in mesoscopic MM kinetics. This
condition, while subsuming the fast equilibrium approximation of the classical
MM kinetics, provides a strict lower bound on the magnitude of the catalytic
rate parameter. Further, from the condition of stationarity of the generating
function solution of the CME, we quantify the duration of the non-classical
regime. Our results show that the violation of detailed balance condition in
the transient regime is inextricably linked to the non-hyperbolic substrate
dependence of the enzymatic velocity. In the steady-state, when an effective
equilibrium between enzyme and complex is asymptotically established, the
condition of detailed balance emerges as a sufficient condition for the
hyperbolic substrate dependence of the MME.


  Spatial organization of chromatin plays a critical role in genome regulation.
Various types of affinity mediators and enzymes have been attributed to
regulate spatial organization of chromatin from a thermodynamics perspective.
However, at the mechanistic level, enzymes act in their unique ways. Here, we
construct a polymer physics model following the mechanistic scheme of
Topoisomerase-II, an enzyme resolving topological constraints of chromatin, and
investigate its role on interphase chromatin organization. Our computer
simulations demonstrate Topoisomerase-II's ability to phase separate chromatin
into eu- and heterochromatic regions with a characteristic wall-like
organization of the euchromatic regions. Exploiting a mean-field framework, we
argue that the ability of the euchromatic regions crossing each other due to
enzymatic activity of Topoisomerase-II induces this phase separation. Motivated
from a recent experimental observation on different structural states of the
eu- and the heterochromatic units, we further extend our model to a bidisperse
setting and show that the characteristic features of the enzymatic activity
driven phase separation survives there. The existence of these characteristic
features, even under the non-localized action of the enzyme, highlights the
critical role of enzymatic activity in chromatin organization, and points out
the importance of further experiments along this line.


  Understanding control mechanisms in biological systems plays a crucial role
in important applications, for instance in cell reprogramming. Boolean modeling
allows the identification of possible efficient strategies, helping to reduce
the usually high and time-consuming experimental efforts. Available approaches
to control strategy identification usually focus either on attractor or
phenotype control, and are unable to deal with more complex control problems,
for instance phenotype avoidance. They also fail to capture, in many
situations, all possible minimal strategies, finding instead only sub-optimal
solutions. In order to fill these gaps, we present a novel approach to control
strategy identification in Boolean networks based on model checking. The method
is guaranteed to identify all minimal control strategies, and provides maximal
flexibility in the definition of the control target. We investigate the
applicability of the approach by considering a range of control problems for
different biological systems, comparing the results, where possible, to those
obtained by alternative control methods.


  Various regimes of a charge motion along a chain in a constant electric field
are investigated. This motion is simulated on the basis of the Holstein model.
Earlier studies demonstrate a possibility of a uniform motion of a charge in a
constant electric field over very long distances. For small values of the
electric field intensity a Holstein polaron can move at a constant velocity. As
the electric field intensity increases, a charge motion acquires oscillatorily
character, performing Bloch oscillations. Since the charge motion depends on
the whole set of the system parameters the character of the motion depends not
only on the value of the electric field intensity. Therefore, the electric
field intensity for which the uniform motion takes place differs for chains
with different parameters. The character of the charge motion and distribution
is considered in chains with different values of the constant of coupling
between the charge and the displacements of the chain. It is shown that the
values of the electric field intensity for which the regime of a charge motion
changes are different in chains with different values of the coupling constant.
It is also demonstrated that for one and the same value of the electric field
intensity, in chains with different values of the coupling constant either a
uniform motion or an oscillatory motion, or a stationary polaron can be
observed.


  The COVID-19 pandemic has severely impacted health systems and economies
worldwide. Significant global efforts are therefore ongoing to improve vaccine
efficacies, optimize vaccine deployment, and develop new antiviral therapies to
combat the pandemic. Mechanistic viral dynamics and quantitative systems
pharmacology models of SARS-CoV-2 infection, vaccines, immunomodulatory agents,
and antiviral therapeutics have played a key role in advancing our
understanding of SARS-CoV-2 pathogenesis and transmission, the interplay
between innate and adaptive immunity to influence the outcomes of infection,
effectiveness of treatments, mechanisms and performance of COVID-19 vaccines,
and the impact of emerging SARS-CoV-2 variants. Here, we review some of the
critical insights provided by these models and discuss the challenges ahead.


  The simplest possible informational heteropolymer requires only a two-letter
alphabet to be able to store information. The evolutionary choice of four
monomers in the informational biomolecules RNA/DNA or their progenitors is
intriguing, given the inherent difficulties in the simultaneous and localized
prebiotic synthesis of all four monomers of progenitors of DNA from common
precursors on early Earth. Excluding the scenario where a two-letter alphabet
genome eventually expanded to include two more letters to code for more amino
acids on teleological grounds, we show here that a heteropolymer sequence in
the RNA-world-like scenario would have had to be composed of at least four
letters in order to predictably fold into a specific secondary structure, and
hence must have outcompeted the two-letter alphabet genomes. Using a model that
we previously used to demonstrate the evolutionary advantages of unidirectional
replication and anti-parallel strand orientation of duplex DNA, we show here
that the competing constraints of maximum replicative potential and predictable
secondary structure formation can be simultaneously satisfied only by
palindromic heteropolymer sequences composed of a minimum of four letters,
within the premise of the presence of sequence-dependent asymmetric
cooperativity in these RNA/DNA progenitors.


  We expose a new excitation transport channel in FMO, consisting of the 7
intramonomeric pigments and their 2nd nearest intermonomeric counterpart. Such
a channel outcompetes the standard alternative, where intramonomeric pigments
partner with their closest intermonomeric pigment. The efficient performance of
the new channel lies at the interplay between incoherent energy transfer and
the capacity of quantum coherence to build optimal acceptor levels. Finally,
the differential fitness and peculiarity of the new over the alternative
channel highlights its possible role within natural selection.


  Phylogenetic diversity indices such as the Fair Proportion (FP) index are
frequently discussed as prioritization criteria in biodiversity conservation.
They rank species according to their contribution to overall diversity by
taking into account the unique and shared evolutionary history of each species
as indicated by its placement in an underlying phylogenetic tree.
Traditionally, phylogenetic trees were inferred from single genes and the
resulting gene trees were assumed to be a valid estimate for the species tree,
i.e., the "true" evolutionary history of the species under consideration.
However, nowadays it is common to sequence whole genomes of hundreds or
thousands of genes, and it is often the case that conflicting genealogical
histories exist in different genes throughout the genome, resulting in
discordance between individual gene trees and the species tree. Here, we
analyze the effects of gene and species tree discordance on prioritization
decisions based on the FP index. In particular, we consider the ranking order
of taxa induced by (i) the FP index on a species tree, and (ii) the expected FP
index across all gene tree histories associated with the species tree. On one
hand, we show that for particular tree shapes, the two rankings always
coincide. On the other hand, we show that for all leaf numbers greater than or
equal to five, there exist species trees for which the two rankings differ.
Finally, we illustrate the variability in the rankings obtained from the FP
index across different gene tree and species tree estimates for an empirical
multilocus mammal data set.


  We propose a decision-making model for joint intentionality by interpreting
it as group-mindedness at the microlevel. We apply this model to give a formal
justification of the first part of the Interdependence Hypothesis due to
Tomasello et al. [Current Anthropology, 2012] which asserts that the emergence
of joint intentionality evolved due to the challenges of difficult
collaborative foraging practices among early humans, and that its evolution led
to robust collaboration and some form of altruism.
  In another application of the microlevel group-mindedness we consider the
problem of establishing cooperation in high-risk-of-defection strategic
conflicts and we show that the emergence of cooperation in such situations can
be explained in the context of cultural group selection as the result of
adaptive learning.


  Motivation: We consider continuous-time Markov chains that describe the
stochastic evolution of a dynamical system by a transition-rate matrix $Q$
which depends on a parameter $\theta$. Computing the probability distribution
over states at time $t$ requires the matrix exponential $\exp(tQ)$, and
inferring $\theta$ from data requires its derivative
$\partial\exp\!(tQ)/\partial\theta$. Both are challenging to compute when the
state space and hence the size of $Q$ is huge. This can happen when the state
space consists of all combinations of the values of several interacting
discrete variables. Often it is even impossible to store $Q$. However, when $Q$
can be written as a sum of tensor products, computing $\exp(tQ)$ becomes
feasible by the uniformization method, which does not require explicit storage
of $Q$.
  Results: Here we provide an analogous algorithm for computing
$\partial\exp\!(tQ)/\partial\theta$, the differentiated uniformization method.
We demonstrate our algorithm for the stochastic SIR model of epidemic spread,
for which we show that $Q$ can be written as a sum of tensor products. We
estimate monthly infection and recovery rates during the first wave of the
COVID-19 pandemic in Austria and quantify their uncertainty in a full Bayesian
analysis.
  Availability: Implementation and data are available at
https://github.com/spang-lab/TenSIR.


  Understanding whether a population will survive and flourish or become
extinct is a central question in population biology. One way of exploring this
question is to study population dynamics using reaction-diffusion equations,
where migration is usually represented as a linear diffusion term, and
birth-death is represented with a bistable source term. While linear diffusion
is most commonly employed to study migration, there are several limitations of
this approach, such as the inability of linear diffusion-based models to
predict a well-defined population front. One way to overcome this is to
generalise the constant diffusivity, $D$, to a nonlinear diffusivity function
$D(C)$, where $C>0$ is the density. While it has been formally established that
the choice of $D(C)$ affects long-term survival or extinction of a bistable
population, working solely in a classical continuum framework makes it
difficult to understand precisely how the choice of $D(C)$ affects survival or
extinction. Here, we address this question by working with a simple discrete
simulation model that is easy to interpret. The continuum limit of the discrete
model is a nonlinear reaction-diffusion equation, where the flux involves a
nonlinear diffusion term and the source term is given by the strong Allee
effect bistable model. We study population extinction/survival using this very
intuitive discrete framework together with numerical solutions of the
reaction-diffusion continuum limit equation. This approach provides clear
insight into how the choice of $D(C)$ either encourages or suppresses
population extinction relative to the classical linear diffusion model.


  Invasive pest establishment is a pervasive threat to global ecosystems,
agriculture, and public health. The recent establishment of the invasive
spotted lanternfly in the northeastern United States has proven devastating to
farms and vineyards, necessitating urgent development of population dynamical
models and effective control practices. In this paper, we propose a stage- and
age-structured system of PDEs to model insect pest populations, in which
underlying dynamics are dictated by ambient temperature through rates of
development, fecundity, and mortality. The model incorporates diapause and
non-diapause pathways, and is calibrated to experimental and field data on the
spotted lanternfly. We develop a novel moving mesh method for capturing
age-advection accurately, even for coarse discretization parameters. We define
a one-year reproductive number ($R_0$) from the spectrum of a one-year solution
operator, and study its sensitivity to variations in the mean and amplitude of
the annual temperature profile. We quantify assumptions sufficient to give rise
to the low-rank structure of the solution operator characteristic of part of
the parameter domain. We discuss establishment potential as it results from the
pairing of a favorable $R_0$ value and transient population survival, and
address implications for pest control strategies.


  Many mechanisms behind the evolution of cooperation, such as reciprocity,
indirect reciprocity, and altruistic punishment, require group knowledge of
individual actions. But what keeps people cooperating when no one is looking?
Conformist norm internalization, the tendency to abide by the behavior of the
majority of the group, even when it is individually harmful, could be the
answer. In this paper, we analyze a world where (1) there is group selection
and punishment by indirect reciprocity but (2) many actions (half) go
unobserved, and therefore unpunished. Can norm internalization fill this
"observation gap" and lead to high levels of cooperation, even when agents may
in principle cooperate only when likely to be caught and punished?
Specifically, we seek to understand whether adding norm internalization to the
strategy space in a public goods game can lead to higher levels of cooperation
when both norm internalization and cooperation start out rare. We found the
answer to be positive, but, interestingly, not because norm internalizers end
up making up a substantial fraction of the population, nor because they
cooperate much more than other agent types. Instead, norm internalizers, by
polarizing, catalyzing, and stabilizing cooperation, can increase levels of
cooperation of other agent types, while only making up a minority of the
population themselves.


  Mixed communities of organisms are found in many environments (from the human
gut to marine ecosystems) and can have profound impact on human health and the
environment. Metagenomics studies the genomic material of such communities
through high-throughput sequencing that yields DNA subsequences for subsequent
analysis. A fundamental problem in the standard workflow, called binning, is to
discover clusters, of genomic subsequences, associated with the unknown
constituent organisms. Inherent noise in the subsequences, various biological
constraints that need to be imposed on them and the skewed cluster size
distribution exacerbate the difficulty of this unsupervised learning problem.
In this paper, we present a new formulation using a graph where the nodes are
subsequences and edges represent homophily information. In addition, we model
biological constraints providing heterophilous signal about nodes that cannot
be clustered together. We solve the binning problem by developing new
algorithms for (i) graph representation learning that preserves both homophily
relations and heterophily constraints (ii) constraint-based graph clustering
method that addresses the problems of skewed cluster size distribution.
Extensive experiments, on real and synthetic datasets, demonstrate that our
approach, called RepBin, outperforms a wide variety of competing methods. Our
constraint-based graph representation learning and clustering methods, that may
be useful in other domains as well, advance the state-of-the-art in both
metagenomics binning and graph representation learning.


  Several epidemiological models have been proposed to study the evolution of
COVID-19 pandemic. In this paper, we propose an extension of the SUIHTER model,
first introduced in [Parolini et al, Proc R. Soc. A., 2021] to analyse the
COVID-19 spreading in Italy, which accounts for the vaccination campaign and
the presence of new variants when they become dominant. In particular, the
specific features of the variants (e.g. their increased transmission rate) and
vaccines (e.g. their efficacy to prevent transmission, hospitalization and
death) are modeled, based on clinical evidence. The new model is validated
comparing its near-future forecast capabilities with other epidemiological
models and exploring different scenario analyses.


  No systematic method exists to derive intra-chromosomal potentials between
nucleosomes along a chromosome consistently across a given genome. Such
potentials can yield information on nucleosomal ordering, thermal as well as
mechanical properties of chromosomes. Thus, indirectly, they shed light on a
possible mechanical genomic code along a chromosome. To develop a method
yielding effective intra-chromosomal potentials between nucleosomes a
generalized Lennard-Jones potential for the parameterization is developed based
on nucleosomal positioning data. This approach eliminates some of the problems
that the underlying nucleosomal positioning data has, rendering the extraction
difficult on the individual nucleosomal level. Furthermore, patterns on which
to base a classification along a chromosome appear on larger domains, such as
hetero- and euchromatin. An intuitive selection strategy for the
noisy-optimization problem is employed to derive effective exponents for the
generalized potential. The method is tested on the Candida albicans genome.
Applying k-means clustering based on potential parameters and thermodynamic
compressibilities, a genome-wide clustering of nucleosome sequences is obtained
for Candida albicans. This clustering shows that a chromosome beyond the
classical dichotomic categories of hetero- and euchromatin, is more
feature-rich.


  We study the formation of trails in populations of self-propelled agents that
make oriented deposits of pheromones and also sense such deposits to which they
then respond with gradual changes of their direction of motion. Based on
extensive off-lattice computer simulations aiming at the scale of insects,
e.g., ants, we identify a number of emerging stationary patterns and obtain
qualitatively the non-equilibrium state diagram of the model, spanned by the
strength of the agent--pheromone interaction and the number density of the
population. In particular, we demonstrate the spontaneous formation of
persistent, macroscopic trails, and highlight some behaviour that is consistent
with a dynamic phase transition. This includes a characterisation of the mass
of system-spanning trails as a potential order parameter. We also propose a
dynamic model for a few macroscopic observables, including the sub-population
size of trail-following agents, which captures the early phase of trail
formation.