  The drug overdose crisis in the United States continues to intensify.
Fatalities have increased five-fold since 1999 reaching a record high of
108,000 deaths in 2021. The epidemic has unfolded through distinct waves of
different drug types, uniquely impacting various age, gender, race and ethnic
groups in specific geographical areas. One major challenge in designing
effective interventions is the forecasting of age-specific overdose patterns at
the local level so that prevention and preparedness can be effectively
delivered. We develop a forecasting method that assimilates observational data
obtained from the CDC WONDER database with an age-structured model of addiction
and overdose mortality. We apply our method nationwide and to three select
areas: Los Angeles County, Cook County and the five boroughs of New York City,
providing forecasts of drug-overdose mortality and estimates of relevant
epidemiological quantities, such as mortality and age-specific addiction rates.


  Behavioral experiments on the trust game have shown that trust and
trustworthiness are universal among human beings, contradicting the prediction
by assuming \emph{Homo economicus} in orthodox Economics. This means some
mechanism must be at work that favors their emergence. Most previous
explanations however need to resort to some factors based upon imitative
learning, a simple version of social learning. Here, we turn to the paradigm of
reinforcement learning, where individuals update their strategies by evaluating
the long-term return through accumulated experience. Specifically, we
investigate the trust game with the Q-learning algorithm, where each
participant is associated with two evolving Q-tables that guide one's decision
making as trustor and trustee respectively. In the pairwise scenario, we reveal
that high levels of trust and trustworthiness emerge when individuals
appreciate both their historical experience and returns in the future.
Mechanistically, the evolution of the Q-tables shows a crossover that resembles
human's psychological changes. We also provide the phase diagram for the game
parameters, where the boundary analysis is conducted. These findings are robust
when the scenario is extended to a latticed population. Our results thus
provide a natural explanation for the emergence of trust and trustworthiness
without external factors involved. More importantly, the proposed paradigm
shows the potential in deciphering many puzzles in human behaviors.


  The time varying reproduction number R is a critical variable for situational
awareness during infectious disease outbreaks, but delays between infection and
reporting hinder its accurate estimation in real time. We propose a nowcasting
method for improving the timeliness and accuracy of R estimates, based on
comparisons of successive versions of surveillance databases. The method was
validated against COVID-19 surveillance data collected in Italy over an
18-month period. Compared to traditional methods, the nowcasted reproduction
number reduced the estimation delay from 13 to 8 days, while maintaining a
better accuracy. Moreover, it allowed anticipating the detection of periods of
epidemic growth by between 6 and 23 days. The method offers a simple and
generally applicable tool to improve situational awareness during an epidemic
outbreak, allowing for informed public health response planning.


  The single-particle cryo-EM field faces the persistent challenge of preferred
orientation, lacking general computational solutions. We introduce cryoPROS, an
AI-based approach designed to address the above issue. By generating the
auxiliary particles with a conditional deep generative model, cryoPROS
addresses the intrinsic bias in orientation estimation for the observed
particles. We effectively employed cryoPROS in the cryo-EM single particle
analysis of the hemagglutinin trimer, showing the ability to restore the
near-atomic resolution structure on non-tilt data. Moreover, the enhanced
version named cryoPROS-MP significantly improves the resolution of the membrane
protein NaX using the no-tilted data that contains the effects of micelles.
Compared to the classical approaches, cryoPROS does not need special
experimental or image acquisition techniques, providing a purely computational
yet effective solution for the preferred orientation problem. Finally, we
conduct extensive experiments that establish the low risk of model bias and the
high robustness of cryoPROS.


  Abrupt changes in behavior can often be associated with changes in underlying
behavioral states. When placed off food, the foraging behavior of C. elegans
can be described as a change between an initial local-search behavior
characterized by a high rate of reorientations, followed by a global-search
behavior characterized by sparse reorientations. This is commonly observed in
individual worms, but when numerous worms are characterized, only about half
appear to exhibit this behavior. We propose an alternative model that predicts
both abrupt and continuous changes to reorientation that does not rely on
behavioral states. This model is inspired by molecular dynamics modeling that
defines the foraging reorientation rate as a decaying parameter. By
stochastically sampling from the probability distribution defined by this rate,
both abrupt and gradual changes to reorientation rates can occur, matching
experimentally observed results. Crucially, this model does not depend on
behavioral states or information accumulation. Even though abrupt behavioral
changes do occur, they may not necessarily be indicative of abrupt changes in
behavioral states, especially when abrupt changes are not universally observed
in the population.


  Suppose that females choose males based on attributes that do not signal any
genetic quality that is not related to the choice itself. Can being choosy
confer selective advantage in this situation? We introduce correlated
strategies, which means that females, when making their choice, may take into
consideration external and independent random factors that are known to be
observable by all. Individual-based simulation is used to show that, in this
case, choosiness can emerge against the cost of over 25\% when pitted against
randomly mating females. Moreover, after being established in the population,
it can sustain costs of over 35\% .
  While such costs are not biologically plausible, they demonstrate
unequivocally that sexual choice is a strong evolutionary force. Thus,
correlated strategies are shown to be an evolutionary tool that channels
randomness from the environment into genetic diversity. In addition, it turns
out that a higher number of attributes in the ornament makes the choice more
advantageous, which may result in a runaway complexity of sexual traits.
Implications for the evolution of (female) cognitive abilities and speciation
are discussed.


  Objective: The study explored therapeutic targets and mechanisms of Gegen
Qinlian Decoction for Helicobacter pylori infection and related gastric cancer
using network pharmacology, molecular docking, and Mendelian randomization.
  Methods: Medicinal components of Gegen Qinlian Decoction were extracted from
TCMSP and HERB databases. Disease treatment targets were sourced from DisGeNET
and PubChem. Interaction networks were constructed via the STRING database and
visualized using Cytoscape 3.9.1. Enrichment analysis of intersected targets
was performed using DAVID and Metascapes. Molecular docking employed Autodock
Tools 1.5.6 and PyMOL 2.5.2. Mendelian randomization was based on the ukb-b-531
sample from UK Biobank.
  Results: 146 active components and 248 targets from Gegen Qinlian Decoction
were identified. 66 targets overlapped with Helicobacter pylori infection
genes. Molecular docking highlighted interactions between primary drug
components like quercetin, wogonin, kaempferol, and target genes PTGS1, PTGS2,
MAPK14. Mendelian randomization pinpointed genes like IGF2, PIK3CG, GJA1, and
PLAU associated with Helicobacter pylori infection.
  Conclusion: Gegen Qinlian Decoction's active components target Helicobacter
pylori infection through diverse targets and pathways, presenting potential
research avenues.


  This paper addresses the problem of statistical inference for latent
continuous-time stochastic processes, which is often intractable, particularly
for discrete state space processes described by Markov jump processes. To
overcome this issue, we propose a new tractable inference scheme based on an
entropic matching framework that can be embedded into the well-known
expectation propagation algorithm. We demonstrate the effectiveness of our
method by providing closed-form results for a simple family of approximate
distributions and apply it to the general class of chemical reaction networks,
which are a crucial tool for modeling in systems biology. Moreover, we derive
closed form expressions for point estimation of the underlying parameters using
an approximate expectation maximization procedure. We evaluate the performance
of our method on various chemical reaction network instantiations, including a
stochastic Lotka-Voltera example, and discuss its limitations and potential for
future improvements. Our proposed approach provides a promising direction for
addressing complex continuous-time Bayesian inference problems.


  The COVID-19 pandemic has left a lasting impact on individuals, with many
experiencing persistent symptoms, including inflammation, in the post-acute
phase of the disease. Detecting and monitoring these inflammatory biomarkers is
critical for timely intervention and improved patient outcomes. This study
employs machine learning techniques to automate the identification of
persistent inflammatory biomarkers in 290 post-COVID-19 patients, based on
medical data collected from hospitals in Iraq. The data encompassed a wide
array of clinical parameters, such as C-reactive protein and interleukin-6
levels, patient demographics, comorbidities, and treatment histories. Rigorous
data preprocessing and feature selection processes were implemented to optimize
the dataset for machine learning analysis. Various machine learning algorithms,
including logistic regression, random forests, support vector machines, and
gradient boosting, were deployed to construct predictive models. These models
exhibited promising results, showcasing high accuracy and precision in the
identification of patients with persistent inflammation. The findings of this
study underscore the potential of machine learning in automating the detection
of persistent inflammatory biomarkers in post-COVID-19 patients. These models
can serve as valuable tools for healthcare providers, facilitating early
diagnosis and personalized treatment strategies for individuals at risk of
persistent inflammation, ultimately contributing to improved post-acute
COVID-19 care and patient well-being. Keywords: COVID-19, post-COVID-19,
inflammation, biomarkers, machine learning, early detection.


  We have devised a data-driven framework for uncovering hidden control
strategies used by an evolutionary system described by an evolutionary
probability distribution. This innovative framework enables deciphering of the
concealed mechanisms that contribute to the progression or mitigation of such
situations as the spread of COVID-19. Novel algorithms are used to estimate the
optimal control in tandem with the parameters for evolution in general
dynamical systems, thereby extending the concept of model predictive control.
This is a significant departure from conventional control methods, which
require knowledge of the system to manipulate its evolution and of the
controller's strategy or parameters. We used a generalized additive model,
supplemented by extensive statistical testing, to identify a set of predictor
covariates closely linked to the control. Using real-world COVID-19 data, we
successfully delineated the descriptive behaviors of the COVID-19 epidemics in
five prefectures in Japan and nine countries. We compared these nine countries
and grouped them on the basis of shared profiles, providing valuable insights
into their pandemic responses. Our findings underscore the potential of our
framework as a powerful tool for understanding and managing complex
evolutionary processes.


  In this paper, we show that an extension of the classical Tikhonov--Fenichel
