  Nutrient-Phytoplankton-Zooplankton-Detritus (NPZD) models, describing the
interactions between phytoplankton, zooplankton systems, and their ecosystem,
are used to predict their ecological and evolutionary population dynamics.
These organisms form the base two trophic levels of aquatic ecosystems. Hence
understanding their population dynamics and how disturbances can affect these
systems is crucial. Here, starting from a base NPZ modeling framework, we
incorporate the harmful effects of phytoplankton overpopulation on zooplankton
- representing a crucial next step in harmful algal bloom (HAB) modeling - and
split the nutrient compartment to formulate an NPZD model. We then
mathematically analyze the NPZ system upon which this new model is based,
including local and global stability of equilibria, Hopf bifurcation condition,
and forward hysteresis, where the bi-stability occurs with multiple attractors.
Finally, we extend the threshold analysis to the NPZD model, which displays
both forward hysteresis with bi-stability and Hopf bifurcation under different
parameter regimes, and examine ecological implications after incorporating
seasonality and ecological disturbances. Ultimately, we quantify ecosystem
health in terms of the relative values of the robust persistence thresholds for
phytoplankton and zooplankton and find (i) ecosystems sufficiently favoring
phytoplankton, as quantified by the relative values of the plankton persistence
numbers, are vulnerable to both HABs and (local) zooplankton extinction (ii)
even healthy ecosystems are extremely sensitive to nutrient depletion over
relatively short time scales.


  Self-supervised learning has recently gained growing interest in molecular
modeling for scientific tasks such as AI-assisted drug discovery. Current
studies consider leveraging both 2D and 3D molecular structures for
representation learning. However, relying on straightforward alignment
strategies that treat each modality separately, these methods fail to exploit
the intrinsic correlation between 2D and 3D representations that reflect the
underlying structural characteristics of molecules, and only perform
coarse-grained molecule-level alignment. To derive fine-grained alignment and
promote structural molecule understanding, we introduce an atomic-relation
level "blend-then-predict" self-supervised learning approach, MoleBLEND, which
first blends atom relations represented by different modalities into one
unified relation matrix for joint encoding, then recovers modality-specific
information for 2D and 3D structures individually. By treating atom
relationships as anchors, MoleBLEND organically aligns and integrates visually
dissimilar 2D and 3D modalities of the same molecule at fine-grained atomic
level, painting a more comprehensive depiction of each molecule. Extensive
experiments show that MoleBLEND achieves state-of-the-art performance across
major 2D/3D molecular benchmarks. We further provide theoretical insights from
the perspective of mutual-information maximization, demonstrating that our
method unifies contrastive, generative (cross-modality prediction) and
mask-then-predict (single-modality prediction) objectives into one single
cohesive framework.


  Adenosine-5'-triphosphate (ATP) plays a crucial role in many biocatalytic
reactions and its regeneration can influence the performance of a related
enzymatic reaction significantly. Here, we established an electrochemically
coupled ATP regeneration by pyruvate oxidase and acetate kinase (ACK) for the
phosphorylation of mevalonate catalyzed by mevalonate kinase. A yield of 84%
for the product mevalonate phosphate was reached and a total turnover number
for ADP of 68. These metrics are promising for the development of an economic
feasible bioprocess and surpass many other enzymatic ATP regeneration systems.
A comparison was made to polyphosphate kinases (PPKs), ACK, pyruvate kinase,
and creatine kinase in terms of the phosphate donor properties and biocatalytic
metrics of exemplary reactions. Furthermore, our system was expanded by a PPK
that enables the phosphorylation of AMP, which can broaden the spectrum of
applications even further.


  Compared with mandatory quarantine, voluntary quarantine offers individuals
the liberty to decide on whether to isolate themselves or not in case of
infection exposure, driven by their own sense of responsibility and concern for
public health. To tackle this problem, we perform agent-based simulations to
study the coupled dynamics of disease spreading and the adaptive behavior of
quarantine. We show that voluntary quarantine can be an effective intervention
measure to mitigate the transmission of the disease. Moreover, we discuss the
roles played by the level of temptation to refrain from quarantine and the
degree of social compassion in reducing the risk of spread and controlling the
outbreak. We find that a low level of temptation or a high degree of social
compassion not only can lead to the amelioration of the disease impact to some
extent but also can result in a complete containment of the epidemic. By
shedding light on the intricate interplay between disease spreading and
voluntary quarantine, our results provide new insights into the fundamental
social dilemma aspect of disease control through nonpharmaceutical
interventions including voluntary quarantine.


  Microbial populations generally evolve in volatile environments, under
conditions fluctuating between harsh and mild, e.g. as the result of sudden
changes in toxin concentration or nutrient abundance. Environmental variability
thus shapes the long-time population dynamics, notably by influencing the
ability of different strains of microorganisms to coexist. Inspired by the
evolution of antimicrobial resistance, we study the dynamics of a community
consisting of two competing strains subject to twofold environmental
variability. The level of toxin varies in time, favouring the growth of one
strain under low drug concentration and the other strain when the toxin level
is high. We also model time-changing resource abundance by a randomly switching
carrying capacity that drives the fluctuating size of the community. While one
strain dominates in a static environment, we show that species coexistence is
possible in the presence of environmental variability. By computational and
analytical means, we determine the environmental conditions under which
long-lived coexistence is possible and when it is almost certain. Notably, we
study the circumstances under which environmental and demographic fluctuations
promote, or hinder, the strains coexistence. We also determine how the make-up
of the coexistence phase and the average abundance of each strain depend on the
environmental variability.


  There is a pressing need to better understand how microbial populations
respond to antimicrobial drugs, and to find mechanisms to possibly eradicate
antimicrobial-resistant cells. The inactivation of antimicrobials by resistant
microbes can often be viewed as a cooperative behavior leading to the
coexistence of resistant and sensitive cells in large populations and static
environments. This picture is however greatly altered by the fluctuations
arising in volatile environments, in which microbial communities commonly
evolve. Here, we study the eco-evolutionary dynamics of a population consisting
of an antimicrobial resistant strain and microbes sensitive to antimicrobial
drugs in a time-fluctuating environment, modeled by a carrying capacity
randomly switching between states of abundance and scarcity. We assume that
antimicrobial resistance is a shared public good when the number of resistant
cells exceeds a certain threshold. Eco-evolutionary dynamics is thus
characterised by demographic noise (birth and death events) coupled to
environmental fluctuations which can cause population bottlenecks. By combining
analytical and computational means, we determine the environmental conditions
for the long-lived coexistence and fixation of both strains, and characterise a
fluctuation-driven antimicrobial resistance eradication mechanism, where
resistant microbes experience bottlenecks leading to extinction. We also
discuss the possible applications of our findings to laboratory-controlled
experiments.


  New automated and high-throughput methods allow the manipulation and
selection of numerous bacterial populations. In this manuscript we are
interested in the neutral diversity patterns that emerge from such a setup in
which many bacterial populations are grown in parallel serial transfers, in
some cases with population-wide extinction and splitting events. We model
bacterial growth by a birth-death process and use the theory of coalescent
point processes. We show that there is a dilution factor that optimises the
expected amount of neutral diversity for a given amount of cycles, and study
the power law behaviour of the mutation frequency spectrum for different
experimental regimes. We also explore how neutral variation diverges between
two recently split populations by establishing a new formula for the expected
number of shared and private mutations. Finally, we show the interest of such a
setup to select a phenotype of interest that requires multiple mutations.


  Symmetry arguments are frequently used -- often implicitly -- in mathematical
modeling of natural selection. Symmetry simplifies the analysis of models and
reduces the number of distinct population states to be considered. Here, I
introduce a formal definition of symmetry in mathematical models of natural
selection. This definition applies to a broad class of models that satisfy a
minimal set of assumptions, using a framework developed in previous works. In
this framework, population structure is represented by a set of sites at which
alleles can live, and transitions occur via replacement of some alleles by
copies of others. A symmetry is defined as a permutation of sites that
preserves probabilities of replacement and mutation. The symmetries of a given
selection process form a group, which acts on population states in a way that
preserves the Markov chain representing selection. Applying classical results
on group actions, I formally characterize the use of symmetry to reduce the
states of this Markov chain, and obtain bounds on the number of states in the
reduced chain.


  In this research, we revisit the paper by Pal et al. [Int. J. Appl. Comput.
Math (2017) 3:1833-1845] and comment on the claim that global stability of the
interior equilibrium point depends on some key parameters. This is not true,
and we have provided detailed proof to this effect. The considered model is a
modified Lotka-Volterra model where square root functional response is
involved. The square root functional response is non-differentiable at the
origin and hence cannot be studied with standard local stability tools.
Furthermore, our numerical simulations indicate that we can classify the phase
portrait into two modes of behavior where some positive initial conditions
converge towards the predator axis in finite time.


  In this study, we address the challenge of using energy-based models to
produce high-quality, label-specific data in complex structured datasets, such
as population genetics, RNA or protein sequences data. Traditional training
methods encounter difficulties due to inefficient Markov chain Monte Carlo
mixing, which affects the diversity of synthetic data and increases generation
times. To address these issues, we use a novel training algorithm that exploits
non-equilibrium effects. This approach, applied on the Restricted Boltzmann
Machine, improves the model's ability to correctly classify samples and
generate high-quality synthetic data in only a few sampling steps. The
effectiveness of this method is demonstrated by its successful application to
four different types of data: handwritten digits, mutations of human genomes
classified by continental origin, functionally characterized sequences of an
enzyme protein family, and homologous RNA sequences from specific taxonomies.


  We consider compartmental models of communicable disease with uncertain
