  The Worldwide spread of the Omicron lineage variants has now been confirmed.
It is crucial to understand the process of cellular life and to discover new
drugs need to identify the important proteins in a protein interaction network
(PPIN). PPINs are often represented by graphs in bioinformatics, which describe
cell processes. There are some proteins that have significant influences on
these tissues, and which play a crucial role in regulating them. The discovery
of new drugs is aided by the study of significant proteins. These significant
proteins can be found by reducing the graph and using graph analysis. Studies
examining protein interactions in the Omicron lineage (B.1.1.529) and its
variants (BA.5, BA.4, BA.3, BA.2, BA.1.1, BA.1) are not yet available. Studying
Omicron has been intended to find a significant protein. 68 nodes represent 68
proteins and 52 edges represent the relationship among the protein in the
network. A few entrality measures are computed namely page rank centrality
(PRC), degree centrality (DC), closeness centrality (CC), and betweenness
centrality (BC) together with node degree and Local Clustering Co-efficient
(LCC). We also discover 18 network clusters using Markov clustering. 8
significant proteins (candidate gene of Omicron lineage variants) were detected
among the 68 proteins, including AHSG, KCNK1, KCNQ1, MAPT, NR1H4, PSMC2, PTPN11
and, UBE21 which scored the highest among the Omicron proteins. It is found
that in the variant of Omicron protein-protein interaction networks, the MAPT
protein's impact is the most significant.


  Hub genes play an essential role in biological systems because of their
interaction with other genes. A vocabulary used in bioinformatics called Gene
Ontology (GO) describes how genes and proteins operate. This flexible ontology
illustrates the operation of molecular, biological, and cellular processes
(Pmol, Pbio, Pcel). There are various methodologies that can be analyzed to
determine semantic similarity. Research in this study, we employ the jack-knife
method by taking into account 4 well-liked Semantic similarity measures namely
Jaccard similarity, Cosine similarity, Pairsewise document similarity, and
Levenshtein distance. Based on these similarity values, the protein-protein
interaction network (PPI) of Malaria GO (Gene Ontology) data is built, which
causes clusters of identical or related protein complexes (Px) to form. The hub
nodes of the network are these necessary proteins. We use a variety of
centrality measures to establish clusters of these networks in order to
determine which node is the most important. The clusters' unique formation
makes it simple to determine which class of Px they are allied to.


  Gene set analysis is a mainstay of functional genomics, but it relies on
manually curated databases of gene functions that are incomplete and unaware of
biological context. Here we evaluate the ability of OpenAI's GPT-4, a Large
Language Model (LLM), to develop hypotheses about common gene functions from
its embedded biomedical knowledge. We created a GPT-4 pipeline to label gene
sets with names that summarize their consensus functions, substantiated by
analysis text and citations. Benchmarking against named gene sets in the Gene
Ontology, GPT-4 generated very similar names in 50% of cases, while in most
remaining cases it recovered the name of a more general concept. In gene sets
discovered in 'omics data, GPT-4 names were more informative than gene set
enrichment, with supporting statements and citations that largely verified in
human review. The ability to rapidly synthesize common gene functions positions
LLMs as valuable functional genomics assistants.


  Biology is perhaps the most complex of the sciences, given the incredible
variety of chemical species that are interconnected in spatial and temporal
pathways that are daunting to understand. Their interconnections lead to
emergent properties such as memory, consciousness, and recognition of self and
non-self. To understand how these interconnected reactions lead to cellular
life characterized by activation, inhibition, regulation, homeostasis, and
adaptation, computational analyses and simulations are essential, a fact
recognized by the biological communities. At the same time, students struggle
to understand and apply binding and kinetic analyses for the simplest reactions
such as the irreversible first-order conversion of a single reactant to a
product. This likely results from cognitive difficulties in combining
structural, chemical, mathematical, and textual descriptions of binding and
catalytic reactions. To help students better understand dynamic reactions and
their analyses, we have introduced two kinds of interactive graphs and
simulations into the online educational resource, Fundamentals of Biochemistry,
a multivolume biochemistry textbook that is part of the LibreText collection.
One type is available for simple binding and kinetic reactions. The other
displays progress curves (concentrations vs time) for both simple reactions and
more complex metabolic and signal transduction pathways, including those
available through databases using systems biology markup language (SBML) files.
Users can move sliders to change dissociation and kinetic constants as well as
initial concentrations and see instantaneous changes in the graphs. They can
also export data into a spreadsheet for further processing, such as producing
derivative Lineweaver-Burk and traditional Michaelis-Menten graphs of initial
velocity (v0) vs substrate concentration.


  Plant electrical signals have been shown to be generated in response to
various environmental stresses, but the relationship between these signals and
stress tolerance is not well understood. In this study, we used the Arabidopsis
STTM165/166 mutant, which exhibits enhanced temperature tolerance, to examine
this relationship. Surface recording techniques were utilized to compare the
generation ratio and duration characteristics of electrical signals in the
STTM165/166 mutant and wild type (WT). Patch-clamp recording was employed to
assess ion channel currents, specifically those of calcium ions. The current
intensity of the mutant was found to be lower than that of the WT. As calcium
ions are involved in the generation of plant electrical signals, we
hypothesized that the reduced calcium channel activity in the mutant increased
its electrical signal threshold. RNA-Seq analysis revealed differential
expression of AHA genes in the STTM165/166 mutant, which may contribute to the
prolonged depolarization phenotype. Gene Ontology enrichment of differentially
expressed genes (DEGs) identified associations between these DEGs and various
stresses, including temperature, salt, and those related to the jasmonic acid
and abscisic acid pathways. These findings provide experimental evidence for
the use of plant electrical signals in characterizing stress tolerance and
explore potential ion mechanisms through patch-clamp recording and DEG Gene
Ontology analysis. They also emphasize the need for further research on the
relationship between plant electrical signals and stress tolerance.


  ds-RNA and standard ds-DNA show specific structural differences in their
di-nucleotide steps, piled along the helical axis. Modeling the helices of
short fragments by a 3D mesoscopic Hamiltonian model, I use path integral
techniques to compute the average helical repeat and show that these structural
features are at the origin of the opposite twist-stretch patterns of the A- and
B- form.


  We propose and investigate an extension of the Caspar-Klug symmetry
principles for viral capsid assembly to the programmable assembly of
size-controlled triply-periodic polyhedra, discrete variants of the Primitive,
Diamond, and Gyroid cubic minimal surfaces. Inspired by a recent class of
programmable DNA origami colloids, we demonstrate that the economy of design in
these crystalline assemblies -- in terms of the growth of the number of
distinct particle species required with the increased size-scale (e.g.
periodicity) -- is comparable to viral shells. We further test the role of
geometric specificity in these assemblies via dynamical assembly simulations,
which show that conditions for simultaneously efficient and high-fidelity
assembly require an intermediate degree of flexibility of local angles and
lengths in programmed assembly. Off-target misassembly occurs via incorporation
of a variant of disclination defects, generalized to the case of hyperbolic
crystals. The possibility of these topological defects is a direct consequence
of the very same symmetry principles that underlie the economical design,
exposing a basic tradeoff between design economy and fidelity of programmable,
size controlled assembly.


  Gene regulatory networks, as a powerful abstraction for describing complex
biological interactions between genes through their expression products within
a cell, are often regarded as virtually deterministic dynamical systems.
However, this view is now being challenged by the fundamentally stochastic,
'bursty' nature of gene expression revealed at the single cell level. We
present a Python package called Harissa which is dedicated to simulation and
inference of such networks, based upon an underlying stochastic dynamical model
driven by the transcriptional bursting phenomenon. As part of this tool,
network inference can be interpreted as a calibration procedure for a
mechanistic model: once calibrated, the model is able to capture the typical
variability of single-cell data without requiring ad hoc external noise, unlike
ordinary or even stochastic differential equations frequently used in this
context. Therefore, Harissa can be used both as an inference tool, to
reconstruct biologically relevant networks from time-course scRNA-seq data, and
as a simulation tool, to generate quantitative gene expression profiles in a
non-trivial way through gene interactions.


  Maxwell defined a true or total current in a way not widely used today. He
said "... true electric current ... is not the same thing as the current of
conduction but that the time-variation of the electric displacement must be
taken into account in estimating the total movement of electricity". We show
that true or total current is a universal property of electrodynamics
independent of properties of matter. We use mathematics without a dielectric
constant. The resulting Maxwell Current Law is a generalization of the
Kirchhoff Law of Current used in circuit analysis, that also includes
displacement current. The generalization is not a long-time low frequency
approximation in contrast to traditional presentation of Kirchhoff's Law. The
Maxwell Current Law does not require currents to be in circuits. It has been
applied to three dimensional systems like the signaling system of nerve and
muscle fibers. The Maxwell Current Law clarifies the flow of electrons,
protons, and ions in mitochondria that generate ATP, the molecule that stores
chemical energy throughout life. The currents are globally coupled because
mitochondria are short. Focusing on Maxwell current reinterprets the classical
chemiosmotic hypothesis of ATP production. The conduction current of protons in
mitochondria is driven by the protonmotive force including its component
electrical potential, just as in the classical chemiosmotic hypothesis. The
electrical potential is now the electrical potential as defined in physical
sciences by Maxwell partial differential equations. The conduction current is
now just a part of the true current analyzed by Maxwell. Details of
accumulation of charges do not have to be considered in analysis of true
current because true current does not accumulate. It is true total current that
provides the energy that generates the ATP, not just the protonmotive force.


  Summary: Raw nanopore signals can be analyzed while they are being generated,
