  Predicting the docking between proteins and ligands is a crucial and
challenging task for drug discovery. However, traditional docking methods
mainly rely on scoring functions, and deep learning-based docking approaches
usually neglect the 3D spatial information of proteins and ligands, as well as
the graph-level features of ligands, which limits their performance. To address
these limitations, we propose an equivariant transformer neural network for
protein-ligand docking pose prediction. Our approach involves the fusion of
ligand graph-level features by feature processing, followed by the learning of
ligand and protein representations using our proposed TAMformer module.
Additionally, we employ an iterative optimization approach based on the
predicted distance matrix to generate refined ligand poses. The experimental
results on real datasets show that our model can achieve state-of-the-art
performance.


  Discovery of novel protein biomarkers for clinical applications is an active
research field across a manifold of diseases. Despite some successes and
progress, the biomarker development pipeline still frequently ends in failure
as biomarker candidates cannot be validated or translated to immunoassays.
Selection of strong disease biomarker candidates that further constitute
suitable targets for antibody binding in immunoassays is thus important. This
essential selection step can be supported and rationalized using bioinformatics
tools such as protein databases. Here, I present a workflow in the form of
decision trees to computationally investigate biomarker candidates and their
available affinity reagents in depth. This analysis can identify the most
promising biomarker candidates for assay development while minimal time and
effort is required.


  We identify the thermodynamic conditions necessary to observe indefinite
growth in homogeneous open chemical reaction networks (CRNs) satisfying mass
action kinetics. We also characterize the thermodynamic efficiency of growth by
considering the fraction of the chemical work supplied from the surroundings
that is converted into CRN free energy. We find that indefinite growth cannot
arise in CRNs chemostatted by fixing the concentration of some species at
constant values, or in continuous-flow stirred tank reactors. Indefinite growth
requires a constant net influx from the surroundings of at least one species.
In this case, unimolecular CRNs always generate equilibrium linear growth,
i.e., a continuous linear accumulation of species with equilibrium
concentrations and efficiency one. Multimolecular CRNs are necessary to
generate nonequilibrium growth, i.e., the continuous accumulation of species
with nonequilibrium concentrations. Pseudo-unimolecular CRNs - a subclass of
multimolecular CRNs - always generate asymptotic linear growth with zero
efficiency. Our findings demonstrate the importance of the CRN topology and the
chemostatting procedure in determining the dynamics and thermodynamics of
growth.


  The thermodynamic and dynamical conditions necessary to observe indefinite
growth in homogeneous open chemical reaction networks (CRNs) satisfying mass
action kinetics were presented in Srinivas et al. (2023): Unimolecular CRNs can
only accumulate equilibrium concentrations of species while multimolecular CRNs
are needed to produce indefinite growth with nonequilibrium concentrations.
Within multimolecular CRNs, pseudo-unimolecular CRNs produce nonequilibrium
concentrations with zero efficiencies. Nonequilibrium growth with finite
efficiencies requires dynamically nonlinear CRNs. In this paper, we provide a
detailed analysis supporting these results. Mathematical proofs are provided
for growth in unimolecular and pseudo-unimolecular CRNs. For multimolecular
CRNs, four models displaying very distinctive topological properties are
extensively studied, both numerically and partly analytically.


  Zoonotic disease transmission between animals and humans is a growing risk
and the agricultural context acts as a likely point of transition, with
individual heterogeneity acting as an important contributor. Thus,
understanding the dynamics of disease spread in the wildlife-livestock
interface is crucial for mitigating these risks of transmission. Specifically,
the interactions between pigeons and in-door cows at dairy farms can lead to
significant disease transmission and economic losses for farmers; putting
livestock, adjacent human populations, and other wildlife species at risk. In
this paper, we propose a novel spatio-temporal multi-pathogen model with
continuous spatial movement. The model expands on the
Susceptible-Exposed-Infected-Recovered-Dead (SEIRD) framework and accounts for
both within-species and cross-species transmission of pathogens, as well as the
exploration-exploitation movement dynamics of pigeons, which play a critical
role in the spread of infection agents. In addition to model formulation, we
also implement it as an agent-based simulation approach and use empirical field
data to investigate different biologically realistic scenarios, evaluating the
effect of various parameters on the epidemic spread. Namely, in agreement with
theoretical expectations, the model predicts that the heterogeneity of the
pigeons' movement dynamics can drastically affect both the magnitude and
stability of outbreaks. In addition, joint infection by multiple pathogens can
have an interactive effect unobservable in single-pathogen SIR models,
reflecting a non-intuitive inhibition of the outbreak. Our findings highlight
the impact of heterogeneity in host behavior on their pathogens and allow
realistic predictions of outbreak dynamics in the multi-pathogen
wildlife-livestock interface with consequences to zoonotic diseases in various
systems.


  In complex ecosystems such as microbial communities, there is constant
ecological and evolutionary feedback between the residing species and the
environment occurring on concurrent timescales. Species respond and adapt to
their surroundings by modifying their phenotypic traits, which in turn alters
their environment and the resources available. To study this interplay between
ecological and evolutionary mechanisms, we develop a consumer-resource model
that incorporates phenotypic mutations. In the absence of noise, we find that
phase transitions require finely-tuned interaction kernels. Additionally, we
quantify the effects of noise on frequency dependent selection by defining a
time-integrated mutation current, which accounts for the rate at which
mutations and speciation occurs. We find three distinct phases: homogeneous,
patterned, and patterned traveling waves. The last phase represents one way in
which co-evolution of species can happen in a fluctuating environment. Our
results highlight the principal roles that noise and non-reciprocal
interactions between resources and consumers play in phase transitions within
eco-evolutionary systems.


  In the face of rapidly accumulating genomic data, our understanding of the
RNA regulatory code remains incomplete. Recent self-supervised methods in other
domains have demonstrated the ability to learn rules underlying the
data-generating process such as sentence structure in language. Inspired by
this, we extend contrastive learning techniques to genomic data by utilizing
functional similarities between sequences generated through alternative
splicing and gene duplication. Our novel dataset and contrastive objective
enable the learning of generalized RNA isoform representations. We validate
their utility on downstream tasks such as RNA half-life and mean ribosome load
prediction. Our pre-training strategy yields competitive results using linear
probing on both tasks, along with up to a two-fold increase in Pearson
correlation in low-data conditions. Importantly, our exploration of the learned
latent space reveals that our contrastive objective yields semantically
meaningful representations, underscoring its potential as a valuable
initialization technique for RNA property prediction.


  Phylogenetics is a branch of computational biology that studies the
evolutionary relationships among biological entities. Its long history and
numerous applications notwithstanding, inference of phylogenetic trees from
sequence data remains challenging: the high complexity of tree space poses a
significant obstacle for the current combinatorial and probabilistic
techniques. In this paper, we adopt the framework of generative flow networks
(GFlowNets) to tackle two core problems in phylogenetics: parsimony-based and
Bayesian phylogenetic inference. Because GFlowNets are well-suited for sampling
complex combinatorial structures, they are a natural choice for exploring and
sampling from the multimodal posterior distribution over tree topologies and
evolutionary distances. We demonstrate that our amortized posterior sampler,
PhyloGFN, produces diverse and high-quality evolutionary hypotheses on real
benchmark datasets. PhyloGFN is competitive with prior works in marginal
likelihood estimation and achieves a closer fit to the target distribution than
state-of-the-art variational inference methods.


  Background: Despite the consensus that vaccines play an important role in
combating the global spread of infectious diseases, vaccine inequity is still
rampant with deep-seated mentality of self-priority. This study aims to
evaluate the existence and possible outcomes of a more equitable global vaccine
distribution and explore a concrete incentive mechanism that promotes vaccine
equity. Methods: We design a metapopulation epidemiological model that
simultaneously considers global vaccine distribution and human mobility, which
is then calibrated by the number of infections and real-world vaccination
records during COVID-19 pandemic from March 2020 to July 2021. We explore the
possibility of the enlightened self-interest incentive mechanism, i.e.,
improving one's own epidemic outcomes by sharing vaccines with other countries,
by evaluating the number of infections and deaths under various vaccine sharing
strategies using the proposed model. To understand how these strategies affect
the national interests, we distinguish the imported and local cases for further
cost-benefit analyses that rationalize the enlightened self-interest incentive
mechanism behind vaccine sharing. ...


  An important aspect in the development of small molecules as drugs or
agro-chemicals is their systemic availability after intravenous and oral
administration. The prediction of the systemic availability from the chemical
structure of a potential candidate is highly desirable, as it allows to focus
the drug or agrochemical development on compounds with a favorable kinetic
profile. However, such pre-dictions are challenging as the availability is the
result of the complex interplay between molecular properties, biology and
physiology and training data is rare. In this work we improve the hybrid model
developed earlier [1]. We reduce the median fold change error for the total
oral exposure from 2.85 to 2.35 and for intravenous administration from 1.95 to
1.62. This is achieved by training on a larger data set, improving the neural
network architecture as well as the parametrization of mechanistic model.
Further, we extend our approach to predict additional endpoints and to handle
different covariates, like sex and dosage form. In contrast to a pure machine
learning model, our model is able to predict new end points on which it has not
been trained. We demonstrate this feature by predicting the exposure over the
first 24h, while the model has only been trained on the total exposure.


  Culture involves the origination and transmission of ideas, but the
