  We propose and analyze an epidemiological model for vector borne diseases
that integrates a multi-stage vector population and several host
sub-populations which may be characterized by a variety of compartmental model
types: subpopulations all include Susceptible and Infected compartments, but
may or may not include Exposed and/or Recovered compartments. The model was
originally designed to evaluate the effectiveness of various prophylactic
measures in malaria-endemic areas, but can be applied as well to other
vector-borne diseases. This model is expressed as a system of several
differential equations, where the number of equations depends on the particular
assumptions of the model. We compute the basic reproduction number $\mathcal
R_0$, and show that if $\mathcal R_0\leqslant 1$, the disease free equilibrium
(DFE) is globally asymptotically stable (GAS) on the nonnegative orthant. If
$\mathcal R_0>1$, the system admits a unique endemic equilibrium (EE) that is
GAS. We analyze the sensitivity of $R_0$ and the EE to different system
parameters, and based on this analysis we discuss the relative effectiveness of
different control measures.


  Decoding the linguistic intricacies of the genome is a crucial problem in
biology, and pre-trained foundational models such as DNABERT and Nucleotide
Transformer have made significant strides in this area. Existing works have
largely hinged on k-mer, fixed-length permutations of A, T, C, and G, as the
token of the genome language due to its simplicity. However, we argue that the
computation and sample inefficiencies introduced by k-mer tokenization are
primary obstacles in developing large genome foundational models. We provide
conceptual and empirical insights into genome tokenization, building on which
we propose to replace k-mer tokenization with Byte Pair Encoding (BPE), a
statistics-based data compression algorithm that constructs tokens by
iteratively merging the most frequent co-occurring genome segment in the
corpus. We demonstrate that BPE not only overcomes the limitations of k-mer
tokenization but also benefits from the computational efficiency of
non-overlapping tokenization. Based on these insights, we introduce DNABERT-2,
a refined genome foundation model that adapts an efficient tokenizer and
employs multiple strategies to overcome input length constraints, reduce time
and memory expenditure, and enhance model capability. Furthermore, we identify
the absence of a comprehensive and standardized benchmark for genome
understanding as another significant impediment to fair comparative analysis.
In response, we propose the Genome Understanding Evaluation (GUE), a
comprehensive multi-species genome classification dataset that amalgamates $28$
distinct datasets across $7$ tasks, with input lengths ranging from $70$ to
$1000$. Through comprehensive experiments on the GUE benchmark, we demonstrate
that DNABERT-2 achieves comparable performance to the state-of-the-art model
with $21 \times$ fewer parameters and approximately $56 \times$ less GPU time
in pre-training.


  We introduce a two-strain model with asymmetric temporary immunity periods
and partial cross-immunity. We derive explicit conditions for competitive
exclusion and coexistence of the strains depending on the strain-specific basic
reproduction numbers, temporary immunity periods, and degree of cross-immunity.
The results of our bifurcation analysis suggest that, even when two strains
share similar basic reproduction numbers and other epidemiological parameters,
a disparity in temporary immunity periods and partial or complete
cross-immunity can provide a significant competitive advantage. To analyze the
dynamics, we introduce a quasi-steady state reduced model which assumes the
original strain remains at its endemic steady state. We completely analyze the
resulting reduced planar hybrid switching system using linear stability
analysis, planar phase-plane analysis, and the Bendixson-Dulac criterion. We
validate both the full and reduced models with COVID-19 incidence data,
focusing on the Delta (B.1.617.2), Omicron (B.1.1.529), and Kraken (XBB.1.5)
variants. These numerical studies suggest that, while early novel strains of
COVID-19 had a tendency toward dramatic takeovers and extinction of ancestral
strains, more recent strains have the capacity for co-existence.


  Optically transparent wound dressings offer a range of potential applications
in the biomedical field, as they allow for the monitoring of wound healing
progress without having to replace the dressing. These dressings must be
impermeable to water and bacteria, yet permeable to moisture vapor and
atmospheric gases in order to maintain a moist environment at the wound site.
This review article provides a comprehensive overview of the types of wound
dressings, novel wound dressing materials, advanced fabrication techniques for
transparent wound dressing materials, and the key features and applications of
transparent dressings for the healing process, as well as how it can improve
healing outcomes. This review mainly focuses on representing specifications of
transparent polymeric wound dressing materials, such as transparent electrospun
nanofibers, transparent crosslinked hydrogels, and transparent composite
films/membranes. Due to the advance properties of electrospun nanofiber such as
large surface area, enable efficient incorporation of antibacterial molecules,
a structure similar to the extracellular matrix, and high mechanical stability,
is often used in wound dressing applications. We also highlight the hydrogels
or films for wound healing applications, it's promote the healing process,
provide a moisture environment, and offer pain relief with their cool,
high-water content, excellent biocompatibility, and bio-biodegradability.


  We consider a moving boundary mathematical model of biological invasion. The
model describes the spatiotemporal evolution of two adjacent populations: each
population undergoes linear diffusion and logistic growth, and the boundary
between the two populations evolves according to a two--phase Stefan condition.
This mathematical model describes situations where one population invades into
regions occupied by the other population, such as the spreading of a malignant
tumour into surrounding tissues. Full time--dependent numerical solutions are
obtained using a level--set numerical method. We use these numerical solutions
to explore several properties of the model including: (i) survival and
extinction of one population initially surrounded by the other; and (ii) linear
stability of the moving front boundary in the context of a travelling wave
solution subjected to transverse perturbations. Overall, we show that many
features of the well--studied one--phase single population analogue of this
model can be very different in the more realistic two--phase setting. These
results are important because realistic examples of biological invasion involve
interactions between multiple populations and so great care should be taken
when extrapolating predictions from a one--phase single population model to
cases for which multiple populations are present. Open source Julia--based
software is available on GitHub to replicate all results in this study.


  Ecological networks describe the interactions between different species,
informing us of how they rely on one another for food, pollination and
survival. If a species in an ecosystem is under threat of extinction, it can
affect other species in the system and possibly result in their secondary
extinction as well. Understanding how (primary) extinctions cause secondary
extinctions on ecological networks has been considered previously using
computational methods. However, these methods do not provide an explanation for
the properties which make ecological networks robust, and can be
computationally expensive. We develop a new analytic model for predicting
secondary extinctions which requires no non-deterministic computational
simulation. Our model can predict secondary extinctions when primary
extinctions occur at random or due to some targeting based on the number of
links per species or risk of extinction, and can be applied to an ecological
network of any number of layers. Using our model, we consider how false
positives and negatives in network data affect predictions for network
robustness. We have also extended the model to predict scenarios in which
secondary extinctions occur once species lose a certain percentage of
interaction strength, and to model the loss of interactions as opposed to just
species extinction. From our model, it is possible to derive new analytic
results such as how ecological networks are most robust when secondary species
degree variance is minimised. Additionally, we show that both specialisation
and generalisation in distribution of interaction strength can be advantageous
for network robustness, depending upon the extinction scenario being
considered.


  Genomic (DNA) sequences encode an enormous amount of information for gene
regulation and protein synthesis. Similar to natural language models,
researchers have proposed foundation models in genomics to learn generalizable
features from unlabeled genome data that can then be fine-tuned for downstream
tasks such as identifying regulatory elements. Due to the quadratic scaling of
attention, previous Transformer-based genomic models have used 512 to 4k tokens
as context (<0.001% of the human genome), significantly limiting the modeling
of long-range interactions in DNA. In addition, these methods rely on
tokenizers or fixed k-mers to aggregate meaningful DNA units, losing single
nucleotide resolution where subtle genetic variations can completely alter
protein function via single nucleotide polymorphisms (SNPs). Recently, Hyena, a
large language model based on implicit convolutions was shown to match
attention in quality while allowing longer context lengths and lower time
complexity. Leveraging Hyena's new long-range capabilities, we present
HyenaDNA, a genomic foundation model pretrained on the human reference genome
with context lengths of up to 1 million tokens at the single nucleotide-level -
an up to 500x increase over previous dense attention-based models. HyenaDNA
scales sub-quadratically in sequence length (training up to 160x faster than
Transformer), uses single nucleotide tokens, and has full global context at
each layer. We explore what longer context enables - including the first use of
in-context learning in genomics. On fine-tuned benchmarks from the Nucleotide
Transformer, HyenaDNA reaches state-of-the-art (SotA) on 12 of 18 datasets
using a model with orders of magnitude less parameters and pretraining data. On
the GenomicBenchmarks, HyenaDNA surpasses SotA on 7 of 8 datasets on average by
+10 accuracy points. Code at https://github.com/HazyResearch/hyena-dna.


  Protein-DNA interaction is critical for life activities such as replication,