contact rates. Stochastic fluctuations are often added to the contact rate to
account for uncertainties. White noise, which is the typical choice for the
fluctuations, leads to significant underestimation of the disease severity.
Here, starting from reasonable assumptions on the social behavior of
individuals, we model the contacts as a Markov process which takes into account
the temporal correlations present in human social activities. Consequently, we
show that the mean-reverting Ornstein-Uhlenbeck (OU) process is the correct
model for the stochastic contact rate. We demonstrate the implication of our
model on two examples: a Susceptibles-Infected-Susceptibles (SIS) model and a
Susceptibles-Exposed-Infected-Removed (SEIR) model of the COVID-19 pandemic. In
particular, we observe that both compartmental models with white noise
uncertainties undergo transitions that lead to the systematic underestimation
of the spread of the disease. In contrast, modeling the contact rate with the
OU process significantly hinders such unrealistic noise-induced transitions.
For the SIS model, we derive its stationary probability density analytically,
for both white and correlated noise. This allows us to give a complete
description of the model's asymptotic behavior as a function of its bifurcation
parameters, i.e., the basic reproduction number, noise intensity, and
correlation time. For the SEIR model, where the probability density is not
available in closed form, we study the transitions using Monte Carlo
simulations. Our study underscores the necessity of temporal correlations in
stochastic compartmental models and the need for more empirical studies that
would systematically quantify such correlations.


  Over recent years, super-resolution and single-molecule imaging methods have
delivered unprecedented details on the nanoscale organization and dynamics of
individual molecules in different contexts. Yet, visualizing single-molecule
processes in living cells with the required spatial and temporal resolution
remains highly challenging. Here, we report on an analytical approach that
extracts such information from live-cell single-molecule imaging at
high-labeling densities using standard fluorescence probes. Our
high-density-mapping (HiDenMap) methodology provides single-molecule nanometric
localization accuracy together with millisecond temporal resolution over
extended observation times, delivering multi-scale spatiotemporal data that
report on the interaction of individual molecules with their dynamic
environment. We validated HiDenMaps by simulations of Brownian trajectories in
the presence of patterns that restrict free diffusion with different
probabilities. We further generated and analyzed HiDenMaps from single-molecule
images of transmembrane proteins having different interaction strengths to
cortical actin, including the transmembrane receptor CD44. HiDenMaps uncovered
a highly heterogenous and multi-scale spatiotemporal organization for all the
proteins that interact with the actin cytoskeleton. Notably, CD44 alternated
between periods of random diffusion and transient trapping, likely resulting
from actin-dependent CD44 nanoclustering. Whereas receptor trapping was dynamic
and lasted for hundreds of milliseconds, actin remodeling occurred at the
timescale of tens of seconds, coordinating the assembly and disassembly of CD44
nanoclusters rich regions. Together, our data demonstrate the power of
HiDenMaps to explore how individual molecules interact with and are organized
by their environment in a dynamic fashion.


  Overcoming chemoresistance is a challenge for multiple chemotherapeutics
agents like cisplatin. ABC transporters such as MDR1 or MRPs and PI3K/AKT
pathway have been proposed as actors of chemoresistance in several cancers. In
this review we describe two downstream targets of Akt: c-myc and p53 in the
chemoresistance. We suggest a potential link between p53, c-myc and ABC
transporters expression. Consequently a link between Akt and ABC
transporters-mediated chemoresistance may exist. Finally we show that Akt
activation may be Orai-dependent and/or TRPC-dependent, suggesting that these
ion channels could constitute a therapeutic target in cancer.


  Perfect adaptation is a phenomenon whereby the output variables of a system
can maintain certain values despite external disturbances. Robust perfect
adaptation (RPA) refers to an adaptation property that does not require
fine-tuning of system parameters. RPA plays a vital role for the survival of
living systems in unpredictable environments. However, complex interaction
patterns in biochemical systems pose a significant challenge in identifying RPA
and associated regulatory mechanisms. The goal of this paper is to present a
novel approach for finding all RPA properties that are realized for a generic
choice of kinetics for general deterministic chemical reaction systems. This is
accomplished by proving that an RPA property is represented by a subnetwork
with certain topological features. This connection is exploited to show that
these structures generate all kinetics-independent RPA properties, allowing us
to systematically identify all RPA properties by enumerating these subnetworks.
An efficient method is developed for this enumeration, and we provide a
computational package for this purpose. We pinpoint the integral feedback
controllers that work in concert to realize each RPA property, casting our
results into the familiar control-theoretic paradigm of the Internal Model
Principle. We further generalize the regulation problem to the multi-output
scenario where the target values belong to a manifold of nonzero dimension, and
provide a sufficient condition for this. The present work significantly
advances our understanding of regulatory mechanisms that lead to RPA in
endogenous biochemical systems, and it also provides rational design principles
for synthetic controllers. The present results indicate that an RPA property is
essentially equivalent to the existence of a "topological invariant", which is
an instance of what we call the "Robust Adaptation is Topological"(RAT)
principle.


  Biological networks often encapsulate promotion/inhibition as signed
edge-weights of a graph. Nodes may correspond to genes assigned expression
levels (mass) of respective proteins. The promotion/inhibition nature of
co-expression between nodes is encoded in the sign of the corresponding entry
of a sign-indefinite adjacency matrix, though the strength of such
co-expression (i.e., the precise value of edge weights) cannot typically be
directly measured. Herein we address the inverse problem to determine network
edge-weights based on a sign-indefinite adjacency and expression levels at the
nodes. While our motivation originates in gene networks, the framework applies
to networks where promotion/inhibition dictates a stationary mass distribution
at the nodes. In order to identify suitable edge-weights we adopt a framework
of ``negative probabilities,'' advocated by P.\ Dirac and R.\ Feynman, and we
set up a likelihood formalism to obtain values for the sought edge-weights. The
proposed optimization problem can be solved via a generalization of the
well-known Sinkhorn algorithm; in our setting the Sinkhorn-type ``diagonal
scalings'' are multiplicative or inverse-multiplicative, depending on the sign
of the respective entries in the adjacency matrix, with value computed as the
positive root of a quadratic polynomial.


  Chemical formula annotation for tandem mass spectrometry (MS/MS) data is the
first step toward structurally elucidating unknown metabolites. While great
strides have been made toward solving this problem, the current
state-of-the-art method depends on time-intensive, proprietary, and
expert-parameterized fragmentation tree construction and scoring. In this work
we extend our previous spectrum Transformer methodology into an energy based
modeling framework, MIST-CF, for learning to rank chemical formula and adduct
assignments given an unannotated MS/MS spectrum. Importantly, MIST-CF learns in
a data dependent fashion using a Formula Transformer neural network
architecture and circumvents the need for fragmentation tree construction. We
train and evaluate our model on a large open-access database, showing an
absolute improvement of 10% top 1 accuracy over other neural network
architectures. We further validate our approach on the CASMI2022 challenge
dataset, achieving nearly equivalent performance to the winning entry within
the positive mode category without any manual curation or post-processing of
our results. These results demonstrate an exciting strategy to more powerfully
leverage MS2 fragment peaks for predicting MS1 precursor chemical formula with
data driven learning.


  In this paper we investigate a structured population model with distributed
delay. Our model incorporates two different types of nonlinearities.
Specifically we assume that individual growth and mortality are affected by
scramble competition, while fertility is affected by contest competition. In
particular, we assume that there is a hierarchical structure in the population,
which affects mating success. The dynamical behavior of the model is analysed
via linearisation by means of semigroup and spectral methods. In particular, we
introduce a reproduction function and use it to derive linear stability
criteria for our model. Further we present numerical simulations to underpin
the stability results we obtained.


  Prokaryotic organisms usually possess compact genomes, which are particularly
suitable to complete sequencing with existing technologies, which led to an
escalating accumulation of available genome data. In response to this
ever-expanding repository of information, we introduce ProSt, a computational
system designed for the batch computation, storage, and interactive
visualization of the values of attributes of prokaryotic genomes. The system
allows for parallel attribute value batch computation, dynamically designed to
incrementally integrate new attribute values as additional genomes become
available.
  ProSt is flexible permitting the definition of attributes by implementing
attribute value computation plugins, supporting several languages (Python, Nim,
Rust and Bash). This allows the system to continually evolve in accordance with
changing research needs and developments. Additionally, our computation and
storage systems maintain comprehensive metadata, thereby enabling data
provenance tracking for the computed attribute values.


  We adopt a maximum-likelihood framework to estimate parameters of a
stochastic susceptible-infected-recovered (SIR) model with contact tracing on a
rooted random tree. Given the number of detectees per index case, our estimator
allows to determine the degree distribution of the random tree as well as the
tracing probability. Since we do not discover all infectees via contact
tracing, this estimation is non-trivial. To keep things simple and stable, we
develop an approximation suited for realistic situations (contract tracing
probability small, or the probability for the detection of index cases small).
In this approximation, the only epidemiological parameter entering the
estimator is $R_0$.
  The estimator is tested in a simulation study and is furthermore applied to
covid-19 contact tracing data from India. The simulation study underlines the
efficiency of the method. For the empirical covid-19 data, we compare different
degree distributions and perform a sensitivity analysis. We find that
particularly a power-law and a negative binomial degree distribution fit the
data well and that the tracing probability is rather large. The sensitivity
analysis shows no strong dependency of the estimates on the reproduction
number. Finally, we discuss the relevance of our findings.


  In the long run, the eventual extinction of any biological population is an
inevitable outcome. While extensive research has focused on the average time it
takes for a population to go extinct under various circumstances, there has
been limited exploration of the distributions of extinction times and the
likelihood of significant fluctuations. Recently, Hathcock and Strogatz [PRL
128, 218301 (2022)] identified Gumbel statistics as a universal asymptotic
distribution for extinction-prone dynamics in a stable environment. In this
study, we aim to provide a comprehensive survey of this problem by examining a
range of plausible scenarios, including extinction-prone, marginal (neutral),
and stable dynamics. We consider the influence of demographic stochasticity,
which arises from the inherent randomness of the birth-death process, as well
as cases where stochasticity originates from the more pronounced effect of
random environmental variations. Our work proposes several generic criteria
that can be used for the classification of experimental and empirical systems,
thereby enhancing our ability to discern the mechanisms governing extinction
dynamics. By employing these criteria, we can improve our understanding of the
underlying mechanisms driving extinction processes.


  Advances in bioengineering and biomedicine demand a deep understanding of the