conditions in which culture can emerge and evolve are unclear. We constructed
and studied a highly simplified neural-network model of these processes. In
this model ideas originate by individual learning from the environment and are
transmitted by communication between individuals. Individuals (or "agents")
comprise a single neuron which receives structured data from the environment
via plastic synaptic connections. The data are generated in the simplest
possible way: linear mixing of independently fluctuating sources and the goal
of learning is to unmix the data. To make this problem tractable we assume that
at least one of the sources fluctuates in a nonGaussian manner. Linear mixing
creates structure in the data, and agents attempt to learn (from the data and
possibly from other individuals) synaptic weights that will unmix, i.e., to
"understand" the agent's world. For a variety of reasons even this goal can be
difficult for a single agent to achieve; we studied one particular type of
difficulty (created by imperfection in synaptic plasticity), though our
conclusions should carry over to many other types of difficulty. We previously
studied whether a small population of communicating agents, learning from each
other, could more easily learn unmixing coefficients than isolated individuals,
learning only from their environment. We found, unsurprisingly, that if agents
learn indiscriminately from any other agent (whether or not they have learned
good solutions), communication does not enhance understanding. Here we extend
the model slightly, by allowing successful learners to be more effective
teachers, and find that now a population of agents can learn more effectively
than isolated individuals. We suggest that a key factor in the onset of culture
might be the development of selective learning.


  The relationship between intelligence and evolution is bidirectional: while
evolution can help evolve intelligences, the degree of intelligence itself can
impact evolution (Baldwin, 1896). In the field of Evolutionary Computation, the
inverse relationship (impact of intelligence on evolution) is approached from
the perspective of organism level behaviour (Hinton, 1996). We extend these
ideas to the developmental (cellular morphogenetic) level in the context of an
expanded view of intelligence as not only the ability of a system to navigate
the three-dimensional world, but also as the ability to navigate other
arbitrary spaces (transcriptional, anatomical, physiological, etc.). Here, we
specifically focus on the intelligence of a minimal model of a system
navigating anatomical morphospace, and assess how the degree and manner of
problem solving competency during morphogenesis effects evolutionary dynamics.
To this end, we evolve populations of artificial embryos using a standard
genetic algorithm in silico. Artificial embryos were cellular collectives given
the capacity to undergo morphogenetic rearrangement (e.g., regulative
development) prior to selection within an evolutionary cycle. Results from our
model indicates that morphogenetic competency significantly alters evolutionary
dynamics, with evolution preferring to improve anatomical intelligence rather
than perfect the structural genes. These observations hint that evolution in
the natural world may be leveraging the problem solving competencies of cells
at multiple scales to boost evolvability and robustness to novel conditions. We
discuss implications of our results for the Developmental Biology and
Artificial Life communities.


  Designing flexible probabilistic models over tree topologies is important for
developing efficient phylogenetic inference methods. To do that, previous works
often leverage the similarity of tree topologies via hand-engineered heuristic
features which would require pre-sampled tree topologies and may suffer from
limited approximation capability. In this paper, we propose a deep
autoregressive model for phylogenetic inference based on graph neural networks
(GNNs), called ARTree. By decomposing a tree topology into a sequence of leaf
node addition operations and modeling the involved conditional distributions
based on learnable topological features via GNNs, ARTree can provide a rich
family of distributions over the entire tree topology space that have simple
sampling algorithms and density estimation procedures, without using heuristic
features. We demonstrate the effectiveness and efficiency of our method on a
benchmark of challenging real data tree topology density estimation and
variational Bayesian phylogenetic inference problems.


  Graph Neural Networks (GNN) are reshaping our understanding of biomedicine
and diseases by revealing the deep connections among genes and cells. As both
algorithmic and biomedical technologies have advanced significantly, we're
entering a transformative phase of personalized medicine. While pioneering
tools like Graph Attention Networks (GAT) and Graph Convolutional Neural
Networks (Graph CNN) are advancing graph-based learning, the rise of
single-cell sequencing techniques is reshaping our insights on cellular
diversity and function. Numerous studies have combined GNNs with single-cell
data, showing promising results. In this work, we highlight the GNN
methodologies tailored for single-cell data over the recent years. We outline
the diverse range of graph deep learning architectures that center on GAT
methodologies. Furthermore, we underscore the several objectives of GNN
strategies in single-cell data contexts, ranging from cell-type annotation,
data integration and imputation, gene regulatory network reconstruction,
clustering and many others. This review anticipates a future where GNNs become
central to single-cell analysis efforts, particularly as vast omics datasets
are continuously generated and the interconnectedness of cells and genes
enhances our depth of knowledge in biomedicine.


  The diameter distribution of a given species of deciduous trees in mature,
temperate zone forests is well approximated by a Gamma distribution. Here we
give new experimental evidence for this conjecture by analyzing deciduous tree
size data in mature semi-natural forest and ancient, traditionally managed
wood-pasture from Central Europe. These distribution functions collapse on a
universal shape if the tree sizes are normalized to the mean value in the
considered sample. A novel evolutionary master equation is used to model the
observed distribution. The model incorporates three probabilistic processes:
tree growth, mortality and diversification. By using simple, and realistic
state dependent kernel functions for the growth and reset rates together with
an assumed multiplicative dilution due to diversification, the stationary
solution of the master equation yields the experimentally observed Gamma
distribution. The model as it is formulated allows analytically compact
solution and has only two fitting parameters whose values are consistent with
the experimental data for the growth and reset processes. Our results suggest
also that tree size statistics can be used to infer woodland naturalness.


  Engineering new molecules with desirable functions and properties has the
potential to extend our ability to engineer proteins beyond what nature has so
far evolved. Advances in the so-called "de novo" design problem have recently
been brought forward by developments in artificial intelligence. Generative
architectures, such as language models and diffusion processes, seem adept at
generating novel, yet realistic proteins that display desirable properties and
perform specified functions. State-of-the-art design protocols now achieve
experimental success rates nearing 20%, thus widening the access to de novo
designed proteins. Despite extensive progress, there are clear field-wide
challenges, for example in determining the best in silico metrics to prioritise
designs for experimental testing, and in designing proteins that can undergo
large conformational changes or be regulated by post-translational
modifications and other cellular processes. With an increase in the number of
models being developed, this review provides a framework to understand how
these tools fit into the overall process of de novo protein design. Throughout,
we highlight the power of incorporating biochemical knowledge to improve
performance and interpretability.


  We consider a deterministic mathematical model to study the role of
Information on the dynamics of Typhoid Fever. We analyses the model to study
its boundedness and compute the threshold value known as basic reproduction
number for determination of number of secondary cases and establishment of the
the condition for local and global asymptotic stability of the stationary
points. The numerical simulation is used to depict dynamical behavior of
Typhoid Fever in the considered population. The results indicate a clear role
of information in influencing a behavior change that may in a way lead to an
increase in the transmission typhoid fever. The result further show that the
increase of the number of individuals with Typhoid fever is greatly influenced
by failure of the people to follow health precaution that reduce the spread of
the disease. The result necessitate the importance of the government to
continue educating and/or giving information to her people on the behavior that
in one way or another may lead to the increase of the transmission of typhoid
fever which should also be an importance topic to be discussed when planning
for any control strategies against the disease.


  To investigate1,1-Diphenyl-2-picrylhydrazyl (DPPH) and superoxide radical
(SOR) 17 scavenging activities of 2-oxo-1,2,3,4-tetrahydropyrimidines
derivatives. Free radicals are 18 highly unstable and reactive molecules/atoms.
In the body, free radicals form during 19 normal and abnormal metabolism in the
body and cause serious damage to other 20 biomolecules through generating
oxidative stress (OS). If free radicals induced OS does 21 not neutralized
properly, host multiple pathologies including several types of cancers.


  In the current manuscript, a first two-patch model with Allee effect and
nonlinear dispersal is presented. We study both the ODE case and the PDE case
here. In the ODE model, the stability of the equilibrium points and the
existence of saddle-node bifurcation are discussed. The phase diagram and
bifurcation curve of our model are also given by numerical simulation. Besides,
the corresponding linear dispersal case is also presented. We show that when
the Allee effect is large, high intensity of linear dispersal is not favorable
to the persistence of the species. We further show when the Allee effect is
large, nonlinear diffusion is more favorable to the survival of the population
than linear diffusion. Moreover, the results of the PDE model extends our
findings from discrete patches to continuous patches.


  Through evolution, nature has presented a set of remarkable protein
materials, including elastins, silks, keratins and collagens with superior
mechanical performances that play crucial roles in mechanobiology. However,
going beyond natural designs to discover proteins that meet specified
mechanical properties remains challenging. Here we report a generative model
that predicts protein designs to meet complex nonlinear mechanical
property-design objectives. Our model leverages deep knowledge on protein
sequences from a pre-trained protein language model and maps mechanical
unfolding responses to create novel proteins. Via full-atom molecular
simulations for direct validation, we demonstrate that the designed proteins
are novel, and fulfill the targeted mechanical properties, including unfolding
energy and mechanical strength, as well as the detailed unfolding
force-separation curves. Our model offers rapid pathways to explore the
enormous mechanobiological protein sequence space unconstrained by biological
synthesis, using mechanical features as target to enable the discovery of
protein materials with superior mechanical properties.


  Traditional molecular string representations, such as SMILES, often pose
challenges for AI-driven molecular design due to their non-sequential depiction
of molecular substructures. To address this issue, we introduce Sequential
Attachment-based Fragment Embedding (SAFE), a novel line notation for chemical
structures. SAFE reimagines SMILES strings as an unordered sequence of
interconnected fragment blocks while maintaining compatibility with existing
SMILES parsers. It streamlines complex generative tasks, including scaffold
decoration, fragment linking, polymer generation, and scaffold hopping, while
facilitating autoregressive generation for fragment-constrained design, thereby
eliminating the need for intricate decoding or graph-based models. We
demonstrate the effectiveness of SAFE by training an 87-million-parameter
GPT2-like model on a dataset containing 1.1 billion SAFE representations.
Through targeted experimentation, we show that our SAFE-GPT model exhibits
versatile and robust optimization performance. SAFE opens up new avenues for
the rapid exploration of chemical space under various constraints, promising
breakthroughs in AI-driven molecular design.


  Symmetry principles have proven important in physics, deep learning and