transcription, and splicing. Identifying protein-DNA binding residues is
essential for modeling their interaction and downstream studies. However,
developing accurate and efficient computational methods for this task remains
challenging. Improvements in this area have the potential to drive novel
applications in biotechnology and drug design. In this study, we propose a
novel approach called CLAPE, which combines a pre-trained protein language
model and the contrastive learning method to predict DNA binding residues. We
trained the CLAPE-DB model on the protein-DNA binding sites dataset and
evaluated the model performance and generalization ability through various
experiments. The results showed that the AUC values of the CLAPE-DB model on
the two benchmark datasets reached 0.871 and 0.881, respectively, indicating
superior performance compared to other existing models. CLAPE-DB showed better
generalization ability and was specific to DNA-binding sites. In addition, we
trained CLAPE on different protein-ligand binding sites datasets, demonstrating
that CLAPE is a general framework for binding sites prediction. To facilitate
the scientific community, the benchmark datasets and codes are freely available
at https://github.com/YAndrewL/clape.


  Chronic Hepatitis B (CHB) is an independent risk factor for hepatocellular
carcinoma (HCC) initiation without cirrhosis occurrence. Apart from the
favorable effects of some antiviral drugs following tumor resection on the
survival of HCC patients, the use of these agents is essential lifelong. Thus,
designing the target-oriented therapeutic strategies to increase life
expectancy in HCC patients would be very important. The present study aimed to
identify downregulated hub genes and enriched pathways in HB-related HCC using
a systems biology-based approach. Microarray data of GSE121248 were downloaded
from gene expression omnibus (GEO) database. The differentially expressed genes
(DEGs) with the cut-off criteria of adjusted p < 0.05 and Log Fold-change (FC)
< -1.5 were selected. Then, the genes with the highest centrality were
detected. Finally, the prognostic values of the hub genes were assessed. Six
under-expressed hub genes with the highest interaction degree, Betweenness and
Eigenvector centrality were including IGF-1, PTGS2, PLG, HGF, ESR-1 and CYP2B6.
Among genes with high centrality, several genes including CYP2C9, ESR-1, CXCL2,
CYP2C8, IGF-1, CYP3A4, CYP2E1, CERPINE-1 and PXR were prognostic in HCC. The
important repressed pathways were including metabolic pathways and PI3K-Akt and
chemokine signaling pathways. The under-expression of several genes implicated
in metabolism, differentiation and chemotaxis might be a hallmark of the
progression of HCC that can be considered as diagnostic and therapeutic
targets.


  There is an unmet need to develop medications or drug combinations which can
stop advancement of NASH to liver cirrhosis and HCC. Therefore, identifying key
biomarkers based on overall survival and the exploring biological processes
involved in the pathogenesis and progression of NASH toward cirrhosis and HCC
to improve therapeutic interventions is necessary. The microarray dataset from
the GPL13667 platform was downloaded from the Gene Expression Omnibus (GEO).
The inclusion criteria for the DEGs included an adjusted p-value <0.05 and a
log(2) fold change >1. In the first step, protein-protein interaction (PPI)
network of differentially expressed genes (DEGs) in every three groups (NASH,
Cirrhosis and HCC) was constructed using STRING online database. In the second
step, the DEGs of each group were imported to Cytoscape software separately,
and the genes with Degree >3 were selected. In the third step, the genes with
Degree >3 were imported to Gephi software (0.9.2 version) and the genes with
betweenness centrality >0 were selected. Venn diagram of these genes was
depicted for the NASH, cirrhotic and HCC groups. According to venn diagram, 96
(NASH), 30 (cirrhotic) and 213 (HCC) genes were specifically upregulated (based
on inclusion criteria). Among specific genes, 22 (NASH), 5 (cirrhotic) and 82
(HCC) genes were with poor overall survival. The overlap among the 3 groups
(NASH, Cirrhosis and HCC) contained 4 upregulated genes HLA-F, HLA-DPA1, TPM1
and YWHAZ. From 4 upregulated genes, only YWHAZ gene was with poor overall
survival. The present study detected new candidate genes and key biological
processes in NASH, cirrhosis and HCC based on overall survival and using in
silico analysis. Therefore, performing in vitro and in vivo researches to
verify the results is necessary.


  We develop a new structured compartmental model for the coevolutionary
dynamics between susceptible and infectious individuals in heterogeneous SI
epidemiological systems. In this model, the susceptible compartment is
structured by a continuous variable that represents the level of resistance to
infection of susceptible individuals, while the infectious compartment is
structured by a continuous variable that represents the viral load of
infectious individuals. We first formulate an individual-based model wherein
the dynamics of single individuals is described through stochastic processes,
which permits a fine-grain representation of individual dynamics and captures
stochastic variability in evolutionary trajectories amongst individuals. Next
we formally derive the mesoscopic counterpart of this model, which consists of
a system of coupled integro-differential equations for the population density
functions of susceptible and infectious individuals. Then we consider an
appropriately rescaled version of this system and we carry out formal
asymptotic analysis to derive the corresponding macroscopic model, which
comprises a system of coupled ordinary differential equations for the
proportions of susceptible and infectious individuals, the mean level of
resistance to infection of susceptible individuals, and the mean viral load of
infectious individuals. Overall, this leads to a coherent mathematical
representation of the coevolutionary dynamics between susceptible and
infectious individuals across scales. We provide well-posedness results for the
mesoscopic and macroscopic models, and we show that there is excellent
agreement between analytical results on the long-time behaviour of the
components of the solution to the macroscopic model, the results of Monte Carlo
simulations of the individual-based model, and numerical solutions of the
macroscopic model.


  Molecular property is usually observed with a limited number of samples, and
researchers have considered property prediction as a few-shot problem. One
important fact that has been ignored by prior works is that each molecule can
be recorded with several different properties simultaneously. To effectively
utilize many-to-many correlations of molecules and properties, we propose a
Graph Sampling-based Meta-learning (GS-Meta) framework for few-shot molecular
property prediction. First, we construct a Molecule-Property relation Graph
(MPG): molecule and properties are nodes, while property labels decide edges.
Then, to utilize the topological information of MPG, we reformulate an episode
in meta-learning as a subgraph of the MPG, containing a target property node,
molecule nodes, and auxiliary property nodes. Third, as episodes in the form of
subgraphs are no longer independent of each other, we propose to schedule the
subgraph sampling process with a contrastive loss function, which considers the
consistency and discrimination of subgraphs. Extensive experiments on 5
commonly-used benchmarks show GS-Meta consistently outperforms state-of-the-art
methods by 5.71%-6.93% in ROC-AUC and verify the effectiveness of each proposed
module. Our code is available at https://github.com/HICAI-ZJU/GS-Meta.


  Machine learning for molecules holds great potential for efficiently
exploring the vast chemical space and thus streamlining the drug discovery
process by facilitating the design of new therapeutic molecules. Deep
generative models have shown promising results for molecule generation, but the
benefits of specific inductive biases for learning distributions over small
graphs are unclear. Our study aims to investigate the impact of subgraph
structures and vocabulary design on distribution learning, using small drug
molecules as a case study. To this end, we introduce Subcover, a new
subgraph-based fragmentation scheme, and evaluate it through a two-step
variational auto-encoder. Our results show that Subcover's improved
identification of chemically meaningful subgraphs leads to a relative
improvement of the FCD score by 30%, outperforming previous methods. Our
findings highlight the potential of Subcover to enhance the performance and
scalability of existing methods, contributing to the advancement of drug
discovery.


  Contact matrices have become a key ingredient of modern epidemic models. They
account for the stratification of contacts for the age of individuals and, in
some cases, the context of their interactions. However, age and context are not
the only factors shaping contact structures and affecting the spreading of
infectious diseases. Socio-economic status (SES) variables such as wealth,
ethnicity, and education play a major role as well. Here, we introduce
generalized contact matrices capable of stratifying contacts across any number
of dimensions including any SES variable. We derive an analytical expression
for the basic reproductive number of an infectious disease unfolding on a
population characterized by such generalized contact matrices. Our results, on
both synthetic and real data, show that disregarding higher levels of
stratification might lead to the under-estimation of the reproductive number
and to a mis-estimation of the global epidemic dynamics. Furthermore, including
generalized contact matrices allows for more expressive epidemic models able to
capture heterogeneities in behaviours such as different levels of adoption of
non-pharmaceutical interventions across different groups. Overall, our work
contributes to the literature attempting to bring socio-economic, as well as
other dimensions, to the forefront of epidemic modeling. Tackling this issue is
crucial for developing more precise descriptions of epidemics, and thus to
design better strategies to contain them.


  Identifying and understanding the large-scale biodiversity patterns in time
and space is vital for conservation and addressing fundamental ecological and
evolutionary questions. Network-based methods have proven useful for
simplifying and highlighting important structures in species distribution data.
However, current network-based biogeography approaches cannot exploit the
evolutionary information available in phylogenetic data. We introduce a method
for incorporating evolutionary relationships into species occurrence networks
to produce more biologically informative and robust bioregions. To keep the
bipartite network structure where bioregions are grid cells indirectly
connected through shared species, we incorporate the phylogenetic tree by
connecting ancestral nodes to the grid cells where their descendant species
occur. To incorporate the whole tree without destroying the spatial signal of
narrowly distributed species or ancestral nodes, we weigh tree nodes by the
geographic information they provide. For a more detailed analysis, we enable
integration of the evolutionary relationships at a specific time in the tree.
By sweeping through the phylogenetic tree in time, our method interpolates
between finding bioregions based only on distributional data and finding
spatially segregated clades, uncovering evolutionarily distinct bioregions at
different time slices. We also introduce a way to segregate the connections
between evolutionary branches at a selected time to enable exploration of
overlapping evolutionarily distinct regions. We have implemented these methods
in Infomap Bioregions, an interactive web application that makes it easy to
explore the possibly hierarchical and fuzzy patterns of biodiversity on
different scales in time and space.


  Even in large systems, the effect of noise arising from when populations are