dynamic behavior of biological systems, ranging from protein pathways to
complex cellular processes. Biological networks like gene regulatory networks
and protein pathways are key drivers of embryogenesis and physiological
processes. Comprehending their diverse behaviors is essential for tackling
diseases, including cancer, as well as for engineering novel biological
constructs. Despite the availability of extensive mathematical models
represented in Systems Biology Markup Language (SBML), researchers face
significant challenges in exploring the full spectrum of behaviors and
optimizing interventions to efficiently shape those behaviors. Existing tools
designed for simulation of biological network models are not tailored to
facilitate interventions on network dynamics nor to facilitate automated
discovery. Leveraging recent developments in machine learning (ML), this paper
introduces SBMLtoODEjax, a lightweight library designed to seamlessly integrate
SBML models with ML-supported pipelines, powered by JAX. SBMLtoODEjax
facilitates the reuse and customization of SBML-based models, harnessing JAX's
capabilities for efficient parallel simulations and optimization, with the aim
to accelerate research in biological network analysis.


  Gene regulatory networks are powerful tools for modeling interactions among
genes to regulate their expression for homeostasis and differentiation.
Single-cell sequencing offers a unique opportunity to build these networks with
high-resolution data. There are many proposed computational methods to build
these networks using single-cell data and different approaches are followed to
benchmark these methods. In this review, we lay the basic terminology in the
field and define the success metrics. Next, we present an overview of
approaches for benchmarking computational gene regulatory network approaches
for building gene regulatory networks and point out gaps and future directions
in this regard.


  The concept of sustainable intensification in agriculture necessitates the
implementation of management practices that prioritize sustainability without
compromising productivity. However, the effects of such practices are known to
depend on environmental conditions, and are therefore expected to change as a
result of a changing climate. We study the impact of crop diversification on
productivity in the context of climate change. We leverage heterogeneous Earth
Observation data and contribute a data-driven approach based on causal machine
learning for understanding how crop diversification impacts may change in the
future. We apply this method to the country of Cyprus throughout a 4-year
period. We find that, on average, crop diversification significantly benefited
the net primary productivity of crops, increasing it by 2.8%. The effect
generally synergized well with higher maximum temperatures and lower soil
moistures. In a warmer and more drought-prone climate, we conclude that crop
diversification exhibits promising adaptation potential and is thus a sensible
policy choice with regards to agricultural productivity for present and future.


  Phylogenetic networks play an important role in evolutionary biology as,
other than phylogenetic trees, they can be used to accommodate reticulate
evolutionary events such as horizontal gene transfer and hybridization. Recent
research has provided a lot of progress concerning the reconstruction of such
networks from data as well as insight into their graph theoretical properties.
However, methods and tools to quantify structural properties of networks or
differences between them are still very limited. For example, for phylogenetic
trees, it is common to use balance indices to draw conclusions concerning the
underlying evolutionary model, and more than twenty such indices have been
proposed and are used for different purposes. One of the most frequently used
balance index for trees is the so-called total cophenetic index, which has
several mathematically and biologically desirable properties. For networks, on
the other hand, balance indices are to-date still scarce.
  In this contribution, we introduce the \textit{weighted} total cophenetic
index as a generalization of the total cophenetic index for trees to make it
applicable to general phylogenetic networks. As we shall see, this index can be
determined efficiently and behaves in a mathematical sound way, i.e., it
satisfies so-called locality and recursiveness conditions. In addition, we
analyze its extremal properties and, in particular, we investigate its maxima
and minima as well as the structure of networks that achieve these values
within the space of so-called level-$1$ networks. We finally briefly compare
this novel index to the two other network balance indices available so-far.


  The regulAS software package is a bioinformatics tool designed to support
computational biology researchers in investigating regulatory mechanisms of
splicing alterations through integrative analysis of large-scale RNA-Seq data
from cancer and healthy human donors, characterized by TCGA and GTEx projects.
This technical report provides a comprehensive overview of regulAS, focusing on
its core functionality, basic modules, experiment configuration, further
extensibility and customisation.
  The core functionality of regulAS enables the automation of computational
experiments, efficient results storage and processing, and streamlined workflow
management. Integrated basic modules extend regulAS with features such as
RNA-Seq data retrieval from the public multi-omics UCSC Xena data repository,
predictive modeling and feature ranking capabilities using the scikit-learn
package, and flexible reporting generation for analysing gene expression
profiles and relevant modulations of alternative splicing aberrations across
tissues and cancer types. Experiment configuration is handled through YAML
files with the Hydra and OmegaConf libraries, offering a user-friendly
approach. Additionally, regulAS allows for the development and integration of
custom modules to handle specialized tasks.
  In conclusion, regulAS provides an automated solution for alternative
splicing and cancer biology studies, enhancing efficiency, reproducibility, and
customization of experimental design, while the extensibility of the pipeline
enables researchers to further tailor the software package to their specific
needs. Source code is available under the MIT license at
https://github.com/slipnitskaya/regulAS.


  Experimental screening and selection pipelines for the discovery of novel
riboswitches are expensive, time-consuming, and inefficient. Using
computational methods to reduce the number of candidates for the screen could
drastically decrease these costs. However, existing computational approaches do
not fully satisfy all requirements for the design of such initial screening
libraries. In this work, we present a new method, libLEARNA, capable of
providing RNA focus libraries of diverse variable-length qualified candidates.
Our novel structure-based design approach considers global properties as well
as desired sequence and structure features. We demonstrate the benefits of our
method by designing theophylline riboswitch libraries, following a previously
published protocol, and yielding 30% more unique high-quality candidates.


  High-dimensional clinical data have become invaluable resources for genetic
studies, due to their accessibility in biobank-scale datasets and the
development of high performance modeling techniques especially using deep
learning. Recent work has shown that low dimensional embeddings of these
clinical data learned by variational autoencoders (VAE) can be used for
genome-wide association studies and polygenic risk prediction. In this work, we
consider multiple unsupervised learning methods for learning disentangled
representations, namely autoencoders, VAE, beta-VAE, and FactorVAE, in the
context of genetic association studies. Using spirograms from UK Biobank as a
running example, we observed improvements in the number of genome-wide
significant loci, heritability, and performance of polygenic risk scores for
asthma and chronic obstructive pulmonary disease by using FactorVAE or
beta-VAE, compared to standard VAE or non-variational autoencoders. FactorVAEs
performed effectively across multiple values of the regularization
hyperparameter, while beta-VAEs were much more sensitive to the hyperparameter
values.


  In recent years, there has been an explosion of research on the application
of deep learning to the prediction of various peptide properties, due to the
significant development and market potential of peptides. Molecular dynamics
has enabled the efficient collection of large peptide datasets, providing
reliable training data for deep learning. However, the lack of systematic
analysis of the peptide encoding, which is essential for AI-assisted
peptide-related tasks, makes it an urgent problem to be solved for the
improvement of prediction accuracy. To address this issue, we first collect a
high-quality, colossal simulation dataset of peptide self-assembly containing
over 62,000 samples generated by coarse-grained molecular dynamics (CGMD).
Then, we systematically investigate the effect of peptide encoding of amino
acids into sequences and molecular graphs using state-of-the-art sequential
(i.e., RNN, LSTM, and Transformer) and structural deep learning models (i.e.,
GCN, GAT, and GraphSAGE), on the accuracy of peptide self-assembly prediction,
an essential physiochemical process prior to any peptide-related applications.
Extensive benchmarking studies have proven Transformer to be the most powerful
sequence-encoding-based deep learning model, pushing the limit of peptide
self-assembly prediction to decapeptides. In summary, this work provides a
comprehensive benchmark analysis of peptide encoding with advanced deep
learning models, serving as a guide for a wide range of peptide-related
predictions such as isoelectric points, hydration free energy, etc.


  Molecular knowledge resides within three different modalities of information
sources: molecular structures, biomedical documents, and knowledge bases.
Effective incorporation of molecular knowledge from these modalities holds
paramount significance in facilitating biomedical research. However, existing
multimodal molecular foundation models exhibit limitations in capturing
intricate connections between molecular structures and texts, and more
importantly, none of them attempt to leverage a wealth of molecular expertise
derived from knowledge graphs. In this study, we introduce MolFM, a multimodal
molecular foundation model designed to facilitate joint representation learning
from molecular structures, biomedical texts, and knowledge graphs. We propose
cross-modal attention between atoms of molecular structures, neighbors of
molecule entities and semantically related texts to facilitate cross-modal
comprehension. We provide theoretical analysis that our cross-modal
pre-training captures local and global molecular knowledge by minimizing the
distance in the feature space between different modalities of the same
molecule, as well as molecules sharing similar structures or functions. MolFM
achieves state-of-the-art performance on various downstream tasks. On
cross-modal retrieval, MolFM outperforms existing models with 12.13% and 5.04%
absolute gains under the zero-shot and fine-tuning settings, respectively.
Furthermore, qualitative analysis showcases MolFM's implicit ability to provide
grounding from molecular substructures and knowledge graphs. Code and models
are available on https://github.com/BioFM/OpenBioMed.


  The classical Sankoff algorithm for the simultaneous folding and alignment of