geometry, allowing for the reduction of complicated systems to simpler, more
comprehensible models that preserve the system's features of interest.
Biological systems often show a high level of complexity and consist of a high
number of interacting parts. Using symmetry fibrations, the relevant symmetries
for biological 'message-passing' networks, we reduced the gene regulatory
networks of E. coli and B. subtilis bacteria in a way that preserves
information flow and highlights the computational capabilities of the network.
Nodes that share isomorphic input trees are grouped into equivalence classes
called fibers, whereby genes that receive signals with the same 'history'
belong to one fiber and synchronize. We further reduce the networks to its
computational core by removing 'dangling ends' via k-core decomposition. The
computational core of the network consists of a few strongly connected
components in which signals can cycle while signals are transmitted between
these 'information vortices' in a linear feed-forward manner. These components
are in charge of decision making in the bacterial cell by employing a series of
genetic toggle-switch circuits that store memory, and oscillator circuits.
These circuits act as the central computation machine of the network, whose
output signals then spread to the rest of the network.


  We introduce a multi-scale model for affinity maturation, which aims to
capture the intra-clonal, inter-clonal and epitope-specific organization of the
B cell population in a germinal center. We apply the Eigen's model of
quasispecies evolution in a random fitness landscape, where the desired
multi-scale structure is reflected in the mutational connectivity of the B cell
receptor (BCR) space, and on statistical characteristics of the fitness
ensembles. Within this mathematical framework, we study the competition among
classes of BCRs targeting different antigen epitopes, and construct an
effective immunogenic space where epitope immunodominance relations can be
universally characterized. We study how varying the composition of a mixture of
antigens with variable and conserved domains allows for a parametric
exploration of this space, and identify general principles for the rational
design of two-antigen cocktails. We find that optimization of the cocktail
composition must take into account possible statistical asymmetries between the
two classes of BCRs targeting the two variants of the mutated epitope, as well
as the evolutionary regime in which affinity maturation occurs.


  In the realm of pandemic dynamics, understanding the intricate interplay
between disease transmission, interventions, and immunity is pivotal for
effective control strategies. Through a rigorous agent-based simulation, we
embarked on a comprehensive exploration, traversing unmitigated spread,
lockdown scenarios, and the transformative potential of vaccination. we unveil
a paradoxical trend: while quarantine unquestionably delays the pandemic peak,
it does not act as an impenetrable barrier to halt the progression of
infectious diseases. Vaccination factor revealed a potent weapon against
outbreaks higher vaccination percentage not only delayed infection peaks but
also substantially curtailed their impact. Our investigation into bond dilution
below the percolation threshold presents an additional dimension to pandemic
control. We observed that localized isolation through bond dilution offers a
targeted control strategy that can be more resource-efficient compared to
blanket lockdowns or large-scale vaccination campaigns.


  When biological populations expand into new territory, the evolutionary
outcomes can be strongly influenced by genetic drift, the random fluctuations
in allele frequencies. Meanwhile, spatial variability in the environment can
also significantly influence the competition between sub-populations vying for
space. Little is known about the interplay of these intrinsic and extrinsic
sources of noise in population dynamics: When does environmental heterogeneity
dominate over genetic drift or vice versa, and what distinguishes their
population genetics signatures? Here, in the context of neutral evolution, we
examine the interplay between a population's intrinsic, demographic noise and
an extrinsic, quenched-random noise provided by a heterogeneous environment.
Using a multi-species Eden model, we simulate a population expanding over a
landscape with random variations in local growth rates and measure how this
variability affects genealogical tree structure, and thus genetic diversity. We
find that, when the heterogeneity is sufficiently strong, the population front
is dominated by genealogical lineages that are pinned to a small number of
optimal paths. The landscape-dependent statistics of these optimal paths then
supersede those of the population's intrinsic noise as the main determinant of
evolutionary dynamics. Remarkably, the statistics for coalescence of
genealogical lineages, derived from those deterministic paths, strongly
resemble the statistics emerging from demographic noise alone in uniform
landscapes. This cautions interpretations of coalescence statistics and raises
new challenges for inferring past population dynamics.


  Innovations like protein diffusion have enabled significant progress in de
novo protein design, which is a vital topic in life science. These methods
typically depend on protein structure encoders to model residue backbone
frames, where atoms do not exist. Most prior encoders rely on atom-wise
features, such as angles and distances between atoms, which are not available
in this context. Thus far, only several simple encoders, such as IPA, have been
proposed for this scenario, exposing the frame modeling as a bottleneck. In
this work, we proffer the Vector Field Network (VFN), which enables network
layers to perform learnable vector computations between coordinates of
frame-anchored virtual atoms, thus achieving a higher capability for modeling
frames. The vector computation operates in a manner similar to a linear layer,
with each input channel receiving 3D virtual atom coordinates instead of scalar
values. The multiple feature vectors output by the vector computation are then
used to update the residue representations and virtual atom coordinates via
attention aggregation. Remarkably, VFN also excels in modeling both frames and
atoms, as the real atoms can be treated as the virtual atoms for modeling,
positioning VFN as a potential universal encoder. In protein diffusion (frame
modeling), VFN exhibits an impressive performance advantage over IPA, excelling
in terms of both designability (67.04% vs. 53.58%) and diversity (66.54% vs.
51.98%). In inverse folding (frame and atom modeling), VFN outperforms the
previous SoTA model, PiFold (54.7% vs. 51.66%), on sequence recovery rate. We
also propose a method of equipping VFN with the ESM model, which significantly
surpasses the previous ESM-based SoTA (62.67% vs. 55.65%), LM-Design, by a
substantial margin.


  In a predator-prey metapopulation, the two traits are adversely related:
synchronization and persistence. A decrease in synchrony apparently leads to an
increase in persistence and, therefore, necessitates the study of desynchrony
in a metapopulation. In this article, we study predator-prey patches that
communicate with one another while being interconnected through distinct
dispersal structures in the layers of a three-layer multiplex network. We
investigate the synchronization phenomenon among the patches of the outer
layers by introducing higher-order interactions (specifically three-body
interactions) in the middle layer. We observe a decrease in the synchronous
behavior or, alternatively, an increase in desynchrony due to the inclusion of
group interactions among the patches of the middle layer. The advancement of
desynchrony becomes more prominent with increasing strength and numbers of
three-way interactions in the middle layer. We analytically validated our
numerical results by performing the stability analysis of the referred
synchronous solution using the master stability function approach.
Additionally, we verify our findings by taking into account two distinct
predator-prey models and dispersal topologies, which ultimately assert that the
findings are generalizable across various models and dispersal structures.


  In the presence of recombination, the evolutionary relationships between a
set of sampled genomes cannot be described by a single genealogical tree.
Instead, the genomes are related by a complex, interwoven collection of
genealogies formalized in a structure called an ancestral recombination graph
(ARG). An ARG extensively encodes the ancestry of the genome(s) and thus is
replete with valuable information for addressing diverse questions in
evolutionary biology. Despite its potential utility, technological and
methodological limitations, along with a lack of approachable literature, have
severely restricted awareness and application of ARGs in empirical evolution
research. Excitingly, recent progress in ARG reconstruction and simulation have
made ARG-based approaches feasible for many questions and systems. In this
review, we provide an accessible introduction and exploration of ARGs, survey
recent methodological breakthroughs, and describe the potential for ARGs to
further existing goals and open avenues of inquiry that were previously
inaccessible in evolutionary genomics. Through this discussion, we aim to more
widely disseminate the promise of ARGs in evolutionary genomics and encourage
the broader development and adoption of ARG-based inference.


  This study embarks on a comprehensive exploration of the decadal trends and
future trajectories of twenty distinct infectious diseases in China from 1998
to 2021. A refined Hybrid Bayesian Structural Time Series (BSTS)-Markov Chain
Monte Carlo (MCMC) model is employed, intertwining with Long Short-Term Memory
(LSTM) networks to dissect intricate relationships amidst population
demographics, economic indices, and the evolution of infectious diseases. The
findings reveal the persistent prevalence of high incidence diseases in future
10 years, like AIDS, Gonorrhea, and Syphilis, and stable occurrences of middle
incidence rate diseases such as Brucellosis and Scarlet Fever, while also
foretelling the potential disappearance of lower incidence rate diseases like
Cholera, Encephalitis B, and Measles. The study particularly underscores the
transformative impact of the COVID-19 pandemic, showcasing its extensive
implications on the incidences and management of a plethora of diseases, urging
a deeper probe into the nuanced alterations in disease transmission, testing,
and reporting modalities amidst global health crises. This research accentuates
the critical role of advanced predictive analytics in fostering global
preparedness and response mechanisms, and in fortifying the resilience and
adaptability of China public health framework against burgeoning infectious
disease threats.


  The correlation among the gene genealogies at different loci is crucial in
biology, yet challenging to understand because such correlation depends on many
factors including genetic linkage, recombination, natural selection and
population structure. Based on a diploid Wright-Fisher model with a single
mating type and partial selfing for a constant large population with size $N$,
we quantify the combined effect of genetic drift and two competing factors,
recombination and selfing, on the correlation of coalescence times at two
linked loci for samples of size two. Recombination decouples the genealogies at
different loci and decreases the correlation while selfing increases the
correlation. We obtain explicit asymptotic formulas for the correlation for
four scaling scenarios that depend on whether the selfing probability and the
recombination probability are of order $O(1/N)$ or $O(1)$ as $N$ tends to
infinity. Our analytical results confirm that the asymptotic lower bound in
[King, Wakeley, Carmi (Theor. Popul. Biol. 2018)] is sharp when the loci are
unlinked and when there is no selfing, and provide a number of new formulas for
other scaling scenarios that have not been considered before. We present
asymptotic results for the variance of Tajima's estimator of the population
mutation rate for infinitely many loci as $N$ tends to infinity. When the
selfing probability is of order $O(1)$ and is equal to a positive constant $s$
for all $N$ and if the samples at both loci are in the same individual, then
the variance of the Tajima's estimator tends to $s/2$ (hence remains positive)
even when the recombination rate, the number of loci and the population size
all tend to infinity.


  Indirect reciprocity is a mechanism that explains large-scale cooperation in