initially small can persist to be measurable on the macroscale. A deterministic
approximation to a stochastic model will fail to capture this effect, but it
can be accurately approximated by including an additional random time-shift to
the initial conditions. We present a efficient numerical method to compute this
time-shift distribution for a large class of stochastic models. The method
relies on differentiation of certain functional equations, which we show can be
effectively automated by deriving rules for different types of model rates that
arise commonly when mass-action mixing is assumed. Explicit computation of the
time-shift distribution can be used to build a practical tool for the efficient
generation of macroscopic trajectories of stochastic population models, without
the need for costly stochastic simulations. Full code is provided to implement
this and we demonstrate our method on an epidemic model and a model of
within-host viral dynamics.


  The framework of iterated Prisoner's Dilemma (IPD) is commonly used to study
direct reciprocity and cooperation, with a focus on the assessment of the
generosity and reciprocal fairness of an IPD strategy in one-on-one settings.
In order to understand the persistence and resilience of reciprocal
cooperation, here we study long-term population dynamics of IPD strategies
using the Moran process where stochastic dynamics of strategy competition can
lead to the rise and fall of cooperation. Although prior work has included a
handful of typical IPD strategies in the consideration, it remains largely
unclear which type of IPD strategies is pivotal in steering the population away
from defection and providing an escape hatch for establishing cooperation. We
use a network-based approach to analyze and characterize networks of
evolutionary pathways that bridge transient episodes of evolution dominated by
depressing defection and ultimately catalyze the evolution of reciprocal
cooperation in the long run. We group IPD strategies into three types according
to their stationary cooperativity with an unconditional cooperator: the good
(fully cooperative), the bad (fully exploitive), and the ugly (in between the
former two types). We consider the mutation-selection equilibrium with rare
mutations and quantify the impact of the presence versus absence of any given
IPD strategy on the resulting population equilibrium. We identify catalysts
(certain IPD strategies) as well as bridges (particular evolutionary pathways)
that are most crucial for boosting the abundance of good types and suppressing
that of bad types or having the highest betweenness centrality. Our work has
practical implications and broad applicability to real-world cooperation
problems by leveraging catalysts and bridges that are capable of strengthening
persistence and resilience.


  In this report we present the work carried out during the Complexity72h
workshop, held at IFISC in Palma de Mallorca, Spain, 26-30 June 2023. We
describe a temporal network-theoretic approach to study fluid flows with
applications to marine ecology. The network representation is derived from the
Lagrangian fluid dynamics and represents fluid transportation between patches
of the sea. It is a directed, weighted and time-dependent network. This
approach enables us to use advanced network-theoretic tools for analysis and
modeling. A common approximation adopted in the literature consists in using an
aggregated time-independent network representation of the fluid flow. In this
report we focus in particular on the role played by the temporal component and
to the information loss related to neglecting that dimension and inspect the
role played by seasonal or long time-period variations. We conduct an analysis
of basic network features of the aggregated and temporal graphs, we analyze
their community structure and we model population dynamics of marine lives
driven by the flow. Ultimately, we determine that time-independent
approximations can effectively represent long-term transportation evolution
spanning multiple years. However, for an accurate depiction of transportation
within a single year, it is necessary to incorporate explicit time-dependence
in the transport matrix to account for seasonality.


  When we apply comparative phylogenetic analyses to genome data, it is a
well-known problem and challenge that some of given species (or taxa) often
have missing genes. In such a case, we have to impute a missing part of a gene
tree from a sample of gene trees. In this short paper we propose a novel method
to infer a missing part of a phylogenetic tree using an analogue of a classical
linear regression in the setting of tropical geometry. In our approach, we
consider a tropical polytope, a convex hull with respect to the tropical metric
closest to the data points. We show a condition that we can guarantee that an
estimated tree from our method has at most four Robinson-Foulds (RF) distance
from the ground truth and computational experiments with simulated data show
our method works well.


  Diffusion models have been successful on a range of conditional generation
tasks including molecular design and text-to-image generation. However, these
achievements have primarily depended on task-specific conditional training or
error-prone heuristic approximations. Ideally, a conditional generation method
should provide exact samples for a broad range of conditional distributions
without requiring task-specific training. To this end, we introduce the Twisted
Diffusion Sampler, or TDS. TDS is a sequential Monte Carlo (SMC) algorithm that
targets the conditional distributions of diffusion models. The main idea is to
use twisting, an SMC technique that enjoys good computational efficiency, to
incorporate heuristic approximations without compromising asymptotic exactness.
We first find in simulation and on MNIST image inpainting and class-conditional
generation tasks that TDS provides a computational statistical trade-off,
yielding more accurate approximations with many particles but with empirical
improvements over heuristics with as few as two particles. We then turn to
motif-scaffolding, a core task in protein design, using a TDS extension to
Riemannian diffusion models. On benchmark test cases, TDS allows flexible
conditioning criteria and often outperforms the state of the art.


  Colon cancer is a prevalent gastrointestinal malignancy arising in the colon.
Ulcerative colitis(UC) is one of the risk factors of colorectal cancer. The
detection of under-expressed biomarkers and molecular mechanisms in UC and
colon cancer can lead to effective management of colitis-associated cancer. A
total of two mRNA expression datasets (GSE87473 and GSE44076) were downloaded
from the Gene Expression Omnibus (GEO) database. GEO2R was used to screen
differentially expressed genes (DEGs) between extensive ulcerative colitis
samples and healthy samples, limited ulcerative colitis samples and healthy
samples, and colon cancer samples and healthy samples. In extensive ulcerative
colitis, limited ulcerative colitis and colon cancer groups, 95,69 and 635
under-expressed genes with adjusted p-value<0.05 and log(2) fold change<-2 were
detected respectively. Using Cytoscape software, the genes with degree> 15
including CLCA1, SLC26A3, SI, KIT, HPGDS, NR1H4, ADIPOQ, PPARGC1A, GCG, MS4A12,
GUCA2A and FABP1 were screened as hub under-expressed genes in colon cancer. In
extensive ulcerative colitis, the genes with degree>5 including ABCB1, ABCG2,
UGT1A6, CYP2B6 and AQP8 were identified as hub genes. Moreover, the genes
including NR1H4, CYP2B6, ABCB1, ABCG2, UGT2A3 and PLA2G12B were detected as hub
genes with degree>5 in limited ulcerative colitis. According to inclusion
criteria and venn diagram, the downregulated gene NR1H4 was common gene in
limited ulcerative colitis and colon cancer. The current in silico study showed
that downregulation of CLCA1, PPARGC1A and AQP8 genes may increase cancer cell
invasion and metastasis ability.


  When an infectious disease propagates throughout society, the incidence
function may rise at first due to an increase in pathogenicity and then
decrease due to inhibitory effects until it reaches saturation. Effective
vaccination and treatment are very helpful for controlling the effects of such
infectious diseases. To analyze the impacts of these diseases, we proposed a
new compartmental model with a generalized non-linear incidence function,
vaccination function, and treatment function, along with time delays in the
respective functions, which show how its monotonic features influence the
stability of the model. Fundamental properties of a model, such as positivity,
boundedness, and the existence of equilibria, are examined in this work. The
basic reproduction number has been computed, and correlative studies for local
stability in view of the basic reproduction number have been examined at the
disease-free and endemic equilibrium points. A delay-independent global
stability result has been established, and to be more precise, we explicitly
derived the result on global stability by restricting delay parameters within a
very specific range. Furthermore, numerical simulations and some examples based
on COVID-19 real-time data are pointed out to emphasize the significance of how
the disease's dynamical behavior is characterized by various functions for
controlling the spread of disease in a population and to justify the
mathematical conclusions.


  The worldwide COVID-19 pandemic has led to a significant growth of interest
in the development of mathematical models that allow to describe effects such
as social distancing measures, the development of vaccines, and mutations.
Several of these models are based on concepts from soft matter theory.
Considerably less well investigated is the reverse direction, i.e., how results
from epidemiological research can be of interest for the physics of colloids
and polymers. In this work, we consider the SIR-DDFT model, a combination of
the susceptible-infected-recovered (SIR) model from epidemiology with dynamical
density functional theory (DDFT) from nonequilibrium soft matter physics, which
allows for an explicit modeling of social distancing. We extend the SIR-DDFT
model both from an epidemiological perspective by incorporating vaccines,
asymptomaticity, reinfections, and mutations, and from a soft matter
perspective by incorporating noise and self-propulsion and by deriving a phase
field crystal (PFC) model that allows for a simplified description. On this
basis, we investigate via computer simulations how epidemiological models are
affected by the presence of non-reciprocal interactions. This is done in a
numerical study of a zombie outbreak.


  The ability to design novel proteins with higher fitness on a given task