homologous RNA sequences is highly influential, but it suffers from two major
limitations in efficiency and modeling power. First, it takes $O(n^6)$ for two
sequences where n is the average sequence length. Most implementations and
variations reduce the runtime to $O(n^3)$ by restricting the alignment search
space, but this is still too slow for long sequences such as full-length viral
genomes. On the other hand, the Sankoff algorithm and all its existing
implementations use a rather simplistic alignment model, which can result in
poor alignment accuracy. To address these problems, we propose LinearSankoff,
which seamlessly integrates the original Sankoff algorithm with a powerful
Hidden Markov Model-based alignment module. This extension substantially
improves alignment quality, which in turn benefits secondary structure
prediction quality, confirmed over a diverse set of RNA families. LinearSankoff
also applies beam search heuristics and the A$^\star$-like algorithm to achieve
that runtime scales linearly with sequence length. LinearSankoff is the first
linear-time algorithm for simultaneous folding and alignment, and the first
such algorithm to scale to coronavirus genomes (n $\approx$ 30,000nt). It only
takes 10 minutes for a pair of SARS-CoV-2 and SARS-related genomes, and
outperforms previous work at identifying crucial conserved structures between
the two genomes.


  Molecular Dynamics (MD) simulations are fundamental computational tools for
the study of proteins and their free energy landscapes. However, sampling
protein conformational changes through MD simulations is challenging due to the
relatively long timescales of these processes. Many enhanced sampling
approaches have emerged to tackle this problem, including biased and
path-sampling methods. In this perspective, we focus on adaptive sampling
algorithms. These techniques differ from other approaches because the
thermodynamic ensemble is preserved and the sampling is enhanced solely by
restarting MD trajectories at particularly chosen seeds, rather than
introducing biasing forces. We begin our treatment with an overview of
theoretically transparent methods where we discuss principles and guidelines
for adaptive sampling. Then, we present a brief summary of select methods that
have been applied to realistic systems in the past. Finally, we discuss recent
advances in adaptive sampling methodology powered by machine learning
techniques as well as their shortcomings.


  The field of RNA secondary structure prediction has made significant progress
with the adoption of deep learning techniques. In this work, we present the
RNAformer, a lean deep learning model using axial attention and recycling in
the latent space. We gain performance improvements by designing the
architecture for modeling the adjacency matrix directly in the latent space and
by scaling the size of the model. Our approach achieves state-of-the-art
performance on the popular TS0 benchmark dataset and even outperforms methods
that use external information. Further, we show experimentally that the
RNAformer can learn a biophysical model of the RNA folding process.


  Gromov-Wasserstein distance has found many applications in machine learning
due to its ability to compare measures across metric spaces and its invariance
to isometric transformations. However, in certain applications, this invariance
property can be too flexible, thus undesirable. Moreover, the
Gromov-Wasserstein distance solely considers pairwise sample similarities in
input datasets, disregarding the raw feature representations. We propose a new
optimal transport-based distance, called Augmented Gromov-Wasserstein, that
allows for some control over the level of rigidity to transformations. It also
incorporates feature alignments, enabling us to better leverage prior knowledge
on the input data for improved performance. We present theoretical insights
into the proposed metric. We then demonstrate its usefulness for single-cell
multi-omic alignment tasks and a transfer learning scenario in machine
learning.


  This is the second part of the previous review. In the previous review we
suspected that Orai3 channels were involved in lung cancer and more precisely
in several cancers. Here we confirm that calcium dysregulation is important for
cancer development. in this paper we show that Orai3 is an upstream activator
of AKT and we prove that AKT is involved in chemoresistance in NSCLC.


  Compartmental models are valuable tools for investigating infectious
diseases. Researchers building such models typically begin with a simple
structure where compartments correspond to individuals with different
epidemiological statuses, e.g., the classic SIR model which splits the
population into susceptible, infected, and recovered compartments. However, as
more information about a specific pathogen is discovered, or as a means to
investigate the effects of heterogeneities, it becomes useful to stratify
models further -- for example by age, geographic location, or pathogen strain.
The operation of constructing stratified compartmental models from a pair of
simpler models resembles the Cartesian product used in graph theory, but
several key differences complicate matters. In this article we give explicit
mathematical definitions for several so-called ``model products'' and provide
examples where each is suitable. We also provide examples of model
stratification where no existing model product will generate the desired
result.


  Prokaryotic gene prediction plays an important role in understanding the
biology of organisms and their function with applications in medicine and
biotechnology. Although the current gene finders are highly sensitive in
finding long genes, their sensitivity decreases noticeably in finding shorter
genes (<180 nts). The culprit is insufficient annotated gene data to identify
distinguishing features in short open reading frames (ORFs). We develop a deep
learning-based method called ProtiGeno, specifically targeting short
prokaryotic genes using a protein language model trained on millions of evolved
proteins. In systematic large-scale experiments on 4,288 prokaryotic genomes,
we demonstrate that ProtiGeno predicts short coding and noncoding genes with
higher accuracy and recall than the current state-of-the-art gene finders. We
discuss the predictive features of ProtiGeno and possible limitations by
visualizing the three-dimensional structure of the predicted short genes. Data,
codes, and models are available at https://github.com/tonytu16/protigeno.


  Graph Neural Networks (GNNs) have emerged as a prominent class of data-driven
methods for molecular property prediction. However, a key limitation of typical
GNN models is their inability to quantify uncertainties in the predictions.
This capability is crucial for ensuring the trustworthy use and deployment of
models in downstream tasks. To that end, we introduce AutoGNNUQ, an automated
uncertainty quantification (UQ) approach for molecular property prediction.
AutoGNNUQ leverages architecture search to generate an ensemble of
high-performing GNNs, enabling the estimation of predictive uncertainties. Our
approach employs variance decomposition to separate data (aleatoric) and model
(epistemic) uncertainties, providing valuable insights for reducing them. In
our computational experiments, we demonstrate that AutoGNNUQ outperforms
existing UQ methods in terms of both prediction accuracy and UQ performance on
multiple benchmark datasets. Additionally, we utilize t-SNE visualization to
explore correlations between molecular features and uncertainty, offering
insight for dataset improvement. AutoGNNUQ has broad applicability in domains
such as drug discovery and materials science, where accurate uncertainty
quantification is crucial for decision-making.


  Language models, primarily transformer-based ones, obtained colossal success
in NLP. To be more precise, studies like BERT in NLU and works such as GPT-3
for NLG are very crucial. DNA sequences are very close to natural language in
terms of structure, so if the DNA-related bioinformatics domain is concerned,
discriminative models, like DNABert, exist. Yet, the generative side of the
coin is mainly unexplored to the best of our knowledge. Consequently, we
focused on developing an autoregressive generative language model like GPT-3
for DNA sequences. Because working with whole DNA sequences is challenging
without substantial computational resources, we decided to carry out our study
on a smaller scale, focusing on nucleotide sequences of human genes, unique
parts in DNA with specific functionalities, instead of the whole DNA. This
decision did not change the problem structure a lot due to the fact that both
DNA and genes can be seen as 1D sequences consisting of four different
nucleotides without losing much information and making too much simplification.
First of all, we systematically examined an almost entirely unexplored problem
and observed that RNNs performed the best while simple techniques like N-grams
were also promising. Another beneficial point was learning how to work with
generative models on languages we do not understand, unlike natural language.
How essential using real-life tasks beyond the classical metrics such as
perplexity is observed. Furthermore, checking whether the data-hungry nature of
these models can be changed through selecting a language with minimal
vocabulary size, four owing to four different types of nucleotides, is
examined. The reason for reviewing this was that choosing such a language might
make the problem easier. However, what we observed in this study was it did not
provide that much of a change in the amount of data needed.


  This review article is the second part of the project ``Selected Topics of
Social Physics". The first part has been devoted to equilibrium systems. The
present part considers nonequilibrium systems. The style of the paper combines
the features of a tutorial and a review, which, from one side, makes it easy to
read for nonspecialists aiming at grasping the basics of social physics, and
from the other side, describes several rather recent original models containing
new ideas that could be of interest to experienced researchers in the field.
The present material is based on the lectures that the author had been giving
during several years at the Swiss Federal Institute of Technology in Zurich
(ETH Zurich).


  While extensive literature exists on the COVID-19 pandemic at regional and
national levels, understanding its dynamics and consequences at the city level
remains limited. This study investigates the pandemic in Maring\'a, a
medium-sized city in Brazil's South Region, using data obtained by actively
monitoring the disease from March 2020 to June 2022. Despite prompt and robust
interventions, COVID-19 cases increased exponentially during the early spread
of COVID-19, with a reproduction number lower than that observed during the
initial outbreak in Wuhan. Our research demonstrates the remarkable impact of
non-pharmaceutical interventions on both mobility and pandemic indicators,
particularly during the onset and the most severe phases of the emergency.
However, our results suggest that the city's measures were primarily reactive
rather than proactive. Maring\'a faced six waves of cases, with the third and
fourth waves being the deadliest, responsible for over two-thirds of all deaths
and overwhelming the local healthcare system. Excess mortality during this
period exceeded deaths attributed to COVID-19, indicating that the burdened
healthcare system may have contributed to increased mortality from other
causes. By the end of the fourth wave, nearly three-quarters of the city's
population had received two vaccine doses, significantly decreasing deaths
despite the surge caused by the Omicron variant. Finally, we compare these
findings with the national context and other similarly sized cities,
highlighting substantial heterogeneities in the spread and impact of the
disease.


  Laboratory experiments of bacterial colonies (e.g., \emph{Escherichia coli})
under well-controlled conditions often lead to evolutionary diversification in
which (at least) two ecotypes, each one specialized in the consumption of a
different set of metabolic resources, branch out from an initially monomorphic
population. Empirical evidence suggests that, even under fixed and stable
conditions, such an ``evolutionary branching'' occurs in a stochastic way,
meaning that: (i) it is observed in a significant fraction, but not all, of the
experimental repetitions, (ii) it may emerge at broadly diverse times, and
(iii) the relative abundances of the resulting subpopulations are variable
across experiments. Theoretical approaches shedding light on the possible
emergence of evolutionary branching in this type of conditions have been
previously developed within the theory of ``adaptive dynamics''. Such
approaches are typically deterministic -- or incorporate at most demographic or
finite-size fluctuations which become negligible for the extremely large
populations of these experiments -- and, thus, do not permit to reproduce the
empirically observed large degree of variability. Here, we make further
progress and shed new light on the stochastic nature of evolutionary outcomes
by introducing the idea of ``soft'' trade-offs (as opposed to ``hard'' ones).
This introduces a natural new source of stochasticity which allows one to
account for the empirically observed variability as well as to make predictions
for the likelihood of evolutionary branching to be observed, thus helping to
bridge the gap between theory and experiments.


  Molecular dynamics is a powerful tool for studying the thermodynamics and