asymptotic procedure applied to multiscale models of vector-borne diseases,
with time scales determined by the dynamics of human and vector populations,
yields a simplified model approximating the original one in a consistent, and
uniform for large times, way. Furthermore, we construct a higher-order
approximation based on the classical Chapman-Enskog procedure of kinetic theory
and show, in particular, that it is equivalent to the dynamics on the
first-order approximation of the slow manifold in the Fenichel theory.


  The global SARS-CoV-2 (COVID-19) pandemic highlighted the challenge of
equitable vaccine distribution between high- and low-income countries. Many
high-income countries were reluctant or slow to distribute extra doses of the
vaccine to lower-income countries via the COVID-19 Vaccines Global Access
(COVAX) collaboration. In addition to moral objections to such vaccine
nationalism, vaccine inequity during a pandemic could contribute to the
evolution of new variants of the virus and possibly increase total deaths,
including in the high-income countries. Using the COVID-19 pandemic as a case
study, we use the epidemiological model of Holleran et al. that incorporates
virus mutation. We identify realistic scenarios under which a donor country
prefers to donate vaccines before distributing them locally in order to
minimize local deaths during a pandemic. We demonstrate that a nondonor-first
vaccination policy can delay, sometimes dramatically, the emergence of
more-contagious variants. Even more surprising, donating all vaccines is
sometimes better for the donor country than a sharing policy in which half of
the vaccines are donated and half are retained because of the impact donation
can have on delaying the emergence of a more contagious virus. Nondonor-first
vaccine allocation is optimal in scenarios in which the local health impact of
the vaccine is limited or when delaying emergence of a variant is especially
valuable. In all cases, we find that vaccine distribution is not a zero-sum
game between donor and nondonor countries. Thus, in addition to moral reasons
to avoid vaccine nationalism, donor nations can also realize local health
benefits from donating vaccines. The insights yielded by this framework can be
used to guide equitable vaccine distribution in future pandemics.


  The majority of viruses are organised according to the structural blueprints
of the seminal Caspar-Klug theory. However, there are a number of notable
exceptions to this geometric design principle. Prominent examples are the
cancer-causing papilloma viridae and the \textit{de novo} designed AaLS cages
that exhibit non-quasiequivalent capsid structures with protein numbers
excluded by Caspar-Klug theory. The biophysical properties of these
geometrically distinct architectures and the fitness advantages driving their
evolution are currently unclear. We investigate here the resilience to
fragmentation and disassembly behaviour of these capsid geometries by
introducing a percolation theory on weighted graphs. We show that these cage
architectures follow one of two distinct disassembly pathways, preferring
either hole formation or capsid fragmentation. This suggests that preference
for specific disassembly scenarios could be a driving force for the evolution
of the non Caspar-Klug protein container architectures.


  In this study, we leveraged machine learning techniques to identify risk
factors associated with post-COVID-19 mental health disorders. Our analysis,
based on data collected from 669 patients across various provinces in Iraq,
yielded valuable insights. We found that age, gender, and geographical region
of residence were significant demographic factors influencing the likelihood of
developing mental health disorders in post-COVID-19 patients. Additionally,
comorbidities and the severity of COVID-19 illness were important clinical
predictors. Psychosocial factors, such as social support, coping strategies,
and perceived stress levels, also played a substantial role. Our findings
emphasize the complex interplay of multiple factors in the development of
mental health disorders following COVID-19 recovery. Healthcare providers and
policymakers should consider these risk factors when designing targeted
interventions and support systems for individuals at risk. Machine
learning-based approaches can provide a valuable tool for predicting and
preventing adverse mental health outcomes in post-COVID-19 patients. Further
research and prospective studies are needed to validate these findings and
enhance our understanding of the long-term psychological impact of the COVID-19
pandemic. This study contributes to the growing body of knowledge regarding the
mental health consequences of the COVID-19 pandemic and underscores the
importance of a multidisciplinary approach to address the diverse needs of
individuals on the path to recovery. Keywords: COVID-19, mental health, risk
factors, machine learning, Iraq


  The COVID-19 pandemic has globally posed numerous health challenges, notably
the emergence of post-COVID-19 cardiovascular complications. This study
addresses this by utilizing data-driven machine learning models to predict such
complications in 352 post-COVID-19 patients from Iraq. Clinical data, including
demographics, comorbidities, lab results, and imaging, were collected and used
to construct predictive models. These models, leveraging various machine
learning algorithms, demonstrated commendable performance in identifying
patients at risk. Early detection through these models promises timely
interventions and improved outcomes. In conclusion, this research underscores
the potential of data-driven machine learning for predicting post-COVID-19
cardiovascular complications, emphasizing the need for continued validation and
research in diverse clinical settings.


  The success of language models, especially transformer-based architectures,
has trickled into other domains giving rise to "scientific language models"
that operate on small molecules, proteins or polymers. In chemistry, language
models contribute to accelerating the molecule discovery cycle as evidenced by
promising recent findings in early-stage drug discovery. Here, we review the
role of language models in molecular discovery, underlining their strength in
de novo drug design, property prediction and reaction chemistry. We highlight
valuable open-source software assets thus lowering the entry barrier to the
field of scientific language modeling. Last, we sketch a vision for future
molecular design that combines a chatbot interface with access to computational
chemistry tools. Our contribution serves as a valuable resource for
researchers, chemists, and AI enthusiasts interested in understanding how
language models can and will be used to accelerate chemical discovery.


  Does an ecological community allow stable coexistence? Identifying the
general principles that determine the answer to this question is a central
problem of theoretical ecology. Random matrix theory approaches have uncovered
the general trends of the effect of competitive, mutualistic, and predator-prey
interactions between species on stability of coexistence. However, an
ecological community is determined not only by the counts of these different
interaction types, but also by their network arrangement. This cannot be
accounted for in a direct statistical description that would enable random
matrix theory approaches. Here, we therefore develop a different approach, of
exhaustive analysis of small ecological communities, to show that this
arrangement of interactions can influence stability of coexistence more than
these general trends. We analyse all interaction networks of $N\leqslant 5$
species with Lotka-Volterra dynamics by combining exact results for $N\leqslant
3$ species and numerical exploration. Surprisingly, we find that a very small
subset of these networks are "impossible ecologies", in which stable
coexistence is non-trivially impossible. We prove that the possibility of
stable coexistence in general ecologies is determined by similarly rare
"irreducible ecologies". By random sampling of interaction strengths, we then
show that the probability of stable coexistence varies over many orders of
magnitude even in ecologies that differ only in the network arrangement of
identical ecological interactions. Finally, we demonstrate that our approach
can reveal the effect of evolutionary or environmental perturbations of the
interaction network. Overall, this work reveals the importance of the full
structure of the network of interactions for stability of coexistence in
ecological communities.


  Airborne infection risk analysis is usually performed for enclosed spaces
where susceptible individuals are exposed to infectious airborne respiratory
droplets by inhalation. It is usually based on exponential, dose-response
models of which a widely used variant is the Wells-Riley (WR) model. We revisit
this infection-risk estimate and extend it to the population level. We use an
epidemiological model where the mode of pathogen transmission, either airborne
or contact, is explicitly considered. We illustrate the link between
epidemiological models and the WR model. We argue that airborne infection
quanta are, up to an overall density, airborne infectious respiratory droplets
modified by a parameter that depends on biological properties of the pathogen,
physical properties of the droplet, and behavioural parameters of the
individual. We calculate the time-dependent risk to be infected during the
epidemic for two scenarios. We show how the epidemic infection risk depends on
the viral latent period and the event time, the time infection occurs. The
infection risk follows the dynamics of the infected population. As the latency
period decreases, infection risk increases. The longer a susceptible is present
in the epidemic, the higher is its risk of infection by equal exposure time to
the mode of transmission.


  Recent advancements in Cryo-EM and protein structure prediction algorithms
have made large-scale protein structures accessible, paving the way for machine
learning-based functional annotations.The field of geometric deep learning
focuses on creating methods working on geometric data. An essential aspect of
learning from protein structures is representing these structures as a
geometric object (be it a grid, graph, or surface) and applying a learning
method tailored to this representation. The performance of a given approach
will then depend on both the representation and its corresponding learning
method.
  In this paper, we investigate representing proteins as $\textit{3D mesh
surfaces}$ and incorporate them into an established representation benchmark.
Our first finding is that despite promising preliminary results, the surface
representation alone does not seem competitive with 3D grids. Building on this,
we introduce a synergistic approach, combining surface representations with
graph-based methods, resulting in a general framework that incorporates both
representations in learning. We show that using this combination, we are able
to obtain state-of-the-art results across $\textit{all tested tasks}$. Our code
and data can be found online: https://github.com/Vincentx15/atom2D .


  Ligand molecule conformation generation is a critical challenge in drug
discovery. Deep learning models have been developed to tackle this problem,
particularly through the use of generative models in recent years. However,
these models often generate conformations that lack meaningful structure and
randomness due to the absence of essential side information. Examples of such
side information include the chemical and geometric features of the target
protein, ligand-target compound interactions, and ligand chemical properties.
Without these constraints, the generated conformations may not be suitable for
further selection and design of new drugs. To address this limitation, we
propose a novel method for generating ligand conformations that leverage side
information and incorporate flexible constraints into standard diffusion
models. Drawing inspiration from the concept of message passing, we introduce
ligand-target massage passing block, a mechanism that facilitates the exchange
of information between target nodes and ligand nodes, thereby incorporating
target node features. To capture non-covalent interactions, we introduce
ligand-target compound inter and intra edges. To further improve the biological
relevance of the generated conformations, we train energy models using scalar
chemical features. These models guide the progress of the standard Denoising
Diffusion Probabilistic Models, resulting in more biologically meaningful
conformations. We evaluate the performance of SIDEGEN using the PDBBind-2020
dataset, comparing it against other methods. The results demonstrate
improvements in both Aligned RMSD and Ligand RMSD evaluations. Specifically,
our model outperforms GeoDiff (trained on PDBBind-2020) by 20% in terms of the
median aligned RMSD metric.


  Without knowledge of specific pockets, generating ligands based on the global