would be revolutionary for many fields of medicine. However, brute-force search
through the combinatorially large space of sequences is infeasible. Prior
methods constrain search to a small mutational radius from a reference
sequence, but such heuristics drastically limit the design space. Our work
seeks to remove the restriction on mutational distance while enabling efficient
exploration. We propose Gibbs sampling with Graph-based Smoothing (GGS) which
iteratively applies Gibbs with gradients to propose advantageous mutations
using graph-based smoothing to remove noisy gradients that lead to false
positives. Our method is state-of-the-art in discovering high-fitness proteins
with up to 8 mutations from the training set. We study the GFP and AAV design
problems, ablations, and baselines to elucidate the results. Code:
https://github.com/kirjner/GGS


  In data-driven drug discovery, designing molecular descriptors is a very
important task. Deep generative models such as variational autoencoders (VAEs)
offer a potential solution by designing descriptors as probabilistic latent
vectors derived from molecular structures. These models can be trained on large
datasets, which have only molecular structures, and applied to transfer
learning. Nevertheless, the approximate posterior distribution of the latent
vectors of the usual VAE assumes a simple multivariate Gaussian distribution
with zero covariance, which may limit the performance of representing the
latent features. To overcome this limitation, we propose a novel molecular deep
generative model that incorporates a hierarchical structure into the
probabilistic latent vectors. We achieve this by a denoising diffusion
probabilistic model (DDPM). We demonstrate that our model can design effective
molecular latent vectors for molecular property prediction from some
experiments by small datasets on physical properties and activity. The results
highlight the superior prediction performance and robustness of our model
compared to existing approaches.


  In this study, we analyze discordance between the transcriptome and proteome
using paired scRNA-Seq and multiplexed spatial proteomics data from HuBMAP. Our
findings highlight persistent transcripts in key immune markers, including
CD45-RO, Ki67, CD45, CD20, and HLA-DR. CD45-RO is consistently expressed in
memory T cells, while Ki67, associated with cell proliferation, also displays
sustained expression. Furthermore, HLA-DR, part of the MHC class II molecules,
demonstrates continuous expression, possibly crucial for APCs to trigger an
effective immune response. This investigation provides novel insights into the
complexity of gene expression regulation and protein function.


  The cooperative binding of molecular agents onto a substrate is pervasive in
living systems. To study whether a system shows cooperativity, one can rely on
a fluctuation analysis of quantities such as the number of substrate-bound
units and the residence time in an occupancy state. Since the relative standard
deviation from the statistical mean monotonically decreases with the number of
binding sites, these techniques are only suitable for small enough systems,
such as those implicated in stochastic processes inside cells. Here, we present
a general-purpose grand canonical Hamiltonian description of a small
one-dimensional (1D) lattice gas with either nearest-neighbor or long-range
interactions as prototypical examples of cooperativity-influenced adsorption
processes. First, we elucidate how the strength and sign of the interaction
potential between neighboring bound particles on the lattice determine the
intensity of the fluctuations of the mean occupancy. We then employ this
relationship to compare the theoretical predictions of our model to data from
single molecule experiments on bacterial flagellar motors (BFM) of E. coli. In
this way, we find evidence that cooperativity controls the mechano-sensitive
dynamical assembly of the torque-generating units, the so-called stator units,
onto the BFM. Finally, we estimate the stator-stator interaction potential and
attempt to quantify the adaptability of the BFM.


  The interior of a living cell is an active, fluctuating, and crowded
environment. Yet, it maintains a high level of coherent organization, which is
readily apparent in the intracellular transport network. Membrane-bound
compartments called endosomes play a key role in carrying cargo, in conjunction
with myriad components including cargo adaptor proteins, membrane sculptors,
motor proteins, and the cytoskeleton. These components coordinate to
effectively navigate the crowded cell interior and transport cargo to specific
intracellular locations, even though the underlying protein interactions and
enzymatic reactions exhibit stochastic behavior. A major challenge is to
measure, analyze, and understand how, despite the inherent stochasticity of the
constituent processes, the collective outcomes show an emergent spatiotemporal
order that is precise and robust. This review focuses on this intriguing
dichotomy, providing insights into the known mechanisms of noise suppression
and noise utilization in intracellular transport processes, and also identifies
opportunities for future inquiry.


  The COVID-19 pandemic has created unprecedented challenges for governments
and healthcare systems worldwide, highlighting the critical importance of
understanding the factors that contribute to virus transmission. This study
aimed to identify the most influential age groups in COVID-19 infection rates
at the US county level using the Modified Morris Method and deep learning for
time series. Our approach involved training the state-of-the-art time-series
model Temporal Fusion Transformer on different age groups as a static feature
and the population vaccination status as the dynamic feature. We analyzed the
impact of those age groups on COVID-19 infection rates by perturbing individual
input features and ranked them based on their Morris sensitivity scores, which
quantify their contribution to COVID-19 transmission rates. The findings are
verified using ground truth data from the CDC and US Census, which provide the
true infection rates for each age group. The results suggest that young adults
were the most influential age group in COVID-19 transmission at the county
level between March 1, 2020, and November 27, 2021. Using these results can
inform public health policies and interventions, such as targeted vaccination
strategies, to better control the spread of the virus. Our approach
demonstrates the utility of feature sensitivity analysis in identifying
critical factors contributing to COVID-19 transmission and can be applied in
other public health domains.


  In intracellular transports, motor proteins transport macromolecules as
cargos to desired locations by moving on biopolymers such as microtubules.
Recent experiments suggest that cargos that can associate motor proteins during
their translocation have larger run-length, association time and can overcome
the motor traffic on microtubule tracks. Here, we model the dynamics of a cargo
that can associate at the most m free motors present on the track as obstacles
to its motion. The proposed models display competing effects of association and
crowding, leading to a peak in the run-length with the free motor density. This
result is consistent with past experimental observations. For m=2 and 3, we
show that this feature is governed by the largest eigenvalue of the transition
matrix describing the cargo dynamics. In all the above cases, free motors are
assumed to be present as stalled obstacles. We finally compare simulation
results for the run-length for general scenarios where the free motors undergo
processive motion in addition to binding and unbinding to or from the
microtubule.


  Hydrophobic domains provide specific microenvironment for essential
functional activities in life. Herein, we studied how the coassembling of
peptides with catecholamines regulate the hydrophobic domain-containing
nanostructures for fluorescence enhancement. By peptide encoding and
coassembling with catecholamines of different hydrophilicities, a series of
hierarchical assembling systems were constructed. In combination with molecular
dynamics simulation, we experimentally discovered the hydrophobic domain of
chromophore microenvironment regulates the fluorescence of coassembled
nanostructures. Our results shed light on the rational design of fluorescent
bio-coassembled nanoprobes for biomedical applications.


  Prediction of protein-ligand interactions (PLI) plays a crucial role in drug
discovery as it guides the identification and optimization of molecules that
effectively bind to target proteins. Despite remarkable advances in deep
learning-based PLI prediction, the development of a versatile model capable of
accurately scoring binding affinity and conducting efficient virtual screening
remains a challenge. The main obstacle in achieving this lies in the scarcity
of experimental structure-affinity data, which limits the generalization
ability of existing models. Here, we propose a viable solution to address this
challenge by introducing a novel data augmentation strategy combined with a
physics-informed graph neural network. The model showed significant
improvements in both scoring and screening, outperforming task-specific deep
learning models in various tests including derivative benchmarks, and notably
achieving results comparable to the state-of-the-art performance based on
distance likelihood learning. This demonstrates the potential of this approach
to drug discovery.


  PROTACs, as a highly promising new. therapeutic paradigm, have attracted
widespread attention from the academic and pharmaceutical communities in recent
years. To date, the design and validation of PROTACs molecule's druggability
primarily rely on experimental approaches, making the development process of
this kind of drug molecule time-consuming. Computer-aided tools for PROTACs
design may offer a potential solution to expedite the design process and
enhance its efficiency. This mini review briefly summarizes the in silico tools
for PROTACs drug molecule design reported recently.


  Since they became observable, neuron morphologies have been informally
compared with biological trees but they are studied by distinct communities,
neuroscientists, and ecologists. The apparent structural similarity suggests
there may be common quantitative rules and constraints. However, there are also
reasons to believe they should be different. For example, while the
environments of trees may be relatively simple, neurons are constructed by a
complex iterative program where synapses are made and pruned. This complexity
may make neurons less self-similar than trees. Here we test this hypothesis by
comparing the features of segmented sub-trees with those of the whole tree. We
indeed find more self-similarity within actual trees than neurons. At the same
time, we find that many other features are somewhat comparable across the two.
Investigation of shapes and behaviors promises new ways of conceptualizing the
form-function link.


  Relation-aware graph structure embedding is promising for predicting
