  Multidimensional scaling (MDS) is an unsupervised learning technique that
preserves pairwise distances between observations and is commonly used for
analyzing multivariate biological datasets. Recent advances in MDS have
achieved successful classification results, but the configurations heavily
depend on the choice of hyperparameters, limiting its broader application.
Here, we present a self-supervised MDS approach informed by the dispersions of
observations that share a common binary label ($F$-ratio). Our visualization
accurately configures the $F$-ratio while consistently preserving the global
structure with a low data distortion compared to existing dimensionality
reduction tools. Using an algal microbiome dataset, we show that this new
method better illustrates the community's response to the host, suggesting its
potential impact on microbiology and ecology data analysis.


  Summary: PDBImages is an innovative, open-source Node.js package that
harnesses the power of the popular macromolecule structure visualization
software Mol*. Designed for use by the scientific community, PDBImages provides
a means to generate high-quality images for PDB and AlphaFold DB models. Its
unique ability to render and save images directly to files in a browserless
mode sets it apart, offering users a streamlined, automated process for
macromolecular structure visualization. Here, we detail the implementation of
PDBImages, enumerating its diverse image types and elaborating on its
user-friendly setup. This powerful tool opens a new gateway for researchers to
visualize, analyse, and share their work, fostering a deeper understanding of
bioinformatics. Availability and Implementation: PDBImages is available as an
npm package from https://www.npmjs.com/package/pdb-images. The source code is
available from https://github.com/PDBeurope/pdb-images. Contact:
mvaradi@ebi.ac.uk, sameer@ebi.ac.uk


  We present QCAM, a quantum analogue of Content-Addressable Memory (CAM),
useful for finding matches in two sequences of bit-strings. Our QCAM
implementation takes advantage of Grover's search algorithm and proposes a
highly-optimized quantum circuit implementation of the QCAM oracle. Our circuit
construction uses the parallel uniformly controlled rotation gates, which were
used in previous work to generate QBArt encodings. These circuits have a high
degree of quantum parallelism which reduces their critical depth. The optimal
number of repetitions of the Grover iterator used in QCAM depends on the number
of true matches and hence is input dependent. We additionally propose a
hardware-efficient implementation of the quantum counting algorithm (HEQC) that
can infer the optimal number of Grover iterations from the measurement of a
single observable. We demonstrate the QCAM application for computing the
Jaccard similarity between two sets of k-mers obtained from two DNA sequences.


  In order to understand the role of a potentially high infectiousness of
symptomatic individuals or, on the contrary, social distancing, quarantine,
etc. in the course of an infectious disease, this paper introduces a natural
variant of the popular contact process that distinguishes between asymptomatic
and symptomatic individuals. Infected individuals all recover at rate one but
infect nearby individuals at a rate that depends on whether they show the
symptoms of the disease or not. Newly infected individuals are always
asymptomatic and may or may not show the symptoms before they recover. The
analysis of the corresponding mean-field model reveals that, in the absence of
local interactions, regardless of the rate at which asymptomatic individuals
become symptomatic, there is an epidemic whenever at least one of the infection
rates is sufficiently large. In contrast, our analysis of the interacting
particle system shows that, when the rate at which asymptomatic individuals
become symptomatic and the rate at which asymptomatic individuals infect nearby
individuals are small, there cannot be any epidemic even when symptomatic
individuals are highly infectious.


  Susceptible-Infected-Recovered (SIR) models have been used for decades to
understand epidemic outbreak dynamics. We develop an SIR model specifically
designed to study the effects of population behavior with respect to health and
vaccination protocols in a generic epidemic. Through a collection of
parameters, our model includes the traditional SIR components: population
birth, death, infection, recovery and vaccination rates, as well as limited
immunity.
  We first use this simple setup to compare the effects of two vaccination
schemes, one in which people are vaccinated at a rate proportional with the
population, and one in which vaccines are administered to a fraction of the
susceptible people (both of which are know strategies in real life epidemics).
We then expand on the model and the analysis by investigating how these two
vaccination schemes hold under two scenarios of population behavior: one in
which people abide by health protocols and work towards diminishing
transmission when infection is high; one in which people relax health protocols
when infection is high.
  We illustrate these two aspects (vaccination and adherence to health
protocols) act together to control the epidemic outbreak. While it is ideal
that the tow components act jointly, we also show that tight observance of
health protocols may diminish the need for vaccination in the effort to clear
or mitigate the outbreak. Conversely, an efficient vaccination strategy can
compensate for some degree of laxity in people's behavior.


  Population-size dependent branching processes (PSDBP) and controlled
branching processes (CBP) are two classes of branching processes widely used to
model biological populations that exhibit logistic growth. In this paper we
develop connections between the two, with the ultimate goal of determining when
a population is more appropriately modelled with a PSDBP or a CBP. In
particular, we state conditions for the existence of equivalent PSDBPs and
CBPs, we then consider the subclass of CBPs with deterministic control
functions (DCBPs), stating a necessary and sufficient condition for DCBP-PSDBP
equivalence. Finally, we derive an upper bound on the total variation distance
between non-equivalent DCBPs and PSDBPs with matching first and second moments
and equal initial population size, and show that under certain conditions this
bound tends to zero as the initial population size becomes large.


  Gene network information is believed to be beneficial for disease module and
pathway identification, but has not been explicitly utilized in the standard
random forest (RF) algorithm for gene expression data analysis. We investigate
the performance of a network-guided RF where the network information is
summarized into a sampling probability of predictor variables which is further
used in the construction of the RF. Our results suggest that network-guided RF
does not provide better disease prediction than the standard RF. In terms of
disease gene discovery, if disease genes form module(s), network-guided RF
identifies them more accurately. In addition, when disease status is
independent from genes in the given network, spurious gene selection results
can occur when using network information, especially on hub genes. Our
empirical analysis on two balanced microarray and RNA-Seq breast cancer
datasets from The Cancer Genome Atlas (TCGA) for classification of progesterone
receptor (PR) status also demonstrates that network-guided RF can identify
genes from PGR-related pathways, which leads to a better connected module of
identified genes.


  Many populations of animals or plants, exhibit a metapopulation structure
with close, spatially-separated subpopulations. The field of metapopulation
theory has made significant advancements since the influential Levins model.
Various modeling approaches have provided valuable insights to theoretical
Ecology. Despite extensive research on metapopulation models, there are still
challenging questions that are difficult to answer from ecological
metapopulational data or multi-patch models. Low-dimension mathematical models
offer a promising avenue to address these questions, especially for global
dynamics which have been scarcely investigated. In this study, we investigate a
two-patch metapopulation model with logistic growth and diffusion between
patches. By using analytical and numerical methods, we thoroughly analyze the
impact of diffusion on the dynamics of the metapopulation. We identify the
equilibrium points and assess their local and global stability. Furthermore, we
analytically derive the optimal diffusion rate that leads to the highest
metapopulation values. Our findings demonstrate that increased diffusion plays
a crucial role in the preservation of both subpopulations and the full
metapopulation, especially under the presence of stochastic perturbations.
Specifically, at low diffusion values, the origin is a repeller, causing orbits
starting around it to travel closely parallel to the axes. This configuration
makes the metapopulation less resilient and thus more susceptible to local and
global extinctions. However, as diffusion increases, the repeller transitions
to a saddle point, and orbits starting near the origin rapidly converge to the
unstable manifold of the saddle. This phenomenon reduces the likelihood of
stochastic extinctions and the metapopulation becomes more resilient due to
these changes in the vector field of the phase space.


  While RNA technologies hold immense therapeutic potential in a range of
applications from vaccination to gene editing, the broad implementation of
these technologies is hindered by the challenge of delivering these agents
effectively. Lipid nanoparticles have emerged as one of the most widely used
delivery agents, but their design optimization relies on laborious and costly
experimental methods. We propose to in silico optimize LNP design with machine
learning models. On a curated dataset of 622 LNPs from published studies, we
demonstrate the effectiveness of our model in predicting the transfection
efficiency of unseen LNPs, with the multilayer perceptron achieving a
classification accuracy of 98% on the test set. Our work represents a
pioneering effort in combining ML and LNP design, offering significant
potential for improving screening efficiency by computationally prioritizing
LNP candidates for experimental validation and accelerating the development of
effective mRNA delivery systems.


  Invasive pine trees pose a threat to biodiversity in a variety of Southern
Hemisphere countries, but understanding of the dynamics of invasions and the
factors that retard or accelerate spread is limited. Here, we consider the past
models of wilding pine spread and develop a new model of pine invasion. We show
that many prior models feature parameter estimates which are not biologically
supported and rely on a conjecture to obtain an asymptotic spread speed of
invasive pine populations, the main output of these models. In contrast to
prior approaches, we use partial differential equations to model an invasion.
We show that invasions are almost static for a significant period of time
before rapidly accelerating to spread at a constant rate, matching observed
behaviour in at least some field sites. Our work suggests that prior methods
for estimating invasion speeds may not accurately predict spread and are
sensitive to assumptions about the distribution of parameters. However, we
present alternative estimation methods and suggest directions for further
research.


  Multiomics data fusion integrates diverse data modalities, ranging from
transcriptomics to proteomics, to gain a comprehensive understanding of
biological systems and enhance predictions on outcomes of interest related to
disease phenotypes and treatment responses. Cooperative learning, a recently
proposed method, unifies the commonly-used fusion approaches, including early
and late fusion, and offers a systematic framework for leveraging the shared
underlying relationships across omics to strengthen signals. However, the
challenge of acquiring large-scale labeled data remains, and there are cases
where multiomics data are available but in the absence of annotated labels. To
harness the potential of unlabeled multiomcis data, we introduce
semi-supervised cooperative learning. By utilizing an "agreement penalty", our
method incorporates the additional unlabeled data in the learning process and
achieves consistently superior predictive performance on simulated data and a
real multiomics study of aging. It offers an effective solution to multiomics
data fusion in settings with both labeled and unlabeled data and maximizes the
utility of available data resources, with the potential of significantly
improving predictive models for diagnostics and therapeutics in an increasingly
multiomics world.


  Ecological communities with many species can be classified into dynamical