human societies. In indirect reciprocity, an individual chooses whether or not
to cooperate with another based on reputation information, and others evaluate
the action as good or bad. Under what evaluation rule (called ``social norm'')
cooperation evolves has long been of central interest in the literature. It has
been reported that if individuals can share their evaluations (i.e., public
reputation), social norms called ``leading eight'' can be evolutionarily
stable. On the other hand, when they cannot share their evaluations (i.e.,
private assessment), the evolutionary stability of cooperation is still in
question. To tackle this problem, we create a novel method to analyze the
reputation structure in the population under private assessment. Specifically,
we characterize each individual by two variables, ``goodness'' (what proportion
of the population considers the individual as good) and ``self-reputation''
(whether an individual thinks of him/herself as good or bad), and analyze the
stochastic process of how these two variables change over time. We discuss
evolutionary stability of each of the leading eight social norms by studying
the robustness against invasions of unconditional cooperators and defectors. We
identify key pivots in those social norms for establishing a high level of
cooperation or stable cooperation against mutants. Our finding gives an insight
into how human cooperation is established in a real-world society.


  Generative probabilistic models emerge as a new paradigm in data-driven,
evolution-informed design of biomolecular sequences. This paper introduces a
novel approach, called Edge Activation Direct Coupling Analysis (eaDCA),
tailored to the characteristics of RNA sequences, with a strong emphasis on
simplicity, efficiency, and interpretability. eaDCA explicitly constructs
sparse coevolutionary models for RNA families, achieving performance levels
comparable to more complex methods while utilizing a significantly lower number
of parameters. Our approach demonstrates efficiency in generating artificial
RNA sequences that closely resemble their natural counterparts in both
statistical analyses and SHAPE-MaP experiments, and in predicting the effect of
mutations. Notably, eaDCA provides a unique feature: estimating the number of
potential functional sequences within a given RNA family. For example, in the
case of cyclic di-AMP riboswitches (RF00379), our analysis suggests the
existence of approximately $\mathbf{10^{39}}$ functional nucleotide sequences.
While huge compared to the known $< \mathbf{4,000}$ natural sequences, this
number represents only a tiny fraction of the vast pool of nearly
$\mathbf{10^{82}}$ possible nucleotide sequences of the same length (136
nucleotides). These results underscore the promise of sparse and interpretable
generative models, such as eaDCA, in enhancing our understanding of the
expansive RNA sequence space.


  Interacting biological systems at all organizational levels display emergent
behavior. Modeling these systems is made challenging by the number and variety
of biological components and interactions (from molecules in gene regulatory
networks to species in ecological networks) and the often-incomplete state of
system knowledge (e.g., the unknown values of kinetic parameters for
biochemical reactions). Boolean networks have emerged as a powerful tool for
modeling these systems. We provide a methodological overview of Boolean network
models of biological systems. After a brief introduction, we describe the
process of building, analyzing, and validating a Boolean model. We then present
the use of the model to make predictions about the system's response to
perturbations and about how to control (or at least influence) its behavior. We
emphasize the interplay between structural and dynamical properties of Boolean
networks and illustrate them in three case studies from disparate levels of
biological organization.


  Measuring diversity accurately is important for many scientific fields,
including machine learning (ML), ecology, and chemistry. The Vendi Score was
introduced as a generic similarity-based diversity metric that extends the Hill
number of order q=1 by leveraging ideas from quantum statistical mechanics.
Contrary to many diversity metrics in ecology, the Vendi Score accounts for
similarity and does not require knowledge of the prevalence of the categories
in the collection to be evaluated for diversity. However, the Vendi Score
treats each item in a given collection with a level of sensitivity proportional
to the item's prevalence. This is undesirable in settings where there is a
significant imbalance in item prevalence. In this paper, we extend the other
Hill numbers using similarity to provide flexibility in allocating sensitivity
to rare or common items. This leads to a family of diversity metrics -- Vendi
scores with different levels of sensitivity -- that can be used in a variety of
applications. We study the properties of the scores in a synthetic controlled
setting where the ground truth diversity is known. We then test their utility
in improving molecular simulations via Vendi Sampling. Finally, we use the
Vendi scores to better understand the behavior of image generative models in
terms of memorization, duplication, diversity, and sample quality.


  Stabilizing proteins is a foundational step in protein engineering. However,
the evolutionary pressure of all extant proteins makes identifying the scarce
number of mutations that will improve thermodynamic stability challenging. Deep
learning has recently emerged as a powerful tool for identifying promising
mutations. Existing approaches, however, are computationally expensive, as the
number of model inferences scales with the number of mutations queried. Our
main contribution is a simple, parallel decoding algorithm. Our Mutate
Everything is capable of predicting the effect of all single and double
mutations in one forward pass. It is even versatile enough to predict
higher-order mutations with minimal computational overhead. We build Mutate
Everything on top of ESM2 and AlphaFold, neither of which were trained to
predict thermodynamic stability. We trained on the Mega-Scale cDNA proteolysis
dataset and achieved state-of-the-art performance on single and higher-order
mutations on S669, ProTherm, and ProteinGym datasets. Code is available at
https://github.com/jozhang97/MutateEverything


  Conventional deep learning methods typically employ supervised learning for
drug response prediction (DRP). This entails dependence on labeled response
data from drugs for model training. However, practical applications in the
preclinical drug screening phase demand that DRP models predict responses for
novel compounds, often with unknown drug responses. This presents a challenge,
rendering supervised deep learning methods unsuitable for such scenarios. In
this paper, we propose a zero-shot learning solution for the DRP task in
preclinical drug screening. Specifically, we propose a Multi-branch
Multi-Source Domain Adaptation Test Enhancement Plug-in, called MSDA. MSDA can
be seamlessly integrated with conventional DRP methods, learning invariant
features from the prior response data of similar drugs to enhance real-time
predictions of unlabeled compounds. We conducted experiments using the GDSCv2
and CellMiner datasets. The results demonstrate that MSDA efficiently predicts
drug responses for novel compounds, leading to a general performance
improvement of 5-10\% in the preclinical drug screening phase. The significance
of this solution resides in its potential to accelerate the drug discovery
process, improve drug candidate assessment, and facilitate the success of drug
discovery.


  Mathematical models play a crucial role in understanding the spread of
infectious disease outbreaks and influencing policy decisions. These models aid
pandemic preparedness by predicting outcomes under hypothetical scenarios and
identifying weaknesses in existing frameworks. However, their accuracy,
utility, and comparability are being scrutinized. Agent-based models (ABMs)
have emerged as a valuable tool, capturing population heterogeneity and spatial
effects, particularly when assessing intervention strategies. Here we present
EpiGeoPop, a user-friendly tool for rapidly preparing spatially accurate
population configurations of entire countries. EpiGeoPop helps to address the
problem of complex and time-consuming model set up in ABMs, specifically
improving the integration of spatial detail. We subsequently demonstrate the
importance of accurate spatial detail in ABM simulations of disease outbreaks
using Epiabm, an ABM based on Imperial College London's CovidSim with improved
modularity, documentation and testing. Our investigation involves the interplay
between population density, the implementation of spatial transmission, and
realistic interventions implemented in Epiabm.


  Cell-free transcription and translation (TXTL) systems have emerged as a
powerful tool for testing genetic regulatory elements and circuits. Cell-free
prototyping can dramatically accelerate the design-build-test cycle of new
functions in synthetic biology, in particular when linear DNA templates are
used. Here we describe a Golden Gate assisted workflow to rapidly produce
linear DNA templates for TXTL reactions by assembling transcriptional units
from basic genetic parts of a modular cloning toolbox. Functional DNA templates
composed of multiple parts such as promoter, ribosomal binding site (RBS),
coding sequence, and terminator are produced in vitro in a one-pot Golden Gate
assembly reaction followed by PCR amplification. By eliminating lengthy
transformation and cloning steps in cells and by taking advantage of modular
cloning toolboxes, our cell-free prototyping workflow can produce data for
large numbers of new constructs within a single day.


  The microtubule cytoskeleton is responsible for sustained, long-range
intracellular transport of mRNAs, proteins, and organelles in neurons. Neuronal
microtubules must be stable enough to ensure reliable transport, but they also
undergo dynamic instability, as their plus and minus ends continuously switch
between growth and shrinking. This process allows for continuous rebuilding of
the cytoskeleton and for flexibility in injury settings. Motivated by
\textit{in vivo} experimental data on microtubule behavior in
\textit{Drosophila} neurons, we propose a mathematical model of dendritic
microtubule dynamics, with a focus on understanding microtubule length,
velocity, and state-duration distributions. We find that limitations on
microtubule growth phases are needed for realistic dynamics, but the type of
limiting mechanism leads to qualitatively different responses to plausible
experimental perturbations. We therefore propose and investigate two
minimally-complex length-limiting factors: limitation due to resource (tubulin)
constraints and limitation due to catastrophe of large-length microtubules. We
combine simulations of a detailed stochastic model with steady-state analysis
of a mean-field ordinary differential equations model to map out qualitatively
distinct parameter regimes. This provides a basis for predicting changes in
microtubule dynamics, tubulin allocation, and the turnover rate of tubulin
within microtubules in different experimental environments. Ultimately, this
work provides a tunable and statistically identifiable framework for studying
the emergent properties of dynamic instability of microtubules.


  The processes of gene expression are inherently stochastic, even for