multi-relational drug-drug interactions (DDIs). Typically, most existing
methods begin by constructing a multi-relational DDI graph and then learning
relation-aware graph structure embeddings (RaGSEs) of drugs from the DDI graph.
Nevertheless, most existing approaches are usually limited in learning RaGSEs
of new drugs, leading to serious over-fitting when the test DDIs involve such
drugs. To alleviate this issue, we propose a novel DDI prediction method based
on relation-aware graph structure embedding with co-contrastive learning,
RaGSECo. The proposed RaGSECo constructs two heterogeneous drug graphs: a
multi-relational DDI graph and a multi-attribute drug-drug similarity (DDS)
graph. The two graphs are used respectively for learning and propagating the
RaGSEs of drugs, aiming to ensure all drugs, including new ones, can possess
effective RaGSEs. Additionally, we present a novel co-contrastive learning
module to learn drug-pairs (DPs) representations. This mechanism learns DP
representations from two distinct views (interaction and similarity views) and
encourages these views to supervise each other collaboratively to obtain more
discriminative DP representations. We evaluate the effectiveness of our RaGSECo
on three different tasks using two real datasets. The experimental results
demonstrate that RaGSECo outperforms existing state-of-the-art prediction
methods.


  In environments that vary frequently and unpredictably, bet-hedgers can
overtake the population. Diversifying bet-hedgers have a diverse set of
offspring so that, no matter the conditions they find themselves in, at least
some offspring will have high fitness. In contrast, conservative bet-hedgers
have a set of offspring that all have an in-between phenotype compared to the
specialists. Here, we use an evolutionary algorithm of gene regulatory networks
to de novo evolve the two strategies and investigate their relative success in
different parameter settings. We found that diversifying bet-hedgers almost
always evolved first, but then eventually got outcompeted by conservative
bet-hedgers. We argue that even though similar selection pressures apply to the
two bet-hedger strategies, conservative bet-hedgers could win due to the
robustness of their evolved networks, in contrast to the sensitive networks of
the diversifying bet-hedgers. These results reveal an unexplored aspect of the
evolution of bet-hedging that could shed more light on the principles of
biological adaptation in variable environmental conditions.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  Within the living cell, protein molecules perform specific functions,
typically by interacting with other proteins, DNA, RNA or small molecules. They
take on a specific three dimensional structure, encoded by its amino acid
sequence, which allows them to function within the cell. Hence, the
understanding of a protein's function is tightly coupled to its sequence and
its three dimensional structure. Before going into protein structure analysis
and prediction, and protein folding and dynamics, here, we give a short and
concise introduction into the basics of protein structures.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  The Protein DataBank (PDB) contains a wealth of structural information. In
order to investigate the similarity between different proteins in this
database, one can compare the primary sequence through pairwise alignment and
calculate the sequence identity (or similarity) over the two sequences. This
strategy will work particularly well if the proteins you want to compare are
close homologs. However, in this chapter we will explain that a structural
comparison through structural alignment will give you much more valuable
information, that allows you to investigate similarities between proteins that
cannot be discovered by comparing the sequences alone.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  Structural bioinformatics involves a variety of computational methods, all of
which require input data. Typical inputs include protein structures and
sequences, which are usually retrieved from a public or private database. This
chapter introduces several key resources that make such data available, as well
as a handful of tools that derive additional information from experimentally
determined or computationally predicted protein structures and sequences.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  Some structural properties of proteins that are closely linked to their
function may be easier (or much faster) to predict from sequence than the
complete tertiary structure; for example, secondary structure, surface
accessibility, flexibility, disorder, interface regions or hydrophobic patches.
Serving as building blocks for the native protein fold, these structural
properties also contain important structural and functional information not
apparent from the amino acid sequence. Here, we will first give an introduction
into the application of machine learning for structural property prediction,
and explain the concepts of cross-validation and benchmarking. Next, we will
review various methods that incorporate knowledge of these concepts to predict
those structural properties, such as secondary structure, surface
accessibility, disorder and flexibility, and aggregation.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  There are still huge gaps in understanding the molecular function of
proteins. This raises the question on how we may predict protein function, when
little to no knowledge from direct experiments is available. Protein function
is a broad concept which spans different scales: from quantum scale effects for
catalyzing enzymatic reactions, to phenotypes that manifest at the organism
level. In fact, many of these functional scales are entirely different research
areas. Here, we will consider prediction of a smaller range of functions,
roughly spanning the protein residue-level up to the pathway level. We will
give a conceptual overview of which functional aspects of proteins we can
predict, which methods are currently available, and how well they work in
practice.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  In this chapter we explore basic physical and chemical concepts required to
understand protein folding. We introduce major (de)stabilising factors of
folded protein structures such as the hydrophobic effect and backbone entropy.
In addition, we consider different states along the folding pathway, as well as
natively disordered proteins and aggregated protein states. In this chapter, an
intuitive understanding is provided about the protein folding process, to
prepare for the next chapter on the thermodynamics of protein folding. In
particular, it is emphasized that protein folding is a stochastic process and
that proteins unfold and refold in a dynamic equilibrium. The effect of
temperature on the stability of the folded and unfolded states is also
explained.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  In the previous chapter, "Introduction to Protein Folding", we introduced the
concept of free energy and the protein folding landscape. Here, we provide a
deeper, more formal underpinning of free energy in terms of the entropy and
enthalpy; to this end, we will first need to better define the meaning of
equilibrium, entropy and enthalpy. When we understand these concepts, we will
come back for a more quantitative explanation of protein folding and dynamics.
We will discuss the influence of temperature on the free energy landscape, and
the difference between microstates and macrostates.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  We know that many proteins have functional motions, and in Chapter "Structure
Determination" we already introduced the famous example of the allosteric
cooperative binding of oxygen to the haem group in hemoglobin. However,
experimentally, such motions are hard to observe. Here, we will introduce MD
simulations to investigate the dynamic behaviour of proteins. In a simulation
the forces and interactions between particles are used to numerically derive
the resulting three-dimensional movement of these particles over a certain
time-scale. We will also highlight some applications, and will see how
simulation results may be interpreted.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  In the previous chapter "Molecular Dynamics" we have considered protein
simulations from a dynamical point of view, using Newton's laws. In the current
Chapter, we first take a step back and return to the bare minimum needed to
simulate proteins, and show that proteins may be simulated in a more simple
fashion, using the partition function directly. This means we do not have to
calculate explicit forces, velocities, moments and do not even consider time
explicitly. Instead, we will rely on the fact that for most systems we will
want to simulate, the system is in a dynamic equilibrium; and that we want to
find the most stable states in such systems by determining the relative
stabilities between those states.


  Based on the normal distribution and its properties, i.e., average and
variance, Fisher works have provided a conceptual framework to identify
genotype-phenotype associations. While Fisher intuition has proved fruitful
over the past century, the current demands for higher mapping precisions have
led to the formulation of a new genotype-phenotype association method a.k.a.
GIFT (Genomic Informational Field Theory). Not only is the method more powerful
in extracting information from genotype and phenotype datasets, GIFT can also
deal with any phenotype distribution density function. Here we apply GIFT to a
hypothetical Cauchy-distributed phenotype. As opposed to the normal
distribution that restricts fluctuations to a finite variance defined by the
bulk of the distribution, Cauchy distribution embraces large phenotypic
fluctuations and as a result, averages and variances from Cauchy-distributed
phenotypes cannot be defined mathematically. While classic genotype-phenotype
association methods (GWAS) are unable to function without proper average and
variance, it is demonstrated here that GIFT can associate genotype to phenotype
in this case. As phenotypic plasticity, i.e., phenotypic fluctuation, is
central to surviving sudden environmental changes, by applying GIFT the unique
characteristic of the genotype permitting evolution of biallelic organisms to
take place is determined in this case.


  The introduction of the spatial Lambda-Fleming-Viot model (LV) in population
genetics was mainly driven by the pioneering work of Alison Etheridge, in
collaboration with Nick Barton and Amandine V\'eber about ten years ago (1,2).
The LV model provides a sound mathematical framework for describing the
evolution of a population of related individuals along a spatial continuum. It
alleviates the "pain in the torus" issue with Wright and Mal\'ecot's isolation
by distance model and is sampling consistent, making it a tool of choice for
statistical inference. Yet, little is known about the potential connections
between the LV and other stochastic processes generating trees and the spatial
coordinates along the corresponding lineages. This work focuses on a version of
the LV whereby lineages move infinitely rapidly over infinitely small
distances. Using simulations, we show that the induced LV tree-generating
process is well approximated by a birth-death model. Our results also indicate
that Brownian motions modelling the movements of lineages along birth-death
trees do not generally provide a good approximation of the LV due to habitat
boundaries effects that play an increasingly important role in the long run.
Finally, we describe efficient algorithms for fast simulation of the backward
and forward in time versions of the LV model.


  Avian prestin is sensitive to membrane thickness as much as mammalian
prestin, which undergoes conformational transitions in membrane area and
thereby drives length changes of the cylindrical cell body of outer hair cells.
The membrane thickness dependence of mammalian prestin stems from changes in
hydrophobic profile in conformational states, accompanied by changes in their
membrane area. Even though such area changes are not detected for avian
prestin, it nonetheless bends hair bundles of avian short hair cells. Here it
is suggested that the motile function of avian prestin can be based on
conformational transitions involving shearing deformation of the membrane
protein, which also leads to membrane thickness sensitivity.


  Comparisons of single-cell RNA sequencing (scRNA-seq) data across species can