structure of a protein target plays a crucial role in drug discovery as it
helps reduce the search space for potential drug-like candidates in the
pipeline. However, contemporary methods require optimizing tailored networks
for each protein, which is arduous and costly. To address this issue, we
introduce TargetVAE, a target-aware variational auto-encoder that generates
ligands with high binding affinities to arbitrary protein targets, guided by a
novel multimodal deep neural network built based on graph Transformers as the
prior for the generative model. This is the first effort to unify different
representations of proteins (e.g., sequence of amino-acids, 3D structure) into
a single model that we name as Protein Multimodal Network (PMN). Our multimodal
architecture learns from the entire protein structures and is able to capture
their sequential, topological and geometrical information. We showcase the
superiority of our approach by conducting extensive experiments and
evaluations, including the assessment of generative model quality, ligand
generation for unseen targets, docking score computation, and binding affinity
prediction. Empirical results demonstrate the promising performance of our
proposed approach. Our software package is publicly available at
https://github.com/HySonLab/Ligand_Generation


  This study explores the application and performance of Transformational
Machine Learning (TML) in drug discovery. TML, a meta learning algorithm,
excels in exploiting common attributes across various domains, thus developing
composite models that outperform conventional models. The drug discovery
process, which is complex and time-consuming, can benefit greatly from the
enhanced prediction accuracy, improved interpretability and greater
generalizability provided by TML. We explore the efficacy of different machine
learning classifiers, where no individual classifier exhibits distinct
superiority, leading to the consideration of ensemble classifiers such as the
Random Forest.
  Our findings show that TML outperforms base Machine Learning (ML) as the
number of training datasets increases, due to its capacity to better
approximate the correct hypothesis, overcome local optima, and expand the space
of representable functions by combining separate classifiers capabilities.
However, this superiority is relative to the resampling methods applied, with
Near Miss demonstrating poorer performance due to noisy data, overlapping
classes, and nonlinear class boundaries. Conversely, Random Over Sampling (ROS)
provides a more robust performance given its resistance to noise and outliers,
improved class overlap management, and suitability for nonlinear class
boundaries.


  In light of the continuous transmission and evolution of SARS-CoV-2 coupled
with a significant decline in clinical testing, there is a pressing need for
scalable, cost-effective, long-term, passive surveillance tools to effectively
monitor viral variants circulating in the population. Wastewater genomic
surveillance of SARS-CoV-2 has arrived as an alternative to clinical genomic
surveillance, allowing to continuously monitor the prevalence of viral lineages
in communities of various size at a fraction of the time, cost, and logistic
effort and serving as an early warning system for emerging variants, critical
for developed communities and especially for underserved ones. Importantly,
lineage prevalence estimates obtained with this approach aren't distorted by
biases related to clinical testing accessibility and participation. However,
the relative performance of bioinformatics methods used to measure relative
lineage abundances from wastewater sequencing data is unknown, preventing both
the research community and public health authorities from making informed
decisions regarding computational tool selection. Here, we perform
comprehensive benchmarking of 18 bioinformatics methods for estimating the
relative abundance of SARS-CoV-2 (sub)lineages in wastewater by using data from
36 in vitro mixtures of synthetic lineage and sublineage genomes. In addition,
we use simulated data from 78 mixtures of lineages and sublineages co-occurring
in the clinical setting with proportions mirroring their prevalence ratios
observed in real data. Importantly, we investigate how the accuracy of the
evaluated methods is impacted by the sequencing technology used, the associated
error rate, the read length, read depth, but also by the exposure of the
synthetic RNA mixtures to wastewater, with the goal of capturing the effects
induced by the wastewater matrix, including RNA fragmentation and degradation.


  Motivation: Read mapping is a computationally expensive process and a major
bottleneck in genomics analyses. The performance of read mapping is mainly
limited by the performance of three key computational steps: Index Querying,
Seed Chaining, and Sequence Alignment. The first step is dominated by how fast
and frequent it accesses the main memory (i.e., memory-bound), while the latter
two steps are dominated by how fast the CPU can compute their
computationally-costly dynamic programming algorithms (i.e., compute-bound).
Accelerating these three steps by exploiting new algorithms and new hardware
devices is essential to accelerate most genome analysis pipelines that widely
use read mapping. Given the large body of work on accelerating Sequence
Alignment, this work focuses on significantly improving the remaining steps.
Results: We introduce GateSeeder, the first CPU-FPGA-based near-memory
acceleration of both short and long read mapping. GateSeeder exploits
near-memory computation capability provided by modern FPGAs that couple a
reconfigurable compute fabric with high-bandwidth memory (HBM) to overcome the
memory-bound and compute-bound bottlenecks. GateSeeder also introduces a new
lightweight algorithm for finding the potential matching segment pairs. Using
real ONT, HiFi, and Illumina sequences, we experimentally demonstrate that
GateSeeder outperforms Minimap2, without performing sequence alignment, by up
to 40.3x, 4.8x, and 2.3x, respectively. When performing read mapping with
sequence alignment, GateSeeder outperforms Minimap2 by 1.15-4.33x (using KSW2)
and by 1.97-13.63x (using WFA-GPU). Availability:
https://github.com/CMU-SAFARI/GateSeeder


  The chronological age used in demography describes the linear evolution of
the life of a living being. The chronological age cannot give precise
information about the exact developmental stage or aging processes an organism
has reached. On the contrary, the biological age (or epigenetic age) represents
the true evolution of the tissues and organs of the living being. Biological
age is not always linear and sometimes proceeds by discontinuous jumps. These
jumps can be positive (we then speak of rejuvenation) or negative (in the event
of premature aging), and they can be dependent on endogenous events such as
pregnancy (negative jump) or stroke (positive jump) or exogenous ones such as
surgical treatment (negative jump) or infectious disease (positive jump). The
article proposes a mathematical model of the biological age by defining a valid
model for the two types of jumps (positive and negative). The existence and
uniqueness of the solution are solved, and its temporal dynamic is analyzed
using a moments equation. We also provide some individual-based stochastic
simulations.


  Predicting drug properties is key in drug discovery to enable de-risking of
assets before expensive clinical trials, and to find highly active compounds
faster. Interest from the Machine Learning community has led to the release of
a variety of benchmark datasets and proposed methods. However, it remains
unclear for practitioners which method or approach is most suitable, as
different papers benchmark on different datasets and methods, leading to
varying conclusions that are not easily compared. Our large-scale empirical
study links together numerous earlier works on different datasets and methods;
thus offering a comprehensive overview of the existing property classes,
datasets, and their interactions with different methods. We emphasise the
importance of uncertainty quantification and the time and therefore cost of
applying these methods in the drug development decision-making cycle. We
discover that the best method depends on the dataset, and that engineered
features with classical ML methods often outperform deep learning.
Specifically, QSAR datasets are typically best analysed with classical methods
such as Gaussian Processes while ADMET datasets are sometimes better described
by Trees or Deep Learning methods such as Graph Neural Networks or language
models. Our work highlights that practitioners do not yet have a
straightforward, black-box procedure to rely on, and sets the precedent for
creating practitioner-relevant benchmarks. Deep learning approaches must be
proven on these benchmarks to become the practical method of choice in drug
property prediction.


  Hemoglycin, a space polymer of glycine and iron, has been identified in the
carbonaceous chondritic meteorites Allende, Acfer 086, Kaba, Sutters Mill and
Orgueil. Its core form has a mass of 1494Da and is basically an antiparallel
pair of polyglycine strands linked at each end by an iron atom. The polymer
forms two and three dimensional lattices with an inter vertex distance of
4.9nm. Here the extraction technique for meteorites is applied to a 2.1Gya
fossil stromatolite to reveal the presence of hemoglycin by mass spectrometry.
Intact ooids from a recent 3,000Year stromatolite exhibited the same visible
hemoglycin fluorescence in response to x-rays as an intact crystal from the
Orgueil meteorite. X-ray analysis of these ooids at wavelengths above and below
the iron K absorption edge yielded a set of high order diffraction rings that
confirmed the existence and nature of a three dimensional lattice of 4.9nm
inter-vertex spacing. The lattice is filled by micro crystals of the aragonite
and calcite forms of calcium carbonate. It seems probable that the copious in
fall of carbonaceous meteoritic material, from Archaean times onward, has left
traces of hemoglycin in sedimentary carbonates and potentially has influenced
ooid formation.


  Molecular property prediction, crucial for early drug candidate screening and
optimization, has seen advancements with deep learning-based methods. While
deep learning-based methods have advanced considerably, they often fall short
in fully leveraging 3D spatial information. Specifically, current molecular
encoding techniques tend to inadequately extract spatial information, leading
to ambiguous representations where a single one might represent multiple
distinct molecules. Moreover, existing molecular modeling methods focus
predominantly on the most stable 3D conformations, neglecting other viable
conformations present in reality. To address these issues, we propose 3D-Mol, a
novel approach designed for more accurate spatial structure representation. It
deconstructs molecules into three hierarchical graphs to better extract
geometric information. Additionally, 3D-Mol leverages contrastive learning for
pretraining on 20 million unlabeled data, treating their conformations with
identical topological structures as weighted positive pairs and contrasting
ones as negatives, based on the similarity of their 3D conformation descriptors
and fingerprints. We compare 3D-Mol with various state-of-the-art baselines on
7 benchmarks and demonstrate our outstanding performance.


  We propose that spontaneous folding and molecular evolution of biopolymers