a process known as real-time analysis. Real-time analysis of raw signals is
essential to utilize the unique features that nanopore sequencing provides,
enabling the early stopping of the sequencing of a read or the entire
sequencing run based on the analysis. The state-of-the-art mechanism, RawHash,
offers the first hash-based efficient and accurate similarity identification
between raw signals and a reference genome by quickly matching their hash
values. In this work, we introduce RawHash2, which provides major improvements
over RawHash, including a more sensitive chaining implementation, weighted
mapping decisions, frequency filters to reduce ambiguous seed hits, minimizers
for hash-based sketching, and support for the R10.4 flow cell version and
various data formats such as POD5. Compared to RawHash, RawHash2 provides
better F1 accuracy (on average by 3.44% and up to 10.32%) and better throughput
(on average by 2.3x and up to 5.4x) than RawHash.
  Availability and Implementation: RawHash2 is available at
https://github.com/CMU-SAFARI/RawHash. We also provide the scripts to fully
reproduce our results on our GitHub page.


  The incredible capabilities of generative artificial intelligence models have
inevitably led to their application in the domain of drug discovery. Within
this domain, the vastness of chemical space motivates the development of more
efficient methods for identifying regions with molecules that exhibit desired
characteristics. In this work, we present a computationally efficient active
learning methodology that requires evaluation of only a subset of the generated
data in the constructed sample space to successfully align a generative model
with respect to a specified objective. We demonstrate the applicability of this
methodology to targeted molecular generation by fine-tuning a GPT-based
molecular generator toward a protein with FDA-approved small-molecule
inhibitors, c-Abl kinase. Remarkably, the model learns to generate molecules
similar to the inhibitors without prior knowledge of their existence, and even
reproduces two of them exactly. We also show that the methodology is effective
for a protein without any commercially available small-molecule inhibitors, the
HNH domain of the CRISPR-associated protein 9 (Cas9) enzyme. We believe that
the inherent generality of this method ensures that it will remain applicable
as the exciting field of in silico molecular generation evolves. To facilitate
implementation and reproducibility, we have made all of our software available
through the open-source ChemSpaceAL Python package.


  The actomyosin cortex is an active material that provides animal cells with a
strong but flexible exterior, whose mechanics, including non-Gaussian
fluctuations and occasional large displacements or cytoquakes, have defied
explanation. We study the active nanoscale fluctuations of the cortex using
high-performance tracking of an array of flexible microposts adhered to
multiple cultured cell types. When the confounding effects of static
heterogeneity and tracking error are removed, the fluctuations are found to be
heavy-tailed and well-described by a truncated L\'evy stable distribution over
a wide range of timescales and multiple cell types. Notably, cytoquakes appear
to correspond to the largest random displacements, unifying all cortical
fluctuations into a single spectrum. These findings reinforce the cortex's
previously noted similarity to soft glassy materials such as foams, while the
form of the fluctuation distribution will constrain future models of the
cytoskeleton.


  In diseases with long-term immunity, vaccination is known to increase the
average age at infection as a result of the decrease in the pathogen
circulation. This implies that a vaccination campaign can have negative effects
when a disease is more costly (financial or health-related costs) for higher
ages. This work considers an age-structured population transmission model with
imperfect vaccination. Our aim is to compare the social and individual costs of
vaccination, assuming that disease costs are age-dependent. A model coupling
pathogen deterministic dynamics for a population consisting of juveniles and
adults, both assumed to be rational agents, is introduced. The parameter region
for which vaccination has a positive social impact is fully characterized and
the Nash equilibrium of the vaccination game is obtained. Finally, collective
strategies designed to promote voluntary vaccination, without compromising
social welfare, are discussed.


  Genetic mutations are footprints of cancer evolution and reveal critical
dynamic parameters of tumour growth, which otherwise are hard to measure in
vivo. The mutation accumulation in tumour cell populations has been described
by various statistics, such as site frequency spectra (SFS) from bulk or
single-cell data, as well as single-cell division distributions (DD) and
mutational burden distributions (MBD). An integrated understanding of these
distributions obtained from different sequencing information is important to
illuminate the ecological and evolutionary dynamics of tumours, and requires
novel mathematical and computational tools. We introduce dynamical matrices to
analyse and unite the SFS, DD and MBD based on a birth-death process. Using the
Markov nature of the model, we derive recurrence relations for the expectations
of these three distributions. While recovering classic exact results in
pure-birth cases for the SFS and the DD through our new framework, we also
derive a new expression for the MBD as well as approximations for all three
distributions when death is introduced, confirming our results with stochastic
simulations. Moreover, we demonstrate a natural link between the SFS and the
single-cell MBD, and show that the MBD can be regenerated through the DD.
Surprisingly, the single-cell MBD is mainly driven by the stochasticity arising
in the DD, rather than the extra stochasticity in the number of mutations at
each cell division.


  We consider the reconciliation problem, in which the task is to find a
mapping of a gene tree into a species tree, so as to maximize the likelihood of
such fitting, given the available data. We describe a model for the evolution
of the species tree, a subfunctionalisation model for the evolution of the gene
tree, and provide an algorithm to compute the likelihood of the reconciliation.
We derive our results using the theory of matrix-analytic methods and describe
efficient algorithms for the computation of a range of useful metrics. We
illustrate the theory with examples and provide the physical interpretations of
the discussed quantities, with a focus on the practical applications of the
theory to incomplete data.


  In this paper, we use the framework of infinitesimal homeostasis to study
general design principles for the occurrence of homeostasis in gene regulatory
networks. We assume that the dynamics of the genes explicitly includes both
transcription and translation, keeping track of both mRNA and protein
concentrations. Given a GRN we construct an associated Protein-mRNA Network
(PRN), where each individual (mRNA and protein) concentration corresponds to a
node and the edges are defined in such a way that the PRN becomes a bipartite
directed graph. By simultaneously working with the GRN and the PRN we are able
to apply our previous results about the classification of homeostasis types
(i.e., topologically defined homeostasis generating mechanism) and their
corresponding homeostasis patterns. Given an arbitrarily large and complex GRN
$\mathcal{G}$ and its associated PRN $\mathcal{R}$, we obtain a correspondence
between all the homeostasis types (and homeostasis patterns) of $\mathcal{G}$
and a subset the homeostasis types (and homeostasis patterns) of $\mathcal{R}$.
Moreover, we completely characterize the homeostasis types of the PRN that do
not have GRN counterparts.


  Pharmacodynamic (PD) models are mathematical models of cellular reaction
networks that include drug mechanisms of action. These models are useful for
studying predictive therapeutic outcomes of novel drug therapies in silico.
However, PD models are known to possess significant uncertainty with respect to
constituent parameter data, leading to uncertainty in the model predictions.
Furthermore, experimental data to calibrate these models is often limited or
unavailable for novel pathways. In this study, we present a Bayesian optimal
experimental design approach for improving PD model prediction accuracy. We
then apply our method using simulated experimental data to account for
uncertainty in hypothetical laboratory measurements. This leads to a
probabilistic prediction of drug performance and a quantitative measure of
which prospective laboratory experiment will optimally reduce prediction
uncertainty in the PD model. The methods proposed here provide a way forward
for uncertainty quantification and guided experimental design for models of
novel biological pathways.


  The trust game, derived from an economics experiment, has recently attracted
interest in the field of evolutionary dynamics. In a recent version of the
evolutionary trust game, players adopt one of three strategies: investor,
trustworthy trustee, or untrustworthy trustee. Trustworthy trustees enhance and
share the investment with the investor, whereas untrustworthy trustees retain
the full amount, betraying the investor. Following this setup, we investigate a
two-player trust game, which is analytically feasible under weak selection. We
explore the evolution of trust in structured populations, factoring in four
strategy updating rules: pairwise comparison (PC), birth-death (BD), imitation
(IM), and death-birth (DB). Comparing structured populations with well-mixed
populations, we arrive at two main conclusions. First, in the absence of
untrustworthy trustees, there is a saddle point between investors and
trustworthy trustees, with collaboration thriving best in well-mixed
populations. The collaboration diminishes sequentially from DB to IM to PC/BD
updating rules in structured populations. Second, an invasion of untrustworthy
trustees makes this saddle point unstable and leads to the extinction of
investors. The 3-strategy system stabilizes at an equilibrium line where the
trustworthy and untrustworthy trustees coexist. The stability span of
trustworthy trustees is maximally extended under the PC and BD updating rules
in structured populations, while it decreases in a sequence from IM to DB
updating rules, with the well-mixed population being the least favorable. This
research thus adding an analytical lens to the evolution of trust in structured
populations.


  Models of invasive species spread often assume that landscapes are spatially
homogeneous; thus simplifying analysis but potentially reducing accuracy. We
extend a recently developed partial differential equation model for invasive
conifer spread to account for spatial heterogeneity in parameter values and
introduce a method to obtain key outputs (e.g. spread rates) from computational
simulations. Simulations produce patterns of spatial spread remarkably similar
to observed patterns in grassland ecosystems invaded by exotic conifers,
validating our spatially explicit strategy. We find that incorporating spatial
variation in different parameters does not significantly affect the evolution
of invasions (which are characterised by a long quiescent period followed by
rapid evolution towards to a constant rate of invasion) but that distributional
assumptions can have a significant impact on the spread rate of invasions. Our
work demonstrates that spatial variation in site-suitability or other
parameters can have a significant impact on invasions


  The current theory of evolution is almost the one Darwin and Wallace proposed
two centuries ago and the following discoveries e.g., Mendelian genetics and
neutral mutation theory have not made significant modifications. The current
evolution theory relies mostly on heritable variations within species
population, natural selection and genetic drift. The inability of the current
theory to explain and predict biological observations, especially the emergence
of evolutionary novelties, highlights the need to incorporate recent
evolutionary, developmental and genetics findings in order to achieve a more
comprehensive explanation of species evolution. The present paper provides
significant body of evidence to substantiate a new theory to account for
species evolution. The main axes of the proposed theory include: First,
mutations leading to genetic novelties in a given species during evolution
should be guided by the environment surrounding that species. Second,
environment and germline are connected with each other via soma-germline
messengers e.g., non-coding RNAs (ncRNAs). Third, based on the information that
germline continuously receives in terms of epigenetic messengers three stages
of heritable changes may occur in germline genome to produce more adaptable
offspring: i, Epigenetic modifications, ii, Genetic mutations in the sequence
of pre-existing genes in order to improve their potency, and iii, The
production of new genes with distinct gene-coding regions, which may have their
own regulatory regions.


  Plasma membrane calcium influx through ion channels is crucial for many