phases. In systems with all-to-all interactions, a phase where a fixed point is
always reached and a dynamically-fluctuating phase have been found. The
dynamics when interactions are sparse, with each species interacting with only
several others, has remained largely unexplored. Here we show that a new type
of phase appears in the phase diagram, where for the same control parameters
different communities may reach either a fixed point or a state where the
abundances of a finite subset of species fluctuate, and calculate the
probability for each outcome. These fluctuating species are organized around
short cycles in the interaction graph, and their abundances undergo large
non-linear fluctuations. We characterize the approach from this phase to a
phase with extensively many fluctuating species, and show that the probability
of fluctuations grows continuously to one as the transition is approached, and
that the number of fluctuating species diverges. This is qualitatively distinct
from the transition to extensive fluctuations coming from a fixed point phase,
which is marked by a loss of linear stability. The differences are traced back
to the emergent binary character of the dynamics when far away from short
cycles in the local fluctuations phase.


  Single-cell data integration can provide a comprehensive molecular view of
cells, and many algorithms have been developed to remove unwanted technical or
biological variations and integrate heterogeneous single-cell datasets. Despite
their wide usage, existing methods suffer from several fundamental limitations.
In particular, we lack a rigorous statistical test for whether two
high-dimensional single-cell datasets are alignable (and therefore should even
be aligned). Moreover, popular methods can substantially distort the data
during alignment, making the aligned data and downstream analysis difficult to
interpret. To overcome these limitations, we present a spectral manifold
alignment and inference (SMAI) framework, which enables principled and
interpretable alignability testing and structure-preserving integration of
single-cell data. SMAI provides a statistical test to robustly determine the
alignability between datasets to avoid misleading inference, and is justified
by high-dimensional statistical theory. On a diverse range of real and
simulated benchmark datasets, it outperforms commonly used alignment methods.
Moreover, we show that SMAI improves various downstream analyses such as
identification of differentially expressed genes and imputation of single-cell
spatial transcriptomics, providing further biological insights. SMAI's
interpretability also enables quantification and a deeper understanding of the
sources of technical confounders in single-cell data.


  The present work shows that it is possible to analytically solve a general
model to explain the transmission dynamics of SARS-CoV-2. First, the
within-host model is described, and later a between-host model, where the
coupling between them is the viral load of SARS-CoV-2. The within-host model
describes the equations involved in the life cycle of SARS-CoV-2, and also the
immune response; while that the between-Host model analyzes the dynamics of
virus spread from the original source of contagion associated with bats,
subsequently transmitted to a host, and then reaching the reservoir (Huanan
Seafood Wholesale Market in Wuhan ), until finally infecting the human
population.


  Nanobodies (Nb) are monomeric heavy-chain fragments derived from heavy-chain
only antibodies naturally found in Camelids and Sharks. Their considerably
small size (~3-4 nm; 13 kDa) and favorable biophysical properties make them
attractive targets for recombinant production. Furthermore, their unique
ability to bind selectively to specific antigens, such as toxins, chemicals,
bacteria, and viruses, makes them powerful tools in cell biology, structural
biology, medical diagnostics, and future therapeutic agents in treating cancer
and other serious illnesses. However, a critical challenge in nanobodies
production is the unavailability of nanobodies for a majority of antigens.
Although some computational methods have been proposed to screen potential
nanobodies for given target antigens, their practical application is highly
restricted due to their reliance on 3D structures. Moreover, predicting
nanobodyantigen interactions (binding) is a time-consuming and labor-intensive
task. This study aims to develop a machine-learning method to predict
Nanobody-Antigen binding solely based on the sequence data. We curated a
comprehensive dataset of Nanobody-Antigen binding and nonbinding data and
devised an embedding method based on gapped k-mers to predict binding based
only on sequences of nanobody and antigen. Our approach achieves up to 90%
accuracy in binding prediction and is significantly more efficient compared to
the widely-used computational docking technique.


  Fast screening of drug molecules based on the ligand binding affinity is an
important step in the drug discovery pipeline. Graph neural fingerprint is a
promising method for developing molecular docking surrogates with high
throughput and great fidelity. In this study, we built a COVID-19 drug docking
dataset of about 300,000 drug candidates on 23 coronavirus protein targets.
With this dataset, we trained graph neural fingerprint docking models for
high-throughput virtual COVID-19 drug screening. The graph neural fingerprint
models yield high prediction accuracy on docking scores with the mean squared
error lower than $0.21$ kcal/mol for most of the docking targets, showing
significant improvement over conventional circular fingerprint methods. To make
the neural fingerprints transferable for unknown targets, we also propose a
transferable graph neural fingerprint method trained on multiple targets. With
comparable accuracy to target-specific graph neural fingerprint models, the
transferable model exhibits superb training and data efficiency. We highlight
that the impact of this study extends beyond COVID-19 dataset, as our approach
for fast virtual ligand screening can be easily adapted and integrated into a
general machine learning-accelerated pipeline to battle future bio-threats.


  We construct a set of new epidemiological thresholds to address the general
problem of spreading and containment of a disease with influx of infected
individuals when the classic $\mathcal R_0$ is no longer meaningful. We provide
analytical properties of these indices and illustrate their usefulness in a
compartmental model of COVID-19 with data taken from Chile showing a good
predictive potential when contrasted with the recorded disease behaviour. This
approach and the associated analytical and numerical results allow us to
quantify the severity of an immigration of infectious individuals into a
community, and identification of the key parameters that are capable of
changing or reversing the spread of an infectious disease in specific models.


  Drug discovery is a highly complicated process, and it is unfeasible to fully
commit it to the recently developed molecular generation methods. Deep
learning-based lead optimization takes expert knowledge as a starting point,
learning from numerous historical cases about how to modify the structure for
better drug-forming properties. However, compared with the more established de
novo generation schemes, lead optimization is still an area that requires
further exploration. Previously developed models are often limited to resolving
one (or few) certain subtask(s) of lead optimization, and most of them can only
generate the two-dimensional structures of molecules while disregarding the
vital protein-ligand interactions based on the three-dimensional binding poses.
To address these challenges, we present a novel tool for lead optimization,
named Delete (Deep lead optimization enveloped in protein pocket). Our model
can handle all subtasks of lead optimization involving fragment growing,
linking, and replacement through a unified deleting (masking) strategy, and is
aware of the intricate pocket-ligand interactions through the geometric design
of networks. Statistical evaluations and case studies conducted on individual
subtasks demonstrate that Delete has a significant ability to produce molecules
with superior binding affinities to protein targets and reasonable
drug-likeness from given fragments or atoms. This feature may assist medicinal
chemists in developing not only me-too/me-better products from existing drugs
but also hit-to-lead for first-in-class drugs in a highly efficient manner.


  Cell proliferation and migration highly relate to normal tissue self-healing,
therefore it is highly significant for artificial controlling. Recently,
vibrational strong coupling between biomolecules and Mid-infrared (MIR) light
photons has been successfully used to modify in vitro bioreactions, neuronal
signaling and even animal behavior. However, the synergistic effects from
molecules to cells remains unclear, and the regulation of MIR on cells needs to
be explained from the molecular level. Herein, the proliferation rate and
migration capacity of fibroblasts were increased by 156% and 162.5%,
respectively, by vibratory coupling of 5.6 micrometers photons with carbonyl
groups in biomolecules. Through transcriptome sequencing analysis, the
regulatory mechanism of infrared light in 5.6 micrometers was explained from
the level of signal pathway and cell components. 5.6 micrometers optical high
power lasers can regulate cell function through vibrational strong coupling
while minimizing photothermal damage. This work not only sheds light on the
non-thermal effect on MIR light-based on wound healing, but also provides new
evidence to future frequency medicine.


  A fundamental question in complex systems is how to relate interactions
between individual components ("microscopic description") to the global
properties of the system ("macroscopic description"). Another fundamental
question is whether such a macroscopic description exists at all and how well
it describes the large-scale properties. Here, we address these questions using
as a canonical example of a self-organizing complex system - the collective
motion of desert locusts. One of the world's most devastating insect plagues
begins when flightless juvenile locusts form "marching bands". Moving through
semiarid habitats in the search for food, these bands display remarkable
coordinated motion. We investigated how well physical models can describe the
flow of locusts within a band. For this, we filmed locusts within marching
bands during an outbreak in Kenya and automatically tracked all individuals
passing through the camera frame. We first analysed the spatial topology of
nearest neighbors and found individuals to be isotropically distributed.
Despite this apparent randomness, a local order was observed in regions of high
density with a clear second neighbor peak in the radial distribution function,
akin to an ordered fluid. Furthermore, reconstructing individual locust
trajectories revealed a highly-aligned movement, consistent with the
one-dimensional version of the Toner-Tu equations, which are a generalization
of the Navier-Stokes equations for fluids, used to describe the equivalent
macroscopic fluid properties of active particles. Using this effective Toner-Tu
equation, which relates the gradient of the pressure to the acceleration, we
show that the effective "pressure" of locusts increases as a linear function of
density in segments with highest polarization. Our study thus demonstrates an
effective hydrodynamic description of flow dynamics in plague locust swarms.


  Computing the rate of evolution in spatially structured populations can be