are two universal aspects that must concur for life to happen. These aspects
are fundamentally related to the chemical composition of biopolymers and
crucially depend on the solvent in which they are embedded. We show that
molecular information theory and energy landscape theory allow us to explore
the limits that solvents impose on biopolymer existence. We consider 54
solvents, including water, alcohols, hydrocarbons, halogenated solvents,
aromatic solvents, and low molecular weight substances made up of elements
abundant in the universe, which may potentially take part in alternative
biochemistries. We find that along with water, there are many solvents for
which the liquid regime is compatible with biopolymer folding and evolution. We
present a ranking of the solvents in terms of biopolymer compatibility. Many of
these solvents have been found in molecular clouds or may be expected to occur
in extrasolar planets.


  While investigating methods to predict small molecule potencies, we found
random forests or support vector machines paired with extended-connectivity
fingerprints (ECFP) consistently outperformed recently developed methods. A
detailed investigation into regression algorithms and molecular fingerprints
revealed gradient-boosted decision trees, particularly CatBoost, in conjunction
with a combination of ECFP, Avalon, and ErG fingerprints, as well as 200
molecular properties, to be most effective. Incorporating a graph neural
network fingerprint further enhanced performance. We successfully validated our
model across 22 Therapeutics Data Commons ADMET benchmarks. Our findings
underscore the significance of richer molecular representations for accurate
property prediction.


  As the size of accessible compound libraries expands to over 10 billion, the
need for more efficient structure-based virtual screening methods is emerging.
Different pre-screening methods have been developed for rapid screening, but
there is still a lack of structure-based methods applicable to various proteins
that perform protein-ligand binding conformation prediction and scoring in an
extremely short time. Here, we describe for the first time a deep-learning
framework for structure-based pharmacophore modeling to address this challenge.
We frame pharmacophore modeling as an instance segmentation problem to
determine each protein hotspot and the location of corresponding
pharmacophores, and protein-ligand binding pose prediction as a graph-matching
problem. PharmacoNet is significantly faster than state-of-the-art
structure-based approaches, yet reasonably accurate with a simple scoring
function. Furthermore, we show the promising result that PharmacoNet
effectively retains hit candidates even under the high pre-screening filtration
rates. Overall, our study uncovers the hitherto untapped potential of a
pharmacophore modeling approach in deep learning-based drug discovery.


  Objective: To investigate free radical scavenging potential of
thiocarbohydrazones derivatives for the discovery of antioxidant compounds of
synthetic origin. Method: Eighteen thiocarbohydrazones derivatives were
screened for antioxidant activities by using in vitro radical inhibition assay.
Cytotoxicity all derivatives were evaluated using invitro cytotoxicity assay on
mouse fibroblast cell line. Results: various compounds of the have shown
significant radical scavenging activities as compared with the standard drug
butylated hydroxytoluene. In addition in an in vitro cell cytotoxicity assay
selected analogues were found to be non-cytotoxic against fibroblast 3T3 cell
line. Conclusion: Taken together, these findings all compound this series could
be further evaluated in in vivo model of oxidative stress to develop as a
therapeutic agent against oxidative stress related disorders due to exhibiting
excellent free radical scavenging properties and having their non cytotoxic
nature against normal cells.


  $\alpha$-glucosidase and prolylendopeptidase has altered expression and
activity patterns in neurological disease, type 2diabetes respectively and
several cancers. Here we screened a series 1-29 benzophenone semicarbazone
derivatives for in vitro free radical scavenging, alpha-glucosidase and
prolylendopeptidase inhibition activities. Seven derivatives were identified as
potential free radical scavengers, 14 as alpha-glucosidase, and 9 as
prolylendopeptidase inhibitors. Kinetic studies on the most promising
inhibitors were performed. Compounds 23, 27, 25 and 28 were found as inhibitor
of alpha-glucosidase, while compound 26 inhibited both prolylendopeptidase and
alpha-glucosidase. The binding modes and binding free energy of the multi
targeted inhibitor 26 were predicted by molecular docking studies. These
results provide insights on prolylendopeptidase and alpha-glucosidase
inhibition of compound 26 for further development as therapeutic agents for
neoplastic, neurological, and endocrine disorders.


  In a predator-prey interaction, there are many factors present which affect
the growth of the species, either positively or negatively. Fear of predation
is such a factor which causes psychological stress in a prey species so that
their overall growth starts to reduce. In this work, a predator-prey model is
proposed where the prey species is facing such reduction in their growth out of
fear, and also the predator is provided with an alternative food source which
helps the prey to hide themselves in a safe place meanwhile. The dynamics
produces a wide range of interesting results, including the significance of the
presence of a certain amount of fear or even prey refuge for population
coexistence. The analysis is extended later to the nonlocal model to analyze
how the non-equilibrium phenomena influence the dynamical behaviour. It is
observed through numerical simulations that though the scope of pattern
formation reduces with the increase of fear level in the spatio-temporal model,
but the incorporation of nonlocal interaction further increases the chance of
species colonization.


  Alkaline fens are known to be wetlands that host large quantities of carbon,
but also a huge biodiversity. However, anthropic pressures are degrading the
peatlands with drainage, and conversion to agriculture or urbanization. These
pressures induce externalities like a large release of greenhouse gas (GHG) by
peat mineralization, and a loss of biodiversity, since they host numerous
protected or endangered species. One of them is the Liparis loeselii, a small
orchid facing declines in Europe and for which conservation measures are taken.
Nevertheless, if recent studies inferred some factors shaping its population
dynamics, they are still not clearly understood, particularly in fen contexts.
This study aims at disentangle the processes shaping a continental population
of L. loeselii in the Somme valley, among factors related to the hydrology of
the site and the elevation of the individuals. We used Bayesian generalized
linear models to infer the seasonal water level effects on the presence of L.
loeselii. Moreover, we showed that the elevation of L. loeselii population is
around the mean water levels of winter, and that the occurrence can be promoted
by water level, notably in summer. Conversely, the highest water levels
reported seemed detrimental to the L. loeselii population, suggesting a
probably negative effect of flooding acting on local dispersal. Finally, this
study provides insights to take restoration measures of alkaline fens, like a
better hydrological functioning and an optimal water management.


  Mutational signatures connect characteristic mutational patterns in the
genome with biological processes that take place in the tumor tissues. Analysis
of mutational signatures can help elucidate tumor evolution, prognosis, and
therapeutic strategies. Although tools for extracting mutational signatures de
novo have been extensively benchmarked, a similar effort is lacking for tools
that fit known mutational signatures to a given catalog of mutations. We fill
this gap by comprehensively evaluating eleven signature fitting tools
(well-established as well as recent) on synthetic input data. To create
realistic input data, we use empirical signature weights in tumor tissue
samples from the COSMIC database. The study design allows us to assess the
effects of the number of mutations, type of cancer, and the catalog of
reference signatures on the results obtained with various fitting tools. We
find substantial performance differences between the evaluated tools. Averaged
over 120,000 simulated mutational catalogs corresponding to eight different
cancer types, SigProfilerSingleSample and SigProfilerAssignment perform best
for small and large numbers of mutations per sample, respectively. We further
show that ad hoc constraining the list of reference signatures is likely to
produce inferior results and that noisy estimates of signature weights in
samples with as few as 100 mutations can still be useful in downstream
analysis.


  The genetic basis of phenotypic differences between species is among the most
longstanding questions in evolutionary biology. How new genes form and the
processes selection acts to produce differences across species are fundamental
to understand how species persist and evolve in an ever-changing environment.
Adaptation and genetic innovation arise in the genome by a variety of sources.
Functional genomics requires both intrinsic genetic discoveries, as well as
empirical testing to observe adaptation between lineages. Here we explore two
species of Drosophila on the island of Sao Tome and mainland Africa, D.
santomea and D. yakuba. These two species both inhabit the island, but occupy
differing species distributions based on elevation, with D. yakuba also having
populations on mainland Africa. Intrinsic evidence shows genes between species
may have a role in adaptation to higher UV tolerance with DNA repair mechanisms
(PARP) and resistance to humeral stress lethal effects (Victoria). We conducted
empirical assays between island D. santomea, D. yakuba, and mainland D. yakuba.
Flies were shocked with UVB radiation (@ 302 nm) at 1650-1990 mW/cm2 for 30
minutes on a transilluminator apparatus. Custom 5-wall acrylic enclosures were
constructed for viewing and containment of flies. All assays were filmed.
Island groups did show significant differences between fall-time under UV
stress and recovery time post-UV stress test between regions and sex. This
study shows evidence that mainland flies are less resistant to UV radiation
than their island counterparts. Further work exploring the genetic basis for UV
tolerance will be conducted from empirical assays. Understanding the mechanisms
and processes that promote adaptation and testing extrinsic traits within the
context of the genome is crucially important to understand evolutionary
machinery.


  Discovering genes with similar functions across diverse biomedical contexts
poses a significant challenge in gene representation learning due to data
heterogeneity. In this study, we resolve this problem by introducing a novel
model called Multimodal Similarity Learning Graph Neural Network, which
combines Multimodal Machine Learning and Deep Graph Neural Networks to learn
gene representations from single-cell sequencing and spatial transcriptomic
data. Leveraging 82 training datasets from 10 tissues, three sequencing
techniques, and three species, we create informative graph structures for model
training and gene representations generation, while incorporating
regularization with weighted similarity learning and contrastive learning to
learn cross-data gene-gene relationships. This novel design ensures that we can
offer gene representations containing functional similarity across different
contexts in a joint space. Comprehensive benchmarking analysis shows our
model's capacity to effectively capture gene function similarity across
multiple modalities, outperforming state-of-the-art methods in gene
representation learning by up to 97.5%. Moreover, we employ bioinformatics
tools in conjunction with gene representations to uncover pathway enrichment,
regulation causal networks, and functions of disease-associated or
dosage-sensitive genes. Therefore, our model efficiently produces unified gene
representations for the analysis of gene functions, tissue functions, diseases,
and species evolution.


  Protein-Protein Interaction Networks aim to model the interactome, providing