kinetics of complex molecular events. However, these simulations can rarely
sample the required time scales in practice. Transition path sampling overcomes
this limitation by collecting unbiased trajectories capturing the relevant
events. Moreover, the integration of machine learning can boost the sampling
while simultaneously learning a quantitative representation of the mechanism.
Still, the resulting trajectories are by construction
non-Boltzmann-distributed, preventing the calculation of free energies and
rates. We developed an algorithm to approximate the equilibrium path ensemble
from machine learning-guided path sampling data. At the same time, our
algorithm provides efficient sampling, the mechanism, free energy, and rates of
rare molecular events at a very moderate computational cost. We tested the
method on the folding of the mini-protein chignolin. Our algorithm is
straightforward and data-efficient, opening the door to applications on many
challenging molecular systems.


  The interaction between elephants and their environment has profound
implications for both ecology and conservation strategies. This study presents
an analytical approach to decipher the intricate patterns of elephant movement
in Sub-Saharan Africa, concentrating on key ecological drivers such as seasonal
variations and rainfall patterns. Despite the complexities surrounding these
influential factors, our analysis provides a holistic view of elephant
migratory behavior in the context of the dynamic African landscape. Our
comprehensive approach enables us to predict the potential impact of these
ecological determinants on elephant migration, a critical step in establishing
informed conservation strategies. This projection is particularly crucial given
the impacts of global climate change on seasonal and rainfall patterns, which
could substantially influence elephant movements in the future. The findings of
our work aim to not only advance the understanding of movement ecology but also
foster a sustainable coexistence of humans and elephants in Sub-Saharan Africa.
By predicting potential elephant routes, our work can inform strategies to
minimize human-elephant conflict, effectively manage land use, and enhance
anti-poaching efforts. This research underscores the importance of integrating
movement ecology and climatic variables for effective wildlife management and
conservation planning.


  The potential effects of conservation actions on threatened species can be
predicted using ensemble ecosystem models by forecasting populations with and
without intervention. These model ensembles commonly assume stable coexistence
of species in the absence of available data. However, existing
ensemble-generation methods become computationally inefficient as the size of
the ecosystem network increases, preventing larger networks from being studied.
We present a novel sequential Monte Carlo sampling approach for ensemble
generation that is orders of magnitude faster than existing approaches. We
demonstrate that the methods produce equivalent parameter inferences, model
predictions, and tightly constrained parameter combinations using a novel
sensitivity analysis method. For one case study, we demonstrate a speed-up from
108 days to 6 hours, while maintaining equivalent ensembles. Additionally, we
demonstrate how to identify the parameter combinations that strongly drive
feasibility and stability, drawing ecological insight from the ensembles. Now,
for the first time, larger and more realistic networks can be practically
simulated and analysed.


  Predicting synergistic drug combinations can help accelerate discovery of
cancer treatments, particularly therapies personalized to a patient's specific
tumor via biopsied cells. In this paper, we propose a novel setting and models
for in-context drug synergy learning. We are given a small "personalized
dataset" of 10-20 drug synergy relationships in the context of specific cancer
cell targets. Our goal is to predict additional drug synergy relationships in
that context. Inspired by recent work that pre-trains a GPT language model (LM)
to "in-context learn" common function classes, we devise novel pre-training
schemes that enable a GPT model to in-context learn "drug synergy functions".
Our model -- which does not use any textual corpora, molecular fingerprints,
protein interaction or any other domain-specific knowledge -- is able to
achieve competitive results. We further integrate our in-context approach with
a genetic algorithm to optimize model prompts and select synergy candidates to
test after conducting a patient biopsy. Finally, we explore a novel task of
inverse drug design which can potentially enable the design of drugs that
synergize specifically to target a given patient's "personalized dataset". Our
findings can potentially have an important impact on precision cancer medicine,
and also raise intriguing questions on non-textual pre-training for LMs.


  We study the roles of social and individual learning on outcomes of the
Minority Game model of a financial market. Social learning occurs via agents
adopting the strategies of their neighbours within a social network, while
individual learning results in agents changing their strategies without input
from other agents. In particular, we show how social learning can undermine
efficiency of the market due to negative frequency dependent selection and loss
of strategy diversity. The latter of which can lock the population into a
maximally inefficient state. We show how individual learning can rescue a
population engaged in social learning from such inefficiencies.


  Given an allotment of land divided into parcels, government decision-makers,
private developers, and conservation biologists can collaborate to select which
parcels to protect, in order to accomplish sustainable ecological goals with
various constraints. In this paper, we propose a mixed-integer optimization
model that considers the presence of multiple species on these parcels, subject
to predator-prey relationships and crowding effects.


  Background: A recent epidemiological analysis of staggered policy
implementation reported a 29.4% reduction in COVID-19 cases by maintaining
school mask mandates in the greater Boston area during the first half of 2022.
The robustness of their results and the appropriateness of methodology are
explored. Methods: Using data from the Massachusetts Department of Elementary
and Secondary Education and the Centers for Disease Control and Prevention, we
re-analyze differences in COVID-19 incidence in school districts that did and
did not lift mask mandates using the same districts as the original study and
expanded the analysis to the entire state of Massachusetts. We present changes
in case rates and differences in prior immunity in areas with different mask
lifting policies. Results: The Boston and Chelsea districts, which maintained
mask mandates, were outliers in terms of size, demographics, and testing. We
failed to find a notable change in student cases in mask mandate districts
compared with the 70 districts in the original study (-0.08/1000; p=0.98) and
found a slight increase compared with a statewide control group +3.63/1000
(p=0.291). Results were similar for students and staff combined. Districts that
dropped mask mandates first experienced the largest decreases in cases (22%
drop vs 12% in the masked districts). There was a moderate to strong
relationship (R2 = 0.35-0.66; p-values <0.001) between prior community
infection burden and district case rates in Spring 2022, with prior immunity
alone explaining as much as two-thirds of the variation in case rates in Spring
of 2022. Conclusions: We fail to find any consistent notable negative
relationship between school mask mandates and infection rates in the Greater
Boston Area or state of Massachusetts during the 2021-2022 academic year.


  In this paper, we present a model describing the collective motion of birds.
We explore the dynamic relationship between followers and leaders, wherein a
select few agents, known as leaders, can initiate spontaneous changes in
direction without being influenced by external factors like predators. Starting
at the microscopic level, we develop a kinetic model that characterizes the
behaviour of large crowds with transient leadership. One significant challenge
lies in managing topological interactions, as identifying nearest neighbors in
extensive systems can be computationally expensive. To address this, we propose
a novel stochastic particle method to simulate the mesoscopic dynamics and
reduce the computational cost of identifying closer agents from quadratic to
logarithmic complexity using a $k$-nearest neighbours search algorithm with a
binary tree. Lastly, we conduct various numerical experiments for different
scenarios to validate the algorithm's effectiveness and investigate collective
dynamics in both two and three dimensions.


  Given a set $X$ of species, a phylogenetic tree is an unrooted binary tree
whose leaves are bijectively labelled by $X$. Such trees can be used to show
the way species evolve over time. One way of understanding how topologically
different two phylogenetic trees are, is to construct a minimum-size agreement
forest: a partition of $X$ into the smallest number of blocks, such that the
blocks induce homeomorphic, non-overlapping subtrees in both trees. This
comparison yields insight into commonalities and differences in the evolution
of $X$ across the two trees. Computing a smallest agreement forest is NP-hard
(Hein, Jiang, Wang and Zhang, Discrete Applied Mathematics 71(1-3), 1996). In
this work we study the problem on caterpillars, which are path-like
phylogenetic trees. We will demonstrate that, even if we restrict the input to
this highly restricted subclass, the problem remains NP-hard and is in fact
APX-hard. Furthermore we show that for caterpillars two standard reductions
rules well known in the literature yield a tight kernel of size at most $7k$,
compared to $15k$ for general trees (Kelk and Simone, SIAM Journal on Discrete
Mathematics 33(3), 2019). Finally we demonstrate that we can determine if two
caterpillars have an agreement forest with at most $k$ blocks in $O^*(2.49^k)$
time, compared to $O^*(3^k)$ for general trees (Chen, Fan and Sze, Theoretical
Computater Science 562, 2015), where $O^*(.)$ suppresses polynomial factors.


  Danger requires a strong rapid response. Speedy triggers are prone to false
signals. False alarms can be costly, requiring strong negative regulators to
oppose the initial triggers. Strongly opposed forces can easily be perturbed,
leading to imbalance and disease. For example, immunity and fear response
balance strong rapid triggers against widespread slow negative regulators.
Diseases of immunity and behavior arise from imbalance. A different opposition
of forces occurs in mammalian growth, which balances strong paternally
expressed accelerators against maternally expressed suppressors. Diseases of
overgrowth or undergrowth arise from imbalance. Other examples of opposing
forces and disease include control of dopamine expression and male versus
female favored traits.


  Coarse-grained molecular models of proteins permit access to length and time