difficult. A key quantity that describes evolutionary dynamics is the fixation
time of a single mutant with relative reproduction rate $r\ge 1$ who invades a
population of residents. We say that the fixation time is "fast" if it is at
most polynomial in terms of the population size $N$. In this work, we study
fixation times of advantageous mutants ($r>1$) and neutral mutants ($r=1$) on
directed graphs, which are defined as those graphs that have at least some
one-way connections. We obtain three main results. First, we prove that for any
directed graph the fixation time is fast, provided that $r$ is sufficiently
large. Second, we devise an efficient algorithm that gives an upper bound for
the fixation time for any graph and any $r\ge 1$. Third, we identify a broad
class of directed graphs with fast fixation times for any $r\ge 1$. This class
includes previously studied amplifiers of selection, such as Superstars and
Metafunnels. We also show that on some graphs fixation time is not a
monotonically declining function of $r$; in particular, neutral fixation can
occur faster than fixation for small selective advantages. Our results have
important algorithmic consequences and enable efficient computational
exploration of various properties of directed graphs.


  In our previous work, we demonstrated that it is feasible to perform analysis
on mutation signature data without the need for downloads or installations and
analyze individual patient data at scale without compromising privacy. Building
on this foundation, we developed a Software Development Kit (SDK) called
mSigSDK to facilitate the orchestration of distributed data processing
workflows and graphic visualization of mutational signature analysis results.
We strictly adhered to modern web computing standards, particularly the
modularization standards set by the ECMAScript ES6 framework (JavaScript
modules). Our approach allows for computation to be entirely performed by
secure delegation to the computational resources of the user's own machine
(in-browser), without any downloads or installations. The mSigSDK was developed
primarily as a companion library to the mSig Portal resource of the National
Cancer Institute Division of Cancer Epidemiology and Genetics (NIH/NCI/DCEG),
with a focus on its FAIR extensibility as components of other researchers'
computational constructs. Anticipated extensions include the programmatic
operation of other mutation signature API ecosystems such as SIGNAL and COSMIC,
advancing towards a data commons for mutational signature research (Grossman et
al., 2016).


  De novo DNA sequence assembly is based on finding paths in overlap graphs,
which is a NP-hard problem. We developed a quantum algorithm for de novo
assembly based on quantum walks in graphs. The overlap graph is partitioned
repeatedly to smaller graphs that form a hierarchical structure. We use quantum
walks to find paths in low rank graphs and a quantum algorithm that finds
Hamiltonian paths in high hierarchical rank. We tested the partitioning quantum
algorithm, as well as the quantum algorithm that finds Hamiltonian paths in
high hierarchical rank and confirmed its correct operation using Qiskit. We
developed a custom simulation for quantum walks to search for paths in low rank
graphs. The approach described in this paper may serve as a basis for the
development of efficient quantum algorithms that solve the de novo DNA assembly
problem.


  Cryo-electron microscopy (cryo-EM) is a technique for reconstructing the
3-dimensional (3D) structure of biomolecules (especially large protein
complexes and molecular assemblies). As the resolution increases to the
near-atomic scale, building protein structures de novo from cryo-EM maps
becomes possible. Recently, recognition-based de novo building methods have
shown the potential to streamline this process. However, it cannot build a
complete structure due to the low signal-to-noise ratio (SNR) problem. At the
same time, AlphaFold has led to a great breakthrough in predicting protein
structures. This has inspired us to combine fragment recognition and structure
prediction methods to build a complete structure. In this paper, we propose a
new method named FFF that bridges protein structure prediction and protein
structure recognition with flexible fitting. First, a multi-level recognition
network is used to capture various structural features from the input 3D
cryo-EM map. Next, protein structural fragments are generated using pseudo
peptide vectors and a protein sequence alignment method based on these
extracted features. Finally, a complete structural model is constructed using
the predicted protein fragments via flexible fitting. Based on our benchmark
tests, FFF outperforms the baseline methods for building complete protein
structures.


  Recent studies show that newly sampled monkeypox virus (MPXV) genomes exhibit
mutations consistent with Apolipoprotein B mRNA Editing Catalytic
Polypeptide-like3 (APOBEC3)-mediated editing, compared to MPXV genomes
collected earlier. It is unclear whether these single nucleotide polymorphisms
(SNPs) result from APOBEC3-induced editing or are a consequence of genetic
drift within one or more MPXV animal reservoirs. We develop a simple method
based on a generalization of the General-Time-Reversible (GTR) model to show
that the observed SNPs are likely the result of APOBEC3-induced editing. The
statistical features allow us to extract lineage information and estimate
evolutionary events.


  The two most common data-structures for genome indexing, FM-indices and
hash-tables, exhibit a fundamental trade-off between memory footprint and
performance. We present Ranger, a new indexing technique for nucleotide
sequences that is both memory efficient and fast. We observe that nucleotide
sequences can be represented as integer ranges and leverage a range-matching
algorithm based on neural networks to perform the lookup.
  We prototype Ranger in software and integrate it into the popular Minimap2
tool. Ranger achieves almost identical end-to-end performance as the original
Minimap2, while occupying 1.7$\times$ and 1.2$\times$ less memory for short-
and long-reads, respectively. With a limited memory capacity, Ranger achieves
up to 4.3$\times$ speedup for short reads compared to FM-Index, and up to
4.2$\times$ and 1.8$\times$ speedups for short- and long-reads, compared to
hash-tables. Ranger opens up new opportunities in the context of hardware
acceleration by reducing the memory footprint of long-seed indexes used in
state-of-the-art alignment accelerators by up to 23$\times$ which results with
3$\times$ faster alignment and negligible accuracy degradation. Moreover, its
worst case memory bandwidth and latency can be bounded in advance without the
need to inflate DRAM capacity.


  To substitute petroleum-based materials with bio-based alternatives,
microbial fermentation combined with inexpensive biomass is suggested. In this
study Saccharina latissima hydrolysate, candy-factory waste, and digestate from
full-scale biogas plant were explored as substrates for lactic acid production.
The lactic acid bacteria Enterococcus faecium, Lactobacillus plantarum, and
Pediococcus pentosaceus were tested as starter cultures. Sugars released from
seaweed hydrolysate and candy-waste were successfully utilized by the studied
bacterial strains. Additionally, seaweed hydrolysate and digestate served as
nutrient supplements supporting microbial fermentation. According to the
highest achieved relative lactic acid production, a scaled-up co-fermentation
of candy-waste and digestate was performed. Lactic acid reached a concentration
of 65.65 g/L, with 61.69% relative lactic acid production, and 1.37 g/L/hour
productivity. The findings indicate that lactic acid can be successfully
produced from low-cost industrial residues.


  Molecular docking plays a pivotal role in drug discovery and precision
medicine, enabling us to understand protein functions and advance novel
therapeutics. Here, we introduce a potential alternative solution to this
problem, the digitized-counterdiabatic quantum approximate optimization
algorithm (DC-QAOA), which utilizes counterdiabatic driving and QAOA on a
quantum computer. Our method was applied to analyze diverse biological systems,
including the SARS-CoV-2 Mpro complex with PM-2-020B, the DPP-4 complex with
piperidine fused imidazopyridine 34, and the HIV-1 gp120 complex with
JP-III-048. The DC-QAOA exhibits superior performance, providing more accurate
and biologically relevant docking results, especially for larger molecular
docking problems. Moreover, QAOA-based algorithms demonstrate enhanced hardware
compatibility in the noisy intermediate-scale quantum era, indicating their
potential for efficient implementation under practical docking scenarios. Our
findings underscore quantum computing's potential in drug discovery and offer
valuable insights for optimizing protein-ligand docking processes.


  The decline in ecological connections signifies the potential extinction of
species, which can be attributed to disruptions and alterations. The decrease
in interconnections among species reflects their susceptibility to changes. For
example, certain insects and plants that rely on exclusive interactions with a
limited number of species, or even a specific species, face the risk of
extinction if they lose these crucial connections. Currently, mining activities
pose significant harm to natural ecosystems, resulting in various adverse
environmental impacts. In this study, we utilized network science techniques to
analyze the ecosystem in a graph-based structure, aiming to conserve the
ecosystem affected by mining operations in the northern region of Scotland. The
research encompasses identifying the most vital members of the network,
establishing criteria for identifying communities within the network,
comparing, and evaluating them, using models to predict secondary extinctions
that occur when a species is removed from the network, and assessing the extent
of network damage. Our study's novelty is utilizing network science approaches
to investigate the biological data related to interactions between insects and
plants.


  Reference-guided DNA sequencing and alignment is an important process in
computational molecular biology. The amount of DNA data grows very fast, and
many new genomes are waiting to be sequenced while millions of private genomes
need to be re-sequenced. Each human genome has 3.2 B base pairs, and each one
could be stored with 2 bits of information, so one human genome would take 6.4
B bits or about 760 MB of storage (National Institute of General Medical
Sciences). Today most powerful tensor processing units cannot handle the volume
of DNA data necessitating a major leap in computing power. It is, therefore,
important to investigate the usefulness of quantum computers in genomic data
analysis, especially in DNA sequence alignment. Quantum computers are expected
to be involved in DNA sequencing, initially as parts of classical systems,
acting as quantum accelerators. The number of available qubits is increasing
annually, and future quantum computers could conduct DNA sequencing, taking the
place of classical computing systems. We present a novel quantum algorithm for
reference-guided DNA sequence alignment modeled with gate-based quantum
computing. The algorithm is scalable, can be integrated into existing classical
DNA sequencing systems and is intentionally structured to limit computational
errors. The quantum algorithm has been tested using the quantum processing
units and simulators provided by IBM Quantum, and its correctness has been
confirmed.


  The metabolic network plays a crucial role in regulating bacterial metabolism