essential genes required for growth. How does the cell maximize fitness in
light of noise? To answer this question, we build a mathematical model to
explore the trade-off between metabolic load and growth robustness. The model
predicts novel principles of central dogma regulation: Optimal protein
expression levels are vastly overabundant. Essential genes are transcribed
above a lower limit of one message per cell cycle. Gene expression is achieved
by load balancing between transcription and translation. We show that each of
these novel regulatory principles is observed. These results reveal that
robustness and metabolic load determine the global regulatory principles that
govern central dogma function, and these principles have broad implications for
cellular function.


  An amino acid insertion or deletion, or InDel, can have profound and varying
functional impacts on a protein's structure. InDel mutations in the
transmembrane conductor regulator protein for example give rise to cystic
fibrosis. Unfortunately performing InDel mutations on physical proteins and
studying their effects is a time prohibitive process. Consequently, modeling
InDels computationally can supplement and inform wet lab experiments. In this
work, we make use of our data sets of exhaustive double InDel mutations for
three proteins which we computationally generated using a robotics inspired
inverse kinematics approach available in Rosetta. We develop and train a neural
network, RoseNet, on several structural and energetic metrics output by Rosetta
during the mutant generation process. We explore and present how RoseNet is
able to emulate the exhaustive data set using deep learning methods, and show
to what extent it can predict Rosetta metrics for unseen mutant sequences with
two InDels. RoseNet achieves a Pearson correlation coefficient median accuracy
of 0.775 over all Rosetta scores for the largest protein. Furthermore, a
sensitivity analysis is performed to determine the necessary quantity of data
required to accurately emulate the structural scores for computationally
generated mutants. We show that the model can be trained on minimal data (<50%)
and still retain a high level of accuracy.


  This brief note highlights a largely overlooked similarity between the SIR
ordinary differential equations used for epidemics on the configuration model
of a Poisson network and the classical mass-action SIR equations introduced
nearly a century ago by Kermack and McKendrick. We demonstrate that the decline
pattern in susceptibles is identical for both models. This equivalence carries
practical implications: the susceptibles decay curve, often referred to as the
epidemic or incidence curve, is frequently used in empirical studies to
forecast epidemic dynamics. Although the curves for susceptibles align
perfectly, those for infections do differ. Yet, the infection curves tend to
converge and become almost indistinguishable in high-degree networks. In
summary, our analysis suggests that under many practical scenarios, it is
acceptable to use the classical SIR model as a close approximation to the
Poisson SIR network model.


  The applicability of computational and dynamical systems models to organisms
is scrutinized, using examples from developmental biology and cognition.
Developmental morphogenesis is dependent on the inherent material properties of
developing tissues, a non-computational modality, but cell differentiation,
which utilizes chromatin-based revisable memory banks and program-like
function-calling, via the developmental gene co-expression system unique to
metazoans, has a quasi-computational basis. Multi-attractor dynamical models
are argued to be misapplied to global properties of development, and it is
suggested that along with computationalism, dynamicism is similarly unsuitable
to accounting for cognitive phenomena. Proposals are made for treating brains
and other nervous tissues as novel forms of excitable matter with inherent
properties which enable the intensification of cell-based basal cognition
capabilities present throughout the tree of life.


  Protein-ligand structure prediction is an essential task in drug discovery,
predicting the binding interactions between small molecules (ligands) and
target proteins (receptors). Although conventional physics-based docking tools
are widely utilized, their accuracy is compromised by limited conformational
sampling and imprecise scoring functions. Recent advances have incorporated
deep learning techniques to improve the accuracy of structure prediction.
Nevertheless, the experimental validation of docking conformations remains
costly, it raises concerns regarding the generalizability of these deep
learning-based methods due to the limited training data. In this work, we show
that by pre-training a geometry-aware SE(3)-Equivariant neural network on a
large-scale docking conformation generated by traditional physics-based docking
tools and then fine-tuning with a limited set of experimentally validated
receptor-ligand complexes, we can achieve outstanding performance. This process
involved the generation of 100 million docking conformations, consuming roughly
1 million CPU core days. The proposed model, HelixDock, aims to acquire the
physical knowledge encapsulated by the physics-based docking tools during the
pre-training phase. HelixDock has been benchmarked against both physics-based
and deep learning-based baselines, showing that it outperforms its closest
competitor by over 40% for RMSD. HelixDock also exhibits enhanced performance
on a dataset that poses a greater challenge, thereby highlighting its
robustness. Moreover, our investigation reveals the scaling laws governing
pre-trained structure prediction models, indicating a consistent enhancement in
performance with increases in model parameters and pre-training data. This
study illuminates the strategic advantage of leveraging a vast and varied
repository of generated data to advance the frontiers of AI-driven drug
discovery.


  Molecular representation learning is fundamental for many drug related
applications. Most existing molecular pre-training models are limited in using
single molecular modality, either SMILES or graph representation. To
effectively leverage both modalities, we argue that it is critical to capture
the fine-grained 'semantics' between SMILES and graph, because subtle
sequence/graph differences may lead to contrary molecular properties. In this
paper, we propose a universal SMILE-graph representation learning model, namely
UniMAP. Firstly, an embedding layer is employed to obtain the token and
node/edge representation in SMILES and graph, respectively. A multi-layer
Transformer is then utilized to conduct deep cross-modality fusion. Specially,
four kinds of pre-training tasks are designed for UniMAP, including Multi-Level
Cross-Modality Masking (CMM), SMILES-Graph Matching (SGM), Fragment-Level
Alignment (FLA), and Domain Knowledge Learning (DKL). In this way, both global
(i.e. SGM and DKL) and local (i.e. CMM and FLA) alignments are integrated to
achieve comprehensive cross-modality fusion. We evaluate UniMAP on various
downstream tasks, i.e. molecular property prediction, drug-target affinity
prediction and drug-drug interaction. Experimental results show that UniMAP
outperforms current state-of-the-art pre-training methods.We also visualize the
learned representations to demonstrate the effect of multi-modality
integration.


  Species' interactions are shaped by their traits. Thus, we expect traits --
in particular, trait (dis)similarity -- to play a central role in determining
whether a particular set of species coexists. Traits are, in turn, the outcome
of an eco-evolutionary process summarized by a phylogenetic tree. Therefore,
the phylogenetic tree associated with a set of species should carry information
about the dynamics and assembly properties of the community. Many studies have
highlighted the potentially complex ways in which this phylogenetic information
is translated into species' ecological properties. However, much less emphasis
has been placed on developing clear, quantitative expectations for community
properties under a particular hypothesis. To address this gap, we couple a
simple model of trait evolution on a phylogenetic tree with Lotka-Volterra
community dynamics. This allows us to derive properties of a community of
coexisting species as a function of the number of traits, tree topology and the
size of the species pool. Our analysis highlights how phylogenies, through
traits, affect the coexistence of a set of species. Together, these results
provide much-needed baseline expectations for the ways in which evolutionary
history, summarized by phylogeny, is reflected in the size and structure of
ecological communities.


  Folding of proteins into their correct native structure is key to their
function. Simultaneously, the intricate interplay between cell movement and
protein conformation highlights the complex nature of cellular processes. In
this work, we demonstrate the impact of Terahertz (THz) signaling on
controlling protein conformational changes in a random medium. Our system of
interest consists of a communication link that involves a nanoantenna
transmitter, a protein receiver, and a channel composed of moving red blood
cells. Due to the system dynamics, we investigate the influence of both the
fast and slow channel variations on protein folding. Specifically, we analyze
the system's selectivity to asses the effectiveness of the induced THz
interaction in targeting a specific group of proteins under fading conditions.
By optimizing the selectivity metric with respect to the nanoantenna power and
frequency, it is possible to enhance the controllability of protein
interactions. Our probabilistic analysis provides a new perspective regarding
electromagnetically triggered protein molecules, their microenvironment and
their interaction with surrounding particles. It helps elucidate how external
conditions impact the protein folding kinetics and pathways. This results in
not only understanding the mechanisms underlying THz-induced protein
interactions but also engineering these still-emerging tools.


  Predictive contagion models are ubiquitous in epidemiology, social sciences,
engineering, and management. This paper formulates a prescriptive contagion
analytics model where a decision-maker allocates shared resources across
multiple segments of a population, each governed by continuous-time dynamics.
We define four real-world problems under this umbrella: vaccine distribution,
vaccination centers deployment, content promotion, and congestion mitigation.
These problems feature a large-scale mixed-integer non-convex optimization
structure with constraints governed by ordinary differential equations,
combining the challenges of discrete optimization, non-linear optimization, and
continuous-time system dynamics. This paper develops a branch-and-price
methodology for prescriptive contagion analytics based on: (i) a set
partitioning reformulation; (ii) a column generation decomposition; (iii) a
state-clustering algorithm for discrete-decision continuous-state dynamic
programming; and (iv) a tri-partite branching scheme to circumvent
non-linearities. Extensive experiments show that the algorithm scales to very
large and otherwise-intractable instances, outperforming state-of-the-art
benchmarks. Our methodology provides practical benefits in contagion systems;
in particular, it can increase the effectiveness of a vaccination campaign by
an estimated 12-70%, resulting in 7,000 to 12,000 extra saved lives over a
three-month horizon mirroring the COVID-19 pandemic. We provide an open-source
implementation of the methodology in an online repository to enable
replication.


  The Turing completeness of continuous Chemical Reaction Networks (CRNs)