scales unattainable by all-atom models and the simulation of processes that
occur on long-time scales such as aggregation and folding. The reduced
resolution realizes computational accelerations but an atomistic representation
can be vital for a complete understanding of mechanistic details. Backmapping
is the process of restoring all-atom resolution to coarse-grained molecular
models. In this work, we report DiAMoNDBack (Diffusion-denoising Autoregressive
Model for Non-Deterministic Backmapping) as an autoregressive denoising
diffusion probability model to restore all-atom details to coarse-grained
protein representations retaining only C{\alpha} coordinates. The
autoregressive generation process proceeds from the protein N-terminus to
C-terminus in a residue-by-residue fashion conditioned on the C{\alpha} trace
and previously backmapped backbone and side chain atoms within the local
neighborhood. The local and autoregressive nature of our model makes it
transferable between proteins. The stochastic nature of the denoising diffusion
process means that the model generates a realistic ensemble of backbone and
side chain all-atom configurations consistent with the coarse-grained C{\alpha}
trace. We train DiAMoNDBack over 65k+ structures from Protein Data Bank (PDB)
and validate it in applications to a hold-out PDB test set,
intrinsically-disordered protein structures from the Protein Ensemble Database
(PED), molecular dynamics simulations of fast-folding mini-proteins from DE
Shaw Research, and coarse-grained simulation data. We achieve state-of-the-art
reconstruction performance in terms of correct bond formation, avoidance of
side chain clashes, and diversity of the generated side chain configurational
states. We make DiAMoNDBack model publicly available as a free and open source
Python package.


  To learn accurate representations of molecules, it is essential to consider
both chemical and geometric features. To encode geometric information, many
descriptors have been proposed in constrained circumstances for specific types
of molecules and do not have the properties to be ``robust": 1. Invariant to
rotations and translations; 2. Injective when embedding molecular structures.
In this work, we propose a universal and robust Directional Node Pair (DNP)
descriptor based on the graph representations of 3D molecules. Our DNP
descriptor is robust compared to previous ones and can be applied to multiple
molecular types. To combine the DNP descriptor and chemical features in
molecules, we construct the Robust Molecular Graph Convolutional Network
(RoM-GCN) which is capable to take both node and edge features into
consideration when generating molecule representations. We evaluate our model
on protein and small molecule datasets. Our results validate the superiority of
the DNP descriptor in incorporating 3D geometric information of molecules.
RoM-GCN outperforms all compared baselines.


  Designing protein mutants with high stability and activity is a critical yet
challenging task in protein engineering. Here, we introduce PRIME, an
innovative deep learning approach for the zero-shot prediction of both protein
stability and enzymatic activity. PRIME leverages temperature-guided language
modelling, providing robust and precise predictions without relying on prior
experimental mutagenesis data. Tested against 33 protein datasets, PRIME
demonstrated superior predictive performance and generalizability compared to
current state-of-the-art models


  Antibody-antigen interactions play a crucial role in identifying and
neutralizing harmful foreign molecules. In this paper, we investigate the
optimal representation for predicting the binding sites in the two molecules
and emphasize the importance of geometric information. Specifically, we compare
different geometric deep learning methods applied to proteins' inner (I-GEP)
and outer (O-GEP) structures. We incorporate 3D coordinates and spectral
geometric descriptors as input features to fully leverage the geometric
information. Our research suggests that surface-based models are more efficient
than other methods, and our O-GEP experiments have achieved state-of-the-art
results with significant performance improvements.


  In previous work, we have developed an approach for characterizing the
long-term dynamics of classes of Biological Interaction Networks (BINs), based
on "rate-dependent Lyapunov functions". In this work, we show that stronger
notions of convergence can be established by proving structural contractivity
with respect to non-standard norms. We illustrate our theory with examples from
signaling pathways.


  An accurate and timely estimate of the reproduction ratio R of an infectious
disease epidemic is crucial to make projections on its evolution and set up the
appropriate public health response. Estimates of R routinely come from
statistical inference on timelines of cases or their proxies like symptomatic
cases, hospitalizatons, deaths. Here, however, we prove that these estimates of
R may not be accurate if the population is made up of spatially distinct
communities, as the interplay between space and mobility may hide the true
epidemic evolution from surveillance data. This means that surveillance may
underestimate R over long periods, to the point of mistaking a growing epidemic
for a subsiding one, misinforming public health response. To overcome this, we
propose a correction to be applied to surveillance data that removes this bias
and ensures an accurate estimate of R across all epidemic phases. We use
COVID-19 as case study; our results, however, apply to any epidemic where
mobility is a driver of circulation, including major challenges of the next
decades: respiratory infections (influenza, SARS-CoV-2, emerging pathogens),
vector-borne diseases (arboviruses). Our findings will help set up public
health response to these threats, by improving epidemic monitoring and
surveillance.


  The plant light-harvesting pigment-protein complex LHCII is the major antenna
sub-unit of PSII and is generally (though not universally) accepted to play a
role in photoprotective energy dissipation under high light conditions, a
process known Non-Photochemical Quenching (NPQ). The underlying mechanisms of
energy trapping and dissipation within LHCII are still debated. Various
proposed models differ considerably in their molecular and kinetic detail, but
are often based on different interpretations of very similar transient
absorption measurements of isolated complexes. Here we present a simulated
measurement of the fluorescence decay kinetics of quenched LHCII aggregates to
determine whether this relatively simple measurement can discriminate between
different potential NPQ mechanisms. We simulate not just the underlying physics
(excitation, energy migration, quenching and singlet-singlet annihilation) but
also the signal detection and typical experimental data analysis. Comparing
this to a selection of published fluorescence decay kinetics we find that: (1)
Different proposed quenching mechanisms produce noticeably different
fluorescence kinetics even at low (annihilation free) excitation density,
though the degree of difference is dependent on pulse width. (2) Measured decay
kinetics are consistent with most LHCII trimers becoming relatively slow
excitation quenchers. A small sub-population of very fast quenchers produces
kinetics which do not resemble any observed measurement. (3) It is necessary to
consider at least two distinct quenching mechanisms in order to accurately
reproduce experimental kinetics, supporting the idea that NPQ is not a simple
binary switch switch.


  The assembly and persistence of ecological communities can be understood as
the result of the interaction and migration of species. Here we study a single
community subject to migration from a species pool in which inter-specific
interactions are organised according to a bipartite network. Considering the
dynamics of species abundances to be governed by generalised Lotka-Volterra
equations, we extend work on unipartite networks to we derive exact results for
the phase diagram of this model. Focusing on antagonistic interactions, we
describe factors that influence the persistence of the two guilds, locate
transitions to multiple-attractor and unbounded phases, as well identify a
region of parameter space in which consumers are essentially absent in the
local community.


  The adaptive immune response relies on T cells that combine phenotypic
specialization with diversity of T cell receptors (TCRs) to recognize a wide
range of pathogens. TCRs are acquired and selected during T cell maturation in
the thymus. Characterizing TCR repertoires across individuals and T cell
maturation stages is important for better understanding adaptive immune
responses and for developing new diagnostics and therapies. Analyzing a dataset
of human TCR repertoires from thymocyte subsets, we find that the variability
between individuals generated during the TCR V(D)J recombination is maintained
through all stages of T cell maturation and differentiation. The
inter-individual variability of repertoires of the same cell type is of
comparable magnitude to the variability across cell types within the same
individual. To zoom in on smaller scales than whole repertoires, we defined a
distance measuring the relative overlap of locally similar sequences in
repertoires. We find that the whole repertoire models correctly predict local
similarity networks, suggesting a lack of forbidden T cell receptor sequences.
The local measure correlates well with distances calculated using whole
repertoire traits and carries information about cell types.


  Epidemiological models are an important tool in coping with epidemics, as
they offer a forecast, even if often simplistic, of the behavior of the disease
in the population. This allows responsible health agencies to organize
themselves and adopt strategies to minimize and postpone the population's
infection peaks. While during an epidemic outbreak, the available model can be
used to describe the behavior of the disease in order to aim for fast and
efficient forecasts of the epidemiological scenario, once the epidemiological
emergency is over, the objective of the subsequent works is to extend models by
integrating new facts and information, thus providing more efficient tools to
face future epidemics. In this sense, we present an epidemiological model that
takes into account on one hand, the creation of a vaccine during an epidemic
outbreak (as we saw happen in the case of COVID-19) and that, on the other
hand, considers the impact of the cooperative behavioral choices of
individuals.


  This paper presents a novel approach for predicting the relative populations
of protein conformations using AlphaFold 2, an AI-powered method that has
revolutionized biology by enabling the accurate prediction of protein
structures. While AlphaFold 2 has shown exceptional accuracy and speed, it is
designed to predict proteins' single ground state conformations and is limited
in its ability to predict fold switching and the effects of mutations on
conformational landscapes. Here, we demonstrate how AlphaFold 2 can directly
predict the relative populations of different conformations of proteins and
even accurately predict changes in those populations induced by mutations by
subsampling multiple sequence alignments. We tested our method against NMR
experiments on two proteins with drastically different amounts of available
sequence data, Abl1 kinase and the granulocyte-macrophage colony-stimulating
factor, and predicted changes in their relative state populations with
accuracies in excess of 80%. Our method offers a fast and cost-effective way to
predict protein conformations and their relative populations at even single
point mutation resolution, making it a useful tool for pharmacology, analyzing
NMR data, and studying the effects of evolution.


  Imbalance in the 3D structure of plants can be an important indicator of