a powerful tool for understanding the complex relationships governing cellular
processes. These networks have numerous applications, including functional
enrichment, discovering cancer driver genes, identifying drug targets, and
more. Various databases make protein-protein networks available for many
species, including Homo sapiens. This work topologically compares four Homo
sapiens networks using a coarse-to-fine approach, comparing global
characteristics, sub-network topology, specific nodes centrality, and
interaction significance. Results show that the four human protein networks
share many common protein-encoding genes and some global measures, but
significantly differ in the interactions and neighbourhood. Small sub-networks
from cancer pathways performed better than the whole networks, indicating an
improved topological consistency in functional pathways. The centrality
analysis shows that the same genes play different roles in different networks.
We discuss how studies and analyses that rely on protein-protein networks for
humans should consider their similarities and distinctions.


  We address several questions in reduced versus extended networks via the
elimination or addition of intermediate complexes in the framework of chemical
reaction networks with mass-action kinetics. We clarify and extend advances in
the literature concerning multistationarity in this context. We establish
general results about MESSI systems, which we use to compute the circuits of
multistationarity for significant biochemical networks.


  Designing novel proteins with desired functions is crucial in biology and
chemistry. However, most existing work focus on protein sequence design,
leaving protein sequence and structure co-design underexplored. In this paper,
we propose GeoPro, a method to design protein backbone structure and sequence
jointly. Our motivation is that protein sequence and its backbone structure
constrain each other, and thus joint design of both can not only avoid
nonfolding and misfolding but also produce more diverse candidates with desired
functions. To this end, GeoPro is powered by an equivariant encoder for
three-dimensional (3D) backbone structure and a protein sequence decoder guided
by 3D geometry. Experimental results on two biologically significant
metalloprotein datasets, including $\beta$-lactamases and myoglobins, show that
our proposed GeoPro outperforms several strong baselines on most metrics.
Remarkably, our method discovers novel $\beta$-lactamases and myoglobins which
are not present in protein data bank (PDB) and UniProt. These proteins exhibit
stable folding and active site environments reminiscent of those of natural
proteins, demonstrating their excellent potential to be biologically
functional.


  The design of \emph{de novo} functional proteins that bind specific ligand