and growth, but it is subject to inherent molecular stochasticity. Previous
studies have utilized flux balance analysis and the maximum entropy method to
predict metabolic fluxes and growth rates, while the underlying principles
governing bacterial metabolism and growth, especially the criticality
hypothesis, remain unclear. In this study, we employ a maximum entropy approach
to investigate the universality in various constraint-based metabolic networks
of Escherichia coli. Our findings reveal the existence of universal scaling
relations across different nutritional environments and metabolic network
models, similar to the universality observed in physics. By analyzing
single-cell data, we confirm that metabolism of Escherichia coli operates close
to the state with maximum Fisher information, which serves as a signature of
criticality. This critical state provides functional advantages such as high
sensitivity and long-range correlation. Moreover, we demonstrate that a
metabolic system operating at criticality takes a compromise solution between
growth and adaptation, thereby serving as a survival strategy in fluctuating
environments.


  Due to the extremely low signal-to-noise ratio (SNR) and unknown poses
(projection angles and image translation) in cryo-EM experiments,
reconstructing 3D structures from 2D images is very challenging. On top of
these challenges, heterogeneous cryo-EM reconstruction also has an additional
requirement: conformation classification. An emerging solution to this problem
is called amortized inference, implemented using the autoencoder architecture
or its variants. Instead of searching for the correct
image-to-pose/conformation mapping for every image in the dataset as in
non-amortized methods, amortized inference only needs to train an encoder that
maps images to appropriate latent spaces representing poses or conformations.
Unfortunately, standard amortized-inference-based methods with entangled latent
spaces have difficulty learning the distribution of conformations and poses
from cryo-EM images. In this paper, we propose an unsupervised deep learning
architecture called "ACE-HetEM" based on amortized inference. To explicitly
enforce the disentanglement of conformation classifications and pose
estimations, we designed two alternating training tasks in our method:
image-to-image task and pose-to-pose task. Results on simulated datasets show
that ACE-HetEM has comparable accuracy in pose estimation and produces even
better reconstruction resolution than non-amortized methods. Furthermore, we
show that ACE-HetEM is also applicable to real experimental datasets.


  We introduce AbDiffuser, an equivariant and physics-informed diffusion model
for the joint generation of antibody 3D structures and sequences. AbDiffuser is
built on top of a new representation of protein structure, relies on a novel
architecture for aligned proteins, and utilizes strong diffusion priors to
improve the denoising process. Our approach improves protein diffusion by
taking advantage of domain knowledge and physics-based constraints; handles
sequence-length changes; and reduces memory complexity by an order of magnitude
enabling backbone and side chain generation. We validate AbDiffuser in silico
and in vitro. Numerical experiments showcase the ability of AbDiffuser to
generate antibodies that closely track the sequence and structural properties
of a reference set. Laboratory experiments confirm that all 16 HER2 antibodies
discovered were expressed at high levels and that 57.1% of selected designs
were tight binders.


  This paper demonstrates the utility of organized numerical representations of
genes in research involving flat string gene formats (i.e., FASTA/FASTQ5).
FASTA/FASTQ files have several current limitations, such as their large file
sizes, slow processing speeds for mapping and alignment, and contextual
dependencies. These challenges significantly hinder investigations and tasks
that involve finding similar sequences. The solution lies in transforming
sequences into an alternative representation that facilitates easier clustering
into similar groups compared to the raw sequences themselves. By assigning a
unique vector embedding to each short sequence, it is possible to more
efficiently cluster and improve upon compression performance for the string
representations of cDNA libraries. Furthermore, through learning alternative
coordinate vector embeddings based on the contexts of codon triplets, we can
demonstrate clustering based on amino acid properties. Finally, using this
sequence embedding method to encode barcodes and cDNA sequences, we can improve
the time complexity of the similarity search by coupling vector embeddings with
an algorithm that determines the proximity of vectors in Euclidean space; this
allows us to perform sequence similarity searches in a quicker and more modular
fashion.


  The difficulty to measure or predict species community composition at fine
spatio-temporal resolution and over large spatial scales severely hampers our
ability to understand species assemblages and take appropriate conservation
measures. Despite the progress in species distribution modeling (SDM) over the
past decades, SDM have just begun to integrate high resolution remote sensing
data and their predictions are still entailed by many biases due to
heterogeneity of the available biodiversity observations, most often
opportunistic presence only data. We designed a European scale dataset covering
around ten thousand plant species to calibrate and evaluate SDM predictions of
species composition in space and time at high spatial resolution (~ten meters),
and their spatial transferability. For model training, we extracted and
harmonized five million heterogeneous presence-only records from selected GBIF
datasets and 6 thousand exhaustive presence-absence surveys both sampled during
2017-2021. We associated species observations to diverse environmental rasters
classically used in SDMs, as well as to 10 m resolution RGB and Near-Infra-Red
satellite images and 20 years-time series of climatic variables and satellite
point values. The evaluation dataset is based on 22 thousand standardized
presence-absence surveys separated from the training set with a spatial block
hold out procedure. The GeoLifeCLEF 2023 dataset is open access and the first
benchmark for researchers aiming to improve the prediction of plant species
composition at a very fine spatial grain and at continental scale. It is a
space to explore new ways of combining massive and diverse species observations
and environmental information at various scales. Innovative AI-based
approaches, in particular, should be among the most interesting methods to
experiment with on the GeoLifeCLEF 2023 dataset.


  The capacity to identify and analyze protein-protein interactions, along with
their internal modular organization, plays a crucial role in comprehending the
intricate mechanisms underlying biological processes at the molecular level. We
can learn a lot about the structure and dynamics of these interactions by using
network analysis. We can improve our understanding of the biological roots of
disease pathogenesis by recognizing network communities. This knowledge, in
turn, holds significant potential for driving advancements in drug discovery
and facilitating personalized medicine approaches for disease treatment. In
this study, we aimed to uncover the communities within the variant B.1.1.529
(Omicron virus) using two proposed novel algorithm (ABCDE and ALCDE) and four
widely recognized algorithms: Girvan-Newman, Louvain, Leiden, and Label
Propagation algorithm. Each of these algorithms has established prominence in
the field and offers unique perspectives on identifying communities within
complex networks. We also compare the networks by the global properties,
statistic summary, subgraph count, graphlet and validate by the modulaity. By
employing these approaches, we sought to gain deeper insights into the
structural organization and interconnections present within the Omicron virus
network.


  This paper proposes pandemic mitigation vaccination policies for Newfoundland
and Labrador (NL) based on two compact mixed integer programming (MIP) models
of the distance-based critical node detection problem (DCNDP). Our main focus
is on two variants of the DCNDP that seek to minimize the number of connections
with lengths of at most one (1-DCNDP) and two (2-DCNDP). A polyhedral study for
the 1-DCNDP is conducted, and new aggregated inequalities are provided for the
2-DCNDP. The computational experiments show that the 2-DCNDP with aggregated
inequalities outperforms the one with disaggregated inequalities for graphs
with a density of at least 0.5%. We also study the strategic vaccine allocation
problem as a real-world application of the DCNDP and conduct a set of
computational experiments on a simulated contact network of NL. Our
computational results demonstrate that the DCNDP-based strategies can have a
better performance in comparison with the real-world strategies implemented
during COVID-19.


  Multiple sequence alignments (MSAs) of proteins encode rich biological
information and have been workhorses in bioinformatic methods for tasks like
protein design and protein structure prediction for decades. Recent
breakthroughs like AlphaFold2 that use transformers to attend directly over
large quantities of raw MSAs have reaffirmed their importance. Generation of
MSAs is highly computationally intensive, however, and no datasets comparable
to those used to train AlphaFold2 have been made available to the research
community, hindering progress in machine learning for proteins. To remedy this
problem, we introduce OpenProteinSet, an open-source corpus of more than 16
million MSAs, associated structural homologs from the Protein Data Bank, and
AlphaFold2 protein structure predictions. We have previously demonstrated the
utility of OpenProteinSet by successfully retraining AlphaFold2 on it. We
expect OpenProteinSet to be broadly useful as training and validation data for
1) diverse tasks focused on protein structure, function, and design and 2)
large-scale multimodal machine learning research.


  Retrosynthesis consists of breaking down a chemical compound recursively
step-by-step into molecular precursors until a set of commercially available
molecules is found with the goal to provide a synthesis route. Its two primary
research directions, single-step retrosynthesis prediction, which models the
chemical reaction logic, and multi-step synthesis planning, which tries to find
the correct sequence of reactions, are inherently intertwined. Still, this
connection is not reflected in contemporary research. In this work, we combine
these two major research directions by applying multiple single-step
retrosynthesis models within multi-step synthesis planning and analyzing their
impact using public and proprietary reaction data. We find a disconnection
between high single-step performance and potential route-finding success,
suggesting that single-step models must be evaluated within synthesis planning
in the future. Furthermore, we show that the commonly used single-step
retrosynthesis benchmark dataset USPTO-50k is insufficient as this evaluation
task does not represent model performance and scalability on larger and more
diverse datasets. For multi-step synthesis planning, we show that the choice of
the single-step model can improve the overall success rate of synthesis
planning by up to +28% compared to the commonly used baseline model. Finally,
we show that each single-step model finds unique synthesis routes, and differs
in aspects such as route-finding success, the number of found synthesis routes,
and chemical validity, making the combination of single-step retrosynthesis
prediction and multi-step synthesis planning a crucial aspect when developing
future methods.


  We present an R-package developed for modeling of biochemical reaction
networks, uncertainty quantification (UQ) and sensitivity analysis (SA).
Estimating parameters and quantifying their uncertainty (and resulting
prediction uncertainty), is required for data-driven systems biology modeling.
Sampling methods need to be efficient when confronted with high-dimensional,
correlated parameter distributions. We have developed the UQSA package to be
fast for this problem class and work well with other tools for modelling. We
aim for simplicity, and part of that is our use of the SBtab format for the
unified storage of model and data. Our tool-set is modular enough, that parts
can be replaced. We use intermediate formats that are not hidden from the user
to make this feasible. UQ is performed through Markov chain Monte Carlo (MCMC)
sampling in an Approximate Bayesian Computation (ABC) setting. This can be
followed by a variance-decomposition based global sensitivity analysis. If
needed, complex parameter distributions can be described, evaluated, and
sampled from, with the help of Vine-copulas that are available in R. This
approach is especially useful when new experimental data become available, and
a previously calibrated model needs to be updated.
  Implementation: R is a high-level language and allows the use of