reveal links between cellular gene expression and the evolution of cell
functions, features, and phenotypes. These comparisons invoke evolutionary
histories, as depicted with phylogenetic trees, that define relationships
between species, genes, and cells. Here we illustrate a tree-based framework
for comparing scRNA-seq data, and contrast this framework with existing
methods. We describe how we can use trees to identify homologous and comparable
groups of genes and cells, based on their predicted relationship to genes and
cells present in the common ancestor. We advocate for mapping data to branches
of phylogenetic trees to test hypotheses about the evolution of cellular gene
expression. We describe the kinds of data that can be compared, and the types
of questions that each comparison has the potential to address. Finally, we
reconcile species phylogenies, gene phylogenies, cell phylogenies, and cell
lineages as different representations of the same concept: the tree of cellular
life. By integrating phylogenetic approaches into scRNA-seq analyses, we can
overcome challenges for building informed comparisons across species, and
robustly test hypotheses about gene and cell evolution.


  The problem of comparing probability distributions is at the heart of many
tasks in statistics and machine learning and the most classical comparison
methods assume that the distributions occur in spaces of the same dimension.
Recently, a new geometric solution has been proposed to address this problem
when the measures live in Euclidean spaces of differing dimensions. Here, we
study the same problem of comparing probability distributions of different
dimensions in the tropical geometric setting, which is becoming increasingly
relevant in computations and applications involving complex, geometric data
structures. Specifically, we construct a Wasserstein distance between measures
on different tropical projective tori - the focal metric spaces in both theory
and applications of tropical geometry - via tropical mappings between
probability measures. We prove equivalence of the directionality of the maps,
whether starting from the lower dimensional space and mapping to the higher
dimensional space or vice versa. As an important practical implication, our
work provides a framework for comparing probability distributions on the spaces
of phylogenetic trees with different leaf sets.


  The dynamics of species communities are typically modelled considering fixed
parameters for species interactions. The problem of over-parameterization that
ensues when considering large communities has been overcome by sampling species
interactions from a probability distribution. However, species interactions are
not fixed in time, but they can change on a time scale comparable to population
dynamics. Here we investigate the impact of time-dependent species interactions
using the generalized Lotka-Volterra model, which serves as a paradigmatic
theoretical framework in several scientific fields. In this work we model
species interactions as stochastic colored noise. Assuming a large number of
species and a steady state, we obtain analytical predictions for the species
abundance distribution, which matches well empirical observations. In
particular, our results suggest the absence of extinctions, unlike scenarios
with fixed species interactions.


  Describing reaction rates in stochastic bio-circuits is commonly done by
directly introducing the deterministically deduced Hill function into the
master equation. However, when fluctuations in enzymatic reaction rates are not
neglectable, the Hill function must be derived, considering all the involved
stochastic reactions. In this work, we derived the stochastic version of the
Hill function from the master equation of the complete set of reactions that,
in the macroscopic limit, lead to the Hill function reaction rate. We performed
a series expansion around the average values of the concentrations, which
allowed us to find corrections for the deterministic Hill function. This
process allowed us to quantify the fluctuations of enzymatic reactions. We
found that the underlying variability in propensity rates of gene regulatory
networks has an important non-linear effect that reduces the intrinsic
fluctuations of the mRNA and protein concentrations.


  Stochastic models of gene expression are typically formulated using the
chemical master equation, which can be solved exactly or approximately using a
repertoire of analytical methods. Here, we provide a tutorial review of an
alternative approach based on queuing theory that has rarely been used in the
literature of gene expression. We discuss the interpretation of six types of
infinite server queues from the angle of stochastic single-cell biology and
provide analytical expressions for the stationary and non-stationary
distributions and/or moments of mRNA/protein numbers, and bounds on the Fano
factor. This approach may enable the solution of complex models which have
hitherto evaded analytical solution.


  The site frequency spectrum (SFS) is a widely used summary statistic of
genomic data, offering a simple means of inferring the evolutionary history of
a population. Motivated by recent evidence for the role of neutral evolution in
cancer, we examine the SFS of neutral mutations in an exponentially growing
population. Whereas recent work has focused on the mean behavior of the SFS in
this scenario, here, we investigate the first-order asymptotics of the
underlying stochastic process. Using branching process techniques, we show that
the SFS of a Galton-Watson process evaluated at a fixed time converges almost
surely to a random limit. We also show that the SFS evaluated at the stochastic
time at which the population first reaches a certain size converges in
probability to a constant. Finally, we illustrate how our results can be used
to construct consistent estimators for the extinction probability and the
effective mutation rate of a birth-death process.


  The time-varying reproduction number $R(t)$ measures the number of new
infections per infectious individual and is closely correlated with the time
series of infection incidence by definition. The timings of actual infections
are rarely known, and analysis of epidemics usually relies on time series data
for other outcomes such as symptom onset. A common implicit assumption, when
estimating $R(t)$ from an epidemic time series, is that $R(t)$ has the same
relationship with these downstream outcomes as it does with the time series of
incidence. However, this assumption is unlikely to be valid given that most
epidemic time series are not perfect proxies of incidence. Rather they
represent convolutions of incidence with uncertain delay distributions. Here we
define the apparent time-varying reproduction number, $R_A(t)$, the
reproduction number calculated from a downstream epidemic time series and
demonstrate how differences between $R_A(t)$ and $R(t)$ depend on the
convolution function. The mean of the convolution function sets a time offset
between the two signals, whilst the variance of the convolution function
introduces a relative distortion between them. We present the convolution
functions of epidemic time series that were available during the SARS-CoV-2
pandemic. Infection prevalence, measured by random sampling studies, presents
fewer biases than other epidemic time series. Here we show that additionally
the mean and variance of its convolution function were similar to that obtained
from traditional surveillance based on mass-testing and could be reduced using
more frequent testing, or by using stricter thresholds for positivity.
Infection prevalence studies continue to be a versatile tool for tracking the
temporal trends of $R(t)$, and with additional refinements to their study
protocol, will be of even greater utility during any future epidemics or
pandemics.


  Several studies have developed dynamical models to understand the underlying
mechanisms of insulin signaling, a signaling cascade that leads to the
translocation of glucose, the human body's main source of energy. Fortunately,
reaction network analysis allows us to extract properties of dynamical systems
without depending on their model parameter values. This study focuses on the
comparison of insulin signaling in healthy state (INSMS or INSulin Metabolic
Signaling) and in type 2 diabetes (INRES or INsulin RESistance) using reaction
network analysis. The analysis uses network decomposition to identify the
different subsystems involved in insulin signaling (e.g., insulin receptor
binding and recycling, GLUT4 translocation, and ERK signaling pathway, among
others). Furthermore, results show that INSMS and INRES are similar with
respect to some network, structo-kinetic, and kinetic properties. Their
differences, however, provide insights into what happens when insulin
resistance occurs. First, the variation in the number of species involved in
INSMS and INRES suggests that when irregularities occur in the insulin
signaling pathway, other complexes (and, hence, other processes) get involved,
characterizing insulin resistance. Second, the loss of concordance exhibited by
INRES suggests less restrictive interplay between the species involved in
insulin signaling, leading to unusual activities in the signaling cascade.
Lastly, GLUT4 losing its absolute concentration robustness in INRES may signify
that the transporter has lost its reliability in shuttling glucose to the cell,
inhibiting efficient cellular energy production. This study also suggests
possible applications of the equilibria parametrization and network
decomposition, resulting from the analysis, to potentially establish absolute
concentration robustness in a species.


  Field and mesocosm studies repeatedly show that higher tier process reduce
the predictive accuracy of toxicity evaluation and consequently their value for
pesticide risk assessment. Therefore, understanding the influence of ecological
complexity on toxicant effects is crucial to improve realism of aquatic risk
assessment. Here we investigate the influence of repeated exposure to
ecologically realistic concentrations of esfenvalerate on the similarly
sensitive species Daphnia magna and Culex pipiens in a food limited and highly
competitive environment. We show that significant perturbations in population
development are only present close to the EC50. In contrast, interspecific
competition between species is already reduced at concentrations 3-4 orders of
magnitude below the acute EC50. We conclude that extremely low, environmentally
relevant concentrations can disrupt species interactions. This toxicant
mediated alteration of competitive balances in ecological communities may be
the underlying mechanism for shifts in species distribution at ultra-low
pesticide concentrations. A realistic risk assessment should therefore consider
these processes in order to predict potential pesticide effects on the
structure of communities.


  Phylogenetic inference, grounded in molecular evolution models, is essential