molecules is paramount in therapeutics and bio-engineering. A critical yet
formidable task in this endeavor is the design of the protein pocket, which is
the cavity region of the protein where the ligand binds. Current methods are
plagued by inefficient generation, inadequate context modeling of the ligand
molecule, and the inability to generate side-chain atoms. Here, we present the
Full-Atom Iterative Refinement (FAIR) method, designed to address these
challenges by facilitating the co-design of protein pocket sequences,
specifically residue types, and their corresponding 3D structures. FAIR
operates in two steps, proceeding in a coarse-to-fine manner (transitioning
from protein backbone to atoms, including side chains) for a full-atom
generation. In each iteration, all residue types and structures are
simultaneously updated, a process termed full-shot refinement. In the initial
stage, the residue types and backbone coordinates are refined using a
hierarchical context encoder, complemented by two structure refinement modules
that capture both inter-residue and pocket-ligand interactions. The subsequent
stage delves deeper, modeling the side-chain atoms of the pockets and updating
residue types to ensure sequence-structure congruence. Concurrently, the
structure of the binding ligand is refined across iterations to accommodate its
inherent flexibility. Comprehensive experiments show that FAIR surpasses
existing methods in designing superior pocket sequences and structures,
producing average improvement exceeding 10\% in AAR and RMSD metrics. FAIR is
available at \url{https://github.com/zaixizhang/FAIR}.


  Single-cell RNA sequencing (scRNA-seq) has made significant strides in
unraveling the intricate cellular diversity within complex tissues. This is
particularly critical in the brain, presenting a greater diversity of cell
types than other tissue types, to gain a deeper understanding of brain function
within various cellular contexts. However, analyzing scRNA-seq data remains a
challenge due to inherent measurement noise stemming from dropout events and
the limited utilization of extensive gene expression information. In this work,
we introduce scHyena, a foundation model designed to address these challenges
and enhance the accuracy of scRNA-seq analysis in the brain. Specifically,
inspired by the recent Hyena operator, we design a novel Transformer
architecture called singe-cell Hyena (scHyena) that is equipped with a linear
adaptor layer, the positional encoding via gene-embedding, and a
{bidirectional} Hyena operator. This enables us to process full-length
scRNA-seq data without losing any information from the raw data. In particular,
our model learns generalizable features of cells and genes through pre-training
scHyena using the full length of scRNA-seq data. We demonstrate the superior
performance of scHyena compared to other benchmark methods in downstream tasks,
including cell type classification and scRNA-seq imputation.


  Extrusion-based 3D printing is a recent and widely used tissue engineering
tool that provides precise 3D control of bioinks to create organ-size
biomaterial-based objects composed of hierarchically organized cellularized
scaffolds. The internal organization of scaffold constituents should mimic the
structural anisotropy of targeted tissue to stimulate proliferative cellular
behavior during 3D cell culture. Both the choice of polymers constituting the
bioink and the topological properties during the extrusion process greatly
influence the structural anisotropy and cellular response of tissue engineering
constructs. The bioink used in our study was a hydrogel made of three
constituents: fibrinogen, alginate and gelatin. These components provide
biocompatibility, printability and conservation of 3D shape after printing. The
topological properties in flowing polymers are dictated by macromolecule
conformation i.e. their orientation and degree of stretch. In this study, we
used the micro-macro approach to describe the orientation state of hydrogel
macromolecules during the extrusion process, offering a two-scale description
of fluid behavior. The goal of our study was to use the Fokker-Planck equation,
which describes the evolution of probability distribution function over time,
to represent the real state of the constituent population in the representative
elementary volume within a hydrogel during extrusion-based 3D printing. Our
data suggest that for a tubular nozzle syringe, constituent orientation is
driven by a high shear rate, which overcomes the fluid rheological behavior.
Also, the interaction coefficient, which represents the microscopic interaction
between fluid particles overcomes hydrogel behavior for constituent orientation
in the prediction model.


  Aphids are damaging insect pests on many crops. Their density can rapidly
build up on a host plant to several thousand over one growing season.
Occasionally, a competition-driven decline in population early in the season,
followed by a build-up later, is observed in the field. Such dynamics cannot be
captured via standard models, such as introduced in Kindlmann and Dixon, 2010.
In Kindlmann et al., 2007, a logistic non-local population model with variable
carrying capacity is proposed to capture these alternative dynamics. The
proposed model has a rich dynamical structure and can predict multiple
population peaks, as observed in the field. We show that additionally, this
model also possesses solutions that can blow-up/explode in finite time. The
blow-up is seen to occur for both large and small initial conditions and can
occur both early and late in the season. We propose a model extension that,
under certain parametric restrictions, has global time-bounded solutions. We
discuss the ecological applications of our findings.


  This paper proposes a latent prompt Transformer model for solving challenging
optimization problems such as molecule design, where the goal is to find
molecules with optimal values of a target chemical or biological property that
can be computed by an existing software. Our proposed model consists of three
components. (1) A latent vector whose prior distribution is modeled by a Unet
transformation of a Gaussian white noise vector. (2) A molecule generation
model that generates the string-based representation of molecule conditional on
the latent vector in (1). We adopt the causal Transformer model that takes the
latent vector in (1) as prompt. (3) A property prediction model that predicts
the value of the target property of a molecule based on a non-linear regression
on the latent vector in (1). We call the proposed model the latent prompt
Transformer model. After initial training of the model on existing molecules
and their property values, we then gradually shift the model distribution
towards the region that supports desired values of the target property for the
purpose of molecule design. Our experiments show that our proposed model
achieves state of the art performances on several benchmark molecule design
tasks.


  Large Language Models (LLMs) have revolutionized the field of natural
language processing, but they fall short in comprehending biological sequences
such as proteins. To address this challenge, we propose InstructProtein, an
innovative LLM that possesses bidirectional generation capabilities in both
human and protein languages: (i) taking a protein sequence as input to predict
its textual function description and (ii) using natural language to prompt
protein sequence generation. To achieve this, we first pre-train an LLM on both
protein and natural language corpora, enabling it to comprehend individual
languages. Then supervised instruction tuning is employed to facilitate the
alignment of these two distinct languages. Herein, we introduce a knowledge
graph-based instruction generation framework to construct a high-quality
instruction dataset, addressing annotation imbalance and instruction deficits
in existing protein-text corpus. In particular, the instructions inherit the
structural relations between proteins and function annotations in knowledge
graphs, which empowers our model to engage in the causal modeling of protein
functions, akin to the chain-of-thought processes in natural languages.
Extensive experiments on bidirectional protein-text generation tasks show that
InstructProtein outperforms state-of-the-art LLMs by large margins. Moreover,
InstructProtein serves as a pioneering step towards text-based protein function
prediction and sequence design, effectively bridging the gap between protein
and human language understanding.


  Quantification of measurement uncertainty is crucial for robust scientific
inference, yet accurate estimates of this uncertainty remain elusive for
ecological measures of diversity. Here, we address this longstanding challenge
by deriving a closed-form unbiased estimator for the sampling variance of
Simpson's diversity index. In numerical tests the estimator consistently
outperforms existing approaches, particularly for scenarios in which species
richness exceeds sample size. We apply the estimator to quantify biodiversity
loss in marine ecosystems and to demonstrate ligand-dependent contributions of
T cell receptor chains to specificity, illustrating its versatility across
fields. The novel estimator provides researchers with a reliable method for
comparing diversity between samples, essential for quantifying biodiversity
trends and making informed conservation decisions.


  Understanding causal relationships within Gene Regulatory Networks (GRNs) is
essential for unraveling the gene interactions in cellular processes. However,
causal discovery in GRNs is a challenging problem for multiple reasons
including the existence of cyclic feedback loops and uncertainty that yields
diverse possible causal structures. Previous works in this area either ignore
cyclic dynamics (assume acyclic structure) or struggle with scalability. We
introduce Swift-DynGFN as a novel framework that enhances causal structure
learning in GRNs while addressing scalability concerns. Specifically,
Swift-DynGFN exploits gene-wise independence to boost parallelization and to
lower computational cost. Experiments on real single-cell RNA velocity and
synthetic GRN datasets showcase the advancement in learning causal structure in
GRNs and scalability in larger systems.


  Target proteins that lack accessible binding pockets and conformational
stability have posed increasing challenges for drug development. Induced
proximity strategies, such as PROTACs and molecular glues, have thus gained
attention as pharmacological alternatives, but still require small molecule
docking at binding pockets for targeted protein degradation (TPD). The
computational design of protein-based binders presents unique opportunities to
access undruggable targets, but have often relied on stable 3D structures or
predictions for effective binder generation. Recently, we have leveraged the
expressive latent spaces of protein language models (pLMs) for the
prioritization of peptide binders from sequence alone, which we have then fused
to E3 ubiquitin ligase domains, creating a CRISPR-analogous TPD system for
target proteins. However, our methods rely on training discriminator models for
ranking heuristically or unconditionally-derived guide peptides for their
target binding capability. In this work, we introduce PepMLM, a purely target
sequence-conditioned de novo generator of linear peptide binders. By employing
a novel masking strategy that uniquely positions cognate peptide sequences at
the terminus of target protein sequences, PepMLM tasks the state-of-the-art
ESM-2 pLM to fully reconstruct the binder region, achieving low perplexities
matching or improving upon previously-validated peptide-protein sequence pairs.
After successful in silico benchmarking with AlphaFold-Multimer, we
experimentally verify PepMLM's efficacy via fusion of model-derived peptides to
E3 ubiquitin ligase domains, demonstrating endogenous degradation of target
substrates in cellular models. In total, PepMLM enables the generative design
of candidate binders to any target protein, without the requirement of target
structure, empowering downstream programmable proteome editing applications.


  Machine learning and deep learning have been celebrating many successes in
the application to biological problems, especially in the domain of protein
folding. Another equally complex and important question has received relatively
little attention by the machine learning community, namely the one of
prediction of complex traits from genetics. Tackling this problem requires
in-depth knowledge of the related genetics literature and awareness of various
subtleties associated with genetic data. In this guide, we provide an overview
for the machine learning community on current state of the art models and
associated subtleties which need to be taken into consideration when developing
new models for phenotype prediction. We use height as an example of a
continuous-valued phenotype and provide an introduction to benchmark datasets,
confounders, feature selection, and common metrics.


  Basecalling, an essential step in many genome analysis studies, relies on
large Deep Neural Networks (DNNs) to achieve high accuracy. Unfortunately,
these DNNs are computationally slow and inefficient, leading to considerable
delays and resource constraints in the sequence analysis process. A
Computation-In-Memory (CIM) architecture using memristors can significantly
accelerate the performance of DNNs. However, inherent device non-idealities and
architectural limitations of such designs can greatly degrade the basecalling
accuracy, which is critical for accurate genome analysis. To facilitate the
adoption of memristor-based CIM designs for basecalling, it is important to (1)
conduct a comprehensive analysis of potential CIM architectures and (2) develop
effective strategies for mitigating the possible adverse effects of inherent
device non-idealities and architectural limitations.
  This paper proposes Swordfish, a novel hardware/software co-design framework
that can effectively address the two aforementioned issues. Swordfish
incorporates seven circuit and device restrictions or non-idealities from
characterized real memristor-based chips. Swordfish leverages various
hardware/software co-design solutions to mitigate the basecalling accuracy loss
due to such non-idealities. To demonstrate the effectiveness of Swordfish, we
take Bonito, the state-of-the-art (i.e., accurate and fast), open-source
basecaller as a case study. Our experimental results using Sword-fish show that
a CIM architecture can realistically accelerate Bonito for a wide range of real
datasets by an average of 25.7x, with an accuracy loss of 6.01%.


  The Gene Regulatory Network (GRN) of biological cells governs a number of key
functionalities that enables them to adapt and survive through different
environmental conditions. Close observation of the GRN shows that the structure
and operational principles resembles an Artificial Neural Network (ANN), which
can pave the way for the development of Biological Artificial Intelligence. In
particular, a gene's transcription and translation process resembles a
sigmoidal-like property based on transcription factor inputs. In this paper, we
develop a mathematical model of gene-perceptron using a dual-layered
transcription-translation chemical reaction model, enabling us to transform a
GRN into a Gene Regulatory Neural Network (GRNN). We perform stability analysis
for each gene-perceptron within the fully-connected GRNN sub network to
determine temporal as well as stable concentration outputs that will result in
reliable computing performance. We focus on a non-linear classifier application
for the GRNN, where we analyzed generic multi-layer GRNNs as well as E.Coli
GRNN that is derived from trans-omic experimental data. Our analysis found that
varying the parameters of the chemical reactions can allow us shift the
boundaries of the classification region, laying the platform for programmable
GRNNs that suit diverse application requirements.


  In the past decade, Artificial Intelligence driven drug design and discovery
has been a hot research topic, where an important branch is molecule generation
by generative models, from GAN-based models and VAE-based models to the latest
diffusion-based models. However, most existing models pursue only the basic
properties like validity and uniqueness of the generated molecules, a few go
further to explicitly optimize one single important molecular property (e.g.
QED or PlogP), which makes most generated molecules little usefulness in
practice. In this paper, we present a novel approach to generating molecules
with desirable properties, which expands the diffusion model framework with
multiple innovative designs. The novelty is two-fold. On the one hand,
considering that the structures of molecules are complex and diverse, and
molecular properties are usually determined by some substructures (e.g.
pharmacophores), we propose to perform diffusion on two structural levels:
molecules and molecular fragments respectively, with which a mixed Gaussian
distribution is obtained for the reverse diffusion process. To get desirable
molecular fragments, we develop a novel electronic effect based fragmentation
method. On the other hand, we introduce two ways to explicitly optimize
multiple molecular properties under the diffusion model framework. First, as
potential drug molecules must be chemically valid, we optimize molecular
validity by an energy-guidance function. Second, since potential drug molecules
should be desirable in various properties, we employ a multi-objective
mechanism to optimize multiple molecular properties simultaneously. Extensive
experiments with two benchmark datasets QM9 and ZINC250k show that the
molecules generated by our proposed method have better validity, uniqueness,
novelty, Fr\'echet ChemNet Distance (FCD), QED, and PlogP than those generated
by current SOTA models.


  Biomolecular condensates formed through the phase separation of proteins and
nucleic acids are widely observed, offering a fundamental means of organizing
intracellular materials in a membrane-less fashion. Traditionally, these
condensates have been regarded as homogeneous isotropic liquids. However, in
analogy with some synthetic copolymer systems, our recent theoretical research
has demonstrated that model biomolecular condensates can exhibit a
microemulsion-like internal structure, contingent upon the specific sequence,
inter-chain site-site interactions, and concentrated phase polymer density. In
this study, we present a microscopic dynamical theory for the self-diffusion
constant and viscosity of concentrated unentangled A/B regular multiblock
copolymer solutions. Our approach integrates static equilibrium local and
microdomain scale structural information obtained from PRISM integral equation
theory and the time evolution of the autocorrelation function of monomer scale
forces at the center-of-mass level that determine the polymer diffusion
constant and viscosity in a weak caging regime far from a glass or gel
transition. We focus on regular multi-block systems both for simplicity and for
its relevance to synthetic macromolecular science. The impact of sequence and
inter-chain attraction strength on the slowing down of copolymer mass transport
and flow due to local clustering enhanced collisional friction and retardation
of motion due to emergent microdomain scale ordering are established. Analytic
analysis and metrics employed in the study of biomolecular condensates are
employed to identify key order parameters that quantity how attractive forces,
packing structure, multiblock sequence, and copolymer density determine
dynamical slowing down above and below the crossover to a fluctuating polymeric
microemulsion state.


  Dominance is usually considered a constant value that describes the relative
difference in fitness or phenotype between heterozygotes and the average of
homozygotes at a focal polymorphic locus. However, the observed dominance can
vary with the genetic background of the focal locus. Here, alleles at other
loci modify the observed phenotype through position effects or dominance
modifiers that are sometimes associated with pathogen resistance, lineage, sex,
or mating type. Theoretical models have illustrated how variable dominance
appears in the context of multi-locus interaction (epistasis). Here, we review
empirical evidence for variable dominance and how the observed patterns may be
captured by proposed epistatic models. We highlight how integrating epistasis
and dominance is crucial for comprehensively understanding adaptation and
speciation.


  Motivated by a general principle governing regulation mechanisms in
biological cells, we investigate a general interaction scheme between different
populations of particles and specific particles, referred to as agents.
Assuming that each particle follows a random path in the medium, when a
particle and an agent meet, they may bind and form a pair which has some
specific functional properties. Such a pair is also subject to random events
and it splits after some random amount of time. In a stochastic context, using
a Markovian model for the vector of the number of paired particles, and by
taking the total number of particles as a scaling parameter, we study the
asymptotic behavior of the time evolution of the number of paired particles.
Two scenarios are investigated: one with a large but fixed number of agents,
and the other one, the dynamic case, when agents are created at a bounded rate
and may die after some time when they are not paired. A first order limit
theorem is established for the time evolution of the system in both cases. The
proof of an averaging principle of the dynamic case is one of the main
contributions of the paper. Limit theorems for fluctuations are obtained in the
case of a fixed number agents. The impact of dynamical arrivals of agents on
the level of pairing of the system is discussed.


  Nanopore-based sequencers generate a series of raw electrical signal
measurements that represent the contents of a biological sequence molecule
passing through the sequencer's nanopore. If the raw signal is analyzed in
real-time, an irrelevant molecule can be ejected from the nanopore before it is
completely sequenced, reducing sequencing time. To meet the low-latency and
high-throughput requirements of the real-time analysis, a number of recent
works propose the direct analysis of raw nanopore signals instead of
traditional basecalling-based analysis approaches. We observe that while
existing proposals for raw signal read mapping typically do well in all metrics
for small reference databases (e.g., viral genomes), they all fail to scale to
large reference databases (e.g., the human genome) in some aspect. Our goal is
to analyze raw nanopore signals with high accuracy, high throughput, low
latency, low memory usage, and needing few bases to be sequenced for a wide
range of reference database sizes. To this end, we propose RawAlign, the first
Seed-Filter-Align mapper for raw nanopore signals. Our evaluation shows that
RawAlign is the only tool that can map raw nanopore signals to large reference
databases $\geq$3117Mbp with high accuracy. Our evaluation shows that RawAlign
generalizes well to a wide range of reference database sizes. In particular,
RawAlign has a similar throughput to the overall prior state-of-the-art RawHash
(between 0.80$\times$-1.08$\times$) while improving accuracy on all datasets
(between 1.02$\times$-1.64$\times$ F-1 score). RawAlign provides a 2.83$\times$
(2.06$\times$) speedup over Uncalled (Sigmap) on average (geo. mean) while
improving accuracy by 1.35$\times$ (1.34$\times$) in terms of F-1 score on
average (geo. mean). Availability: https://github.com/cmu-safari/RawAlign


  We study the bacterial respiration through the numerical solution of the
Fairen-Velarde coupled nonlinear differential equations. The instantaneous
concentrations of the oxygen and the nutrients are computed. The fixed point
solution and the stable limit cycle are found in different parameter ranges as
predicted by the linearized differential equations. In a specified range of
parameters, it is observed that the system spends some time near the stable
limit cycle and eventually reaches the stable fixed point. This metastability
has been investigated systematically. Interestingly, it is observed that the
system exhibits two distinctly different time scales in reaching the stable
fixed points. The slow time scale of the metastable lifetime near the stable
limit cycle and a fast time scale (after leaving the zone of limit cycle) in
rushing towards the stable fixed point. The gross residence time, near the
limit cycle (described by a slow time scale), can be reduced by varying the
concentrations of nutrients. This idea can be used to control the harmful
metastable lifespan of active bacteria.


  West Nile virus (WNV) is prevalent in the United States but it shows
considerable divergence in transmission patterns and spatio-temporal
intensity.It is to be noted that the mechanism that drives the transmission
potential of WNV is described by the abilities of host species to maintain and
disseminate the pathogens pertinent with different eco-epidemiological factors
that have an influence on the contact rates amongst the interacting
species.There is growing evidence that several vectors exhibit strong feeding
preferences towards different host communities.We construct a process based
weather driven ordinary differential equation (ODE) model to understand the
impact of one vector species Culex pipiens, preferred avian and non-preferred
human hosts and compared it surveillance data for the Culex pipiens complex
collected in Cook County, Illinois, USA.In our mechanistic model, we also
demonstrate that adulticide treatments produced significant reductions in the
Culex pipiens population.We take into account the feeding index that can be
described as the ratio between observed frequency of mosquitoes feeding on one
host compared to another host, divided by the expected frequency of mosquitoes
feeding on these two hosts based on the presence of the particular hosts to
develop this transmission model for WNV. Our findings demonstrate that the
interplay between the feeding index and mosquito abatement strategy is rather a
complex phenomenon and it induces a heterogeneous contact rates that should be
included while modelling multi-host, multi-vector transmission model.


  The practice of marriage is an understudied phenomenon in behavioural
sciences despite being ubiquitous across human cultures. This modelling paper
shows that replacing distant direct kin with in-laws increases the
interconnectedness of the family social network graph, which allows more
cooperative and larger groups. In this framing, marriage can be seen as a
social technology that reduces free-riding within collaborative group. This
approach offers a solution to the puzzle of why our species has this particular
form of regulating mating behaviour, uniquely among pair-bonded animals.


  This paper aims to shed light on the evolutionary dynamics of diverse and
social populations by introducing the rich expressiveness of generative models
into the trait expression of social agent-based evolutionary models.
Specifically, we focus on the evolution of personality traits in the context of
a game-theoretic relationship as a situation in which inter-individual
interests exert strong selection pressures. We construct an agent model in
which linguistic descriptions of personality traits related to cooperative
behavior are used as genes. The deterministic strategies extracted from Large
Language Model (LLM) that make behavioral decisions based on these personality
traits are used as behavioral traits. The population is evolved according to
selection based on average payoff and mutation of genes by asking LLM to
slightly modify the parent gene toward cooperative or selfish. Through
preliminary experiments and analyses, we clarify that such a model can indeed
exhibit the evolution of cooperative behavior based on the diverse and
higher-order representation of personality traits. We also observed the
repeated intrusion of cooperative and selfish personality traits through
changes in the expression of personality traits, and found that the emerging
words in the evolved gene well reflected the behavioral tendency of its
personality in terms of their semantics.


  We demonstrate and characterize a first-principles approach to modeling the
mass action dynamics of metabolism. Starting from a basic definition of entropy
expressed as a multinomial probability density using Boltzmann probabilities
with standard chemical potentials, we derive and compare the free energy
dissipation and the entropy production rates. We express the relation between
the entropy production and the chemical master equation for modeling
metabolism, which unifies chemical kinetics and chemical thermodynamics.
Subsequent implementation of an maximum free energy dissipation model for
systems of coupled reactions is accomplished by using an approximation to the
Marcelin equation for mass action kinetics that maximizes the entropy
production. Because prediction uncertainty with respect to parameter
variability is frequently a concern with mass action models utilizing rate
constants, we compare and contrast the maximum entropy production model, which
has its own set of rate parameters, to a population of standard mass action
models in which the rate constants are randomly chosen. We show that a maximum
entropy production model is characterized by a high probability of free energy
dissipation rate, and likewise entropy production rate, relative to other
models. We then characterize the variability of the maximum entropy production
predictions with respect to uncertainties in parameters (standard free energies
of formation) and with respect to ionic strengths typically found in a cell.


  Recombinant protein production is pivotal in molecular biology, enabling
profound insights into cellular processes through biophysical, biochemical, and
structural analyses of the purified samples. The demand for substantial
biomolecule quantities often presents challenges, particularly for eukaryotic
proteins. Escherichia coli expression systems have evolved to address these
issues, offering advanced features such as solubility tags, posttranslational
modification capabilities, and modular plasmid libraries. Nevertheless,
existing tools are often complex, which limits their accessibility and
necessitates streamlined systems for rapid screening under standardized
conditions. Based on the Golden Gate cloning method, we have developed a simple
'one-pot' approach for the generation of expression constructs using
strategically chosen tags like hexahistidine, SUMO, MBP, GST, and GB1 to
enhance solubility and expression. Tags are removable via TEV protease
cleavage, and the system allows visual cloning verification through mScarlet
fluorescence. We provide a comprehensive protocol encompassing oligonucleotide
design, cloning, expression, Ni-NTA affinity chromatography, and size-exclusion
chromatography. This method therefore streamlines prokaryotic and eukaryotic
protein production, rendering it accessible to standard molecular biology
laboratories with basic protein biochemical equipment.


  Infections are known to interact as previous infections may have an effect on
risk of succumbing to a new infection. The co-dynamics can be mediated by
immunosuppression or -modulation, shared environmental or climatic drivers, or
competition for susceptible hosts. Research and statistical methods in
epidemiology often concentrate on large pooled datasets, or high quality data
from cities, leaving rural areas underrepresented in literature. Data
considering rural populations are typically sparse and scarce, especially in
the case of historical data sources, which may introduce considerable
methodological challenges. In order to overcome many obstacles due to such
data, we present a general Bayesian spatio-temporal model for disease
co-dynamics. Applying the proposed model on historical (1820-1850) Finnish
parish register data, we study the spread of infectious diseases in
pre-healthcare Finland. We observe that measles, pertussis and smallpox exhibit
positively correlated dynamics such that any new infection increased mortality
in all three diseases, indicating possibly general immunosuppressive effects at
population level.


  Numerous cellular functions rely on protein$\unicode{x2013}$protein
interactions. Efforts to comprehensively characterize them remain challenged
however by the diversity of molecular recognition mechanisms employed within
the proteome. Deep learning has emerged as a promising approach for tackling
this problem by exploiting both experimental data and basic biophysical
knowledge about protein interactions. Here, we review the growing ecosystem of
deep learning methods for modeling protein interactions, highlighting the
diversity of these biophysically-informed models and their respective
trade-offs. We discuss recent successes in using representation learning to
capture complex features pertinent to predicting protein interactions and
interaction sites, geometric deep learning to reason over protein structures
and predict complex structures, and generative modeling to design de novo
protein assemblies. We also outline some of the outstanding challenges and
promising new directions. Opportunities abound to discover novel interactions,
elucidate their physical mechanisms, and engineer binders to modulate their
functions using deep learning and, ultimately, unravel how protein interactions
orchestrate complex cellular behaviors.


  Modeling the interaction between proteins and ligands and accurately
predicting their binding structures is a critical yet challenging task in drug
discovery. Recent advancements in deep learning have shown promise in
addressing this challenge, with sampling-based and regression-based methods
emerging as two prominent approaches. However, these methods have notable
limitations. Sampling-based methods often suffer from low efficiency due to the
need for generating multiple candidate structures for selection. On the other
hand, regression-based methods offer fast predictions but may experience
decreased accuracy. Additionally, the variation in protein sizes often requires
external modules for selecting suitable binding pockets, further impacting
efficiency. In this work, we propose $\mathbf{FABind}$, an end-to-end model
that combines pocket prediction and docking to achieve accurate and fast
protein-ligand binding. $\mathbf{FABind}$ incorporates a unique ligand-informed
pocket prediction module, which is also leveraged for docking pose estimation.
The model further enhances the docking process by incrementally integrating the
predicted pocket to optimize protein-ligand binding, reducing discrepancies
between training and inference. Through extensive experiments on benchmark
datasets, our proposed $\mathbf{FABind}$ demonstrates strong advantages in
terms of effectiveness and efficiency compared to existing methods. Our code is
available at https://github.com/QizhiPei/FABind


  Score-based generative models (SGMs) have proven to be powerful tools for
designing new proteins. Designing proteins that bind a pre-specified target is
highly relevant to a range of medical and industrial applications. Despite the
flurry of new SGMs in the last year, there has been little systematic
exploration of the impact of design choices in SGMs for protein design. Here we
present LoopGen, a flexible SGM framework for the design of short binding
peptide structures. We apply our framework to design antibody binding loop
structures conditional on a target epitope and evaluate a variety of modelling
choices in SGM-based protein design. We demonstrate that modelling residue
orientations in addition to positions improves not only the quality of the
output structures but also their diversity. Additionally, we identify variance
schedules that result in significant performance improvements and observe
patterns that may motivate the development of better schedules for protein
design. Finally, we develop three novel tests to evaluate whether the model
generates structures that are appropriately conditioned on an epitope,
demonstrating that LoopGen's generated structures are dependent on the
structure, sequence, and position of the epitope. Our findings will help guide
future development and evaluation of generative models for binding proteins.


  Diffusiophoresis, which refers to the movement of entities driven by
gradients in the concentration of attractants, is observable in colloids and
chemotactic bacteria. We suggest that humans also exhibit diffusiophoresis when
they perceive concentration gradients in economic opportunities, social
connections, safety factors, and more. Consequently, we build upon the
Fisher-KPP reaction-diffusion model, which is foundational to human population
dynamics, to incorporate diffusiophoretic advection. Through simulations, we
demonstrate that diffusiophoresis can predict the emergence of population
hotspots from initially dispersed populations and can enable precise control
over inter- and intra-population segregation, a capability not achievable
through Fisher-KPP alone. This framework is particularly pertinent given the
anticipated impacts of climate change, which may result in significant human
displacement influenced by the diffusiophoretic movement of individuals in
response to concentration gradients in safety conditions.


  Intraspecific trait variation has been increasingly recognized as an
important factor in determining species interaction and diversity.
Eco-evolutionary models have studied the distribution of trait values within a
population that changes over the generations as a result of selection and
heritability. Non-heritable traits that can change within the lifetime, such as
behavior, can cause trait-mediated indirect effects, often studied by modeling
the dynamics of a homogeneous trait. Complementary to these two approaches, we
study the distribution of traits within a population and its dynamics on short
timescales due to ecological processes. We consider several mechanisms by which
the trait distribution can shift dynamically: phenotypic plasticity within each
individual, differential growth among individuals, and preferential consumption
by the predator. Through a simple predator-prey model that explicitly tracks
the trait distribution within the prey, we illustrate some unusual consequences
of trait shifts, including the healthy herd effect and the emergent promotion
of the prey by the predator. We analyze the cause of these phenomena by
identifying the density and trait effects from the predator. Our results show
that intraspecific trait variation and the dynamics of its distribution can
lead to the modification of species interactions and result in outcomes that
are otherwise unexpected.


  Malaria is one of the most common mosquito-borne diseases widespread in
tropical and subtropical regions, causing thousands of deaths every year in the
world. In a previous paper, we formulated an age-structured model containing
three structural variables: (i) the chronological age of human and mosquito
populations, (ii) the time since they are infected, and (iii) humans waning
immunity (i.e. the progressive loss of protective antibodies after recovery).
In the present paper, we expand the analysis of this age-structured model and
focus on the derivation of entomological and epidemiological results commonly
used in the literature, following the works of Smith and McKenzie. We
generalize their results to the age-structured case. In order to quantify the
impact of neglecting structuring variables such as chronological age, we
assigned values from the literature to our model parameters. While some
parameters values are readily accessible from the literature, at least those
about the human population, the parameters concerning mosquitoes are less
commonly documented and the values of a number of them (e.g. mosquito survival
in the presence or in absence of infection) can be discussed extensively.


  Recent advancements in biological research leverage the integration of
molecules, proteins, and natural language to enhance drug discovery. However,
current models exhibit several limitations, such as the generation of invalid
molecular SMILES, underutilization of contextual information, and equal
treatment of structured and unstructured knowledge. To address these issues, we
propose $\mathbf{BioT5}$, a comprehensive pre-training framework that enriches
cross-modal integration in biology with chemical knowledge and natural language
associations. $\mathbf{BioT5}$ utilizes SELFIES for $100%$ robust molecular
representations and extracts knowledge from the surrounding context of
bio-entities in unstructured biological literature. Furthermore,
$\mathbf{BioT5}$ distinguishes between structured and unstructured knowledge,
leading to more effective utilization of information. After fine-tuning, BioT5
shows superior performance across a wide range of tasks, demonstrating its
strong capability of capturing underlying relations and properties of
bio-entities. Our code is available at
$\href{https://github.com/QizhiPei/BioT5}{Github}$.


  Increasing cultivated lands, crop homogenization and global food trade have
fostered the spread of crop pests and diseases. Optimizing crop protection is
urgently needed to ensure food safety. One aspect of crop protection is
surveillance, which focuses on the early detection of a parasite, and control,
aiming to fight and possibly eradicate it. Network theory has been widely used
to model the spread of human and animal infectious diseases in systems
described through nodes and edges. It has been successfully used to optimize
monitoring and immunization campaigns. In crop protection, there is a growing
literature using this theory to describe parasites spread and to conceive
protection strategies. Here we review the use of network theory in crop
protection, from the more descriptive to the more applied approaches aimed to
optimize crop protection. We retrace the logical process that has led
epidemiological models to rely on network theory, and we provide examples of
how the spread of crop parasites has been represented via a network
description. We define the objectives of surveillance and control, and we show
how these have been declined in the network-based epidemiological sphere and
then adapted in the agricultural context. We eventually discuss the discrepancy
between the application of network theory in surveillance and control to
identify culprits and solutions. We find that: i) scientists have successfully
interpreted very different modes of parasitic transmission under the lens of
network theory; ii) while network-based surveillance has progressively
clarified its objectives and sound tools have been proposed, network-based
control has been less studied and applied; iii) network-thinking must address
how to properly define edges and nodes at different geographic scale to broad
its application in crop protection.


  Design patterns are generalized solutions to frequently recurring problems.