states that any computable real function can be computed by a continuous CRN on
a finite set of molecular species, possibly restricted to elementary reactions,
i.e. with at most two reactants and mass action law kinetics. In this paper, we
introduce a more stringent notion of robust online analog computation, and
Absolute Functional Robustness (AFR), for the CRNs that stabilize the
concentration values of some output species to the result of one function of
the input species concentrations, in a perfectly robust manner with respect to
perturbations of both intermediate and output species. We prove that the set of
real functions stabilized by a CRN with mass action law kinetics is precisely
the set of real algebraic functions. Based on this result, we present a
compiler which takes as input any algebraic function (defined by one polynomial
and one point for selecting one branch of the algebraic curve defined by the
polynomial) and generates an abstract CRN to stabilize it. Furthermore, we
provide error bounds to estimate and control the error of an unperturbed
system, under the assumption that the environment inputs are driven by
k-Lipschitz functions.


  Although bulk transcriptomic analyses have significantly contributed to an
enhanced comprehension of multifaceted diseases, their exploration capacity is
impeded by the heterogeneous compositions of biological samples. Indeed, by
averaging expression of multiple cell types, RNA-Seq analysis is oblivious to
variations in cellular changes, hindering the identification of the internal
constituents of tissues, involved in disease progression. On the other hand,
single-cell techniques are still time, manpower and resource-consuming
analyses.To address the intrinsic limitations of both bulk and single-cell
methodologies, computational deconvolution techniques have been developed to
estimate the frequencies of cell subtypes within complex tissues. These methods
are especially valuable for dissecting intricate tissue niches, with a
particular focus on tumour microenvironments (TME).In this paper, we offer a
comprehensive overview of deconvolution techniques, classifying them based on
their methodological characteristics, the type of prior knowledge required for
the algorithm, and the statistical constraints they address. Within each
category identified, we delve into the theoretical aspects for implementing the
underlying method, while providing an in-depth discussion of their main
advantages and disadvantages in supplementary materials.Notably, we emphasise
the advantages of cutting-edge deconvolution tools based on probabilistic
models, as they offer robust statistical frameworks that closely align with
biological realities. We anticipate that this review will provide valuable
guidelines for computational bioinformaticians in order to select the
appropriate method in alignment with their statistical and biological
objectives.We ultimately end this review by discussing open challenges that
must be addressed to accurately quantify closely related cell types from RNA
sequencing data, and the complementary role of single-cell RNA-Seq to that
purpose.


  Advances in single-cell omics allow for unprecedented insights into the
transcription profiles of individual cells. When combined with large-scale
perturbation screens, through which specific biological mechanisms can be
targeted, these technologies allow for measuring the effect of targeted
perturbations on the whole transcriptome. These advances provide an opportunity
to better understand the causative role of genes in complex biological
processes such as gene regulation, disease progression or cellular development.
However, the high-dimensional nature of the data, coupled with the intricate
complexity of biological systems renders this task nontrivial. Within the
machine learning community, there has been a recent increase of interest in
causality, with a focus on adapting established causal techniques and
algorithms to handle high-dimensional data. In this perspective, we delineate
the application of these methodologies within the realm of single-cell genomics
and their challenges. We first present the model that underlies most of current
causal approaches to single-cell biology and discuss and challenge the
assumptions it entails from the biological point of view. We then identify open
problems in the application of causal approaches to single-cell data:
generalising to unseen environments, learning interpretable models, and
learning causal models of dynamics. For each problem, we discuss how various
research directions - including the development of computational approaches and
the adaptation of experimental protocols - may offer ways forward, or on the
contrary pose some difficulties. With the advent of single cell atlases and
increasing perturbation data, we expect causal models to become a crucial tool
for informed experimental design.


  Conditional testing via the knockoff framework allows one to identify --
among large number of possible explanatory variables -- those that carry unique
information about an outcome of interest, and also provides a false discovery
rate guarantee on the selection. This approach is particularly well suited to
the analysis of genome wide association studies (GWAS), which have the goal of
identifying genetic variants which influence traits of medical relevance.
  While conditional testing can be both more powerful and precise than
traditional GWAS analysis methods, its vanilla implementation encounters a
difficulty common to all multivariate analysis methods: it is challenging to
distinguish among multiple, highly correlated regressors. This impasse can be
overcome by shifting the object of inference from single variables to groups of
correlated variables. To achieve this, it is necessary to construct "group
knockoffs." While successful examples are already documented in the literature,
this paper substantially expands the set of algorithms and software for group
knockoffs. We focus in particular on second-order knockoffs, for which we
describe correlation matrix approximations that are appropriate for GWAS data
and that result in considerable computational savings. We illustrate the
effectiveness of the proposed methods with simulations and with the analysis of
albuminuria data from the UK Biobank.
  The described algorithms are implemented in an open-source Julia package
Knockoffs.jl, for which both R and Python wrappers are available.


  Prediction of binding sites for transcription factors is important to
understand how they regulate gene expression and how this regulation can be
modulated for therapeutic purposes. Although in the past few years there are
significant works addressing this issue, there is still space for improvement.
In this regard, a transformer based capsule network viz. DNABERT-Cap is
proposed in this work to predict transcription factor binding sites mining
ChIP-seq datasets. DNABERT-Cap is a bidirectional encoder pre-trained with
large number of genomic DNA sequences, empowered with a capsule layer
responsible for the final prediction. The proposed model builds a predictor for
transcription factor binding sites using the joint optimisation of features
encompassing both bidirectional encoder and capsule layer, along with
convolutional and bidirectional long-short term memory layers. To evaluate the
efficiency of the proposed approach, we use a benchmark ChIP-seq datasets of
five cell lines viz. A549, GM12878, Hep-G2, H1-hESC and Hela, available in the
ENCODE repository. The results show that the average area under the receiver
operating characteristic curve score exceeds 0.91 for all such five cell lines.
DNABERT-Cap is also compared with existing state-of-the-art deep learning based
predictors viz. DeepARC, DeepTF, CNN-Zeng and DeepBind, and is seen to
outperform them.


  Motivation: The regulation of gene expression during tissue development is
extremely complex. One of the key regulatory mechanisms of gene expression
involves the recognition of regulatory motifs by various proteins in the
promoter regions of many genes. Localisation of these motifs in proximity to
the transcription start site (TSS) or translation start site (ATG) is critical
for regulating the initiation and rate of transcription. The levels of
transcription of individual genes, regulated by these motifs, can vary signifi
cantly in diff erent tissues and developmental stages, especially during
tightly regulated processes such as sexual reproduction. However, the precise
localisation and visualisation of the regulatory motifs within gene promoters
with respect to gene transcription in specifi c tissues, can be challenging.
Results: Here, we introduce a program called GOLEM (Gene regulatOry eLEMents)
which enables users to precisely locate any motif of interest with respect to
TSS or ATG within the relevant plant genomes across the plant Tree of Life
(Marchantia, Physcomitrium, Amborella, Oryza, Zea, Solanum and Arabidopsis).
The visualisation of the motifs is performed with respect to the transcript
levels of particular genes in leaves and male reproductive tissues, and can be
compared with genome-wide distribution regardless of the transcription level.
Availability and implementation: GOLEM is freely available at
https://golem.ncbr.muni.cz and its source codes are provided under the MIT
licence at GitHub at https://github.com/sb-ncbr/golem.


  Recently, drug repurposing has emerged as an effective and resource-efficient
paradigm for AD drug discovery. Among various methods for drug repurposing,
network-based methods have shown promising results as they are capable of
leveraging complex networks that integrate multiple interaction types, such as
protein-protein interactions, to more effectively identify candidate drugs.
However, existing approaches typically assume paths of the same length in the
network have equal importance in identifying the therapeutic effect of drugs.
Other domains have found that same length paths do not necessarily have the
same importance. Thus, relying on this assumption may be deleterious to drug
repurposing attempts. In this work, we propose MPI (Modeling Path Importance),
a novel network-based method for AD drug repurposing. MPI is unique in that it
prioritizes important paths via learned node embeddings, which can effectively
capture a network's rich structural information. Thus, leveraging learned
embeddings allows MPI to effectively differentiate the importance among paths.
We evaluate MPI against a commonly used baseline method that identifies anti-AD
drug candidates primarily based on the shortest paths between drugs and AD in
the network. We observe that among the top-50 ranked drugs, MPI prioritizes
20.0% more drugs with anti-AD evidence compared to the baseline. Finally, Cox
proportional-hazard models produced from insurance claims data aid us in
identifying the use of etodolac, nicotine, and BBB-crossing ACE-INHs as having
a reduced risk of AD, suggesting such drugs may be viable candidates for
repurposing and should be explored further in future studies.


  A good group reputation often facilitates more efficient synergistic teamwork
in production activities. Here we translate this simple motivation into a
reputation-based synergy and discounting mechanism in the public goods game.
Specifically, the reputation type of a group, either good or bad determined by
a reputation threshold, modifies the nonlinear payoff structure described by a
unified reputation impact factor. Results show that this reputation-based
incentive mechanism could effectively promote cooperation compared with linear
payoffs, despite the coexistence of synergy and discounting effects. Notably,
the complicated interactions between reputation impact and reputation threshold
result in a sharp phase transition from full cooperation to full defection. We
also find that the presence of a few discounting groups could increase the
average payoffs of cooperators, leading to an interesting phenomenon that when
the reputation threshold is raised, the gap between the average payoffs of
cooperations and defectors increases while the overall payoff decreases. Our
work provides important insights into facilitating cooperation in social
groups.


  With the recent move towards sequencing of accurate long reads, finding