for understanding the evolutionary relationships in biological data. Accounting
for the uncertainty of phylogenetic tree variables, which include tree
topologies and evolutionary distances on branches, is crucial for accurately
inferring species relationships from molecular data and tasks requiring
variable marginalization. Variational Bayesian methods are key to developing
scalable, practical models; however, it remains challenging to conduct
phylogenetic inference without restricting the combinatorially vast number of
possible tree topologies. In this work, we introduce a novel, fully
differentiable formulation of phylogenetic inference that leverages a unique
representation of topological distributions in continuous geometric spaces.
Through practical considerations on design spaces and control variates for
gradient estimations, our approach, GeoPhy, enables variational inference
without limiting the topological candidates. In experiments using real
benchmark datasets, GeoPhy significantly outperformed other approximate
Bayesian methods that considered whole topologies.


  The limited availability of annotations in small molecule datasets presents a
challenge to machine learning models. To address this, one common strategy is
to collaborate with additional auxiliary datasets. However, having more data
does not always guarantee improvements. Negative transfer can occur when the
knowledge in the target dataset differs or contradicts that of the auxiliary
molecule datasets. In light of this, identifying the auxiliary molecule
datasets that can benefit the target dataset when jointly trained remains a
critical and unresolved problem. Through an empirical analysis, we observe that
combining graph structure similarity and task similarity can serve as a more
reliable indicator for identifying high-affinity auxiliary datasets. Motivated
by this insight, we propose MolGroup, which separates the dataset affinity into
task and structure affinity to predict the potential benefits of each auxiliary
molecule dataset. MolGroup achieves this by utilizing a routing mechanism
optimized through a bi-level optimization framework. Empowered by the meta
gradient, the routing mechanism is optimized toward maximizing the target
dataset's performance and quantifies the affinity as the gating score. As a
result, MolGroup is capable of predicting the optimal combination of auxiliary
datasets for each target dataset. Our extensive experiments demonstrate the
efficiency and effectiveness of MolGroup, showing an average improvement of
4.41%/3.47% for GIN/Graphormer trained with the group of molecule datasets
selected by MolGroup on 11 target molecule datasets.


  Type 2 diabetes mellitus is a global pandemic, it is a chronic, progressive
and an incompletely understood metabolic condition. The disease is
characterized by higher levels of sugar in blood caused either due to
insufficient production of insulin or because of insulin resistance. Major
drugs used for the treatment of the condition are fraught with side effects.
Hence, it is becoming an obligation to gaze at alternative agents showing
marginal adverse effects. An important source of such agents are the medicinal
plants. Several plants have been positively identified to show anti-diabetic
effects. The species Salacia reticulata Wight., belonging to the family
Celastraceae which is found in the forests of Southern India is one such
promising plant to tackle type 2 diabetics. In this study, numerous volatile
compounds were identified from various tissues through GC-MS analysis. Among
these the compounds possessed suitable ADMET properties, and high binding
affinities were compared with two approved {\alpha}-glucosidase inhibitors,
Acarbose and Miglitol. The analysis indicated that the compounds with PubChem
IDs, CID-240051 and CID-533471 exhibited potential as inhibitors of Human
Maltase-Glucoamylase enzyme.


  Many reaction networks arising in applications are multistationary, that is,
they have the capacity for more than one steady state; while some networks
exhibit absolute concentration robustness (ACR), which means that some species
concentration is the same at all steady states. Both multistationarity and ACR
are significant in biological settings, but only recently has attention focused
on the possibility for these properties to coexist. Our main result states that
such coexistence in at-most-bimolecular networks (which encompass most networks
arising in biology) requires at least $3$ species, $5$ complexes, and $3$
reactions. We prove additional bounds on the number of reactions for general
networks based on the number of linear conservation laws. Finally, we prove
that, outside of a few exceptional cases, ACR is equivalent to
non-multistationarity for bimolecular networks that are small (more precisely,
one-dimensional or up to two species). Our proofs involve analyses of systems
of sparse polynomials, and we also use classical results from chemical reaction
network theory.


  The online estimation of the derivative of an input signal is widespread in
control theory and engineering. In the realm of chemical reaction networks
(CRN), this raises however a number of specific issues on the different ways to
achieve it. A CRN pattern for implementing a derivative block has already been
proposed for the PID control of biochemical processes, and proved correct using
Tikhonov's limit theorem. In this paper, we give a detailed mathematical
analysis of that CRN, thus clarifying the computed quantity and quantifying the
error done as a function of the reaction kinetic parameters. In a synthetic
biology perspective, we show how this can be used to design error correcting
terms to compute online functions involving derivatives with CRNs. In the
systems biology perspective, we give the list of models in BioModels containing
(in the sense of subgraph epimorphisms) the core derivative CRN, most of which
being models of oscillators and control systems in the cell, and discuss in
detail two such examples: one model of the circadian clock and one model of a
bistable switch.


  We investigate location problems whose optimum lies in the tropical convex
hull of the input points. Firstly, we study geodesically star-convex sets under
the asymmetric tropical distance and introduce the class of tropically
quasiconvex functions whose sub-level sets have this shape. The latter are
related to monotonic functions. Then we show that location problems whose
distances are measured by tropically quasiconvex functions as before give an
optimum in the tropical convex hull of the input points. We also show that a
similar result holds if we replace the input points by tropically convex sets.
Finally, we focus on applications to phylogenetics presenting properties of
consensus methods arising from our class of location problems.


  Sequence Alignment is the process of aligning biological sequences in order
to identify similarities between multiple sequences. In this paper, a Quantum
Algorithm for finding the optimal alignment between DNA sequences has been
demonstrated which works by mapping the sequence alignment problem into a
path-searching problem through a 2D graph. The transition, which converges to a
fixed path on the graph, is based on a proposed oracle for profit calculation.
By implementing Grover's search algorithm, our proposed approach is able to
align a pair of sequences and figure out the optimal alignment within linear
time, which hasn't been attained by any classical deterministic algorithm. In
addition to that, the proposed algorithm is capable of quadratic speeding up to
any unstructured search problem by finding out the optimal paths accurately in
a deterministic manner, in contrast to existing randomized algorithms that
frequently sort out the sub-optimal alignments, therefore, don't always
guarantee of finding out the optimal solutions.


  Consistency and reliability are crucial for conducting AI research. Many
famous research fields, such as object detection, have been compared and
validated with solid benchmark frameworks. After AlphaFold2, the protein
folding task has entered a new phase, and many methods are proposed based on
the component of AlphaFold2. The importance of a unified research framework in
protein folding contains implementations and benchmarks to consistently and
fairly compare various approaches. To achieve this, we present Solvent, a
protein folding framework that supports significant components of
state-of-the-art models in the manner of an off-the-shelf interface Solvent
contains different models implemented in a unified codebase and supports
training and evaluation for defined models on the same dataset. We benchmark
well-known algorithms and their components and provide experiments that give
helpful insights into the protein structure modeling field. We hope that
Solvent will increase the reliability and consistency of proposed models and
give efficiency in both speed and costs, resulting in acceleration on protein
folding modeling research. The code is available at
https://github.com/kakaobrain/solvent, and the project will continue to be
developed.


  When a fraction of a population becomes immune to an infectious disease, the
risk of widespread outbreaks decreases non-linearly as a result of the
collective protection known as herd immunity. Typically, immunity to a disease
can be acquired through natural infection or vaccination. It has been argued
that natural infection in a heterogeneous population may induce a stronger herd
immunity effect than homogeneous vaccination, because the early stages of
transmission would primarily affect highly interactive individuals, leading to
more efficient suppression of subsequent infections. In this study, we provide
a comprehensive assessment of the herd immunity effect by analyzing the
behavior of the susceptible-infected-recovered (SIR) model on static contact
networks. We find that the effectiveness of disease-induced herd immunity is
shaped by two competing mechanisms: While it exploits degree heterogeneity in
the contact network to efficiently target highly connected individuals for
immunity, it also introduces mixing heterogeneity between immune and
susceptible individuals, weakening herd immunity. By controlling the degree of
spatial embeddedness of the contact network structure, we find that spatial
structure generally makes herd immunity stronger and also enhances the
advantage of the herd immunity effect induced by random immunization over the
disease-induced one in networks with low degree heterogeneity. These findings
provide valuable insights into the complexity of herd immunity and have
implications for designing effective disease control strategies.


  A very large class of ODE epidemic models (2.2) discussed in this paper
enjoys the property of admitting also an integral renewal formulation, with
respect to an "age of infection kernel" a(t) which has a matrix exponential
form (3.2). We observe first that a very short proof of this fact is available
when there is only one susceptible compartment, and when its associated "new
infections" matrix has rank one. In this case, a(t) normalized to have integral
1, is precisely the probabilistic law which governs the time spent in all the
"infectious states associated to the susceptible compartment", and the
normalization is precisely the basic replacement number. The Laplace transform
(LT) of a(t) is a generalization of the basic replacement number, and its
structure reflects the laws of the times spent in each infectious state.
Subsequently, we show that these facts admit extensions to processes with
several susceptible classes, provided that all of them have a new infections
matrix of rank one. These results reveal that the ODE epidemic models
highlighted below have also interesting probabilistic properties.


  In group foraging situations, the conventional expectation is that increased