They were initially developed by architects and computer scientists to create a
higher level of abstraction for their designs. Here, we extend these concepts
to cell biology in order to lend a new perspective on the evolved designs of
cells' underlying reaction networks. We present a catalog of 21 design patterns
divided into three categories: creational patterns describe processes that
build the cell, structural patterns describe the layouts of reaction networks,
and behavioral patterns describe reaction network function. Applying this
pattern language to the E. coli central metabolic reaction network, the yeast
pheromone response signaling network, and other examples lends new insights
into these systems.


  In many situations, it would be useful to know not just the best phylogenetic
tree for a given data set, but the collection of high-quality trees. This goal
is typically addressed using Bayesian techniques, however, current Bayesian
methods do not scale to large data sets. Furthermore, for large data sets with
relatively low signal one cannot even store every good tree individually,
especially when the trees are required to be bifurcating. In this paper, we
develop a novel object called the "history subpartition directed acyclic graph"
(or "history sDAG" for short) that compactly represents an ensemble of trees
with labels (e.g. ancestral sequences) mapped onto the internal nodes. The
history sDAG can be built efficiently and can also be efficiently trimmed to
only represent maximally parsimonious trees. We show that the history sDAG
allows us to find many additional equally parsimonious trees, extending
combinatorially beyond the ensemble used to construct it. We argue that this
object could be useful as the "skeleton" of a more complete uncertainty
quantification.


  Background: Missing data is a common challenge in mass spectrometry-based