insufficient light or nutrient supply, as well as excessive wind, (formerly
present) physical barriers, neighbor or storm damage. It can also be a simple
means to detect certain illnesses, since some diseases like the apple
proliferation disease, an infection with the barley yellow dwarf virus or plant
canker can cause abnormal growth, like \enquote{witches' brooms} or burls,
resulting in a deviating 3D plant architecture. However, quantifying imbalance
of plant growth is not an easy task, and it requires a mathematically sound 3D
model of plants to which imbalance indices can be applied. Current models of
plants are often based on stacked cylinders or voxel matrices and do not allow
for measuring the degree of 3D imbalance in the branching structure of the
whole plant.
  On the other hand, various imbalance indices are readily available for
so-called graph-theoretical trees and are frequently used in areas like
phylogenetics and computer science. While only some basic ideas of these
indices can be transferred to the 3D setting, graph-theoretical trees are a
logical foundation for 3D plant models that allow for elegant and natural
imbalance measures.
  In this manuscript, our aim is thus threefold: We first present a new
graph-theoretical 3D model of plants and discuss desirable properties of
imbalance measures in the 3D setting. We then introduce and analyze eight
different 3D imbalance indices and their properties. Thirdly, we illustrate all
our findings using a data set of 63 bush beans. Moreover, we implemented all
our indices in the publicly available \textsf{R}-software package
\textsf{treeDbalance} accompanying this manuscript.


  We use fitness graphs, or directed cube graphs, for analyzing evolutionary
reversibility. The main application is antimicrobial drug resistance.
Reversible drug resistance has been observed both clinically and
experimentally. If drug resistance depends on a single point mutation, then a
possible scenario is that the mutation reverts back to the wild-type codon
after the drug has been discontinued, so that susceptibility is fully restored.
In general, a drug pause does not automatically imply fast elimination of drug
resistance. Also if drug resistance is reversible, the threshold concentration
for reverse evolution may be lower than for forward evolution. For a
theoretical understanding of evolutionary reversibility, including threshold
asymmetries, it is necessary to analyze obstacles in fitness landscapes. We
compare local and global obstacles, obstacles for forward and reverse
evolution, and conjecture that favorable landscapes for forward evolution
correlate with evolution being reversible. Both suboptimal peaks and plateaus
are analyzed with some observations on the impact of redundancy and
dimensionality. Our findings are compared with laboratory studies on
irreversible malarial drug resistance.


  Protein engineering is an emerging field in biotechnology that has the
potential to revolutionize various areas, such as antibody design, drug
discovery, food security, ecology, and more. However, the mutational space
involved is too vast to be handled through experimental means alone. Leveraging
accumulative protein databases, machine learning (ML) models, particularly
those based on natural language processing (NLP), have considerably expedited
protein engineering. Moreover, advances in topological data analysis (TDA) and
artificial intelligence-based protein structure prediction, such as AlphaFold2,
have made more powerful structure-based ML-assisted protein engineering
strategies possible. This review aims to offer a comprehensive, systematic, and
indispensable set of methodological components, including TDA and NLP, for
protein engineering and to facilitate their future development.


  Certain periodical cicadas exhibit life cycles with durations of 13 or 17
years, and it is now generally accepted that such large prime numbers arise
evolutionarily to avoid synchrony with predators. Less well explored is the
question of how, in the face of intrinsic biological and environmental noise,
insects within a brood emerge together in large successive swarms from
underground during springtime warming. Here we consider the decision-making
process of underground cicada nymphs experiencing random but
spatially-correlated thermal microclimates like those in nature. Introducing
short-range communication between insects leads to a model of consensus
building that maps on to the statistical physics of an Ising model with a
quenched, spatially correlated random magnetic field and annealed site
dilution, which displays the kinds of collective swarms seen in nature.


  Small extracellular vesicles (sEVs) are known to play an important role in
the communication between distant cells and to deliver biological information
throughout the body. To date, many studies have focused on the role of sEVs
characteristics such as cell origin, surface composition, and molecular cargo
on the resulting uptake by the recipient cell. Yet, a full understanding of the
sEV fusion process with recipient cells and in particular the role of cell
membrane physical properties on the uptake are still lacking. Here we explore
this problem using sEVs from a cellular model of triple-negative breast cancer
fusing to a range of synthetic planar lipid bilayers both with and without
cholesterol, and designed to mimic the formation of raft-like nanodomains in
cell membranes. Using time-resolved Atomic Force Microscopy we were able to
track the sEVs interaction with the different model membranes, showing the
process to be strongly dependent on the local membrane fluidity. The strongest
interaction and fusion is observed over the less fluid regions, with sEVs even
able to disrupt ordered domains at sufficiently high cholesterol concentration.
Our findings suggest the biophysical characteristics of recipient cell
membranes to be crucial for sEVs uptake regulation.


  Accelerating the discovery of novel and more effective therapeutics is an
important pharmaceutical problem in which deep learning is playing an
increasingly significant role. However, real-world drug discovery tasks are
often characterized by a scarcity of labeled data and significant covariate
shift$\unicode{x2013}\unicode{x2013}$a setting that poses a challenge to
standard deep learning methods. In this paper, we present Q-SAVI, a
probabilistic model able to address these challenges by encoding explicit prior
knowledge of the data-generating process into a prior distribution over
functions, presenting researchers with a transparent and probabilistically
principled way to encode data-driven modeling preferences. Building on a novel,
gold-standard bioactivity dataset that facilitates a meaningful comparison of
models in an extrapolative regime, we explore different approaches to induce
data shift and construct a challenging evaluation setup. We then demonstrate
that using Q-SAVI to integrate contextualized prior knowledge of drug-like
chemical space into the modeling process affords substantial gains in
predictive accuracy and calibration, outperforming a broad range of
state-of-the-art self-supervised pre-training and domain adaptation techniques.


  The recently introduced Theory of the Adjacent Possible (TAP) is a model of
combinatorial innovation aiming to explain the "hockey-stick" upward trend of
human technological evolution, where an explosion in the number of produced
items with increasing complexity suddenly occurs. In addition, the TAP model
was also used to explain the rapidly emerging biological complexity. Inspired
by TAP here we propose a reaction-diffusion system aiming to extend the model
in both space and time. We show that the new model exhibits similar
characteristics to the TAP model, like the sudden increase in the production of
items, after a longer period of slow growth. The new model also exhibits wave
propagation of "innovation", resulting in self-sustained complex interference
patterns.


  Phylogenetic networks represent evolutionary histories of sets of taxa where
horizontal evolution or hybridization has occurred. Placing a Markov model of
evolution on a phylogenetic network gives a model that is particularly amenable
to algebraic study by representing it as an algebraic variety. In this paper,
we give a formula for the dimension of the variety corresponding to a
triangle-free level-1 phylogenetic network under a group-based evolutionary
model. On our way to this, we give a dimension formula for codimension zero
toric fiber products. We conclude by illustrating applications to
identifiability.


  The accumulation of somatic mutations is a driver of cancer and has long been
associated with ageing. Due to limitations in quantifying mutation burden with
age in non-cancerous tissues, the impact of somatic mutations in other ageing
phenotypes is unclear. Recent advances in DNA sequencing technologies have
allowed the large-scale quantification of somatic mutations in ageing. These
studies have revealed a gradual accumulation of mutations in most normal
tissues with age as well as a substantial clonal expansion driven mostly by
cancer-related mutations. Nevertheless, because of the relatively modest burden
of age-related somatic mutations identified so far and their stochastic nature,
it is difficult to envision how somatic mutation accumulation alone can explain
most ageing phenotypes that develop gradually. Studies across species have also
found that longer-lived species have lower somatic mutation rates, though these
could be explained by selective pressures to reduce or postpone cancer as
longevity increases. Overall, with a few exceptions like cancer, results from
recent DNA sequencing studies do not add weight to the idea that somatic
mutations with age drive ageing phenotypes and the phenotypic role, if any, of
somatic mutations in ageing remains unclear. Recent studies in patients with
somatic mutation burden and no signs of accelerated ageing further question the
role of somatic mutations in ageing.


  The emergence of oxygenic photosynthesis was a major event in Earth's
evolutionary history and was facilitated by chlorophylls (a major category of
photopigments). The accurate modelling of photopigments is important to
understand the characteristics of putative extraterrestrial life and its
spectral signatures (detectable by future telescopes). In this paper, we
perform a detailed assessment of various time-dependent density-functional
theory (TD-DFT) methods for predicting the absorption spectra of chlorophyll a,
with particular emphasis on modern low-cost approximations. We also investigate
a potential extraterrestrial photopigment called phot0 and demonstrate that the
electronegativity of the metal ion may exert a direct influence on the
locations of the absorption peaks, with higher electronegativity inducing
blue-shifting and vice-versa. Based on these calculations, we established that
global-hybrid approximations with a moderate percentage of exact exchange --
such as M06 and PW6B95 -- are the most appropriate compromise between cost and
accuracy for the computational characterization of photopigments of
astrobiological interest. We conclude with a brief assessment of the
implications and avenues for future research.


  This paper examines the influence of climate change on the human mind and
culture from the 13th century to the 21st century. By quantitatively analyzing
100,000 paintings and the biological data of over 2,000 artists, an interesting
U-shaped pattern in the lightness of paintings was found, which correlated with
trends in global temperature. Event study analysis revealed that when an artist
is subjected to a high-temperature shock, their paintings become brighter in
later periods. Moreover, the effects are more pronounced in art genres that
rely less on real things and more on the artist's imagination, indicating the
influence of artists' minds. Overall, this study demonstrates the significant
and enduring influence of climate on the human mind and culture over centuries.


  We establish an analogy between the Fokker-Planck equation describing
evolutionary landscape dynamics and the Schr\"{o}dinger equation which
characterizes quantum mechanical particles, showing how a population with
multiple genetic traits evolves analogously to a wavefunction under a
multi-dimensional energy potential in imaginary time. Furthermore, we discover
within this analogy that the stationary population distribution on the
landscape corresponds exactly to the ground-state wavefunction. This
mathematical equivalence grants entry to a wide range of analytical tools
developed by the quantum mechanics community, such as the Rayleigh-Ritz
variational method and the Rayleigh-Schr\"{o}dinger perturbation theory,
allowing us to not only make reasonable quantitative assessments but also
explore fundamental biological inquiries. We demonstrate the effectiveness of
these tools by estimating the population success on landscapes where precise
answers are elusive, and unveiling the ecological consequences of
stress-induced mutagenesis -- a prevalent evolutionary mechanism in pathogenic
and neoplastic systems. We show that, even in a unchanging environment, a sharp
mutational burst resulting from stress can always be advantageous, while a
gradual increase only enhances population size when the number of relevant
evolving traits is limited. Our interdisciplinary approach offers novel
insights, opening up new avenues for deeper understanding and predictive
capability regarding the complex dynamics of evolving populations.


  It is obvious that both epigenetic and non-epigenetic actors contribute to
tumorigenesis in chondrosarcomas and more generally in other cancers. Thus, the
main altered pathways in chondrosarcomas are now well established and include
both epigenetic and non-epigenetic pathways such as the PI3K-AKT signaling,
EGFR overexpression, SPARC overexpression, c-myc overexpression, IHH/GLI1 axis,
loss of Rb function, HIF1-alpha stabilization, IDH1 mutations, hypermethylation
and SIRT1. This review aims to provide a detailed analysis of these pathways
and highlights recurrent interactions between non-epigenetic and epigenetic
actors in chondrosarcomas, raising the intriguing possibility of developing
therapeutics targeting both epigenetic and non-epigenetic actors and supporting
data from previous studies. Finally, we propose some strategies for targeting
chondrosarcomas in vivo based on properties of this tumor.


  Phylogenetic networks are an extension of phylogenetic trees that allow for
the representation of reticulate evolution events. One of the classes of
networks that has gained the attention of the scientific community over the
last years is the class of orchard networks, that generalizes tree-child
networks, one of the most studied classes of networks.
  In this paper we focus on the combinatorial and algorithmic problem of the
generation of orchard networks, and also of tree-child networks. To this end,
we use that these networks are defined as those that can be recovered by a
reversing a certain reduction process. Then, we show how to choose a
``minimum'' reduction process among all that can be applied to a network, and
hence we get a unique representation of the network that, in fact, can be given
in terms of sequences of pairs of integers, whose length is related to the
number of leaves and reticulations of the network. Therefore, the generation of
networks is reduced to the generation of such sequences of pairs. Our main
result is a recursive method for the efficient generation of all minimum
sequences, and hence of all orchard (or tree-child) networks with a given
number of leaves and reticulations.
  An implementation in C of the algorithms described in this paper, along with
some computational experiments, can be downloaded from the public repository
https://github.com/gerardet46/OrchardGenerator. Using this implementation, we
have computed the number of orchard networks with at most 6 leaves and 8
reticulations.


  In this work, we have described the mathematical modeling of COVID-19
transmission using fractional differential equations. The mathematical modeling
of infectious disease goes back to the 1760s when the famous mathematician
Daniel Bernoulli used an elementary version of compartmental modeling to find
the effectiveness of deliberate smallpox inoculation on life expectancy. We
have used the well-known SIR (Susceptible, Infected and Recovered) model of
Kermack & McKendrick to extend the analysis further by including exposure,
quarantining, insusceptibility and deaths in a SEIQRDP model. Further, we have
generalized this model by using the solutions of Fractional Differential
Equations to test the accuracy and validity of the mathematical modeling
techniques against Canadian COVID-19 trends and spread of real-world disease.
Our work also emphasizes the importance of Personal Protection Equipment (PPE)
and impact of social distancing on controlling the spread of COVID-19.


  Genome-wide association studies (GWAS) involving increasing sample sizes have
identified hundreds of genetic variants associated with complex diseases, such
as type 2 diabetes (T2D); however, it is unclear how GWAS hits form unique
topological structures in protein-protein interaction (PPI) networks. Using
persistent homology, this study explores the evolution and persistence of the
topological features of T2D GWAS hits in the PPI network with increasing
P-value thresholds. We define an $n$-dimensional persistent disease module as a
higher-order generalization of the largest connected component (LCC). The
0-dimensional persistent T2D disease module is the LCC of the T2D GWAS hits,
which is significantly detected in the PPI network (196 nodes and 235 edges,
P$<$0.05). In the 1-dimensional homology group analysis, all 18 1-dimensional
holes (loops) of the T2D GWAS hits persist over all P-value thresholds. The
1-dimensional persistent T2D disease module comprising these 18 persistent
1-dimensional holes is significantly larger than that expected by chance (59
nodes and 83 edges, P$<$0.001), indicating a significant topological structure
in the PPI network. Our computational topology framework potentially possesses
broad applicability to other complex phenotypes in identifying topological
features that play an important role in disease pathobiology.


  Measuring NT-proBNP biomarker is recommended for preliminary diagnostics of
the heart failure. Recent studies suggest a possibility of early screening of
biomarkers in saliva for non-invasive identification of cardiac diseases at the
point-of-care. However, NT-proBNP concentrations in saliva can be thousand time
lower than in blood plasma, going down to pg/mL level. To reach this level, we
developed a label-free aptasensor based on a liquid-gated field effect
transistor using a film of reduced graphene oxide monolayer (rGO-FET) with
immobilized NT-proBNP specific aptamer. We found that, depending on ionic
strength of tested solutions, there were different levels of correlation in
responses of electrical parameters of the rGO-FET aptasensor, namely, the Dirac
point shift and transconductance change. The correlation in response to
NT-proBNP was high for 1.6 mM phosphate-buffered saline (PBS) and zero for 16
mM PBS in a wide range of analyte concentrations, varied from 1 fg/mL to 10
ng/mL. The effects of transconductance and Dirac point shift in PBS solutions
of different concentrations are discussed. The biosensor exhibited a high
sensitivity for both transconductance (2 uS/decade) and Dirac point shift (2.3
mV/decade) in diluted PBS with the linear range from 10 fg/mL to 1 pg/mL. The
aptasensor performance has been also demonstrated in undiluted artificial
saliva with the achieved limit of detection down to 41 fg/mL (~4.6 fM).


  The aim of this short note is twofold. We formulate the general
Kermack-McKendrick epidemic model incorporating static heterogeneity and show
how it simplifies to a scalar Renewal Equation (RE) when separable mixing is
assumed. A key feature is that all information about the heterogeneity is
encoded in one nonlinear real valued function of a real variable. Inspired by
work of R. Pastor-Satorras and C. Castellano, we next investigate mask
efficiency and demonstrate that it is straightforward to rederive from the RE
their main conclusion, that the best way to protect the population as a whole
is to protect yourself. Thus we establish that this conclusion is robust, in
the sense that it also holds outside the world of network models.


  Coupled spreading processes on interdependent networks are studied, where the
spreading infection of diseases or malware in one layer can trigger cascading
failures in another layer and lead to secondary disasters, e.g., disrupting
public services, supply chains, or power distribution. We utilize the dynamic
message-passing method to devise efficient algorithms for inferring the system
states, which allows one to investigate systematically the nature of such
complex intertwined spreading processes and evaluate their impact. Based on the
dynamic message-passing framework and optimal control, we further develop an
effective optimization algorithm for mitigating network failures.


  Pickering emulsions have been widely studied due to their good stability and
potential applications. Nanocellulose including cellulose nanocrystals (CNCs),
cellulose nanofibrils (CNFs), and bacterial cellulose nanofibrils (BCNFs) has
emerged as sustainable stabilizers/emulsifiers in food-related Pickering
emulsions due to their favorable properties such as renewability, low toxicity,
amphiphilicity, biocompatibility, and high aspect ratio. Nanocellulose can be
widely obtained from different sources and extraction methods and can
effectively stabilize Pickering emulsions via the irreversible adsorption onto
oil-water interface. The synergistic effects of nanocellulose and other
substances can further enhance the interfacial networks. The
nanocellulose-based Pickering emulsions have potential food-related
applications in delivery systems, food packaging materials, and fat
substitutes. In this review, we highlight key fundamental work and recent
reports on nanocellulose-based Pickering emulsion systems. The sources and
extraction of nanocellulose and the fabrication of nanocellulose-based
Pickering emulsions are briefly summarized. Furthermore, the synergistic
stability and food-related applications of nanocellulose-stabilized Pickering
emulsions are spotlighted.


  Despite recent advancement in 3D molecule conformation generation driven by
diffusion models, its high computational cost in iterative diffusion/denoising
process limits its application. In this paper, an equivariant consistency model
(EC-Conf) was proposed as a fast diffusion method for low-energy conformation
generation. In EC-Conf, a modified SE (3)-equivariant transformer model was
directly used to encode the Cartesian molecular conformations and a highly
efficient consistency diffusion process was carried out to generate molecular
conformations. It was demonstrated that, with only one sampling step, it can
already achieve comparable quality to other diffusion-based models running with
thousands denoising steps. Its performance can be further improved with a few
more sampling iterations. The performance of EC-Conf is evaluated on both
GEOM-QM9 and GEOM-Drugs sets. Our results demonstrate that the efficiency of
EC-Conf for learning the distribution of low energy molecular conformation is
at least two magnitudes higher than current SOTA diffusion models and could
potentially become a useful tool for conformation generation and sampling. We
release our code at https://github.com/zhi520/EcConf.


  Foraging is crucial for animals to survive. Many species forage in groups, as
individuals communicate to share information about the location of available
resources. For example, eusocial foragers, such as honey bees and many ants,
recruit members from their central hive or nest to a known foraging site.
However, the optimal level of communication and recruitment depends on the
overall group size, the distribution of available resources, and the extent of
interference between multiple individuals attempting to forage from a site. In
this paper, we develop a discrete-time Markov chain model of eusocial foragers,
who communicate information with a certain probability. We compare the
stochastic model and its corresponding infinite-population limit. We find that
foraging efficiency tapers off when recruitment probability is too high -- a
phenomenon that does not occur in the infinite-population model, even though it
occurs for any finite population size. The marginal inefficiency at high
recruitment probability increases as the population increases, similar to a
boundary layer. In particular, we prove there is a significant gap between the
foraging efficiency of finite and infinite population models in the extreme
case of complete communication. We also analyze this phenomenon by
approximating the stationary distribution of foragers over sites in terms of
mean escape times from multiple quasi-steady states. We conclude that for any
finite group of foragers, an individual who has found a resource should only
sometimes recruit others to the same resource. We discuss the relationship
between our analysis and multi-agent multi-arm bandit problems.