solutions that support efficient analysis of these reads becomes more
necessary. The long execution time required for sequence alignment of long
reads negatively affects genomic studies relying on sequence alignment.
Although pre-alignment filtering as an extra step before alignment was recently
introduced to mitigate sequence alignment for short reads, these filters do not
work as efficiently for long reads. Moreover, even with efficient pre-alignment
filters, the overall end-to-end (i.e., filtering + original alignment)
execution time of alignment for long reads remains high, while the filtering
step is now a major portion of the end-to-end execution time.
  Our paper makes three contributions. First, it identifies data movement of
sequences between memory units and computing units as the main source of
inefficiency for pre-alignment filters of long reads. This is because although
filters reject many of these long sequencing pairs before they get to the
alignment stage, they still require a huge cost regarding time and energy
consumption for the large data transferred between memory and processor.
Second, this paper introduces an adaptation of a short-read pre-alignment
filtering algorithm suitable for long reads. We call this LongGeneGuardian.
Finally, it presents Filter-Fuse as an architecture that supports
LongGeneGuardian inside the memory. FilterFuse exploits the
Computation-In-Memory computing paradigm, eliminating the cost of data movement
in LongGeneGuardian.
  Our evaluations show that FilterFuse improves the execution time of filtering
by 120.47x for long reads compared to State-of-the-Art (SoTA) filter,
SneakySnake. FilterFuse also improves the end-to-end execution time of sequence
alignment by up to 49.14x and 5207.63x compared to SneakySnake with SoTA
aligner and only SoTA aligner, respectively.


  The minimal 3-state scheme of kinetic cooperativity of monomeric enzymes is
subjected to detailed analysis. The rigorous criteria of positive cooperativity
and its sigmoidal version are established in terms of the system parameters
(rate constants). It is shown that the cooperativity extent is especially
sensitive to the rates and direction of the exchange between conformational
states of the free enzyme. However, no necessity of the "kinetic resonance"
(or, moreover, its generality claimed recently) for enhancing cooperativity is
revealed. Overall, while the minimal 3-state model serves well for qualitative
understanding the origin of kinetic cooperativity, it is hardly suitable for
quantitative describing reactions of real enzymes, as it is shown with the case
of glucokinase.


  Darwinian evolution gives rise not only to the adaptation to the environment
but also to the enhancement of the robustness against mutation. Suppose more
than one phenotypes share the same fitness value. We expect that some
phenotypes are hardly selected as a consequence of the selection bias for
mutational robustness. We investigated this phenotype selection for a model of
gene regulatory networks (GRNs). We constructed the randomly generated set of
GRNs using the multicanonical Monte Carlo method and compared it to the
outcomes of evolutionary simulations. The results suggest that the mutationally
least robust phenotype is suppressed in evolution.


  The schooling behavior of fish can be studied through simulations involving a
large number of interacting particles. In such systems, each individual
particle is guided by behavior rules, which include aggregation towards a
centroid, collision avoidance, and direction alignment. The movement vector of
each particle may be expressed as a linear combination of behaviors, with
unknown parameters that define a trade-off among several behavioral
constraints. A fitness function for collective schooling behavior encompasses
all individual particle parameters.
  For a large number of interacting particles in a complex environment,
heuristic methods, such as evolutionary algorithms, are used to optimize the
fitness function, ensuring that the resulting decision rule preserves
collective behavior. However, these algorithms exhibit slow convergence, making
them inefficient in terms of CPU time cost.
  This paper proposes a CPU-efficient iterative (Cluster, Partition, Refine --
CPR) algorithm for estimating decision rule parameters for a large number of
interacting particles. In the first step, we employ the K-Means (unsupervised
learning) algorithm to cluster candidate solutions. Then, we partition the
search space using Voronoi tessellation over the defined clusters. We assess
the quality of each cluster based on the fitness function, with the centroid of
their Voronoi cells representing the clusters. Subsequently, we refine the
search space by introducing new cells into a number of identified well-fitting
Voronoi cells. This process is repeated until convergence.
  A comparison of the performance of the CPR algorithm with a standard Genetic
Algorithm reveals that the former converges faster than the latter. We also
demonstrate that the application of the CPR algorithm results in a schooling
behavior consistent with empirical observations.


  Proteins have evolved through mutations, amino acid substitutions, since life
appeared on Earth, some 109 years ago. The study of these phenomena has been of
particular significance because of their impact on protein stability, function,
and structure. Three of the most recent findings in these areas deserve to be
highlighted. First, an innovative method has made it feasible to massively
determine the impact of mutations on protein stability. Second, a theoretical
analysis showed how mutations impact the evolution of protein folding rates.
Lastly, it has been shown that native-state structural changes brought on by
mutations can be explained in detail by the amide hydrogen exchange protection
factors. This study offers a new perspective on how those findings can be used
to analyze proteins in the light of mutations. The preliminary results indicate
that: (i) mutations can be viewed as sensitive probes to identify "typos" in
the amino-acid sequence and also to assess the resistance of naturally
occurring proteins to unwanted sequence alterations; (ii) the presence of
"typos" in the amino acid sequence, rather than being an evolutionary obstacle,
could promote faster evolvability and, in turn, increase the likelihood of
higher protein stability; (iii) the mutation site is far more important than
the substituted amino acid in terms of the protein's marginal stability
changes, and (iv) the protein evolution unpredictability at the molecular level
by mutations exists even in the absence of epistasis effects. Finally, the
study results support the Darwinian concept of evolution as "descent with
modification" by demonstrating that some regions of any protein sequence are
susceptible to mutations while others are not.


  Would another origin of life resemble Earth's biochemical use of amino acids?
Here we review current knowledge at three levels: 1) Could other classes of
chemical structure serve as building blocks for biopolymer structure and
catalysis? Amino acids now seem both readily available to, and a plausible
chemical attractor for, life as we don't know it. Amino acids thus remain
important and tractable targets for astrobiological research. 2) If amino acids
are used, would we expect the same L-alpha-structural subclass used by life?
Despite numerous ideas, it is not clear why life favors L-enantiomers. It seems
clearer, however, why life on Earth uses the shortest possible (alpha-) amino
acid backbone, and why each carries only one side chain. However, assertions
that other backbones are physicochemically impossible have relaxed into
arguments that they are disadvantageous. 3) Would we expect a similar set of
side chains to those within the genetic code? Not only do many plausible
alternatives exist and evidence exists for both evolutionary advantage and
physicochemical constraint for those encoded by life. Overall, as focus shifts
from amino acids as a chemical class to specific side chains used by post-LUCA
biology, the probable role of physicochemical constraint diminishes relative to
that of biological evolution. Exciting opportunities now present themselves for
laboratory work and computing to explore how changing the amino acid alphabet
alters the universe of protein folds. Near-term milestones include: a)
expanding evidence about amino acids as attractors within chemical evolution;
b) extending characterization of other backbones relative to biological
proteins; c) merging computing and laboratory explorations of structures and
functions unlocked by xeno peptides.


  Understanding the combined effects of climate change and habitat
fragmentation on the adaptive potential of plant populations is essential for
devising effective conservation strategies. This is particularly important
where mating system variation impacts the evolutionary consequences of habitat
fragmentation. Here we aimed to reveal how habitat fragmentation and climate
adaptation jointly influence the evolutionary trajectories in Primula elatior,
a heterostylous self-incompatible and dispersal-constrained forest herb. We
quantified the genomic variation and degree of herkogamy, a floral trait
reducing self-pollination, across 60 geographically paired populations of
Primula elatior across Europe, each pair featuring contrasting levels of
habitat fragmentation. Our findings revealed a large and unique set of adaptive
outliers in more fragmented landscapes, compared to high-connectivity ones,
despite the geographic proximity of the sampling pairs. This suggests elevated
selective pressures in fragmented habitats, mirrored by a reduced adaptive
potential to cope with climate change. Finally, a minority of genetic variants
associated with herkogamy were influenced by current levels of habitat
fragmentation and population size, potentially signalling early indicators of
evolutionary mating system changes in response to pollinator limitation.
Because evolutionary trajectories and adaptive potential are expected to be
increasingly affected by habitat fragmentation, our findings underscore the
importance of considering both habitat fragmentation and climate adaptation in
conservation research and planning.


  Accurate and representative data is vital for precisely reporting the impact
of influenza in healthcare systems. Northern hemisphere winter 2022/23
experienced the most substantial influenza wave since the COVID-19 pandemic
began in 2020. Simultaneously, new data streams become available within health
services because of the pandemic. Comparing these data, surveillance and
administrative, supports the accurate monitoring of population level disease
trends. We analysed admissions rates per capita from four different collection
mechanisms covering National Health Service hospital Trusts in England over the
winter 2022/23 wave. We adjust for difference in reporting and extracted key
epidemic characteristics including the maximum admission rate, peak timing,
cumulative season admissions and growth rates by fitting generalised additive
models at national and regional levels. By modelling the admission rates per
capita across surveillance and administrative data systems we show that
different data measuring the epidemic produce different estimates of key
quantities. Nationally and in most regions the data correspond well for the
maximum admission rate, date of peak and growth rate, however, in subnational
analysis discrepancies in estimates arose, particularly for the cumulative
admission rate. This research shows that the choice of data used to measure
seasonal influenza epidemics can influence analysis substantially at
sub-national levels. For the admission rate per capita there is comparability
in the sentinel surveillance approach (which has other important functions),
rapid situational reports, operational databases and time lagged administrative
data giving assurance in their combined value. Utilising multiple sources of
data aids understanding of the impact of seasonal influenza epidemics in the
population.


  The architecture of the brain is too complex to be intuitively surveyable