metabolomics, which can lead to biased and incomplete analyses. The integration
of whole-genome sequencing (WGS) data with metabolomics data has emerged as a
promising approach to enhance the accuracy of data imputation in metabolomics
studies. Method: In this study, we propose a novel method that leverages the
information from WGS data and reference metabolites to impute unknown
metabolites. Our approach utilizes a multi-view variational autoencoder to
jointly model the burden score, polygenetic risk score (PGS), and linkage
disequilibrium (LD) pruned single nucleotide polymorphisms (SNPs) for feature
extraction and missing metabolomics data imputation. By learning the latent
representations of both omics data, our method can effectively impute missing
metabolomics values based on genomic information. Results: We evaluate the
performance of our method on empirical metabolomics datasets with missing
values and demonstrate its superiority compared to conventional imputation
techniques. Using 35 template metabolites derived burden scores, PGS and
LD-pruned SNPs, the proposed methods achieved r2-scores > 0.01 for 71.55% of
metabolites. Conclusion: The integration of WGS data in metabolomics imputation
not only improves data completeness but also enhances downstream analyses,
paving the way for more comprehensive and accurate investigations of metabolic
pathways and disease associations. Our findings offer valuable insights into
the potential benefits of utilizing WGS data for metabolomics data imputation
and underscore the importance of leveraging multi-modal data integration in
precision medicine research.