sophisticated statistical methods. The ode solver we used is written in C
(gsl_odeiv2, interface to R is ours). We use the SBtab tabular format for the
model description, as well as the data and an event system to be able to model
inputs frequently encountered in systems biology and neuroscience. The code has
been tested on one node with 256 cores of a computing cluster, but smaller
examples are included in the repository that can be run on a laptop.
  Source code: https://github.com/icpm-kth/uqsa


  Degree correlation plays a crucial role in studying network structures;
however, its varied forms pose challenges to understanding its impact on
network dynamics. This study devised a method that uses eigenvalue
decomposition to characterize degree correlations. Additionally, the
applicability of this method was demonstrated by approximating the basic and
type reproduction numbers in an epidemic network model. The findings elucidate
the interplay between degree correlations and epidemic behavior, thus
contributing to a deeper understanding of complex networks and their dynamics.


  Transcription factor combinations determine gene locus activity and thereby
cell identity. However, the precise link between concentrations of such
activating transcription factors and target-gene activity is ambiguous. Here we
investigate this link for the gap gene dependent activation of the even-skipped
(eve) locus in the Drosophila embryo. We simultaneously measure the
spatiotemporal gap gene concentrations in hemizygous and homozygous gap
mutants, and link these to eve activity. Although changes in expression extend
well beyond the genetically manipulated gene, nearly all expression
alternations approximate the canonical combinations of activating levels in
wild-type, sometimes necessitating pattern shifts. Expression levels that
diverge from the wild-type repertoire still drive locus activation. Specific
stripes in the homozygous mutants show partial penetrance, justifying their
renown variable phenotypes. However, all eve stripes appear at highly
reproducible positions, even though a broader span of gap gene expression
levels activates eve. Our results suggest a correction capacity of the gap gene
network and set constraints on the activity of multi-enhancer gene loci.


  Construction of in vitro vascular models is of great significance to various
biomedical research, such as pharmacokinetics and hemodynamics, thus is an
important direction in tissue engineering. In this work, a standing surface
acoustic wave field was constructed to spatially arrange suspended endothelial
cells into a designated patterning. The cell patterning was maintained after
the acoustic field was withdrawn by the solidified hydrogel. Then, interstitial
flow was provided to activate vessel tube formation. Thus, a functional
vessel-on-a-chip was engineered with specific vessel geometry. Vascular
function, including perfusability and vascular barrier function, was
characterized by beads loading and dextran diffusion, respectively. A
computational atomistic simulation model was proposed to illustrate how solutes
cross vascular lipid bilayer. The reported acoustofluidic methodology is
capable of facile and reproducible fabrication of functional vessel network
with specific geometry. It is promising to facilitate the development of both
fundamental research and regenerative therapy.


  A complete time-parameterized statistical model quantifying the divergent