food availability would enhance consumption, especially when animals prioritize
maximizing their food intake. This paper challenges this conventional wisdom by
conducting an in-depth game-theoretic analysis of a basic producer-scrounger
model, in which animals must choose between intensive food searching as
producers or moderate searching while relying on group members as scroungers.
Surprisingly, our study reveals that, under certain circumstances, increasing
food availability can amplify the inclination to scrounge to such an extent
that it paradoxically leads to a reduction in animals' food consumption
compared to scenarios with limited food availability. We further illustrate a
similar phenomenon in a model capturing free-riding dynamics among workers in a
company. We demonstrate that, under certain reward mechanisms, enhancing
workers' production capacities can inadvertently trigger a surge in free-riding
behavior, leading to both diminished group productivity and reduced individual
payoffs. Our findings underscore the significance of contextual factors when
comprehending and predicting the impact of resource availability on individual
and collective outcomes.


  In the world of epidemics, the mathematical modeling of disease co-infection
is gaining importance due to its contributions to mathematics and public
health. Because the co-infection may have a double burden on families,
countries, and the universe, understanding its dynamics is paramount. We study
a SEIQR (susceptible-exposed-infectious-quarantined-recovered) deterministic
epidemic model with a single host population and multiple strains (-$c$ and
-$i$) to account for two competitive diseases with quarantine effects. To model
the role of quarantine and isolation efficacy in disease dynamics, we utilize a
linear function. Further, we shed light on the standard endemic threshold and
determine the conditions for extinction or coexistence with and without forming
co-infection. Next, we show the dependence of the criticality based on specific
parameters of the different pathogens. We found that the disease-free
equilibrium (DFE) of the single-strain model always exists and is globally
asymptotically stable (GAS) if $\tilde{\mathcal{R}}_k^q\leq 1$, else, a stable
endemic equilibrium. On top of that, the model has forward bifurcation at
$\tilde{\mathcal{R}}_k^q = 1$. In the case of a two-strain model, the strain
with a large reproduction number outcompetes the one with a smaller
reproduction number. Further, if the co-infected quarantine reproduction number
is less than one, the infections of already infected individuals will die out,
and co-infection will persist in the population otherwise. We note that the
quarantine and isolation of exposed and infected individuals will reduce the
number of secondary cases below one, consequently reducing the disease
complications if the total number of people in the quarantine is at most the
critical value.


  Accurately inferring the root causes of disease from sequencing data can
improve the discovery of novel therapeutic targets. However, existing root
causal inference algorithms require perfectly measured continuous random
variables. Single cell RNA sequencing (scRNA-seq) datasets contain large
numbers of cells but non-negative counts measured by an error prone process. We
therefore introduce an algorithm called Root Causal Inference with Negative
Binomials (RCI-NB) that accounts for count-based measurement error by
separating negative binomial distributions into their gamma and Poisson
components; the gamma distributions form a fully identifiable but latent post
non-linear causal model representing the true RNA expression levels, which we
only observe with Poisson corruption. RCI-NB identifies patient-specific root
causal contributions from scRNA-seq datasets by integrating novel sparse
regression and goodness of fit testing procedures that bypass Poisson
measurement error. Experiments demonstrate significant improvements over
existing alternatives.


  The swift relaxation of the zero-COVID policy in December 2022 led to an
unprecedented surge in Omicron variant infections in China. With the suspension
of mandatory testing, tracking this epidemic outbreak was challenging because
infections were often underrepresented in survey and testing results, which
only involved partial populations. We used large-scale mobile phone data to
estimate daily infections in Beijing from November 2022 to January 2023. We
demonstrated that an individual's location records of mobile phone could be
used to infer his or her infectious status. Then, the derived status of
millions of individuals could be summed to reconstruct the citywide
spatiotemporal dynamics of infections. We found that the infection incidence
peaked on 21 December, and 80.1% of populations had been infected by 14 January
2023 in Beijing. Furthermore, infection dynamics exhibited significant
demographic and spatiotemporal disparities. Our work provides a ubiquitous and
high-coverage data source for monitoring epidemic outbreaks.


  We generalize the model proposed in [Adimy, Babin, Pujo-Menjouet, \emph{SIAM
Journal on Applied Dynamical Systems} (2022)] for prion infection to a network
of neurons. We do so by applying a so-called \emph{multigroup approach} to the
system of Delay Differential Equations (DDEs) proposed in the aforementioned
paper. We derive the classical threshold quantity $\mathcal{R}_0$,
\textit{i.e.} the basic reproduction number, exploiting the fact that the DDEs
of our model qualitatively behave like Ordinary Differential Equations (ODEs)
when evaluated at the Disease Free Equilibrium. We prove analytically that the
disease naturally goes extinct when $\mathcal{R}_0<1$, whereas it persists when
$\mathcal{R}_0>1$. We conclude with some selected numerical simulations of the
system, to illustrate our analytical results.


  Ribozymes, RNA molecules with distinct 3D structures and catalytic activity,
have widespread applications in synthetic biology and therapeutics. However,
relatively little research has focused on leveraging deep learning to enhance
our understanding of ribozymes. This study implements Word2Vec, an unsupervised
learning technique for natural language processing, to learn ribozyme
embeddings. Ribo2Vec was trained on over 9,000 diverse ribozymes, learning to
map sequences to 128 and 256-dimensional vector spaces. Using Ribo2Vec,
sequence embeddings for five classes of ribozymes (hatchet, pistol, hairpin,
hovlinc, and twister sister) were calculated. Principal component analysis
demonstrated the ability of these embeddings to distinguish between ribozyme
classes. Furthermore, a simple SVM classifier trained on ribozyme embeddings
showed promising results in accurately classifying ribozyme types. Our results
suggest that the embedding vectors contained meaningful information about
ribozymes. Interestingly, 256-dimensional embeddings behaved similarly to
128-dimensional embeddings, suggesting that a lower dimension vector space is
generally sufficient to capture ribozyme features. This approach demonstrates
the potential of Word2Vec for bioinformatics, opening new avenues for ribozyme
research. Future research includes using a Transformer-based method to learn
RNA embeddings, which can capture long-range interactions between nucleotides.


  Pre-trained large language models demonstrate potential in extracting
information from DNA sequences, yet adapting to a variety of tasks and data
modalities remains a challenge. To address this, we propose DNAGPT, a
generalized DNA pre-training model trained on over 200 billion base pairs from
all mammals. By enhancing the classic GPT model with a binary classification
task (DNA sequence order), a numerical regression task (guanine-cytosine
content prediction), and a comprehensive token language, DNAGPT can handle
versatile DNA analysis tasks while processing both sequence and numerical data.
Our evaluation of genomic signal and region recognition, mRNA abundance
regression, and artificial genomes generation tasks demonstrates DNAGPT's
superior performance compared to existing models designed for specific
downstream tasks, benefiting from pre-training using the newly designed model
structure.


  Pain is a significant global health issue, and the current treatment options
for pain management have limitations in terms of effectiveness, side effects,
and potential for addiction. There is a pressing need for improved pain
treatments and the development of new drugs. Voltage-gated sodium channels,
particularly Nav1.3, Nav1.7, Nav1.8, and Nav1.9, play a crucial role in
neuronal excitability and are predominantly expressed in the peripheral nervous
system. Targeting these channels may provide a means to treat pain while
minimizing central and cardiac adverse effects. In this study, we construct
protein-protein interaction (PPI) networks based on pain-related sodium
channels and develop a corresponding drug-target interaction (DTI) network to
identify potential lead compounds for pain management. To ensure reliable
machine learning predictions, we carefully select 111 inhibitor datasets from a
pool of over 1,000 targets in the PPI network. We employ three distinct machine
learning algorithms combined with advanced natural language processing
(NLP)-based embeddings, specifically pre-trained transformer and autoencoder
representations. Through a systematic screening process, we evaluate the side
effects and repurposing potential of over 150,000 drug candidates targeting
Nav1.7 and Nav1.8 sodium channels. Additionally, we assess the ADMET
(absorption, distribution, metabolism, excretion, and toxicity) properties of
these candidates to identify leads with near-optimal characteristics. Our
strategy provides an innovative platform for the pharmacological development of
pain treatments, offering the potential for improved efficacy and reduced side
effects.


  Multi-type birth-death processes underlie approaches for inferring
evolutionary dynamics from phylogenetic trees across biological scales, ranging
from deep-time species macroevolution to rapid viral evolution and somatic
cellular proliferation. A limitation of current phylogenetic birth-death models
is that they require restrictive linearity assumptions that yield tractable
likelihoods, but that also preclude interactions between individuals. Many
fundamental evolutionary processes -- such as environmental carrying capacity
or frequency-dependent selection -- entail interactions, and may strongly
influence the dynamics in some systems. Here, we introduce a multi-type
birth-death process in mean-field interaction with an ensemble of replicas of
the focal process. We prove that, under quite general conditions, the
ensemble's stochastically evolving interaction field converges to a
deterministic trajectory in the limit of an infinite ensemble. In this limit,
the replicas effectively decouple, and self-consistent interactions appear as
nonlinearities in the infinitesimal generator of the focal process. We
investigate a special case that is amenable to calculations in the context of a
phylogenetic birth-death model, and is rich enough to model both carrying
capacity and frequency-dependent selection.