events in cellular physiology. Cell surface stimuli lead to the production of
inositol 1,4,5-trisphosphate (IP3), which binds to IP3 receptors in the
endoplasmic reticulum (ER) to release calcium pools from the ER lumen. This
leads to depletion of ER calcium pools which has been termed store-depletion.
Store-depletion leads the dissociation of calcium ions from the EF-hand motif
of the ER calcium sensor Stromal Interaction Molecule 1 (STIM1). This leads to
a conformational change in STIM1 which helps it to interact with a plasma
membrane (PM) at ER:PM junctions. At these ER:PM junctions, STIM1 binds to and
activates a calcium channel known as Orai1 to form calcium-release activated
calcium (CRAC) channels. Activation of Orai1 leads to calcium influx, known as
store-operated calcium entry (SOCE). In addition to Orai1 and STIM1, the
homologs of Orai1 and STIM1, such as Orai2/3 and STIM2 also play a crucial role
in calcium homeostasis. The influx of calcium through the Orai channel
activates a calcium current that has been termed CRAC currents. CRAC channels
form multimers and cluster together in large macromolecular assemblies termed
puncta. How these CRAC channels form puncta has been contentious since their
discovery. In this review, we will outline the history of SOCE, the molecular
players involved in this process (Orai and STIM proteins, TRP channels,
SOCE-associated regulatory factor etc.), as well as the models that have been
proposed to explain this important mechanism in cellular physiology.


  We analyzed the auditory-perceptual space across a substantial portion of the
human vocal range (220-1046 Hz) using multidimensional scaling analysis of
cochlea-scaled spectra from 250-ms vowel segments, initially studied in
Friedrichs et al. (2017) J. Acoust. Soc. Am. 142 1025-1033. The dataset
comprised the vowels /i y e {\o} {\epsilon} a o u/ (N=240) produced by three
native German female speakers, encompassing a broad range of their respective
voice frequency ranges. The initial study demonstrated that, during a
closed-set identification task involving 21 listeners, the point vowels /i a u/
were significantly recognized at fundamental frequencies (fo) nearing 1 kHz,
whereas the recognition of other vowels decreased at higher pitches. Building
on these findings, our study revealed systematic spectral shifts associated
with vowel height and frontness as fo increased, with a notable clustering
around /i a u/ above 523 Hz. These observations underscore the pivotal role of
spectral shape in vowel perception, illustrating the reliance on acoustic
anchors at higher pitches. Furthermore, this study sheds light on the quantal
nature of these vowels and their potential impact on language evolution,
offering a plausible explanation for their widespread presence in the world's
languages.


  Recently, we demonstrated the validity of anomalous numbers law, known as
Newcomb-Benford's law, at the mammalian neuromuscular transmission, considering
different extracellular calcium concentrations arXiv:2002.01986. The present
work continues to examine how alterations in extracellular solution modulate
the first digit law in the context of the spontaneous release of acetylcholine
from the neuromuscular junction. We investigated if the intervals of miniature
potentials collected at the neuromuscular junction obey the law in a
hyperkalemic environment. The analysis showed that the interval between the
miniature potentials at high potassium concentrations follows Newcomb-Benford's
law. Also, our data allowed us to uncover a conformity fluctuation as a
function of the number of intervals of the miniature potentials. Finally, we
discuss the biophysical implications of the present findings.


  Comparing graphs by means of optimal transport has recently gained
significant attention, as the distances induced by optimal transport provide
both a principled metric between graphs as well as an interpretable description
of the associated changes between graphs in terms of a transport plan. As the
lack of symmetry introduces challenges in the typically considered
formulations, optimal transport distances for graphs have mostly been developed
for undirected graphs. Here, we propose two distance measures to compare
directed graphs based on variants of optimal transport: (i) an earth movers
distance (Wasserstein) and (ii) a Gromov-Wasserstein (GW) distance. We evaluate
these two distances and discuss their relative performance for both simulated
graph data and real-world directed cell-cell communication graphs, inferred
from single-cell RNA-seq data.


  The COVID-19 pandemic and its multiple outbreaks have challenged governments
around the world. Much of the epidemiological modeling was based on
pre-pandemic contact information of the population, which changed drastically
due to governmental health measures, so called non-pharmaceutical interventions
made to reduce transmission of the virus, like social distancing and complete
lockdown. In this work, we evaluate an ensemble-based data assimilation
framework applied to a meta-population model to infer the transmission of the
disease between different population agegroups. We perform a set of idealized
twin-experiments to investigate the performance of different possible
parameterizations of the transmission matrix. These experiments show that it is
not possible to unambiguously estimate all the independent parameters of the
transmission matrix. However, under certain parameterizations, the transmission
matrix in an age-compartmental model can be estimated. These estimated
parameters lead to an increase of forecast accuracy in agegroups compartments
assimilating age-dependent accumulated cases and deaths observed in Argentina
compared to a single-compartment model, and reliable estimations of the
effective reproduction number. The age-dependent data assimilation and
forecasting of virus transmission may be important for an accurate prediction
and diagnosis of health care demand.


  The Maasai Mara in Kenya, renowned for its biodiversity, is witnessing
ecosystem degradation and species endangerment due to intensified human
activities. Addressing this, we introduce a dynamic system harmonizing
ecological and human priorities. Our agent-based model replicates the Maasai
Mara savanna ecosystem, incorporating 71 animal species, 10 human
classifications, and 2 natural resource types. The model employs the metabolic
rate-mass relationship for animal energy dynamics, logistic curves for animal
growth, individual interactions for food web simulation, and human intervention
impacts. Algorithms like fitness proportional selection and particle swarm
mimic organism preferences for resources. To guide preservation activities, we
formulated 21 management strategies encompassing tourism, transportation,
taxation, environmental conservation, research, diplomacy, and poaching,
employing a game-theoretic framework. Using the TOPSIS method, we prioritized
four key developmental indicators: environmental health, research advancement,
economic growth, and security. The interplay of 16 factors determines these
indicators, each influenced by our policies to varying degrees. By evaluating
the policies' repercussions, we aim to mitigate adverse animal-human
interactions and equitably address human concerns. We classified the policy
impacts into three categories: Environmental Preservation, Economic Prosperity,
and Holistic Development. By applying these policy groupings to our ecosystem
model, we tracked the effects on the intricate animal-human-resource dynamics.
Utilizing the entropy weight method, we assessed the efficacy of these policy
clusters over a decade, identifying the optimal blend emphasizing both
environmental conservation and economic progression.


  Tens of thousands of simultaneous hypothesis tests are routinely performed in
genomic studies to identify differentially expressed genes. However, due to
unmeasured confounders, many standard statistical approaches may be
substantially biased. This paper investigates the large-scale hypothesis
testing problem for multivariate generalized linear models in the presence of
confounding effects. Under arbitrary confounding mechanisms, we propose a
unified statistical estimation and inference framework that harnesses
orthogonal structures and integrates linear projections into three key stages.
It begins by disentangling marginal and uncorrelated confounding effects to
recover the latent coefficients. Subsequently, latent factors and primary
effects are jointly estimated through lasso-type optimization. Finally, we
incorporate projected and weighted bias-correction steps for hypothesis
testing. Theoretically, we establish the identification conditions of various
effects and non-asymptotic error bounds. We show effective Type-I error control
of asymptotic $z$-tests as sample and response sizes approach infinity.
Numerical experiments demonstrate that the proposed method controls the false
discovery rate by the Benjamini-Hochberg procedure and is more powerful than
alternative methods. By comparing single-cell RNA-seq counts from two groups of
samples, we demonstrate the suitability of adjusting confounding effects when
significant covariates are absent from the model.


  Accurately labeling biomedical data presents a challenge. Traditional
semi-supervised learning methods often under-utilize available unlabeled data.
To address this, we propose a novel reliability-based training data cleaning
method employing inductive conformal prediction (ICP). This method capitalizes
on a small set of accurately labeled training data and leverages ICP-calculated
reliability metrics to rectify mislabeled data and outliers within vast
quantities of noisy training data. The efficacy of the method is validated
across three classification tasks within distinct modalities: filtering
drug-induced-liver-injury (DILI) literature with title and abstract, predicting
ICU admission of COVID-19 patients through CT radiomics and electronic health
records, and subtyping breast cancer using RNA-sequencing data. Varying levels
of noise to the training labels were introduced through label permutation.
Results show significant enhancements in classification performance: accuracy
enhancement in 86 out of 96 DILI experiments (up to 11.4%), AUROC and AUPRC
enhancements in all 48 COVID-19 experiments (up to 23.8% and 69.8%), and
accuracy and macro-average F1 score improvements in 47 out of 48 RNA-sequencing
experiments (up to 74.6% and 89.0%). Our method offers the potential to
substantially boost classification performance in multi-modal biomedical
machine learning tasks. Importantly, it accomplishes this without necessitating
an excessive volume of meticulously curated training data.


  A major challenge in imaging genetics and similar fields is to link