without the use of compressed representations that project its variation into a
compact, navigable space. The task is especially challenging with
high-dimensional data, such as gene expression, where the joint complexity of
anatomical and transcriptional patterns demands maximum compression.
Established practice is to use standard principal component analysis (PCA),
whose computational felicity is offset by limited expressivity, especially at
great compression ratios. Employing whole-brain, voxel-wise Allen Brain Atlas
transcription data, here we systematically compare compressed representations
based on the most widely supported linear and non-linear methods-PCA, kernel
PCA, non-negative matrix factorization (NMF), t-stochastic neighbour embedding
(t-SNE), uniform manifold approximation and projection (UMAP), and deep
auto-encoding-quantifying reconstruction fidelity, anatomical coherence, and
predictive utility with respect to signalling, microstructural, and metabolic
targets. We show that deep auto-encoders yield superior representations across
all metrics of performance and target domains, supporting their use as the
reference standard for representing transcription patterns in the human brain.


  Two distinct families of pan-primate endogenous retroviruses, namely HERVL
and HERVH, infected primates germline, colonized host genomes, and evolved into
the global retroviral genomic regulatory dominion (GRD) operating during human
embryogenesis (HE). HE retroviral GRD constitutes 8839 highly conserved fixed
LTR elements linked to 5444 down-stream target genes forged by evolution into a
functionally-consonant constellation of 26 genome-wide multimodular genomic
regulatory networks (GRNs), each of which is defined by significant enrichment
of numerous single gene ontology (GO)-specific traits. Locations of GRNs appear
scattered across chromosomes to occupy from 5.5%-15.09% of human genome. Each
GRN harbors from 529-1486 retroviral LTRs derived from LTR7, MLT2A1, and MLT2A2
sequences that are quantitatively balanced according to their genome-wide
abundance. GRNs integrate activities from 199-805 down-stream target genes,
including transcription factors, chromatin-state remodelers, signal-sensing and
signal-transduction mediators, enzymatic and receptor binding effectors,
intracellular complexes and extracellular matrix elements, and cell-cell
adhesion molecules. GRNs compositions consist of several hundred to thousands
smaller GO enrichment-defined genomic regulatory modules (GRMs) combining from
a dozen to hundreds LTRs and down-stream target genes, which appear to operate
on individuals life-span timescale along specific phenotypic avenues to exert
profound effects on patterns of transcription, protein-protein interactions,
developmental phenotypes, physiological traits, and pathological conditions of
Modern Humans. Overall, this study identifies 69,573 statistically significant
retroviral LTR-linked GRMs (Binominal FDR q-value threshold of 0.001),
including 27,601 GRMs validated by the single GO-specific directed acyclic
graph (DAG) analyses across six GO annotations.


  In order to better understand the impact of environmental stochastic
fluctuations on the evolution of animal behavior, we introduce the concept of a
stochastic Nash equilibrium (SNE) that extends the classical concept of a Nash
equilibrium (NE). Based on a stochastic stability analysis of a linear
evolutionary game with temporally varying payoffs, we address the question of
the existence of a SNE, either weak when the geometric mean payoff against it
is the same for all other strategies or strong when it is strictly smaller for
all other strategies, and its relationship with a stochastically evolutionarily
stable (SES) strategy. While a strong SNE is always SES, this is not
necessarily the case for a weak SNE. We give conditions for a completely mixed
weak SNE not to be SES and to coexist with at least two strong SNE. More
importantly, we show that a pair of two completely mixed strong SNE can emerge
as the noise level increases. This not only indicates that a noise-induced SNE
may possess some properties that a NE cannot possess, such as being completely
mixed and strong, but also illustrates the complexity of evolutionary game
dynamics in a stochastic environment.


  The DisProt database is a significant resource containing manually curated
data on experimentally validated intrinsically disordered proteins (IDPs) and
regions (IDRs) from the literature. Developed in 2005, its primary goal was to
collect structural and functional information into proteins that lack a fixed
three-dimensional (3D) structure. Today, DisProt has evolved into a major
repository that not only collects experimental data but also contributes
significantly to our understanding of the IDPs/IDRs roles in various biological
processes, such as autophagy or the life cycle mechanisms in viruses, or their
involvement in diseases (such as cancer and neurodevelopmental disorders).
DisProt offers detailed information on the structural states of IDPs/IDRs,
including state transitions, interactions, and their functions, all provided as
curated annotations. One of the central activities of DisProt is the meticulous
curation of experimental data from the literature. For this reason, to ensure
that every expert and volunteer curator possesses the requisite knowledge for
data evaluation, collection, and integration, training courses and curation
materials are available. However, biocuration guidelines concur on the
importance of developing robust guidelines that not only provide critical
information about data consistency but also ensure data acquisition.This
guideline aims to provide both biocurators and external users with best
practices for manually curating IDPs and IDRs in DisProt. It describes every
step of the literature curation process and provides use cases of IDP curation
within DisProt.
  Database URL: https://disprot.org/


  The N-methyl-D-aspartate receptor is a prominent player in brain development
and functioning. Perturbations to its functioning through external stimuli like
magnetic fields can potentially affect the brain in numerous ways. Various
studies have shown that magnetic fields of varying strengths affect these
receptors. We propose that the radical pair mechanism, a quantum mechanical
process, could explain some of these field effects. Radicals of the form
$[\mbox{RO}^\bullet \mbox{ Mg($\mbox{H}_2$O$)_n$}^{+\bullet}]$, where R is a
protein residue that can be Serine or Tyrosine, are considered for this study.
The variation in the singlet fractional yield of the radical pairs, as a
function of magnetic field strength, is calculated to understand how the
magnetic field affects the products of the radical pair reactions. Based on the
results, the radical pair mechanism is a likely candidate for explaining the
magnetic field effects observed on the receptor activity. The model predicts
changes in the behaviour of the system as magnetic field strength is varied and
also predicts certain isotope effects. The results further suggest that similar
effects on radical pairs could be a plausible explanation for various magnetic
field effects within the brain.


  Computational complexity is a key limitation of genomic analyses. Thus, over
the last 30 years, researchers have proposed numerous fast heuristic methods
that provide computational relief. Comparing genomic sequences is one of the
most fundamental computational steps in most genomic analyses. Due to its high
computational complexity, optimized exact and heuristic algorithms are still
being developed. We find that these methods are highly sensitive to the
underlying data, its quality, and various hyperparameters. Despite their wide
use, no in-depth analysis has been performed, potentially falsely discarding
genetic sequences from further analysis and unnecessarily inflating
computational costs. We provide the first analysis and benchmark of this
heterogeneity. We deliver an actionable overview of the 11 most widely used
state-of-the-art methods for comparing genomic sequences. We also inform
readers about their advantages and downsides using thorough experimental
evaluation and different real datasets from all major manufacturers (i.e.,
Illumina, ONT, and PacBio). SequenceLab is publicly available at
https://github.com/CMU-SAFARI/SequenceLab.


  A number of top down bio-control models have been proposed where the
predators' efficacy is enhanced via the provision of additional food (AF).
These can lead to infinite time \emph{blow-up} of the predator . Predator
interference, predator competition, refuge, pest defense and cannibalism have
been proposed to achieve pest extinction, while keeping the predator population
bounded, but entail \emph{quadratic} or higher order damping terms. Thus an AF
model with monotone pest dependent functional response, devoid of these
mechanisms, cannot yield the pest extinction state.
  In the current manuscript we propose a novel patch driven bio-control model,
with monotone pest dependent functional response. Here a predator is introduced
into a ``patch" such as a prairie strip with AF, and then disperses into a
neighboring ``patch" such as a crop field, to target a pest. Assuming only
\emph{linear} dispersal in the predator, we show, (i) blow-up in the predator
can be prevented in both the prairie strip and the crop field, (ii) the pest
extinction state is attainable in the crop field, (iii) the AF model with patch
structure can keep pest densities lower than the AF model without patch
structure and/or the classical top-down bio-control model, (iv) the AF model
with patch structure possesses rich dynamics including chaos. In particular,
one observes ``patch specific chaos" - the pest occupying the crop-field can
oscillate chaotically, while the pest in the prairie strip oscillates
periodically. We discuss these results in light of bio-control strategies that
would utilize prairie strips.


  During a public health crisis like COVID-19, individuals' adoption of
protective behaviors, such as self-isolation and wearing masks, can
significantly impact the spread of the disease. In the meanwhile, the spread of
the disease can also influence individuals' behavioral choices. Moreover, when
facing uncertain losses, individuals' decisions tend to be irrational.
Therefore, it is critical to study individuals' irrational behavior choices in
the context of a pandemic. In this paper, we propose an epidemic-behavior
co-evolution model that captures the dynamic interplay between individual
decision-making and disease spread. To account for irrational decision-making,
we incorporate the Prospect Theory in our individual behavior modeling. We
conduct a theoretical analysis of the model, examining the steady states that
emerge from the co-evolutionary process. We use simulations to validate our
theoretical findings and gain further insights. This investigation aims to
enhance our understanding of the complex dynamics between individual behavior
and disease spread during a pandemic.


  Large protein language models are adept at capturing the underlying
evolutionary information in primary structures, offering significant practical
value for protein engineering. Compared to natural language models, protein
amino acid sequences have a smaller data volume and a limited combinatorial
space. Choosing an appropriate vocabulary size to optimize the pre-trained
model is a pivotal issue. Moreover, despite the wealth of benchmarks and
studies in the natural language community, there remains a lack of a
comprehensive benchmark for systematically evaluating protein language model
quality. Given these challenges, PETA trained language models with 14 different
vocabulary sizes under three tokenization methods. It conducted thousands of
tests on 33 diverse downstream datasets to assess the models' transfer learning
capabilities, incorporating two classification heads and three random seeds to
mitigate potential biases. Extensive experiments indicate that vocabulary sizes
between 50 and 200 optimize the model, whereas sizes exceeding 800
detrimentally affect the model's representational performance. Our code, model
weights and datasets are available at
https://github.com/ginnm/ProteinPretraining.