evolution of protein structures in terms of the patterns of conservation of
their secondary structures is inferred from a large collection of protein 3D
structure alignments. This provides a better alternative to time-parameterized
sequence-based models of protein relatedness, that have clear limitations
dealing with twilight and midnight zones of sequence relationships. Since
protein structures are far more conserved due to the selection pressure
directly placed on their function, divergence time estimates can be more
accurate when inferred from structures. We use the Bayesian and
information-theoretic framework of Minimum Message Length to infer a
time-parameterized stochastic matrix (accounting for perturbed structural
states of related residues) and associated Dirichlet models (accounting for
insertions and deletions during the evolution of protein domains). These are
used in concert to estimate the Markov time of divergence of tertiary
structures, a task previously only possible using proxies (like RMSD). By
analyzing one million pairs of homologous structures, we yield a relationship
between the Markov divergence time of structures and of sequences. Using these
inferred models and the relationship between the divergence of sequences and
structures, we demonstrate a competitive performance in secondary structure
prediction against neural network architectures commonly employed for this
task. The source code and supplementary information are downloadable from
\url{http://lcb.infotech.monash.edu.au/sstsum}.


  Modeling the properties of chemical mixtures is a difficult but important
part of any modeling process intended to be applicable to the often messy and
impure phenomena of everyday life, including food and environmental safety,
healthcare, etc. Part of this difficulty stems from the increased complexity of
designing suitable model validation schemes for mixture data, a fact which has
been elucidated in previous work only in the case of binary mixture models. We
extend these previously defined validation strategies for QSAR modeling of
binary mixtures to the more complex case of general, $N$-ary mixtures and argue
that these strategies are applicable to many modeling tasks beyond simple
chemical mixtures. Additionally, we propose a method of establishing a baseline
model performance for each mixture dataset to be in used in model selection
comparisons. This baseline is intended to account for the statistical
dependence generically present between the properties of mixtures that share
constituents. We contend that without such a baseline, estimates of model
performance can be dramatically overestimated, and we demonstrate this with
multiple case studies using real and simulated data.


  We present a general method for rendering representations of multi-stranded
DNA complexes from textual descriptions into 2D diagrams. The complexes can be
arbitrarily pseudoknotted, and if a planar rendering is possible, the method
will determine one in time which is essentially linear in the size of the
textual description. (That is, except for a final stochastic fine-tuning step.)
If a planar rendering is not possible, the method will compute a visually
pleasing approximate rendering in quadratic time. Examples of diagrams produced
by the method are presented in the paper.


  We study two non-Markovian gene-expression models in which protein production
is a stochastic process with a fat-tailed non-exponential distribution of
waiting times. For both models, we find two distinct scaling regimes separated
by an exponentially long time, proportional to the mean first passage time
(MFPT) to a ground state of the dynamics, from which the system can only exit
via a non-exponential reaction. At times shorter than the MFPT the dynamics are
stationary and ergodic while at times longer than the MFPT the dynamics are
nonstationary and nonergodic. The MFPT to the ground state is shown to directly
affect the average population sizes and we postulate that the transition to
nonergodicity is a universal feature of models with a non-Markovian state.


  Environments such as shopping malls, airports, or hospital emergency
departments often experience crowding, with many people simultaneously
requesting service. Crowding is highly noisy, with sudden overcrowding
"spikes". Past research has either focused on average behavior or used
context-specific non-generalizable models. Here we show that a stochastic
population model, previously applied to a broad range of natural phenomena, can
aptly describe hospital emergency-department crowding, using data from
five-year minute-by-minute emergency-department records. The model provides
reliable forecasting of the crowding distribution. Overcrowding is highly
sensitive to the patient arrival-flux and length-of-stay: a 10% increase in
arrivals triples the probability of overcrowding events. Expediting patient
exit-rate to shorten the typical length-of-stay by just 20 minutes (8.5%)
reduces severe overcrowding events by 50%. Such forecasting is crucial in
prevention and mitigation of breakdown events. Our results demonstrate that
despite its high volatility, crowding follows a dynamic behavior common to many
natural systems.


  The ancestral maximum-likelihood and phylogeography problems are two
fundamental problems involving evolutionary studies. The ancestral
maximum-likelihood problem involves identifying a rooted tree alongside
internal node sequences that maximizes the probability of observing a given set
of sequences as leaves. The phylogeography problem extends the ancestral
maximum-likelihood problem to incorporate geolocation of leaf and internal
nodes. While a constant factor approximation algorithm has been established for
the ancestral maximum-likelihood problem concerning two-state sequences, no
such algorithm has been devised for any generalized instances of the problem.
In this paper, we focus on a generalization of the two-state model, the time
reversible Markov evolutionary models for sequences and geolocations. Under
this evolutionary model, we present a $2\log_2 k $-approximation algorithm,
where $k$ is the number of input samples, addressing both the ancestral
maximum-likelihood and phylogeography problems. This is the first approximation
algorithm for the phylogeography problem. Furthermore, we show how to apply the
algorithm on popular evolutionary models like generalized time-reversible (GTR)
model and its specialization Jukes and Cantor 69 (JC69).


  In the first part of this paper, we review old and new results about the
influence of host population heterogeneity on (various characteristics of)
epidemic outbreaks. In the second part we highlight a modelling issue that so
far has received little attention: how do contact patterns, and hence
transmission opportunities, depend on the size and the composition of the host
population? Without any claim on completeness, we offer a range of potential
(quasi-mechanistic) submodels. The overall aim of the paper is to describe the
state-of-the-art and to catalyse new work.


  This study presents a thorough analysis of an eco-epidemiological model that
integrates infectious diseases in prey, prey aggregation, and the dual fear
effect induced by predators. We establish criteria for determining the
existence of equilibrium points, which carry substantial biological
significance. We establish the conditions for the occurrence of Hopf,
saddle-node, and transcritical bifurcations by employing fear parameters as key
bifurcation parameters. Furthermore, through numerical simulations, we
demonstrate the occurrence of multiple zero-Hopf (ZH) and saddle-node
transcritical (SNTC) bifurcations around the endemic steady states by varying
specific key parameters across the two-parametric plane. We demonstrate that
the introduction of predator-induced fear, which hinders the growth rate of
susceptible prey, can lead to the finite time extinction of an initially stable
susceptible prey population. Finally, we discuss management strategies aimed at
regulating disease transmission, focusing on fear-based interventions and
selective predation via predator attack rate on infectious prey.


  We present ensemblQueryR, a package providing an R interface to the Ensembl
REST API that facilitates flexible, fast, user-friendly and R workflow
integrable querying of Ensembl REST API linkage disequilibrium (LD) endpoints,
optimised for high-throughput querying. ensemblQueryR achieves this through
functions that are intuitive and amenable to custom code integration, use of
familiar R object types as inputs and outputs, code optimisation and optional
parallelisation functionality. For each LD endpoint, ensemblQueryR provides two
functions, permitting both single-query and multi-query modes of operation. The
multi-query functions are optimised for large query sizes and provide optional
parallelisation to leverage available computational resources and minimise
processing time. We demonstrate that ensemblQueryR has improved performance in
terms of random access memory (RAM) usage and speed, delivering a 10-fold speed
increase over analogous software whilst using a third of the RAM. Finally,
ensemblQueryR is near-agnostic to operating system and computational
architecture through availability of Docker and singularity images, making this
tool widely accessible to the scientific community.


  Large language models have made significant strides in natural language
processing, enabling innovative applications in molecular science by processing
textual representations of molecules. However, most existing language models
cannot capture the rich information with complex molecular structures or
images. In this paper, we introduce GIT-Mol, a multi-modal large language model
that integrates the Graph, Image, and Text information. To facilitate the
integration of multi-modal molecular data, we propose GIT-Former, a novel
architecture that is capable of aligning all modalities into a unified latent
space. We achieve a 5%-10% accuracy increase in properties prediction and a
20.2% boost in molecule generation validity compared to the baselines. With the
any-to-language molecular translation strategy, our model has the potential to
perform more downstream tasks, such as compound name recognition and chemical
reaction prediction.


  The birth of ChatGPT, a cutting-edge language model-based chatbot developed
by OpenAI, ushered in a new era in AI. However, due to potential pitfalls, its
role in rigorous scientific research is not clear yet. This paper vividly
showcases its innovative application within the field of drug discovery.
Focused specifically on developing anti-cocaine addiction drugs, the study
employs GPT-4 as a virtual guide, offering strategic and methodological
insights to researchers working on generative models for drug candidates. The
primary objective is to generate optimal drug-like molecules with desired
properties. By leveraging the capabilities of ChatGPT, the study introduces a
novel approach to the drug discovery process. This symbiotic partnership
between AI and researchers transforms how drug development is approached.
Chatbots become facilitators, steering researchers towards innovative
methodologies and productive paths for creating effective drug candidates. This
research sheds light on the collaborative synergy between human expertise and
AI assistance, wherein ChatGPT's cognitive abilities enhance the design and
development of potential pharmaceutical solutions. This paper not only explores
the integration of advanced AI in drug discovery but also reimagines the
landscape by advocating for AI-powered chatbots as trailblazers in
revolutionizing therapeutic innovation.


  This paper develops low-dimensional moment methods for advective problems on
networks of domains. The evolution of a density function is described by a
linear advection-diffusion-reaction equation on each domain, combined via
advective flux coupling across domains in the network graph. The PDEs'
coefficients vary in time and across domains but they are fixed along each
domain. As a result, the solution on each domain is frequently close to a
Gaussian that moves, decays, and widens. For that reason, this work studies
moment methods that track only three degrees of freedom per domain -- in
contrast to traditional PDE discretization methods that tend to require many
more variables per domain. A simple ODE-based moment method is developed, as
well as an asymptotic-preserving scheme. We apply the methodology to an
application that models the life cycle of forest pests that undergo different
life stages and developmental pathways. The model is calibrated for the spotted
lanternfly, an invasive species present in the Eastern USA. We showcase that
the moment method, despite its significant low-dimensionality, can successfully
reproduce the prediction of the pest's establishment potential, compared to
much higher-dimensional computational approaches.


  The response threshold model explains the emergence of division of labor
(i.e., task specialization) in an unstructured population by assuming that the
individuals have different propensities to work on different tasks. The
incentive to attend to a particular task increases when the task is left
unattended and decreases when individuals work on it. Here we derive mean-field
equations for the stimulus dynamics and show that they exhibit complex
attractors through period-doubling bifurcation cascades when the noise
disrupting the thresholds is small. In addition, we show how the fixed
threshold can be set to ensure specialization in both the transient and
equilibrium regimes of the stimulus dynamics. However, a complete explanation
of the emergence of division of labor requires that we address the question of
where the threshold variation comes from, starting from a homogeneous
population. We then study a structured population scenario, where the
population is divided into a large number of independent groups of equal size,
and the fitness of a group is proportional to the weighted mean work performed
on the tasks during a fixed period of time. Using a winner-take-all strategy to
model group competition and assuming an initial homogeneous metapopulation, we
find that a substantial fraction of workers specialize in each task, without
the need to penalize task switching.


  Predicting which proteins interact together from amino-acid sequences is an
important task. We develop a method to pair interacting protein sequences which
leverages the power of protein language models trained on multiple sequence
alignments, such as MSA Transformer and the EvoFormer module of AlphaFold. We
formulate the problem of pairing interacting partners among the paralogs of two
protein families in a differentiable way. We introduce a method called DiffPALM
that solves it by exploiting the ability of MSA Transformer to fill in masked
amino acids in multiple sequence alignments using the surrounding context. MSA
Transformer encodes coevolution between functionally or structurally coupled
amino acids. We show that it captures inter-chain coevolution, while it was
trained on single-chain data, which means that it can be used
out-of-distribution. Relying on MSA Transformer without fine-tuning, DiffPALM
outperforms existing coevolution-based pairing methods on difficult benchmarks
of shallow multiple sequence alignments extracted from ubiquitous prokaryotic
protein datasets. It also outperforms an alternative method based on a
state-of-the-art protein language model trained on single sequences. Paired
alignments of interacting protein sequences are a crucial ingredient of
supervised deep learning methods to predict the three-dimensional structure of
protein complexes. DiffPALM substantially improves the structure prediction of
some eukaryotic protein complexes by AlphaFold-Multimer, without significantly
deteriorating any of those we tested. It also achieves competitive performance
with using orthology-based pairing.


  In this paper we study invasion probabilities and invasion times of
cooperative parasites spreading in spatially structured host populations. The
spatial structure of the host population is given by a random geometric graph
on $[0,1]^n$, $n\in \mathbb{N}$, with a Poisson($N$)-distributed number of
vertices and in which vertices are connected over an edge when they have a
distance of at most $r_N\in \Theta\left(N^{\frac{\beta-1}{n}}\right)$ for some
$0<\beta<1$ and $N\rightarrow \infty$. At a host infection many parasites are
generated and parasites move along edges to neighbouring hosts. We assume that
parasites have to cooperate to infect hosts, in the sense that at least two
parasites need to attack a host simultaneously. We find lower and upper bounds
on the invasion probability of the parasites in terms of survival probabilities
of branching processes with cooperation. Furthermore, we characterize the
asymptotic invasion time.
  An important ingredient of the proofs is a comparison with infection dynamics
of cooperative parasites in host populations structured according to a complete
graph, i.e. in well-mixed host populations. For these infection processes we
can show that invasion probabilities are asymptotically equal to survival
probabilities of branching processes with cooperation.
  Furthermore, we build in the proofs on techniques developed in [BP22], where
an analogous invasion process has been studied for host populations structured
according to a configuration model.
  We substantiate our results with simulations.


  Deep generative models for structure-based drug design (SBDD), where molecule
generation is conditioned on a 3D protein pocket, have received considerable
interest in recent years. These methods offer the promise of higher-quality
molecule generation by explicitly modelling the 3D interaction between a
potential drug and a protein receptor. However, previous work has primarily
focused on the quality of the generated molecules themselves, with limited
evaluation of the 3D molecule \emph{poses} that these methods produce, with
most work simply discarding the generated pose and only reporting a "corrected"
pose after redocking with traditional methods. Little is known about whether
generated molecules satisfy known physical constraints for binding and the
extent to which redocking alters the generated interactions. We introduce
PoseCheck, an extensive analysis of multiple state-of-the-art methods and find
that generated molecules have significantly more physical violations and fewer
key interactions compared to baselines, calling into question the implicit
assumption that providing rich 3D structure information improves molecule
complementarity. We make recommendations for future research tackling
identified failure modes and hope our benchmark can serve as a springboard for
future SBDD generative modelling work to have a real-world impact.


  Scaffold hopping is a drug discovery strategy to generate new chemical
entities by modifying the core structure, the \emph{scaffold}, of a known
active compound. This approach preserves the essential molecular features of
the original scaffold while introducing novel chemical elements or structural
features to enhance potency, selectivity, or bioavailability. However, there is
currently a lack of generative models specifically tailored for this task,
especially in the pocket-conditioned context. In this work, we present
DiffHopp, a conditional E(3)-equivariant graph diffusion model tailored for
scaffold hopping given a known protein-ligand complex.


  This paper proposes a sign gradient descent (SGD) algorithm for predicting
the three-dimensional folded protein molecule structures under the kinetostatic
compliance method (KCM). In the KCM framework, which can be used to simulate
the range of motion of peptide-based nanorobots/nanomachines, protein molecules
are modeled as a large number of rigid nano-linkages that form a kinematic
mechanism under motion constraints imposed by chemical bonds while folding
under the kinetostatic effect of nonlinear interatomic force fields. In a
departure from the conventional successive kinetostatic fold compliance
framework, the proposed SGD-based iterative algorithm in this paper results in
convergence to the local minima of the free energy of protein molecules
corresponding to their final folded conformations in a faster and more robust
manner. KCMbased folding dynamics simulations of the backbone chains of protein
molecules demonstrate the effectiveness of the proposed algorithm.


  Biological functions stem from coordinated interactions among proteins,
nucleic acids and small molecules. Mass spectrometry technologies for reliable,
high throughput single-cell proteomics will add a new modality to genomics and
enable data-driven modeling of the molecular mechanisms coordinating proteins
and nucleic acids at single-cell resolution. This promising potential requires
estimating the reliability of measurements and computational analysis so that
models can distinguish biological regulation from technical artifacts. We
highlight different measurement modes that can support single-cell
proteogenomic analysis and how to estimate their reliability. We then discuss
approaches for developing both abstract and mechanistic models that aim to
biologically interpret the measured differences across modalities, including
specific applications to directed stem cell differentiation and to inferring
protein interactions in cancer cells from the buffing of DNA copy-number
variations. Single-cell proteogenomic data will support mechanistic models of
direct molecular interactions that will provide generalizable and predictive
representations of biological systems.


  An increasingly common viewpoint is that protein dynamics data sets reside in
a non-linear subspace of low conformational energy. Ideal data analysis tools
for such data sets should therefore account for such non-linear geometry. The
Riemannian geometry setting can be suitable for a variety of reasons. First, it
comes with a rich structure to account for a wide range of geometries that can
be modelled after an energy landscape. Second, many standard data analysis
tools initially developed for data in Euclidean space can also be generalised
to data on a Riemannian manifold. In the context of protein dynamics, a
conceptual challenge comes from the lack of a suitable smooth manifold and the
lack of guidelines for constructing a smooth Riemannian structure based on an
energy landscape. In addition, computational feasibility in computing geodesics
and related mappings poses a major challenge. This work considers these
challenges. The first part of the paper develops a novel local approximation
technique for computing geodesics and related mappings on Riemannian manifolds
in a computationally feasible manner. The second part constructs a smooth
manifold of point clouds modulo rigid body group actions and a Riemannian
structure that is based on an energy landscape for protein conformations. The
resulting Riemannian geometry is tested on several data analysis tasks relevant
for protein dynamics data. It performs exceptionally well on coarse-grained
molecular dynamics simulated data. In particular, the geodesics with given
start- and end-points approximately recover corresponding molecular dynamics
trajectories for proteins that undergo relatively ordered transitions with
medium sized deformations. The Riemannian protein geometry also gives
physically realistic summary statistics and retrieves the underlying dimension
even for large-sized deformations within seconds on a laptop.


  De novo assembly is the process of reconstructing the genome sequence of an
organism from sequencing reads. Genome sequences are essential to biology, and
assembly has been a central problem in bioinformatics for four decades. Until
recently, genomes were typically assembled into fragments of a few megabases at
best but technological advances in long-read sequencing now enable near
complete chromosome-level assembly, also known as telomere-to-telomere
assembly, for many organisms. Here we review recent progress on assembly
algorithms and protocols. We focus on how to derive near telomere-to-telomere
assemblies and discuss potential future developments.


  The prediction of protein 3D structure from amino acid sequence is a
computational grand challenge in biophysics, and plays a key role in robust
protein structure prediction algorithms, from drug discovery to genome
interpretation. The advent of AI models, such as AlphaFold, is revolutionizing
applications that depend on robust protein structure prediction algorithms. To
maximize the impact, and ease the usability, of these novel AI tools we
introduce APACE, AlphaFold2 and advanced computing as a service, a novel
computational framework that effectively handles this AI model and its TB-size
database to conduct accelerated protein structure prediction analyses in modern
supercomputing environments. We deployed APACE in the Delta supercomputer, and
quantified its performance for accurate protein structure predictions using
four exemplar proteins: 6AWO, 6OAN, 7MEZ, and 6D6U. Using up to 200 ensembles,
distributed across 50 nodes in Delta, equivalent to 200 A100 NVIDIA GPUs, we
found that APACE is up to two orders of magnitude faster than off-the-shelf
AlphaFold2 implementations, reducing time-to-solution from weeks to minutes.
This computational approach may be readily linked with robotics laboratories to
automate and accelerate scientific discovery.


  Tensor factorizations (TF) are powerful tools for the efficient
representation and analysis of multidimensional data. However, classic TF
methods based on maximum likelihood estimation underperform when applied to
zero-inflated count data, such as single-cell RNA sequencing (scRNA-seq) data.
Additionally, the stochasticity inherent in TFs results in factors that vary
across repeated runs, making interpretation and reproducibility of the results
challenging. In this paper, we introduce Zero Inflated Poisson Tensor
Factorization (ZIPTF), a novel approach for the factorization of
high-dimensional count data with excess zeros. To address the challenge of
stochasticity, we introduce Consensus Zero Inflated Poisson Tensor
Factorization (C-ZIPTF), which combines ZIPTF with a consensus-based
meta-analysis. We evaluate our proposed ZIPTF and C-ZIPTF on synthetic
zero-inflated count data and synthetic and real scRNA-seq data. ZIPTF
consistently outperforms baseline matrix and tensor factorization methods in
terms of reconstruction accuracy for zero-inflated data. When the probability
of excess zeros is high, ZIPTF achieves up to $2.4\times$ better accuracy.
Additionally, C-ZIPTF significantly improves the consistency and accuracy of
the factorization. When tested on both synthetic and real scRNA-seq data, ZIPTF
and C-ZIPTF consistently recover known and biologically meaningful gene
expression programs.


  The migratory dynamics of cells can be influenced by the complex
micro-environment through which they move. It remains unclear how the motility
machinery of confined cells responds and adapts to their micro-environment.
Here, we propose a biophysical mechanism for a geometry-dependent coupling
between cellular protrusions and the nucleus that leads to directed migration.
We apply our model to geometry-guided cell migration to obtain insights into
the origin of directed migration on asymmetric adhesive micro-patterns and the
polarization enhancement of cells observed under strong confinement.
Remarkably, for cells that can choose between channels of different size, our
model predicts an intricate dependence for cellular decision making as a
function of the two channel widths, which we confirm experimentally.


  Distributive phosphorylation is recurrent motif in intracellular signaling
and control. It is either sequential or cyclic. Sequential distributive
phosphorylation has been extensively studied and it has been shown that the
catalytic constants of kinase and phosphatase enable multistationarity and
hence bistability. Here we show for cyclic distributive phosphorylation that if
its catalytic constants satisfy the very same inequality as in the sequential
case, then Hopf bifurcations and hence sustained oscillations can occur. Hence
we conclude that in distributive phosphorylation (sequential and distributive)
it is the catalytic constants that enable nontrivial dynamics.
  In fact, if the rate constants in a network of cyclic distributive
phosphorylation are such that Hopf bifurcations and sustained oscillations
occur, then a network of sequential distributive phosphorylation with the same
rate constants will show multistationarity. For cyclic distributive
phosphorylation we further describe a procedure to generate rate constants
where Hopf bifurcations and hence sustained oscillations occur. This may allow
an efficient sampling of oscillatory regions in parameter space.
  Our analysis is greatly simplified by the fact that it is possible to reduce
the network of cyclic distributive phosphorylation to what we call a network
with a single extreme ray. We summarize key properties of these networks.


  The Research & Development (R&D) phase of drug development is a lengthy and
costly process. To revolutionize this process, we introduce our new concept
QMLS to shorten the whole R&D phase to three to six months and decrease the
cost to merely fifty to eighty thousand USD. For Hit Generation, Machine
Learning Molecule Generation (MLMG) generates possible hits according to the
molecular structure of the target protein while the Quantum Simulation (QS)
filters molecules from the primary essay based on the reaction and binding
effectiveness with the target protein. Then, For Lead Optimization, the
resultant molecules generated and filtered from MLMG and QS are compared, and
molecules that appear as a result of both processes will be made into dozens of
molecular variations through Machine Learning Molecule Variation (MLMV), while
others will only be made into a few variations. Lastly, all optimized molecules
would undergo multiple rounds of QS filtering with a high standard for reaction
effectiveness and safety, creating a few dozen pre-clinical-trail-ready drugs.
This paper is based on our first paper, where we pitched the concept of machine
learning combined with quantum simulations. In this paper we will go over the
detailed design and framework of QMLS, including MLMG, MLMV, and QS.


  The Kelvin-Planck statement of the Second Law of Thermodynamics is a
stricture on the nature of heat receipt by any body suffering a cyclic process.
It makes no mention of temperature or of entropy. Beginning with a
Kelvin-Planck statement of the Second Law, we show that entropy and temperature
-- in particular, existence of functions that relate the local specific entropy
and thermodynamic temperature to the local state in a material body -- emerge
immediately and simultaneously as consequences of the Hahn-Banach Theorem.
Existence of such functions of state requires no stipulation that their domains
be restricted to equilibrium states. Further properties, including uniqueness,
are addressed in a companion paper.


  The recovery of many groundfish stocks throughout the Northwest Atlantic has
been impeded by elevated natural (i.e., non-fishing) mortality (M) among
older/larger individuals. The causes of elevated mortality are not well known,
though predation by rapidly growing grey seal herds and unreported fishing are
thought to be possible drivers of mortality for Atlantic Cod (Gadus morhua) on
the Western Scotian Shelf and in the Bay of Fundy (known as "4X5Y cod") and in
nearby ecosystems. We developed a statistical catch-at-age model for 4X5Y cod
that accounted for both grey seal predation and estimated bycatch/discards to
evaluate the degree to which either of these factors may influence cod
mortality. The model was fit over a range of predation and discarding scenarios
to account for uncertainties and a lack of data for these processes. We found
that most cod M remained unaccounted for unless cod comprised a large
proportion (>0.45) of the grey seal diet by weight. If the reported bycatch
estimates are taken as accurate, then the magnitude of cod discards from
non-directed fisheries was minor, though these estimates are highly uncertain.


  The typically rugged nature of molecular free energy landscapes can frustrate
efficient sampling of the thermodynamically relevant phase space due to the
presence of high free energy barriers. Enhanced sampling techniques can improve
phase space exploration by accelerating sampling along particular collective
variables (CVs). A number of techniques exist for data-driven discovery of CVs
parameterizing the important large scale motions of the system. A challenge to
CV discovery is learning CVs invariant to symmetries of the molecular system,
frequently rigid translation, rigid rotation, and permutational relabeling of
identical particles. Of these, permutational invariance have proved a
persistent challenge in frustrating the the data-driven discovery of
multi-molecular CVs in systems of self-assembling particles and
solvent-inclusive CVs for solvated systems. In this work, we integrate
Permutation Invariant Vector (PIV) featurizations with autoencoding neural
networks to learn nonlinear CVs invariant to translation, rotation, and
permutation, and perform interleaved rounds of CV discovery and enhanced
sampling to iteratively expand sampling of configurational phase space and
obtain converged CVs and free energy landscapes. We demonstrate the
Permutationally Invariant Network for Enhanced Sampling (PINES) approach in
applications to the self-assembly of a 13-atom Argon cluster,
association/dissociation of a NaCl ion pair in water, and hydrophobic collapse
of a C45H92 n-pentatetracontane polymer chain. We make the approach freely
available as a new module within the PLUMED2 enhanced sampling libraries.


  Undocumented species represent one of the largest hurdles for conservation
efforts due to the uncertainty they introduce into conservation planning. Until
the distribution of earth's biodiversity is better understood, substantial
conjecture will continue to be required for protecting species from
anthropogenic extinction. Therefore, we developed a novel approach for
identifying regions with promising biodiscovery prospects, linked to
integrative conservation priorities, which we illustrate using amphibians. Our
approach builds on previous estimates of biodiscovery priorities by
simultaneously (1) considering linkages between spatio-environmental variables
and biodiversity, (2) accounting for the negative relationship between past
sampling intensity and future biodiscovery potential, (3) incorporating a
priori knowledge about global species distribution patterns, (4) addressing
spatial autocorrelation in community composition, and (5) weighting theoretical
undocumented species by their predicted levels of conservation need. Using
boosted regression trees and 50km^2 map pixels spread across the global range
of amphibians, we identified several regions likely to contain many
undocumented amphibian species and conservation needs, including the Southeast
Asian Archipelago, humid portions of sub-Saharan Africa, and undersampled
portions of the Amazon, Andes Mountains, and Central America. We also ranked
top-scoring ecoregions by their mean and maximum biodiscovery potential and
found that the top-20 ranked ecoregions were most concentrated in the Southeast
Asian Archipelago and tropical Africa for undocumented species richness, and in
tropical Africa and tropical South America for integrative undocumented
amphibian conservation needs. However, high-scoring pixels tended to be widely
distributed across different ecoregions for both biodiscovery scoring
approaches.


  This paper examines a discrete predator-prey model that incorporates prey
refuge and its detrimental impact on the growth of the prey population. Age
structure is taken into account for predator species. Furthermore, juvenile
hunting as well as prey counter-attack are also considered. This paper provides
a comprehensive analysis of the existence and stability conditions pertaining
to all possible fixed points. The analytical and numerical investigation into
the occurrence of different bifurcations, such as the Neimark-Sacker
bifurcation and period-doubling bifurcation, in relation to various parameters
is discussed. The impact of the parameters reflecting prey growth and prey
refuge is thoroughly addressed. Numerous numerical simulations are presented in
order to validate the theoretical findings.


  Predicting the bioactivity of a ligand is one of the hardest and most
important challenges in computer-aided drug discovery. Despite years of data
collection and curation efforts by research organizations worldwide,
bioactivity data remains sparse and heterogeneous, thus hampering efforts to
build predictive models that are accurate, transferable and robust. The
intrinsic variability of the experimental data is further compounded by data
aggregation practices that neglect heterogeneity to overcome sparsity. Here we
discuss the limitations of these practices and present a hierarchical
meta-learning framework that exploits the information synergy across disparate
assays by successfully accounting for assay heterogeneity. We show that the
model achieves a drastic improvement in affinity prediction across diverse
protein targets and assay types compared to conventional baselines. It can
quickly adapt to new target contexts using very few observations, thus enabling
large-scale virtual screening in early-phase drug discovery.


  Protein language models learn powerful representations directly from
sequences of amino acids. However, they are constrained to generate proteins
with only the set of amino acids represented in their vocabulary. In contrast,
chemical language models learn atom-level representations of smaller molecules
that include every atom, bond, and ring. In this work, we show that chemical
language models can learn atom-level representations of proteins enabling
protein generation unconstrained to the standard genetic code and far beyond
it. In doing so, we show that language models can generate entire proteins atom
by atom -- effectively learning the multiple hierarchical layers of molecular
information that define proteins from their primary sequence to their
secondary, and tertiary structure. We demonstrate language models are able to
explore beyond protein space -- generating proteins with modified sidechains
that form unnatural amino acids. Even further, we find that language models can
explore chemical space and protein space simultaneously and generate novel
examples of protein-drug conjugates. The results demonstrate the potential for
biomolecular design at the atom level using language models.


  In biomedical research, validation of a new scientific discovery is tied to
the reproducibility of its experimental results. However, in genomics, the
definition and implementation of reproducibility still remain imprecise. Here,
we argue that genomic reproducibility, defined as the ability of bioinformatics
tools to maintain consistent genomics results across technical replicates, is
key to generating scientific knowledge and enabling medical applications. We
first discuss different concepts of reproducibility and then focus on
reproducibility in the context of genomics, aiming to establish clear
definitions of relevant terms. We then focus on the role of bioinformatics
tools and their impact on genomic reproducibility and assess methods of
evaluating bioinformatics tools in terms of genomic reproducibility. Lastly, we
suggest best practices for enhancing genomic reproducibility, with an emphasis
on assessing the performance of bioinformatics tools through rigorous testing
across multiple technical replicates.


  Precision medicine fundamentally aims to establish causality between
dysregulated biochemical mechanisms and cancer subtypes. Omics-based cancer
subtyping has emerged as a revolutionary approach, as different level of omics
records the biochemical products of multistep processes in cancers. This paper
focuses on fully exploiting the potential of multi-omics data to improve cancer
subtyping outcomes, and hence developed MoCLIM, a representation learning
framework. MoCLIM independently extracts the informative features from distinct
omics modalities. Using a unified representation informed by contrastive
learning of different omics modalities, we can well-cluster the subtypes, given
cancer, into a lower latent space. This contrast can be interpreted as a
projection of inter-omics inference observed in biological networks.
Experimental results on six cancer datasets demonstrate that our approach
significantly improves data fit and subtyping performance in fewer
high-dimensional cancer instances. Moreover, our framework incorporates various
medical evaluations as the final component, providing high interpretability in
medical analysis.


  We develop a modeling framework for cancer progression that distinguishes the
order of two possible mutations. Recent observations and information on
myeloproliferative neoplasms are analyzed within our framework. In some
patients with myeloproliferative neoplasms, two genetic mutations can be found,
JAK2 V617F and TET2. Whether or not one mutation is present will influence how
the other subsequent mutation affects the regulation of gene expression. When
both mutations are present, the order of their occurrence has been shown to
influence disease progression and prognosis. In this paper, we propose a
nonlinear ordinary differential equation (ODE) and Markov chain models to
explain the non-additive and non-commutative clinical observations with respect
to different orders of mutations: gene expression patterns, proportions of
cells with different mutations, and ages at diagnosis. We also propose
potential experiments measurements that can be used to verify our models.


  In a companion article it was shown in a certain precise sense that, for any
thermodynamical theory that respects the Kelvin-Planck Second Law, the
Hahn-Banach Theorem immediately ensures the existence of a pair of continuous
functions of the local material state -- a specific entropy (entropy per mass)
and a thermodynamic temperature -- that together satisfy the Clausius-Duhem
inequality for every process. There was no requirement that the local states
considered be states of equilibrium. This article addresses questions about
properties of the entropy and thermodynamic temperature functions so obtained:
To what extent do such temperature functions provide a faithful reflection of
``hotness"? In precisely which Kelvin-Planck theories is such a temperature
function essentially unique, and, among those theories, for which is the
entropy function also essentially unique? What is a thermometer for a
Kelvin-Planck theory, and, for the theory, what properties does the existence
of a thermometer confer? In all of these questions, the Hahn-Banach Theorem
again plays a crucial role.


  Antimicrobial peptides (AMPs) emerge as promising agents against
antimicrobial resistance, providing an alternative to conventional antibiotics.
Artificial intelligence (AI) revolutionized AMP discovery through both
discrimination and generation approaches. The discriminators aid the
identification of promising candidates by predicting key peptide properties
such as activity and toxicity, while the generators learn the distribution over
peptides and enable sampling novel AMP candidates, either de novo, or as
analogues of a prototype peptide. Moreover, the controlled generation of AMPs
with desired properties is achieved by discriminator-guided filtering,
positive-only learning, latent space sampling, as well as conditional and
optimized generation. Here we review recent achievements in AI-driven AMP
discovery, highlighting the most exciting directions.


  In evolutionary biology, networks are becoming increasingly used to represent
evolutionary histories for species that have undergone non-treelike or
reticulate evolution. Such networks are essentially directed acyclic graphs
with a leaf set that corresponds to a collection of species, and in which
non-leaf vertices with indegree 1 correspond to speciation events and vertices
with indegree greater than 1 correspond to reticulate events such as gene
transfer. Recently forest-based networks have been introduced, which are
essentially (multi-rooted) networks that can be formed by adding some arcs to a
collection of phylogenetic trees (or phylogenetic forest), where each arc is
added in such a way that its ends always lie in two different trees in the
forest. In this paper, we consider the complexity of deciding whether or not a
given network is proper forest-based, that is, whether it can be formed by
adding arcs to some underlying phylogenetic forest which contains the same
number of trees as there are roots in the network. More specifically, we show
that it can be decided in polynomial time whether or not a binary, tree-child
network with $m \ge 2$ roots is proper forest-based in case $m=2$, but that
this problem is NP-complete for $m\ge 3$. We also give a fixed parameter
tractable (FPT) algorithm for deciding whether or not a network in which every
vertex has indegree at most 2 is proper forest-based. A key element in proving
our results is a new characterization for when a network with $m$ roots is
proper forest-based which is given in terms of the existence of certain
$m$-colorings of the vertices of the network.