high-dimensional data in one domain, e.g., genetic data, to high dimensional
data in a second domain, e.g., brain imaging data. The standard approach in the
area are mass univariate analyses across genetic factors and imaging
phenotypes. That entails executing one genome-wide association study (GWAS) for
each pre-defined imaging measure. Although this approach has been tremendously
successful, one shortcoming is that phenotypes must be pre-defined.
Consequently, effects that are not confined to pre-selected regions of interest
or that reflect larger brain-wide patterns can easily be missed. In this work
we introduce a Partial Least Squares (PLS)-based framework, which we term
Cluster-Bootstrap PLS (CLUB-PLS), that can work with large input dimensions in
both domains as well as with large sample sizes. One key factor of the
framework is to use cluster bootstrap to provide robust statistics for single
input features in both domains. We applied CLUB-PLS to investigating the
genetic basis of surface area and cortical thickness in a sample of 33,000
subjects from the UK Biobank. We found 107 genome-wide significant
locus-phenotype pairs that are linked to 386 different genes. We found that a
vast majority of these loci could be technically validated at a high rate:
using classic GWAS or Genome-Wide Inferred Statistics (GWIS) we found that 85
locus-phenotype pairs exceeded the genome-wide suggestive (P<1e-05) threshold.


  Stochastic Chemical Reaction Networks are continuous time Markov chain models
that describe the time evolution of the molecular counts of species interacting
stochastically via discrete reactions. Such models are ubiquitous in systems
and synthetic biology, but often have a large or infinite number of states, and
thus are not amenable to computation and analysis. Due to this, approximations
that rely on the molecular counts and the volume being large are commonly used,
with the most common being the Reaction Rate Equations and the Linear Noise
Approximation. For finite time intervals, Kurtz established the validity of the
Reaction Rate Equations and Linear Noise Approximation, by proving law of large
numbers and central limit theorem results respectively. However, the analogous
question for the stationary distribution of the Markov chain model has remained
mostly unanswered, except for chemical reaction networks with special
structures or bounded molecular counts. In this work, we use Stein's Method to
obtain sufficient conditions for the stationary distribution of an
appropriately scaled Stochastic Chemical Reaction Network to converge to the
Linear Noise Approximation as the system size goes to infinity. Our results
give non asymptotic bounds on the error in the Reaction Rate Equations and in
the Linear Noise Approximation as applied to the stationary distribution. As a
special case, we give conditions under which the global exponential stability
of an equilibrium point of the Reaction Rate Equations is sufficient to obtain
our error bounds, thus permitting one to obtain conclusions about the Markov
chain by analyzing the deterministic Reaction Rate Equations.


  African swine fever (ASF) is a highly virulent viral disease that affects
both domestic pigs and wild boar. Current ASF transmission in Europe is in part
driven by wild boar populations, which act as a disease reservoir. Wild boar
are abundant throughout Europe and are highly social animals with complex
social organisation. Despite the known importance of wild boar in ASF spread
and persistence, there remain knowledge gaps surrounding wild boar
transmission. To investigate the influence of density-contact functions and
wild boar social structure on disease dynamics, we developed a wild boar
modelling framework. The framework included an ordinary differential equation
model, a homogeneous stochastic model, and various network-based stochastic
models that explicitly included wild boar social grouping. We found that power
law functions (transmission $\propto$ density$^{0.5}$) and frequency-based
density-contact functions were best able to reproduce recent Baltic outbreaks;
however, power law function models predicted considerable carcass transmission,
while frequency-based models had negligible carcass transmission. Furthermore,
increased model heterogeneity caused a decrease in the relative importance of
carcass-based transmission. The different dominant transmission pathways
predicted by each model type affected the efficacy of potential interventions,
which highlights the importance of evaluating model type and structure when
modelling systems with uncertainties.


  Information transmission via communication between agents is ubiquitous on
Earth, and is a vital facet of living systems. In this paper, we aim to
quantify this rate of information transmission associated with Earth's
biosphere and technosphere (i.e., a measure of global information flow) by
means of a heuristic order-of-magnitude model. By adopting ostensibly
conservative values for the salient parameters, we estimate that the global
information transmission rate for the biosphere might be $\sim 10^{24}$ bits/s,
and that it may perhaps exceed the corresponding rate for the current
technosphere by $\sim 9$ orders of magnitude. However, under the equivocal
assumption of sustained exponential growth, we find that information
transmission in the technosphere can potentially surpass that of the biosphere
$\sim 90$ years in the future, reflecting its increasing dominance.


  Single-cell RNA sequencing (scRNA-seq) enables dissecting cellular
heterogeneity in tissues, resulting in numerous biological discoveries. Various
computational methods have been devised to delineate cell types by clustering
scRNA-seq data where the clusters are often annotated using prior knowledge of
marker genes. In addition to identifying pure cell types, several methods have
been developed to identify cells undergoing state transitions which often rely
on prior clustering results. Present computational approaches predominantly
investigate the local and first-order structures of scRNA-seq data using graph
representations, while scRNA-seq data frequently displays complex
high-dimensional structures. Here, we present a tool, scGeom for exploiting the
multiscale and multidimensional structures in scRNA-seq data by inspecting the
geometry via graph curvature and topology via persistent homology of both cell
networks and gene networks. We demonstrate the utility of these structural
features for reflecting biological properties and functions in several
applications where we show that curvatures and topological signatures of cell
and gene networks can help indicate transition cells and developmental potency
of cells. We additionally illustrate that the structural characteristics can
improve the classification of cell types.


  In this study, we investigate a stochastic SIR model with media effects. The
uniqueness and the existence of a global positive solution are studied. The
sufficient conditions of extinction and persistence of the disease are
established. We obtain the basic reproduction number $R_0^S$ for stochastic
system, which can act as the threshold given small environmental noise. Note
that large noise can induce the disease extinction with probability of 1,
suggesting that environmental noises can not be ignored when investigating
threshold dynamics. Further, inclusion of media induced behaviour changes does
not affect the threshold itself, which is similar to the conclusion of the
deterministic models. However, numerical simulations suggest that media impacts
induce the disease infection decline.


  This paper is the first to propose an allelopathic phytoplankton competition
ODE model influenced by a fear effect based on natural biological phenomena. It
is shown that the interplay of this fear effect and the allelopathic term cause
rich dynamics in the proposed competition model, such as global stability,
transcritical bifurcation, pitchfork bifurcation, and saddle-node bifurcation.
We also consider the spatially explicit version of the model and prove
analogous results. Numerical simulations verify the feasibility of the
theoretical analysis. The results demonstrate that the primary cause of the
extinction of non-toxic species is the fear of toxic species compared to
toxins. Allelopathy only affects the density of non-toxic species. The
discussion provides guidance for the conservation of species and the
maintenance of biodiversity.


  Network biology, an interdisciplinary field at the intersection of
computational and biological sciences, is critical for deepening understanding
of cellular functioning and disease. While the field has existed for about two
decades now, it is still relatively young. There have been rapid changes to it
and new computational challenges have arisen. This is caused by many factors,
including increasing data complexity, such as multiple types of data becoming
available at different levels of biological organization, as well as growing
data size. This means that the research directions in the field need to evolve
as well. Hence, a workshop on Future Directions in Network Biology was
organized and held at the University of Notre Dame in 2022, which brought
together active researchers in various computational and in particular
algorithmic aspects of network biology to identify pressing challenges in this
field. Topics that were discussed during the workshop include: inference and
comparison of biological networks, multimodal data integration and
heterogeneous networks, higher-order network analysis, machine learning on
networks, and network-based personalized medicine. Video recordings of the
workshop presentations are publicly available on YouTube. For even broader
impact of the workshop, this paper, co-authored mostly by the workshop
participants, summarizes the discussion from the workshop. As such, it is
expected to help shape short- and long-term vision for future computational and
algorithmic research in network biology.


  To sample rare events, dissipation-corrected targeted molecular dynamics
(dcTMD) applies a constant velocity constraint along a one-dimensional reaction
coordinate $s$, which drives an atomistic system from an initial state into a
target state. Employing a cumulant approximation of Jarzynski's identity, the
free energy $\Delta G (s)$ is calculated from the mean external work and
dissipated work of the process. By calculating the friction coefficient $\Gamma
(s)$ from the dissipated work, in a second step the equilibrium dynamics of the
process can be studied by propagating a Langevin equation. While so far dcTMD
has been mostly applied to study the unbinding of protein-ligand complexes,
here its applicability to rare conformational transitions within a protein and
the prediction of their kinetics is investigated. As this typically requires
the introduction of multiple collective variables $\{x_j\}= \vec{x}$, a
theoretical framework is outlined to calculate the associated free energy
$\Delta G (\vec{x})$ and friction $\matrix{\Gamma}(\vec{x})$ from dcTMD
simulations along coordinate $s$. Adopting the $\alpha$-$\beta$ transition of
alanine dipeptide as well as the open-closed transition of T4 lysozyme as
representative examples, the virtues and shortcomings of dcTMD to predict
protein conformational transitions and the related kinetics are studied.


  We study the Susceptible-Infected-Recovered (SIR) and the
Susceptible-Exposed-Infected-Recovered (SEIR) models of epidemics, with
possibly time-varying rates, on a class of networks that are locally tree-like,
which includes sparse Erd\H{o}s-R\`enyi random graphs, random regular graphs,
and other configuration models. We identify tractable systems of ODEs that
exactly describe the dynamics of the SIR and SEIR processes in a suitable
asymptotic regime in which the population size goes to infinity. Moreover, in
the case of constant recovery and infection rates, we characterize the outbreak
size as the unique zero of an explicit functional. We use this to show that a
(suitably defined) mean-field prediction always overestimates the outbreak
size, and that the outbreak sizes for SIR and SEIR processes with the same
initial condition and constant infection and recovery rates coincide. In
contrast, we show that the outbreak sizes for SIR and SEIR processes with the
same time-varying infection and recovery rates can in general be quite
different. We also demonstrate via simulations the efficacy of our
approximations for populations of moderate size.


  Free-ranging dogs (FRDs) in human-dominated areas encounter obstacles such as
noise, pollution, limited food sources, and anthropogenic disturbance while
resting. Since FRDs have survived as a population in India, as in many other
parts of the Global South for centuries, they provide a unique opportunity to
study adaptation of animals to the human-dominated urban landscape. We
documented factors impacting resting behaviour and site preferences in three
states of India, for 284 dogs, leading to 6047 observations over 3 years. 7
physical parameters of the resting sites, along with the biological factors
like mating and pup-rearing and time of day affected their choice of resting
sites. The frequency-rank distribution of the unique combinations in which the
parameters were selected followed a Power law distribution, which suggests
underlying biological reasons for the observed preferences. Further, 3 of these
parameters showed maximum consistency of choice in terms of the sub-parameters
selected, explaining 30% of the observations. FRDs prefer to rest close to
their resource sites within the territory, at a place that enabled maximum
visibility of the surroundings. They chose such sites in the core of the
territory for sleeping. At other times, they chose such sites away from the
core, and were less restive, thus allowing for immediate response in case of
intrusion or threat. They generally avoided anthropogenic disturbance for
sleeping, and preferred areas with shade.Incorporating these aspects into urban
management plans can promote human-dog cooperation and reduce situations of
conflict. We envisage more inclusive urban areas in the future, that can allow
for co-existence of the humans and their oldest companions in the commensal
relationship that has been maintained for hundreds of generations of dogs in
this part of the world.


  Previous studies have shown that the interiors of proteins are densely
packed, reaching packing fractions that are as large as those found for static
packings of individual amino-acid-shaped particles. How can the interiors of
proteins take on such high packing fractions given that amino acids are
connected by peptide bonds and many amino acids are hydrophobic with attractive
interactions? We investigate this question by comparing the structural and
mechanical properties of collapsed attractive disk-shaped bead-spring polymers
to those of three reference systems: static packings of repulsive disks, of
attractive disks, and of repulsive disk-shaped bead-spring polymers. We show
that attractive systems quenched to temperatures below the glass transition $T
\ll T_g$ and static packings of both repulsive disks and bead-spring polymers
possess similar interior packing fractions. Previous studies have shown that
static packings of repulsive disks are isostatic at jamming onset, i.e. the
number of contacts $N_c$ matches the number of degrees of freedom, which
strongly influences their mechanical properties. We find that repulsive
polymers are hypostatic at jamming onset, but effectively isostatic when
including quartic modes. While attractive disk and polymer packings are
hyperstatic, we identify a definition for interparticle contacts for which they
can also be considered as effectively isostatic. As a result, we show that the
mechanical properties (e.g. scaling of the potential energy with excess contact
number and low-frequency contribution to the density of vibrational modes) of
weakly attractive disk and polymer packings are similar to those of ${\it
isostatic}$ repulsive disk and polymer packings. Our results demonstrate that
static packings generated via attractive collapse or compression of repulsive
particles possess similar structural and mechanical properties.


  Protein folding is the intricate process by which a linear sequence of amino
acids self-assembles into a unique three-dimensional structure. Protein folding
kinetics is the study of pathways and time-dependent mechanisms a protein
undergoes when it folds. Understanding protein kinetics is essential as a
protein needs to fold correctly for it to perform its biological functions
optimally, and a misfolded protein can sometimes be contorted into shapes that
are not ideal for a cellular environment giving rise to many degenerative,
neuro-degenerative disorders and amyloid diseases. Monitoring at-risk
individuals and detecting protein discrepancies in a protein's folding kinetics
at the early stages could majorly result in public health benefits, as
preventive measures can be taken. This research proposes an efficient pipeline
for predicting protein folding kinetics with high accuracy and low memory
footprint. The deployed machine learning (ML) model outperformed the
state-of-the-art ML models by 4.8% in terms of accuracy while consuming 327x
lesser memory and being 7.3% faster.


  Severe acute respiratory disease SARS-CoV-2 has had a found impact on public
health systems and healthcare emergency response especially with respect to
making decisions on the most effective measures to be taken at any given time.
As demonstrated throughout the last three years with COVID-19, the prediction
of the number of positive cases can be an effective way to facilitate
decision-making. However, the limited availability of data and the highly
dynamic and uncertain nature of the virus transmissibility makes this task very
challenging. Aiming at investigating these challenges and in order to address
this problem, this work studies data-driven (learning, statistical) methods for
incrementally training models to adapt to these nonstationary conditions. An
extensive empirical study is conducted to examine various characteristics, such
as, performance analysis on a per virus wave basis, feature extraction,
"lookback" window size, memory size, all for next-, 7-, and 14-day forecasting
tasks. We demonstrate that the incremental learning framework can successfully
address the aforementioned challenges and perform well during outbreaks,
providing accurate predictions.


  Understanding how cooperation can evolve in populations despite its cost to
individual cooperators is an important challenge. Models of spatially
structured populations with one individual per node of a graph have shown that
cooperation, modeled via the prisoner's dilemma, can be favored by natural
selection. These results depend on microscopic update rules, which determine
how birth, death and migration on the graph are coupled. Recently, we developed
coarse-grained models of spatially structured populations on graphs, where each
node comprises a well-mixed deme, and where migration is independent from
division and death, thus bypassing the need for update rules. Here, we study
the evolution of cooperation in these models in the rare migration regime,
within the prisoner's dilemma. We find that cooperation is not favored by
natural selection in these coarse-grained models on graphs where overall deme
fitness does not directly impact migration from a deme. This is due to a
separation of scales, whereby cooperation occurs at a local level within demes,
while spatial structure matters between demes.


  We study how migration and interactions can rescue species-rich, spatially
structured ecosystems from extinction due to demographic noise. For a single
species, one finds an extinction and a survival state depending on the strength
of the demographic noise with respect to migration. The transition between
these two regimes is a second-order out-of-equilibrium phase transition, which
falls into the universality class of Directed Percolation. Here we show that
the case of many species with heterogeneous interactions is different and
richer. When the demographic fluctuations are strong enough, the transition is
continuous and analogous to the single-species case. Ecological interactions
play a secondary role in it. At small demographic noise, instead, we find novel
features that are a signature of the complexity of the ecosystem. Thanks to
interactions and migration, the metacommunity is able to thrive even when
demographic noise would drive single species to extinction. This is accompanied
by the development of mutualism between non-extinct species. The transition
becomes discontinuous and a global bistability emerges. We present a way to
predict the catastrophic shift from high diversity to extinction by probing
responses to perturbations as early warning signal. Our analytical framework is
grounded in statistical physics and focuses on a model of metacommunities based
on Lotka-Volterra equations with random interactions.


  The COVID-19 pandemic has brought to light a concerning aspect of long-term
neurological complications in post-recovery patients. This study delved into
the investigation of such neurological sequelae in a cohort of 500
post-COVID-19 patients, encompassing individuals with varying illness severity.
The primary aim was to predict outcomes using a machine learning approach based
on diverse clinical data and neuroimaging parameters. The results revealed that
68% of the post-COVID-19 patients reported experiencing neurological symptoms,
with fatigue, headache, and anosmia being the most common manifestations.
Moreover, 22% of the patients exhibited more severe neurological complications,
including encephalopathy and stroke. The application of machine learning models
showed promising results in predicting long-term neurological outcomes.
Notably, the Random Forest model achieved an accuracy of 85%, sensitivity of
80%, and specificity of 90% in identifying patients at risk of developing
neurological sequelae. These findings underscore the importance of continuous
monitoring and follow-up care for post-COVID-19 patients, particularly in
relation to potential neurological complications. The integration of machine
learning-based outcome prediction offers a valuable tool for early intervention
and personalized treatment strategies, aiming to improve patient care and
clinical decision-making. In conclusion, this study sheds light on the
prevalence of long-term neurological complications in post-COVID-19 patients
and demonstrates the potential of machine learning in predicting outcomes,
thereby contributing to enhanced patient management and better health outcomes.
Further research and larger studies are warranted to validate and refine these
predictive models and to gain deeper insights into the underlying mechanisms of
post-COVID-19 neurological sequelae.


  Gene essentiality refers to the degree to which a gene is necessary for the
survival and reproductive efficacy of a living organism. Although the
essentiality of non-coding genes has been documented, there are still aspects
of non-coding genes' essentiality that are unknown to us. For example, We do
not know the contribution of sequence features and network spatial features to
essentiality. As a consequence, in this work, we propose DeepHEN that could
answer the above question. By buidling a new lncRNA-proteion-protein network
and utilizing both representation learning and graph neural network, we
successfully build our DeepHEN models that could predict the essentiality of
lncRNA genes. Compared to other methods for predicting the essentiality of
lncRNA genes, our DeepHEN model not only tells whether sequence features or
network spatial features have a greater influence on essentiality but also
addresses the overfitting issue of those methods caused by the low number of
essential lncRNA genes, as evidenced by the results of enrichment analysis.


  The exploration of pathways and alternative pathways that have a specific
function is of interest in numerous chemical contexts. A framework for
specifying and searching for pathways has previously been developed, but a
focus on which of the many pathway solutions are realisable, or can be made
realisable, is missing. Realisable here means that there actually exists some
sequencing of the reactions of the pathway that will execute the pathway. We
present a method for analysing the realisability of pathways based on the
reachability question in Petri nets. For realisable pathways, our method also
provides a certificate encoding an order of the reactions which realises the
pathway. We present two extended notions of realisability of pathways, one of
which is related to the concept of network catalysts. We exemplify our findings
on the pentose phosphate pathway. Lastly, we discuss the relevance of our
concepts for elucidating the choices often implicitly made when depicting
pathways.


  We explore the connection between migration patterns and emergent behaviors
of evolving populations in spatially heterogeneous environments. Despite
extensive studies in ecologically and medically important systems, a unifying
framework that clarifies this connection and makes concrete predictions remains
much needed. Using a simple evolutionary model on a network of interconnected
habitats with distinct fitness landscapes, we demonstrate a fundamental
connection between migration feedback, emergent ecotypes, and an unusual form
of discontinuous critical transition. We show how migration feedback generates
spatially non-local niches in which emergent ecotypes can specialize. Rugged
fitness landscapes lead to a complex, yet understandable, phase diagram in
which different ecotypes coexist under different migration patterns. The
discontinuous transitions are distinct from the standard first-order phase
transitions in statistical physics. They arise due to simultaneous
transcritical bifurcations and exhibit a "fine structure" due to symmetry
breaking between intra- and inter-ecotype interactions. We suggest feasible
experiments to test our predictions.


  We consider the stochastic Fisher-Kolmogorov-Petrovsky-Piscunov (FKPP)
equation on the circle $\mathbb{S}$, \begin{equation*}
  \partial_t u(t,x) \,= \frac{\alpha}{2}\Delta u +\beta\,u(1-u) +
\sqrt{\gamma\,u(1-u)}\,\dot{W}, \qquad (t,x)\in(0,\infty)\times \mathbb{S},
\end{equation*} where $\dot{W}$ is space-time white noise. While any solution
will eventually be absorbed at one of two states, the constant 1 and the
constant 0 on the circle, essentially nothing had been established about the
absorption time (also called the fixation time in population genetics), or
about the long-time behavior prior to absorption. We establish the existence
and uniqueness of the quasi-stationary distribution (QSD) for the solution of
the stochastic FKPP. Moreover, we show that the solution conditioned on not
being absorbed at time $t$ converges to this unique QSD as $t\to\infty$, for
any initial distribution, and characterize the leading-order asymptotics for
the tail distribution of the fixation time. We obtain explicit calculations in
the neutral case ($\beta=0$), quantifying the effect of spatial diffusion on
fixation time. We explicitly express the fixation rate in terms of the
migration rate $\alpha$ for all $\alpha\in (0,\infty)$, finding in particular
that the fixation rate is given by
$\gamma[1-\frac{\gamma}{12\alpha}+\mathcal{O}(\frac{\gamma^2}{\alpha^2})]$ for
fast migration and
$\pi^2\alpha[1-\frac{8\alpha}{\gamma}+\mathcal{O}(\frac{\alpha^2}{\gamma^2})]$
for slow migration. Our proof relies on the observation that the absorbed (or
killed) stochastic FKPP is dual to a system of $2$-type branching-coalescing
Brownian motions killed when one type dies off, and on leveraging the
relationship between these two killed processes.


  DNA sequence alignment involves assigning short DNA reads to the most
probable locations on an extensive reference genome. This process is crucial
for various genomic analyses, including variant calling, transcriptomics, and
epigenomics. Conventional methods, refined over decades, tackle this challenge
in two steps: genome indexing followed by efficient search to locate likely
positions for given reads. Building on the success of Large Language Models
(LLM) in encoding text into embeddings, where the distance metric captures
semantic similarity, recent efforts have explored whether the same Transformer
architecture can produce numerical representations for DNA sequences. Such
models have shown early promise in tasks involving classification of short DNA
sequences, such as the detection of coding vs non-coding regions, as well as
the identification of enhancer and promoter sequences. Performance at sequence
classification tasks does not, however, translate to sequence alignment, where
it is necessary to conduct a genome-wide search to successfully align every
read. We address this open problem by framing it as an Embed-Search-Align task.
In this framework, a novel encoder model DNA-ESA generates representations of
reads and fragments of the reference, which are projected into a shared vector
space where the read-fragment distance is used as surrogate for alignment. In
particular, DNA-ESA introduces: (1) Contrastive loss for self-supervised
training of DNA sequence representations, facilitating rich sequence-level
embeddings, and (2) a DNA vector store to enable search across fragments on a
global scale. DNA-ESA is >97% accurate when aligning 250-length reads onto a
human reference genome of 3 gigabases (single-haploid), far exceeds the
performance of 6 recent DNA-Transformer model baselines and shows task transfer
across chromosomes and species.


  Bacteria reside in constantly changing environments and require rapid and
precise adjustments of gene expression to ensure survival. Small regulatory
RNAs (sRNAs) are a crucial element that bacteria utilize to achieve this. sRNAs
are short RNA molecules that modulate gene expression usually through
base-pairing interactions with target RNAs, primarily mRNAs. These interactions
can lead to either negative outcomes such as mRNA degradation or translational
repression or positive outcomes such as mRNA stabilization or translation
enhancement. In recent years, high-throughput approaches such as RIL-seq (RNA
interaction by ligation and sequencing) revolutionized the sRNA field by
enabling the identification of sRNA targets on a global scale, unveiling
intricate sRNA-RNA networks. In this review, we discuss the insights gained
from investigating sRNA-RNA networks in well-studied bacterial species as well
as in under-studied bacterial species. Having a complete understanding of
sRNA-mediated regulation is critical for the development of new strategies for
controlling bacterial growth and combating bacterial infections.


  Relating forest productivity to local variations in forest structure has been
a long-standing challenge. Previous studies often focused on the connection
between forest structure and stand-level photosynthesis (GPP). However, biomass
production (NPP) and net ecosystem exchange (NEE) are also subject to
respiration and other carbon losses, which vary with local conditions and life
history traits. Here, we use a simulation approach to study how these losses
impact forest productivity and reveal themselves in forest structure. We fit
the process-based forest model Formind to a 25ha inventory of an old-growth
temperate forest in China and classify trees as "mature" (full-grown) or
"immature" based on their intrinsic carbon use efficiency. Our results reveal a
strong negative connection between the stand-level carbon use efficiency and
the prevalence of mature trees: GPP increases with the total basal area,
whereas NPP and NEE are driven by the basal area of immature trees.
Accordingly, the basal area entropy - a structural proxy for the prevalence of
immature trees - correlated well with NPP and NEE and had higher predictive
power than other structural characteristics such as Shannon diversity and
height standard deviation. Our results were robust across spatial scales
(0.04-1ha) and yield promising hypotheses field studies and new theoretical
work.


  Virtual screening of large compound libraries to identify potential hit
candidates is one of the earliest steps in drug discovery. As the size of
commercially available compound collections grows exponentially to the scale of
billions, brute-force virtual screening using traditional tools such as docking
becomes infeasible in terms of time and computational resources. Active
learning and Bayesian optimization has recently been proven as effective
methods of narrowing down the search space. An essential component in those
methods is a surrogate machine learning model that is trained with a small
subset of the library to predict the desired properties of compounds. Accurate
model can achieve high sample efficiency by finding the most promising
compounds with only a fraction of the whole library being virtually screened.
In this study, we examined the performance of pretrained transformer-based
language model and graph neural network in Bayesian optimization active
learning framework. The best pretrained models identifies 58.97% of the
top-50000 by docking score after screening only 0.6% of an ultra-large library
containing 99.5 million compounds, improving 8% over previous state-of-the-art
baseline. Through extensive benchmarks, we show that the superior performance
of pretrained models persists in both structure-based and ligand-based drug
discovery. Such model can serve as a boost to the accuracy and sample
efficiency of active learning based molecule virtual screening.


  Motivation: Navigating the high dimensional space of discrete trees for
phylogenetics presents a challenging problem for tree optimisation. To address
this, hyperbolic embeddings of trees offer a promising approach to encoding
trees efficiently in continuous spaces. However, they require a differentiable
tree decoder to optimise the phylogenetic likelihood. We present soft-NJ, a
differentiable version of neighbour-joining that enables gradient-based
optimisation over the space of trees.
  Results: We illustrate the potential for differentiable optimisation over
tree space for maximum likelihood inference. We then perform variational
Bayesian phylogenetics by optimising embedding distributions in hyperbolic
space. We compare the performance of this approximation technique on eight
benchmark datasets to state-of-art methods. However, geometric frustrations of
the embedding locations produce local optima that pose a challenge for
optimisation.
  Availability: Dodonaphy is freely available on the web at
www.https://github.com/mattapow/dodonaphy. It includes an implementation of
soft-NJ.


  One-dependent first passage percolation is a spreading process on a graph
where the transmission time through each edge depends on the direct
surroundings of the edge. In particular, the classical iid transmission time
$L_{xy}$ is multiplied by $(W_xW_y)^\mu$, a polynomial of the expected degrees
$W_x, W_y$ of the endpoints of the edge $xy$, which we call the penalty
function. Beyond the Markov case, we also allow any distribution for $L_{xy}$
with regularly varying distribution near $0$. We then run this process on three
spatial scale-free random graph models: finite and infinite Geometric
Inhomogeneous Random Graphs, and Scale-Free Percolation. In these spatial
models, the connection probability between two vertices depends on their
spatial distance and on their expected degrees.
  We show that as the penalty-function, i.e., $\mu$ increases, the transmission
time between two far away vertices sweeps through four universal phases:
explosive (with tight transmission times), polylogarithmic, polynomial but
strictly sublinear, and linear in the Euclidean distance. The strictly
polynomial growth phase here is a new phenomenon that so far was extremely rare
in spatial graph models. The four growth phases are highly robust in the model
parameters and are not restricted to phase boundaries. Further, the transition
points between the phases depend non-trivially on the main model parameters:
the tail of the degree distribution, a long-range parameter governing the
presence of long edges, and the behaviour of the distribution $L$ near $0$. In
this paper we develop new methods to prove the upper bounds in all
sub-explosive phases. Our companion paper complements these results by
providing matching lower bounds in the polynomial and linear regimes.


  In this paper we study a version of (non-Markovian) first passage percolation
on graphs, where the transmission time between two connected vertices is
non-iid, but increases by a penalty factor polynomial in their expected
degrees. Based on the exponent of the penalty-polynomial, this makes it
increasingly harder to transmit to and from high-degree vertices. This choice
is motivated by awareness or time-limitations. For the iid part of the
transmission times we allow any nonnegative distribution with regularly varying
behaviour at $0$. For the underlying graph models we choose spatial random
graphs that have power-law degree distributions, so that the effect of the
penalisation becomes visible: (finite and infinite) Geometric Inhomogeneous
Random Graphs, and Scale-Free Percolation. In these spatial models, the
connection probability between two vertices depends on their spatial distance
and on their expected degrees. We prove that upon increasing the penalty
exponent, the transmission time between two far away vertices $x,y$ sweeps
through four universal phases even for a single underlying graph: explosive
(tight transmission times), polylogarithmic, polynomial but sublinear
($\Theta(|x-y|^{\eta_0})$ for an explicit $\eta_0<1$), and linear
($\Theta(|x-y|)$) in their Euclidean distance. Further, none of these phases
are restricted to phase boundaries, and those are non-trivial in the main model
parameters: the tail of the degree-distribution, a long-range parameter, and
the exponent of regular variation of the iid part of the transmission times. In
this paper we present proofs of lower bounds for the latter two phases and the
upper bound for the linear phase. These complement the matching upper bounds
for the polynomial regime in our companion paper.


  We aim at understanding the evolution of the genetic composition of cancer
cell populations. To this aim, we consider a branching individual based model
representing a cell population where cells divide, die and mutate along the
edges of a finite directed graph $(V,E)$. The process starts with only one cell
of trait $0$. Following typical parameter values in cancer cell populations we
study the model under large population and power law mutation rates limit, in
the sense that the mutation probabilities are parameterized by negative powers
of $n$ and the typical sizes of the population of our interest are positive
powers of $n$. Under non-increasing growth rate condition (namely the growth
rate of any sub-population is smaller than the growth rate of trait $0$), we
describe the time evolution of the first-order asymptotics of the size of each
sub-population on the $log(n)$ time scale, as well as in the random time scale
at which the initial population, resp. the total population, reaches the size
$n^{t}$. In particular, such results allow to characterize whose mutational
paths along the edges of the graph are actually contributing to the size order
of the sub-populations. Without any condition on the growth rate, we describe
the time evolution of the orders of magnitude of each sub-population. Adapting
techniques from Durrett and Mayberry 2011, we show that these converges to
positive deterministic non-decreasing piecewise linear continuous functions,
whose slopes are given by an algorithm.


  Steered molecular dynamics (SMD) simulations for the calculation of free
energies are well suited for high-throughput molecular simulations on a
distributed infrastructure due to the simplicity of the setup and parallel
granularity of the runs. However, so far, the computational cost limited the
estimation of the free energy typically over just a few pullings, thus impeding
the evaluation of statistical uncertainties involved. In this work, we
performed two thousand pulls for the permeation of a potassium ion in the
gramicidin A pore by all-atom molecular dynamics in order to assess the
bidirectional SMD protocol with a proper amount of sampling. The estimated free
energy profile still shows a statistical error of several kcal/mol, while the
work distributions are estimated to be non-Gaussian at pulling speeds of 10
{\AA}/ns. We discuss the methodology and the confidence intervals in relation
to increasing amounts of computed trajectories and how different permeation
pathways for the potassium ion, knock-on and sideways, affect the sampling and
the free energy estimation.


  Global challenges like climate change may be considered as collective action
problems that require sufficient cooperation with pro-mitigation norms, soon
enough to be effective. Socio-political polarization is a barrier to collective
action. Prior agent-based models of behavioural change on structured networks
in a shared socio-political environment have shown that polarization emerges
naturally in such systems and that the speed of consensus formation is limited
by the rate at which polarized clusters can be dissolved. Here we study how
guided social link rewiring affects the speed of network depolarization. We
investigate rewiring algorithms representing random meetings, introduction by
mutual acquaintances, and bridging between socially distant communities. We
find that building lasting links between polarized individuals and communities
can accelerate consensus formation when the sociopolitical environment is
favourable. This strengthens the evidence that promoting connection between
polarized communities could accelerate collective action on urgent global
challenges.


  We analyze the vegetation growth dynamics with a stochastic cellular automata
model and in real-world data obtained from satellite images. We look for areas
where vegetation breaks down into clusters, comparing it to a percolation
transition that happens in the cellular automata model and is an early warning
signal of land degradation. We use satellite imagery data such as the
Normalized Difference Vegetation Index (NDVI) and Leaf Area Index (LAI). We
consider the periodic effect of seasons as a periodic environmental stress, and
show numerically how the vegetation can be resilient to high stress during
seasonal fluctuations. We qualitatively recognize these effects in real-world
vegetation images. Finally, we qualitatively evaluate the environmental stress
in land images by considering both the vegetation density and its
clusterization.


  This paper presents a mathematical model to investigate co-infection with
HIV/AIDS and zika virus (ZIKV) in Colombia and Brazil, where the first cases
were reported in 2015-2016. The model considers the sexual transmission
dynamics of both viruses and vector-host interactions. We begin by exploring
the qualitative behaviour of each model separately. Then, we analyze the
dynamics of the co-infection model using the thresholds and results defined
separately for each model. The model also considers the impact of intervention
strategies, such as, personal protection, antiretroviral therapy (ART), and
sexual protection (condoms use). Using available parameter values for Colombia
and Brazil, the model is calibrated to predict the potential effect of
implementing combinations of those intervention strategies on the co-infection
spread. According to these findings, transmission through sexual contact is a
determining factor in the long-term behaviour of these two diseases.
Furthermore, it is important to note that co-infection with HIV and ZIKV may
result in higher rates of HIV transmission and an increased risk of severe
congenital disabilities linked to ZIKV infection. As a result, control measures
have been implemented to limit the number of infected individuals and
mosquitoes, with the aim of halting disease transmission. This study provides
novel insights into the dynamics of HIV/ZIKV co-infection and highlights the
importance of integrated intervention strategies in controlling the spread of
these viruses, which may impact public health


  For over 165 million years, dinosaurs reigned on this planet. Their entire
existence saw variations in their body size and mass . Understanding the
relationship between various attributes such as femur length, breadth; humerus
length, breadth; tibia length, breadth and body mass of dinosaurs contributes
to our understanding of the Jurassic era and further provides reasoning for
bone and body size evolution of modern day descendants of those from the
Dinosauria clade. The following work consists of statistical evidence derived
from an encyclopedic data set consisting of a wide variety of measurements
pertaining to discovered fossils of a particular taxa of dinosaur. Our study
establishes linearly regressive correspondence between femur and humerus length
and radii. Furthermore, there is also a comparison with terrestrial bird bone
lengths, to verify the claim of birds being closest alive species to dinosaurs.
An analysis into bone ratios of early humans shows that terrestrial birds are
closer to humans than that of dinosaurs. Not only on one hand it challenges the
closeness of birds with dinosaurs, but on the other hand it makes a case of
convergent evolution between birds and humans, due to their closeness in
regressive fits.
  A correlation between bone ratios of dinosaurs and early humans also advances
understanding in the structural and physical distinctions between the two
species. Overall, the work contains evaluation of dinosaur skeletons and
promotes further exploration and research in the paleontological field to
strengthen the conclusions drawn thus far.


  Recently there has been an increasing interest in studying dynamical
processes on networks exhibiting higher-order structures, such as simplicial
complexes, where the dynamics acts above and beyond dyadic interactions. Using
simulations or heuristically derived epidemic spreading models it was shown
that new phenomena can emerge, such as bi-stability/multistability. Here, we
show that such new emerging phenomena do not require complex contact patterns,
such as community structures, but naturally result from the higher-order
contagion mechanisms. We show this by deriving an exact higher-order SIS model
and its limiting mean-field equivalent for fully connected simplicial
complexes. Going beyond previous results, we also give the global bifurcation
picture for networks with 3- and 4-body interactions, with the latter allowing
for two non-trivial stable endemic steady states. Differently from previous
approaches, we are able to study systems featuring interactions of arbitrary
order. In addition, we characterise the contributions from higher-order
infections to the endemic equilibrium as perturbations of the pairwise
baseline, finding that these diminish as the pairwise rate of infection
increases. Our approach represents a first step towards a principled
understanding of higher-order contagion processes beyond triads and opens up
further directions for analytical investigations.


  Explainable Graph Neural Networks (GNNs) have been developed and applied to
drug-protein binding prediction to identify the key chemical structures in a
drug that have active interactions with the target proteins. However, the key
structures identified by the current explainable GNN models are typically
chemically invalid. Furthermore, a threshold needs to be manually selected to
pinpoint the key structures from the rest. To overcome the limitations of the
current explainable GNN models, we propose our SLGNN, which stands for using
Sparse Learning to Graph Neural Networks. Our SLGNN relies on using a
chemical-substructure-based graph (where nodes are chemical substructures) to
represent a drug molecule. Furthermore, SLGNN incorporates generalized fussed
lasso with message-passing algorithms to identify connected subgraphs that are
critical for the drug-protein binding prediction. Due to the use of the
chemical-substructure-based graph, it is guaranteed that any subgraphs in a
drug identified by our SLGNN are chemically valid structures. These structures
can be further interpreted as the key chemical structures for the drug to bind
to the target protein. We demonstrate the explanatory power of our SLGNN by
first showing all the key structures identified by our SLGNN are chemically
valid. In addition, we illustrate that the key structures identified by our
SLGNN have more predictive power than the key structures identified by the
competing methods. At last, we use known drug-protein binding data to show the
key structures identified by our SLGNN contain most of the binding sites.


  Drug repurposing is an emerging approach for drug discovery involving the
reassignment of existing drugs for novel purposes. An alternative to the
traditional de novo process of drug development, repurposed drugs are faster,
cheaper, and less failure prone than drugs developed from traditional methods.
Recently, drug repurposing has been performed in silico, in which databases of
drugs and chemical information are used to determine interactions between
target proteins and drug molecules to identify potential drug candidates. A
proposed algorithm is NeuroCADR, a novel system for drug repurposing via a
multi-pronged approach consisting of k-nearest neighbor algorithms (KNN),
random forest classification, and decision trees. Data was sourced from several
databases consisting of interactions between diseases, symptoms, genes, and
affiliated drug molecules, which were then compiled into datasets expressed in
binary. The proposed method displayed a high level of accuracy, outperforming
nearly all in silico approaches. NeuroCADR was performed on epilepsy, a
condition characterized by seizures, periods of time with bursts of
uncontrolled electrical activity in brain cells. Existing drugs for epilepsy
can be ineffective and expensive, revealing a need for new antiepileptic drugs.
NeuroCADR identified novel drug candidates for epilepsy that can be further
approved through clinical trials. The algorithm has the potential to determine
possible drug combinations to prescribe a patient based on a patient's prior
medical history. This project examines NeuroCADR, a novel approach to
computational drug repurposing capable of revealing potential drug candidates
in neurological diseases such as epilepsy.


  Over the past decade, antibodies have steadily grown in therapeutic
importance thanks to their high specificity and low risk of adverse effects
compared to other drug modalities. While traditional antibody discovery is
primarily wet lab driven, the rapid improvement of ML-based generative modeling
has made in-silico approaches an increasingly viable route for discovery and
engineering. To this end, we train an antibody-specific language model,
AntiBARTy, based on BART (Bidirectional and Auto-Regressive Transformer) and
use its latent space to train a property-conditional diffusion model for guided
IgG de novo design. As a test case, we show that we can effectively generate
novel antibodies with improved in-silico solubility while maintaining antibody
validity and controlling sequence diversity.


  During the SARS-CoV-2 pandemic, wastewater-based genomic surveillance (WWGS)
emerged as an efficient viral surveillance tool that takes into account
asymptomatic cases and can identify known and novel mutations and offers the
opportunity to assign known virus lineages based on the detected mutations
profiles. WWGS can also hint towards novel or cryptic lineages, but it is
difficult to clearly identify and define novel lineages from wastewater (WW)
alone. While WWGS has significant advantages in monitoring SARS-CoV-2 viral
spread, technical challenges remain, including poor sequencing coverage and
quality due to viral RNA degradation. As a result, the viral RNAs in wastewater
have low concentrations and are often fragmented, making sequencing difficult.
WWGS analysis requires advanced computational tools that are yet to be
developed and benchmarked. The existing bioinformatics tools used to analyze
wastewater sequencing data are often based on previously developed methods for
quantifying the expression of transcripts or viral diversity. Those methods
were not developed for wastewater sequencing data specifically, and are not
optimized to address unique challenges associated with wastewater. While
specialized tools for analysis of wastewater sequencing data have also been
developed recently, it remains to be seen how they will perform given the
ongoing evolution of SARS-CoV-2 and the decline in testing and patient-based
genomic surveillance. Here, we discuss opportunities and challenges associated
with WWGS, including sample preparation, sequencing technology, and
bioinformatics methods.


  In 2022, a group of basic and clinical virologists, bioinformaticians, and
evolutionary and structural biologists met in Oxford, UK, to develop a
consensus on methodologies used to classify viruses. They concluded that virus
taxonomy, which is hierarchical and based on evolution, is only one of many
possible ways to classify viruses. This taxonomy, while satisfying the four
principles they set out, faces difficulties in coordinating with other
classification systems useful to clinicians, infectious disease specialists,
agronomists, etc. One example discussed is the grouping of different viral
strains that cause different diseases into the species Enterovirus C. Here we
show that the use of a previously proposed variant of a natural virus
classification system based on the use of Neural Replicator Analysis can
resolve this contradiction by establishing the fine structure of the
Enterovirus C species, in which strains that cause different diseases are
placed in several different cells of the binomial table of viruses. A key
element in enabling this is the sophisticated preprocessing of the original
viral genomes using neural replicators.


  Protein-protein interactions (PPIs) are critical for various biological
processes, and understanding their dynamics is essential for decoding molecular
mechanisms and advancing fields such as cancer research and drug discovery.
Mutations in PPIs can disrupt protein binding affinity and lead to functional
changes and disease. Predicting the impact of mutations on binding affinity is
valuable but experimentally challenging. Computational methods, including
physics-based and machine learning-based approaches, have been developed to
address this challenge. Machine learning-based methods, fueled by extensive PPI
datasets such as Ab-Bind, PINT, SKEMPI, and others, have shown promise in
predicting binding affinity changes. However, accurate predictions and
generalization of these models across different datasets remain challenging.
Geometric graph learning has emerged as a powerful approach, combining graph
theory and machine learning, to capture structural features of biomolecules. We
present GGL-PPI, a novel method that integrates geometric graph learning and
machine learning to predict mutation-induced binding free energy changes.
GGL-PPI leverages atom-level graph coloring and multi-scale weighted colored
geometric subgraphs to extract informative features, demonstrating superior
performance on three validation datasets, namely AB-Bind, SKEMPI 1.0, and
SKEMPI 2.0 datasets. Evaluation on a blind test set highlights the unbiased
predictions of GGL-PPI for both direct and reverse mutations. The findings
underscore the potential of GGL-PPI in accurately predicting binding free
energy changes, contributing to our understanding of PPIs and aiding drug
design efforts.


  In single-cell RNA sequencing (scRNA-seq) data analysis, a critical challenge
is to infer hidden dynamic cellular processes from measured static cell
snapshots. To tackle this challenge, many computational methods have been
developed from distinct perspectives. Besides the common perspectives of
inferring trajectories (or pseudotime) and RNA velocity, another important
perspective is to estimate the differentiation potential of cells, which is
commonly referred to as "cell potency." In this review, we provide a
comprehensive summary of 11 computational methods that estimate cell potency
from scRNA-seq data under different assumptions, some of which are even
conceptually contradictory. We divide these methods into three categories:
mean-based, entropy-based, and correlation-based methods, depending on how a
method summarizes gene expression levels of a cell or cell type into a potency
measure. Our review focuses on the key similarities and differences of the
methods within each category and between the categories, providing a high-level
intuition of each method. Moreover, we use a unified set of mathematical
notations to detail the 11 methods' methodologies and summarize their usage
complexities, including the number of ad-hoc parameters, the number of required
inputs, and the existence of discrepancies between the method description in
publications and the method implementation in software packages. Realizing the
conceptual contradictions of existing methods and the difficulty of fair
benchmarking without single-cell-level ground truths, we conclude that accurate
estimation of cell potency from scRNA-seq data remains an open challenge.


  We investigate the connections between clusters and least common ancestors
(LCAs) in directed acyclic graphs (DAGs). We focus on the class of DAGs having
unique least common ancestors for certain subsets of their minimal elements
since these are of interest, particularly as models of phylogenetic networks.
Here, we use the close connection between the canonical k-ary transit function
and the closure function on a set system to show that pre-k-ary clustering
systems are exactly those that derive from a class of DAGs with unique LCAs.
Moreover, we show that k-ary T-systems and k-weak hierarchies are associated
with DAGs that satisfy stronger conditions on the existence of unique LCAs for
sets of size at most k.


  Quantifying variable importance is essential for answering high-stakes
questions in fields like genetics, public policy, and medicine. Current methods
generally calculate variable importance for a given model trained on a given
dataset. However, for a given dataset, there may be many models that explain
the target outcome equally well; without accounting for all possible
explanations, different researchers may arrive at many conflicting yet equally
valid conclusions given the same data. Additionally, even when accounting for
all possible explanations for a given dataset, these insights may not
generalize because not all good explanations are stable across reasonable data
perturbations. We propose a new variable importance framework that quantifies
the importance of a variable across the set of all good models and is stable
across the data distribution. Our framework is extremely flexible and can be
integrated with most existing model classes and global variable importance
metrics. We demonstrate through experiments that our framework recovers
variable importance rankings for complex simulation setups where other methods
fail. Further, we show that our framework accurately estimates the true
importance of a variable for the underlying data distribution. We provide
theoretical guarantees on the consistency and finite sample error rates for our
estimator. Finally, we demonstrate its utility with a real-world case study
exploring which genes are important for predicting HIV load in persons with
HIV, highlighting an important gene that has not previously been studied in
connection with HIV. Code is available at
https://github.com/jdonnelly36/Rashomon_Importance_Distribution.


  Principal components computed via PCA (principal component analysis) are
traditionally used to reduce dimensionality in genomic data or to correct for
population stratification. In this paper, we explore the penalized eigenvalue
problem (PEP) which reformulates the computation of the first eigenvector as an
optimization problem and adds an L1 penalty constraint. The contribution of our
article is threefold. First, we extend PEP by applying Nesterov smoothing to
the original LASSO-type L1 penalty. This allows one to compute analytical
gradients which enable faster and more efficient minimization of the objective
function associated with the optimization problem. Second, we demonstrate how
higher order eigenvectors can be calculated with PEP using established results
from singular value decomposition (SVD). Third, using data from the 1000 Genome
Project dataset, we empirically demonstrate that our proposed smoothed PEP
allows one to increase numerical stability and obtain meaningful eigenvectors.
We further investigate the utility of the penalized eigenvector approach over
traditional PCA.


  Generative molecular design has moved from proof-of-concept to real-world
applicability, as marked by the surge in very recent papers reporting
experimental validation. Key challenges in explainability and sample efficiency
present opportunities to enhance generative design to directly optimize
expensive high-fidelity oracles and provide actionable insights to domain
experts. Here, we propose Beam Enumeration to exhaustively enumerate the most
probable sub-sequences from language-based molecular generative models and show
that molecular substructures can be extracted. When coupled with reinforcement
learning, extracted substructures become meaningful, providing a source of
explainability and improving sample efficiency through self-conditioned
generation. Beam Enumeration is generally applicable to any language-based
molecular generative model and notably further improves the performance of the
recently reported Augmented Memory algorithm, which achieved the new
state-of-the-art on the Practical Molecular Optimization benchmark for sample
efficiency. The combined algorithm generates more high reward molecules and
faster, given a fixed oracle budget. Beam Enumeration is the first method to
jointly address explainability and sample efficiency for molecular design.