  Atlantic cod (Gadus morhua) in the southern Gulf of St. Lawrence (sGSL)
declined to low abundance in the early 1990s and have since failed to recover
due to high natural mortality, which has been linked to grey seal (Halichoerus
grypus) predation. Increased grey seal harvests have been suggested to improve
cod survival, however, predicting the response of cod to seal abundance changes
in the sGSL is complicated by a hypothesized triangular food web involving
seals, cod, and small pelagic fishes, wherein the pelagic fishes are prey for
cod and grey seals, but may also prey on young cod. Grey seals may therefore
have an indirect positive effect on prerecruit cod survival via predation on
pelagic fish. Using a multispecies model of intermediate complexity fitted to
various scientific and fisheries data, we found that seal predation accounted
for the majority of recent cod mortality and that cod will likely be extirpated
without a strong and rapid reduction in grey seal abundance. We did not find
evidence that reducing grey seal abundance will result in large increases to
herring biomass that could impair cod recovery.


  Developmental bias plays a major role in phenotypic evolution. Some
researchers have argued that phenotypes, regulated by development, can only
evolve along restricted trajectory under certain scenarios, such as the case
for mammalian molar size ratios. However, this view has been challenged.
Broadly speaking, sources for phenotypic variation remain largely unknown. The
study here presents a generalized Inhibitory Cascade Model and explains that
the original model described only means of phenotypes resulting from selection
when viewed under a higher taxonomic scope. Consequently, I propose the
combined property of development: certainty, when the prior intersegmental
inhibition is strong, and uncertainty, when the opposite holds. This property
potentially not only explains counterintuitively high levels of developmental
instability, but also plays an essential role in generating phenotypic
variation.


  The human genome encodes a family of editing enzymes known as APOBEC3
(apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3). Several
family members, such as APO-BEC3G, APOBEC3F, and APOBEC3H haplotype II, exhibit
activity against viruses such as HIV. These enzymes induce C-to-U mutations in
the negative strand of viral genomes, resulting in multiple G-to-A changes,
commonly referred to as 'hypermutation.' Mutations catalyzed by these enzymes
are sequence context-dependent in the HIV genome; for instance, APOBEC3G
preferen-tially mutates G within GG, TGG, and TGGG contexts, while other
members mutate G within GA, TGA, and TGAA contexts. However, the same sequence
context has not been explored in relation to these enzymes and HBV. In this
study, our objective is to identify the mutational footprint of APOBEC3 enzymes
in the HBV genome. To achieve this, we employ a multivariable data analytics
technique to investigate motif preferences and potential sequence hierarchies
of mutation by APOBEC3 enzymes using full genome HBV sequences from a diverse
range of naturally infected patients. This approach allows us to distinguish
between normal and hypermutated sequences based on the representation of mono-
to tetra-nucleotide motifs. Additionally, we aim to identify motifs associated
with hypermutation induced by different APOBEC3 enzymes in HBV genomes. Our
analyses reveal that either APOBEC3 enzymes are not active against HBV, or the
induction of G-to-A mutations by these enzymes is not sequence
context-dependent in the HBV genome.


  Mathematical models in ecology and epidemiology must be consistent with
observed data in order to generate reliable knowledge and sound policy.
Metapopulation systems, which consist of a collection of sub-populations at
various locations, pose technical challenges in statistical inference due to
nonlinear, stochastic interactions. Difficulties encountered in these
methodological issues can obstruct the core scientific questions concerning the
link between the mathematical models and the data. Progress in statistically
efficient simulation-based inference for partially observed stochastic dynamic
systems has enabled the development of statistically rigorous approaches to the
analysis of nonlinear but low-dimensional systems. Recently, an algorithm has
been developed which enables comparable inference for higher-dimensional models
arising in metapopulation systems. The COVID-19 pandemic provides a situation
where mathematical models and their policy implications were widely visible,
and we revisit an influential metapopulation model used to inform basic
epidemiological understanding early in the pandemic. Our methods support
self-critical data analysis, enabling us to identify and address model
limitations, and leading to a new model with substantially improved statistical
fit and parameter identifiability. Our results suggest that the lockdown
initiated on January 23, 2020 in China was more effective than previously
thought. We proceed to recommend statistical analysis standards for future
metapopulation system modeling.


  Molecular communication (MC) is recently featured as a novel communication
tool to connect individual biological nanorobots in vivo. It is expected that a
large number of nanorobots can form large multi-agent MC systems through MC to
accomplish complex and large-scale tasks that cannot be achieved by a single
nanorobot. However, most previous models for MC systems assume a unidirectional
diffusion communication channel and cannot capture the feedback between each
nanorobot, which is important for multi-agent MC systems. In this paper, we
introduce the system theoretic model for large-scale multi-agent MC systems
using transfer functions, and then propose a method to analyze the stability
for multi-agent MC systems. The proposed method decomposes the multi-agent MC
system into multiple single-input and single-output (SISO) systems, which
facilitates to analyze the stability of the large-scale multi-agent MC system.
Finally, we demonstrate the proposed method by analyzing the stability of a
specific large-scale multi-agent MC system.


  Cancer heterogeneity arises from complex molecular interactions. Elucidating
systems-level properties of gene interaction networks distinguishing cancer
from normal cells is critical for understanding disease mechanisms and
developing targeted therapies. Previous works focused only on identifying
differences in network structures. In this study, we used graph frequency
analysis of cancer genetic signals defined on a co-expression network to
describe the spectral properties of underlying cancer systems. We demonstrated
that cancer cells exhibit distinctive signatures in the graph frequency content
of their gene expression signals. Applying graph frequency filtering, graph
Fourier transforms, and its inverse to gene expression from different cancer
stages resulted in significant improvements in average F-statistics of the
genes compared to using their unfiltered expression levels. We propose graph
spectral properties of cancer genetic signals defined on gene co-expression
networks as cancer hallmarks with potential application for differential
co-expression analysis.


  The emergence of a spatially-organized population distribution depends on the
dynamics of the population and mediators of interaction (activators and
inhibitors). Two broad classes of models have been used to investigate when and
how self-organization is triggered, namely, reaction-diffusion and spatially
nonlocal models. Nevertheless, these models implicitly assume smooth
propagation scenarios, neglecting that individuals many times interact by
exchanging short and abrupt pulses of the mediating substance. A recently
proposed framework advances in the direction of properly accounting for these
short-scale fluctuations by applying a coarse-graining procedure on the pulse
dynamics. In this paper, we generalize the coarse-graining procedure and apply
the extended formalism to new scenarios in which mediators influence
individuals' reproductive success or their motility. We show that, in the slow-
and fast-mediator limits, pulsed interactions recover, respectively, the
reaction-diffusion and nonlocal models, providing a mechanistic connection
between them. Furthermore, at each limit, the spatial stability condition is
qualitatively different, leading to a timescale-induced transition where
spatial patterns emerge as mediator dynamics becomes sufficiently fast.


  Background: the stability of serum specimen during time storage is importance
in clinical and medical science researches in addition of diagnosis. Lipids are
organic molecules that classified into 8 classes: fatty acids, phospholipids,
glycerolipids, saccharolipids, polyketides, prenol lipids, sterol, and
sphingolipids. Lipid analysis (cholesterol, triglyceride, high density
lipoprotein, low density lipoprotein, very low-density lipoprotein) is a vital
tool for diagnosis of numerous disease such as cardiovascular disease.
follow-up the lipid profile progression is essential of numerous diseases.
Methods: a 72 apparently healthy individual were participated in current study
the serum sample was taken after 12-14 hours of fasting. The serum cholesterol,
triglyceride, high density cholesterol, low density cholesterol, and very
low-density cholesterol were determined. Results: there are no significant
value found in all groups when compared to control except triglyceride that was
significant after 24 hours of freezing.


  A survey was conducted from February 2022 to May 2022 on the usage of
antibiotics at a poultry farm in different areas of Multan, Punjab Pakistan. A
well-organized questionnaire was used for the collection of data. Sixty poultry
farms were surveyed randomly in the Multan district. All of these Farms were
using antibiotics. Antibiotics are commonly used for the treatment of diseases.
Some are used as preventive medicine and a few are used as growth promotors.
neomycin, erythromycin, oxytetracycline, streptomycin, and colistin are the
broad-spectrum antibiotics that are being used commercially. Enrofloxacin and
Furazolidone are the common antibiotics that are being used in Studies these
days. The class of Fluoroquinolones is commonly used in poultry farms.
Thirty-three patterns of antibiotic usage were observed at poultry farms.
multi-drug practices were also observed on various farms. In this study, 25% of
antibiotics are prescribed by the veterans while more than 90 % were acquired
from the veterinary store. This study provides information about the
antibiotics which are commonly being used in the study location district
Multan. It is expected that the finding of this survey will be helpful in the
development of new strategies against the misuse of antibiotics on farms.


  Viruses and their hosts are involved in an 'arms race' where they continually
evolve mechanisms to overcome each other. It has long been proposed that
intrinsic disorder provides a substrate for the evolution of viral hijack
functions and that short linear motifs (SLiMs) are important players in this
process. Here, we review evidence in support of this tenet from two model
systems: the papillomavirus E7 protein and the adenovirus E1A protein.
Phylogenetic reconstructions reveal that SLiMs appear and disappear multiple
times across evolution, providing evidence of convergent evolution within
individual viral phylogenies. Multiple functionally related SLiMs show strong
co-evolution signals that persist across long distances in the primary sequence
and occur in unrelated viral proteins. Moreover, changes in SLiMs are
associated with changes in phenotypic traits such as host range and tropism.
Tracking viral evolutionary events reveals that host switch events are
associated with the loss of several SLiMs, suggesting that SLiMs are under
functional selection and that changes in SLiMs support viral adaptation.
Fine-tuning of viral SLiM sequences can improve affinity, allowing them to
outcompete host counterparts. However, viral SLiMs are not always competitive
by themselves, and tethering of two suboptimal SLiMs by a disordered linker may
instead enable viral hijack. Coevolution between the SLiMs and the linker
indicates that the evolution of disordered regions may be more constrained than
previously thought. In summary, experimental and computational studies support
a role for SLiMs and intrinsic disorder in viral hijack functions and in viral
adaptive evolution.


  Metaproteomics has become a crucial omics technology for studying
microbiomes. In this area, the Unipept ecosystem, accessible at
https://unipept.ugent.be, has emerged as an invaluable resource for analyzing
metaproteomic data. It offers in-depth insights into both taxonomic
distributions and functional characteristics of complex ecosystems. This
tutorial explains essential concepts like Lowest Common Ancestor (LCA)
determination and the handling of peptides with missed cleavages. It also
provides a detailed, step-by-step guide on using the Unipept Web application
and Unipept Desktop for thorough metaproteomics analyses. By integrating
theoretical principles with practical methodologies, this tutorial empowers
researchers with the essential knowledge and tools needed to fully utilize
metaproteomics in their microbiome studies.


  In the rapidly evolving landscape of genomics, deep learning has emerged as a
useful tool for tackling complex computational challenges. This review focuses
on the transformative role of Large Language Models (LLMs), which are mostly
based on the transformer architecture, in genomics. Building on the foundation
of traditional convolutional neural networks and recurrent neural networks, we
explore both the strengths and limitations of transformers and other LLMs for
genomics. Additionally, we contemplate the future of genomic modeling beyond
the transformer architecture based on current trends in research. The paper
aims to serve as a guide for computational biologists and computer scientists
interested in LLMs for genomic data. We hope the paper can also serve as an
educational introduction and discussion for biologists to a fundamental shift
in how we will be analyzing genomic data in the future.


  In harsh environments, organisms may self-organize into spatially patterned
systems in various ways. So far, studies of ecosystem spatial self-organization
have primarily focused on apparent orders reflected by regular patterns.
However, self-organized ecosystems may also have cryptic orders that can be
unveiled only through certain quantitative analyses. Here we show that
disordered hyperuniformity as a striking class of hidden orders can exist in
spatially self-organized vegetation landscapes. By analyzing the
high-resolution remotely sensed images across the American drylands, we
demonstrate that it is not uncommon to find disordered hyperuniform vegetation
states characterized by suppressed density fluctuations at long range. Such
long-range hyperuniformity has been documented in a wide range of microscopic
systems. Our finding contributes to expanding this domain to accommodate
natural landscape ecological systems. We use theoretical modeling to propose
that disordered hyperuniform vegetation patterning can arise from three
generalized mechanisms prevalent in dryland ecosystems, including (1) critical
absorbing states driven by an ecological legacy effect, (2) scale-dependent
feedbacks driven by plant-plant facilitation and competition, and (3)
density-dependent aggregation driven by plant-sediment feedbacks. Our modeling
results also show that disordered hyperuniform patterns can help ecosystems
cope with arid conditions with enhanced functioning of soil moisture
acquisition. However, this advantage may come at the cost of slower recovery of
ecosystem structure upon perturbations. Our work highlights that disordered
hyperuniformity as a distinguishable but underexplored ecosystem
self-organization state merits systematic studies to better understand its
underlying mechanisms, functioning, and resilience.


  Biomedical information graphs are crucial for interaction discovering of
biomedical information in modern age, such as identification of multifarious
molecular interactions and drug discovery, which attracts increasing interests
in biomedicine, bioinformatics, and human healthcare communities. Nowadays,
more and more graph neural networks have been proposed to learn the entities of
biomedical information and precisely reveal biomedical molecule interactions
with state-of-the-art results. These methods remedy the fading of features from
a far distance but suffer from remedying such problem at the expensive cost of
redundant memory and time. In our paper, we propose a novel Residual Message
Graph Convolution Network (ResMGCN) for fast and precise biomedical interaction
prediction in a different idea. Specifically, instead of enhancing the message
from far nodes, ResMGCN aggregates lower-order information with the next round
higher information to guide the node update to obtain a more meaningful node
representation. ResMGCN is able to perceive and preserve various messages from
the previous layer and high-order information in the current layer with least
memory and time cost to obtain informative representations of biomedical
entities. We conduct experiments on four biomedical interaction network
datasets, including protein-protein, drug-drug, drug-target, and gene-disease
interactions, which demonstrates that ResMGCN outperforms previous
state-of-the-art models while achieving superb effectiveness on both storage
and time.


  Human genetic diseases often arise from point mutations, emphasizing the
critical need for precise genome editing techniques. Among these, base editing
stands out as it allows targeted alterations at the single nucleotide level.
However, its clinical application is hindered by low editing efficiency and
unintended mutations, necessitating extensive trial-and-error experimentation
in the laboratory. To speed up this process, we present an attention-based
two-stage machine learning model that learns to predict the likelihood of all
possible editing outcomes for a given genomic target sequence. We further
propose a multi-task learning schema to jointly learn multiple base editors
(i.e. variants) at once. Our model's predictions consistently demonstrated a
strong correlation with the actual experimental results on multiple datasets
and base editor variants. These results provide further validation for the
models' capacity to enhance and accelerate the process of refining base editing
designs.


  We introduce oxNA, a new model for the simulation of DNA-RNA hybrids which is
based on two previously developed coarse-grained models$\unicode{x2014}$oxDNA
and oxRNA. The model naturally reproduces the physical properties of hybrid
duplexes including their structure, persistence length and force-extension
characteristics. By parameterising the DNA-RNA hydrogen bonding interaction we
fit the model's thermodynamic properties to experimental data using both
average-sequence and sequence-dependent parameters. To demonstrate the model's
applicability we provide three examples of its use$\unicode{x2014}$calculating
the free energy profiles of hybrid strand displacement reactions, studying the
resolution of a short R-loop and simulating RNA-scaffolded wireframe origami.


  Financial incentives are provided by governments to encourage the plantation
of new woodland to increase habitat, biodiversity, carbon sequestration, and
other economic benefits for landowners. Whilst these are largely positive
effects, it is worth considering that greater biodiversity and presence of
wildlife species in proximity to agricultural holdings may pose a risk of
disease transmission between wildlife and livestock. Wildlife transmission and
the provision of a reservoir for infectious disease is particularly important
in the transmission dynamics of bovine tuberculosis.
  In this paper we develop an economic model for changing land use due to
forestry subsidies. We use this asses the impact on wild deer populations in
the newly created woodland areas and the emergent infectious disease risk
arising from the proximity of new and existing wild deer populations and
existing cattle holdings.
  We consider an area in the South-West of Scotland, having existing woodland,
deer populations, and extensive and diverse cattle farm holdings. In this area
we find that, with a varying level of subsidy and plausible new woodland
creation, the contact risk between areas of wild deer and cattle increases
between 26% and 35% over the contact risk present with zero subsidy.
  This model provides a foundation for extending to larger regions and for
examining potential risk mitigation strategies, for example the targeting of
subsidy in low risk areas or provisioning for buffer zones between woodland and
agricultural holdings.


  Molecular dynamics simulations provide a wealth of data whose in-depth
analysis can be computationally demanding and, sometimes, even unnecessary.
Dimensionality reduction techniques are thus routinely employed to simplify and
improve the interpretation of trajectories focusing on specific subsets of the
system's atoms; a key issue, in this context, is to determine the optimal
resolution level, i.e. the smallest number of atoms needed to preserve the
largest information content from the full atomistic trajectory. Here, we
introduce the protein optimal resolution identification method (PROPRE), an
unsupervised approach built on information theory principles that determines
the smallest number of atoms that need to be retained to attain a synthetic yet
informative description of a protein. By applying the method to a protein
dataset and two particular case studies, we show that this number is typically
between 1.5 and 2 times the number of residues in a protein; nonetheless, the
degree of conformational variability of the system influences the specific
number importantly, in that a broader range of large-scale conformations
correlates with fewer retained sites. The PROPRE method is implemented in
efficient and user-friendly python scripts, which are made available for
download on a github repository.


  The flow of ions through permeable channels causes voltage drop in
physiological nanodomains such as synapses, dendrites and dendritic spines, and
other protrusions. How the voltage changes around channels in these nanodomains
has remained poorly studied. We focus this book chapter on summarizing recent
efforts in computing the steady-state current, voltage and ionic concentration
distributions based on the Poisson-Nernst-Planck equations as a model of
electro-diffusion. We first consider the spatial distribution of an uncharged
particle density and derive asymptotic formulas for the concentration
difference by solving the Laplace's equation with mixed boundary conditions. We
study a constant particles injection rate modeled by a Neumann flux condition
at a channel represented by a small boundary target, while the injected
particles can exit at one or several narrow patches. We then discuss the case
of two species (positive and negative charges) and take into account motions
due to both concentration and electrochemical gradients. The voltage resulting
from charge interactions is calculated by solving the Poisson's equation. We
show how deep an influx diffusion propagates inside a nanodomain, for
populations of both uncharged and charged particles. We estimate the
concentration and voltage changes in relations with geometrical parameters and
quantify the impact of membrane curvature.


  We present a review of a series of contact maps for the determination of
native interactions in proteins and nucleic acids based on a
distance-threshold. Such contact maps are mostly based on physical and chemical
construction, and yet they are sensitive to some parameters (e.g. distances or
atomic radii) and can neglect some key interactions. Furthermore, we also
comment on a new class of contact maps that only requires geometric arguments.
The contact map is a necessary ingredient to build a robust G\=oMartini model
for proteins and their complexes in the Martini 3 force field. We present the
extension of a popular structure-based G\=o-like approach for the study of
protein-sugar complexes, and also limitations of this approach are discussed.
The G\=oMartini approach was first introduced by Poma et al. J. Chem. Theory
Comput. 2017, 13(3), 1366-1374 in Martini 2 force field and recently, it has
gained the status of gold-standard for protein simulation undergoing
conformational changes in Martini 3 force field. We discuss several studies
that have provided support to this approach in the context of the biophysical
community.


  In this work we define a kinetic model for understanding the impact of
heterogeneous opinion formation dynamics on epidemics. The considered
many-agent system is characterized by nonsymmetric interactions which define a
coupled system of kinetic equations for the evolution of the opinion density in
each compartment. In the quasi-invariant limit we may show positivity and
uniqueness of the solution of the problem together with its convergence towards
an equilibrium distribution exhibiting bimodal shape. The tendency of the
system towards opinion clusters is further analyzed by means of numerical
methods, which confirm the consistency of the kinetic model with its moment
system whose evolution is approximated in several regimes of parameters.


  RNA design is the search for a sequence or set of sequences that will fold
into predefined structures, also known as the inverse problem of RNA folding.
While numerous RNA design methods have been invented to find sequences capable
of folding into a target structure, little attention has been given to the
identification of undesignable structures according to the minimum free energy
(MFE) criterion under the Turner model. In this paper, we address this gap by
first introducing mathematical theorems outlining sufficient conditions for
recognizing undesignable structures, then proposing efficient algorithms,
guided by these theorems, to verify the undesignability of RNA structures.
Through the application of these theorems and algorithms to the Eterna100
puzzles, we demonstrate the ability to efficiently establish that 15 of the
puzzles indeed fall within the category of undesignable structures. In
addition, we provide specific insights from the study of undesignability, in
the hope that it will enable more understanding of RNA folding and RNA design.


  What is the long term dynamics of the Covid-19 pandemic? How will it end?
Here we constructed an infectious disease model with testing and analyzed the
existence and stability of its endemic states. For a large parameter set,
including those relevant to the SARS-CoV-2 virus, we demonstrated the existence
of one endemic equilibrium without testing and one endemic equilibrium with
testing and proved their local and global stabilities for some cases. Our
results suggest that the pandemic is to end with a testing-less endemic state
through a novel and surprising mechanism called stochastic trapping.


  Contraction assay based on surface measurement have been widely used to
evaluate cell contractility in 3D models. This method is straightforward and
requires no specific equipment, but it does not provide quantitative data about
contraction forces generated by cells. We expanded this method with a new
biomechanical model, based on the work-energy theorem, to provide
non-destructive longitudinal monitoring of contraction forces generated by
cells in 3D.We applied this method on hydrogels seeded with either fibroblasts
or osteoblasts. Hydrogel mechanical characteristics were modulated to enhance
(condition HCA$_{High}$: hydrogel contraction assay high contraction) or limit
(condition HCA$_{Low}$: hydrogel contraction assay low contraction) cell
contractile behaviors. Macroscopic measures were further correlated with cell
contractile behavior and descriptive analysis of their physiology in response
to different mechanical environments. Fibroblasts and osteoblasts contracted
their matrix up to 47% and 77% respectively. Contraction stress peaked at day 5
with 1.1 10$^{-14}$Pa for fibroblasts and 3.5 10$^{-14}$Pa for osteoblasts,
which correlated with cell attachment and spreading. Negligible contraction was
seen in HCA$_{Low}$. Both fibroblasts and osteoblasts expressed $\alpha$-SMA
contractile fibers in HCA$_{High}$ and HCA$_{Low}$. Failure to contract
HCA$_{Low}$ was attributed to increased cross-linking and resistance to
proteolytic degradation of the hydrogel.


  Temporal environmental noise (EN) is a prevalent natural phenomenon that
controls population and community dynamics, shaping the destiny of biological
species and genetic types. Conventional theoretical models often depict EN as a
Markovian process with an exponential distribution of correlation times,
resulting in two distinct qualitative dynamical categories: quenched
(pertaining to short demographic timescales) and annealed (pertaining to long
timescales). However, numerous empirical studies demonstrate a fat-tailed decay
of correlation times. Here, we study the consequences of power-law correlated
EN on the dynamics of isolated and competing populations. We reveal the
emergence of a novel intermediate phase that lies between the quenched and
annealed regimes. Within this phase, dynamics are primarily driven by rare, yet
not exceedingly rare, long periods of almost-steady environmental conditions.
For an isolated population, the time to extinction in this phase exhibits a
novel scaling with the abundance, and also a non-monotonic dependence on the
spectral exponent.


  The illness-death model for chronic conditions is combined with a renewal
equation for the number of newborns taking into account possibly different
fertility rates in the healthy and diseased parts of the population. The
resulting boundary value problem consists of a system of partial differential
equations with an integral boundary condition. As an application, the boundary
value problem is applied to an example about type 2 diabetes.


  In an ecosystem, environmental changes as a result of natural and human
processes can cause some key parameters of the system to change with time.
Depending on how fast such a parameter changes, a tipping point can occur.
Existing works on rate-induced tipping, or R-tipping, offered a theoretical way
to study this phenomenon but from a local dynamical point of view, revealing,
e.g., the existence of a critical rate for some specific initial condition
above which a tipping point will occur. As ecosystems are subject to constant
disturbances and can drift away from their equilibrium point, it is necessary
to study R-tipping from a global perspective in terms of the initial conditions
in the entire relevant phase space region. In particular, we introduce the
notion of the probability of R-tipping defined for initial conditions taken
from the whole relevant phase space. Using a number of real-world, complex
mutualistic networks as a paradigm, we discover a scaling law between this
probability and the rate of parameter change and provide a geometric theory to
explain the law. The real-world implication is that even a slow parameter
change can lead to a system collapse with catastrophic consequences. In fact,
to mitigate the environmental changes by merely slowing down the parameter
drift may not always be effective: only when the rate of parameter change is
reduced to practically zero would the tipping be avoided. Our global dynamics
approach offers a more complete and physically meaningful way to understand the
important phenomenon of R-tipping.


  In recent years, the development of Artificial Intelligence (AI) has offered
the possibility to tackle many interdisciplinary problems, and the field of
chemistry is not an exception. Drug analysis is crucial in drug discovery,
playing an important role in human life. However, this task encounters many
difficulties due to the wide range of computational chemistry methods. Drug
analysis also involves a massive amount of work, including determining taste.
Thus, applying deep learning to predict a molecule's bitterness is inevitable
to accelerate innovation in drug analysis by reducing the time spent. This
paper proposes an artificial neural network (ANN) based approach (EC-ANN) for
the molecule's bitter prediction. Our approach took the SMILE (Simplified
molecular-input line-entry system) string of a molecule as the input data for
the prediction, and the 256-bit ECFP descriptor is the input vector for our
network. It showed impressive results compared to state-of-the-art, with a
higher performance on two out of three test sets according to the experiences
on three popular test sets: Phyto-Dictionary, Unimi, and Bitter-new set [1].
For the Phyto-Dictionary test set, our model recorded 0.95 and 0.983 in
F1-score and AUPR, respectively, depicted as the highest score in F1-score. For
the Unimi test set, our model achieved 0.88 in F1-score and 0.88 in AUPR, which
is roughly 12.3% higher than the peak of previous models [1, 2, 3, 4, 5].


  Accurately modeling protein 3D structure is essential for the design of
functional proteins. An important sub-task of structure modeling is protein
side-chain packing: predicting the conformation of side-chains (rotamers) given
the protein's backbone structure and amino-acid sequence. Conventional
approaches for this task rely on expensive sampling procedures over
hand-crafted energy functions and rotamer libraries. Recently, several deep
learning methods have been developed to tackle the problem in a data-driven
way, albeit with vastly different formulations (from image-to-image translation
to directly predicting atomic coordinates). Here, we frame the problem as a
joint regression over the side-chains' true degrees of freedom: the dihedral
$\chi$ angles. We carefully study possible objective functions for this task,
while accounting for the underlying symmetries of the task. We propose
Holographic Packer (H-Packer), a novel two-stage algorithm for side-chain
packing built on top of two light-weight rotationally equivariant neural
networks. We evaluate our method on CASP13 and CASP14 targets. H-Packer is
computationally efficient and shows favorable performance against conventional
physics-based algorithms and is competitive against alternative deep learning
solutions.


  In evolutionary biology, phylogenetic trees are commonly inferred from a set
of characters (partitions) of a collection of biological entities (e.g.,
species or individuals in a population). Such characters naturally arise from
molecular sequences or morphological data. Interestingly, it has been known for
some time that any binary phylogenetic tree can be (convexly) defined by a set
of at most four characters, and that there are binary phylogenetic trees for
which three characters are not enough. Thus, it is of interest to characterise
those phylogenetic trees that are defined by a set of at most three characters.
In this paper, we provide such a characterisation, in particular proving that a
binary phylogenetic tree $T$ is defined by a set of at most three characters
precisely if $T$ has no internal subtree isomorphic to a certain tree.


  Human aging is a process controlled by both genetics and environment. Many
studies have been conducted to identify a subset of genes related to aging from
the human genome. Biologists implicitly categorize age-related genes into genes
that cause aging and genes that are influenced by aging, which resulted in both
causal inference and inference of associations studies. While inference of
association is better explored, causal inference and computational causal
inference, remains less explored. In this work, we are primarily motivated to
tackle the problem of identifying genes associated with aging, while having a
brief look into genes with probable causal relations, both from a computational
perspective. Specifically, we form a set of hypotheses and accordingly,
introduce a data-tailored framework for inference. First we perform linear
modeling on the expression values of age-related genes, and then examine the
presence of heteroskedastic properties in the residual of the model. We
evaluate this framework and our results suggest that, 1) presence of
heteroskedasticity in these residuals is a potential signature of association
for age-related genes, and 2) consistent heteroskedasticity along the human
life span could imply some sort of causality. To our knowledge, along with
identifying age-associated genes, this is the first work to propose a framework
for computational causal inference on age-related genes, using a dataset of
human dermal fibroblast gene expression data. Hence the results of our simple,
yet effective approach can be used not only to assess future age-related genes,
but also as a possible criterion to select new associative or potential causal
genes with respect to aging.


  The Lotka-Volterra (LV) model is a simple, robust, and versatile model used
to describe large interacting systems such as food webs or microbiomes. The
model consists of $n$ coupled differential equations linking the abundances of
$n$ different species. We consider a large random interaction matrix with
independent entries and a block variance profile. The $i$th diagonal block
represents the intra-community interaction in community $i$, while the
off-diagonal blocks represent the inter-community interactions. The variance
remains constant within each block, but may vary across blocks.
  We investigate the important case of two communities of interacting species,
study how interactions affect their respective equilibrium. We also describe
equilibrium with feasibility (i.e., whether there exists an equilibrium with
all species at non-zero abundances) and the existence of an attrition
phenomenon (some species may vanish) within each community.
  Information about the general case of $b$ communities ($b> 2$) is provided in
the appendix.


  Proteins are intricate molecular machines whose complexity arises from the
heterogeneity of the amino acid building blocks and their dynamic network of
many-body interactions. These nanomachines gain function when put in the
context of a whole organism through interaction with other inhabitants of the
biological realm. And this functionality shapes their evolutionary histories
through intertwined paths of selection and adaptation. Recent advances in
machine learning have solved the decades-old problem of how protein sequence
determines their structure. However, the ultimate question regarding the basic
logic of protein machines remains open: How does the collective physics of
proteins lead to their functionality? and how does a sequence encode the full
range of dynamics and chemical interactions that facilitate function? Here, we
explore these questions within a physical approach that treats proteins as
mechano-chemical machines, which are adapted to function via concerted
evolution of structure, motion, and chemical interactions.


  The drug development pipeline for a new compound can last 10-20 years and
cost over 10 billion. Drug repurposing offers a more time- and cost-effective
alternative. Computational approaches based on biomedical knowledge graph
representations have recently yielded new drug repurposing hypotheses. In this
study, we present a novel, disease-specific hypergraph representation learning
technique to derive contextual embeddings of biological pathways of various
lengths but that all start at any given drug and all end at the disease of
interest. Further, we extend this method to multi-disease hypergraphs. To
determine the repurposing potential of each of the 1,522 drugs, we derive
drug-specific distributions of cosine similarity values and ultimately consider
the median for ranking. Cosine similarity values are computed between (1) all
biological pathways starting at the considered drug and ending at the disease
of interest and (2) all biological pathways starting at drugs currently
prescribed against that disease and ending at the disease of interest. We
illustrate our approach with Alzheimer's disease (AD) and two of its risk
factors: hypertension (HTN) and type 2 diabetes (T2D). We compare each drug's
rank across four hypergraph settings (single- or multi-disease): AD only, AD +
HTN, AD + T2D, and AD + HTN + T2D. Notably, our framework led to the
identification of two promising drugs whose repurposing potential was
significantly higher in hypergraphs combining two diseases: dapagliflozin
(antidiabetic; moved up, from top 32$\%$ to top 7$\%$, across all considered
drugs) and debrisoquine (antihypertensive; moved up, from top 76$\%$ to top
23$\%$). Our approach serves as a hypothesis generation tool, to be paired with
a validation pipeline relying on laboratory experiments and semi-automated
parsing of the biomedical literature.


  We build and study an individual based model of the telomere length's
evolution in a population across multiple generations. This model is a
continuous time typed branching process, where the type of an individual
includes its gamete mean telomere length and its age. We study its Malthusian's
behaviour and provide numerical simulations to understand the influence of
biologically relevant parameters.


  In a study conducted in Singapore, a country prone to dengue outbreaks due to
its climate and urban population, researchers examined the effectiveness of
releasing male Aedes aegypti mosquitoes infected with Wolbachia (wAlbB strain)
to reduce dengue transmission. These infected males, when mating with wild-type
females, produced non-viable eggs, leading to vector suppression. Extensive
field trials involving over 600,000 residents in four townships were conducted
from 2018 to 2022. The results showed a 57% decline in total dengue incidence
and a 64% decline in clustered dengue incidence. This approach offers promise
for large-scale dengue control in regions facing rising dengue cases, providing
a critical solution in combating the disease.


  Inference of the reproduction number through time is of vital importance
during an epidemic outbreak. Typically, epidemiologists tackle this using
observed prevalence or incidence data. However, prevalence and incidence data
alone is often noisy or partial. Models can also have identifiability issues
with determining whether a large amount of a small epidemic or a small amount
of a large epidemic has been observed. Sequencing data however is becoming more
abundant, so approaches which can incorporate genetic data are an active area
of research. We propose using particle MCMC methods to infer the time-varying
reproduction number from a combination of prevalence data reported at a set of
discrete times and a dated phylogeny reconstructed from sequences. We validate
our approach on simulated epidemics with a variety of scenarios. We then apply
the method to a real data set of HIV-1 in North Carolina, USA, between 1957 and
2019.


  In multistable dynamical systems driven by weak Gaussian noise, transitions
between competing states are often assumed to pass via a saddle on the
separating basin boundary. In contrast, we show that timescale separation can
cause saddle avoidance in non-gradient systems. Using toy models from
neuroscience and ecology, we study cases where sample transitions deviate
strongly from the instanton predicted by Large Deviation Theory, even for weak
finite noise. We attribute this to a flat quasipotential and propose an
approach based on the Onsager-Machlup action to aptly predict transition paths.


  Modeling environmental ecosystems is critical for the sustainability of our
planet, but is extremely challenging due to the complex underlying processes
driven by interactions amongst a large number of physical variables. As many
variables are difficult to measure at large scales, existing works often
utilize a combination of observable features and locally available measurements
or modeled values as input to build models for a specific study region and time
period. This raises a fundamental question in advancing the modeling of
environmental ecosystems: how to build a general framework for modeling the
complex relationships amongst various environmental data over space and time?
In this paper, we introduce a new framework, FREE, which maps available
environmental data into a text space and then converts the traditional
predictive modeling task in environmental science to the semantic recognition
problem. The proposed FREE framework leverages recent advances in Large
Language Models (LLMs) to supplement the original input features with natural
language descriptions. This facilitates capturing the data semantics and also
allows harnessing the irregularities of input features. When used for long-term
prediction, FREE has the flexibility to incorporate newly collected
observations to enhance future prediction. The efficacy of FREE is evaluated in
the context of two societally important real-world applications, predicting
stream water temperature in the Delaware River Basin and predicting annual corn
yield in Illinois and Iowa. Beyond the superior predictive performance over
multiple baseline methods, FREE is shown to be more data- and
computation-efficient as it can be pre-trained on simulated data generated by
physics-based models.


  Networks or graphs are widely used across the sciences to represent
relationships of many kinds. igraph (https://igraph.org) is a general-purpose
software library for graph construction, analysis, and visualisation, combining
fast and robust performance with a low entry barrier. igraph pairs a fast core
written in C with beginner-friendly interfaces in Python, R, and Mathematica.
Over the last two decades, igraph has expanded substantially. It now scales to
billions of edges, supports Mathematica and interactive plotting, integrates
with Jupyter notebooks and other network libraries, includes new graph layouts
and community detection algorithms, and has streamlined the documentation with
examples and Spanish translations. Modern testing features such as continuous
integration, address sanitizers, stricter typing, and memory-managed vectors
have also increased robustness. Hundreds of bug reports have been fixed and a
community forum has been opened to connect users and developers. Specific
effort has been made to broaden use and community participation by women,
non-binary people, and other demographic groups typically underrepresented in
open source software.


  Given a flow network, the Minimum Flow Decomposition (MFD) problem is finding
the smallest possible set of weighted paths whose superposition equals the
flow. It is a classical, strongly NP-hard problem that is proven to be useful
in RNA transcript assembly and applications outside of Bioinformatics. We
improve an existing ILP (Integer Linear Programming) model by Dias et al.
[RECOMB 2022] for DAGs by decreasing the solver's search space using solution
safety and several other optimizations. This results in a significant speedup
compared to the original ILP, of up to 55-90x on average on the hardest
instances. Moreover, we show that our optimizations apply also to MFD problem
variants, resulting in similar speedups, going up to 123x on the hardest
instances. We also developed an ILP model of reduced dimensionality for an MFD
variant in which the solution path weights are restricted to a given set. This
model can find an optimal MFD solution for most instances, and overall, its
accuracy significantly outperforms that of previous greedy algorithms while
being up to an order of magnitude faster than our optimized ILP.


  Logarithmic growth-rates are fundamental observables for describing
ecological systems and the characterization of their distributions with
analytical techniques can greatly improve their comprehension. Here a neutral
model based on a stochastic differential equation with demographic noise, which
presents a closed form for these distributions, is used to describe the
population dynamics of microbiota. Results show that this model can
successfully reproduce the log-growth rate distribution of the considered
abundance time-series. More significantly, it predicts its temporal dependence,
by reproducing its kurtosis evolution when the time lag $\tau$ is increased.
Furthermore, its typical shape for large $\tau$ is assessed, verifying that the
distribution variance does not diverge with $\tau$. The simulated processes
generated by the calibrated stochastic equation and the analysis of each
time-series, taken one by one, provided additional support for our approach.
Alternatively, we tried to describe our dataset by using a logistic model with
an environmental stochastic term. Analytical and numerical results show that
this model is not suited for describing the leptokurtic log-growth rates
distribution found in our data. These results effectively support a neutral
model with demographic stochasticity for describing the growth-rate dynamics
and the stationary abundance distribution of the considered microbiota. This
suggests that there are no significant parametric demographic differences among
the species, which can be statistically characterized by the same vital rates.


  The U.S. Covid-19 data exhibit a high-frequency oscillation along a
low-frequency wave for outbreaks. There is no model to account for it. A
modified SIR model is proposed to explain this spiking phenomenon. It is also
used to best-fit the data and to make forecast. For the simulated duration of
590 days, the model is capable of achieving a 0.5 percent mean squared relative
error (MSRE) fit to the seven-day average of the daily case number. The
outright 28-day prediction by the model generates a 20 percent MSRE for the
cumulative case total due to a persistent underestimation of the data by the
model. With the proposed correction to the aberration, the model is able to
keep the 28-day cumulative case total forecast within 10 percent MSRE of the
data.


  In this study, we present a novel molecular fingerprint generation method
based on multiparameter persistent homology. This approach reveals the latent
structures and relationships within molecular geometry, and detects topological
features that exhibit persistence across multiple scales along multiple
parameters, such as atomic mass, partial charge, and bond type, and can be
further enhanced by incorporating additional parameters like ionization energy,
electron affinity, chirality and orbital hybridization. The proposed
fingerprinting method provides fresh perspectives on molecular structure that
are not easily discernible from single-parameter or single-scale analysis.
Besides, in comparison with traditional graph neural networks, multiparameter
persistent homology has the advantage of providing a more comprehensive and
interpretable characterization of the topology of the molecular data. We have
established theoretical stability guarantees for multiparameter persistent
homology, and have conducted extensive experiments on the Lipophilicity,
FreeSolv, and ESOL datasets to demonstrate its effectiveness in predicting
molecular properties.


  Why do phylogenetic algorithms fail when they return incorrect answers? This
simple question has not been answered in detail, even for maximum parsimony
(MP), the simplest phylogenetic criterion. Understanding MP has recently gained
relevance in the regime of extremely dense sampling, where each virus sample
commonly differs by zero or one mutation from another previously sampled virus.
Although recent research shows that evolutionary histories in this regime are
close to being maximally parsimonious, the structure of their deviations from
MP is not yet understood. In this paper, we develop algorithms to understand
how the correct tree deviates from being MP in the densely sampled case. By
applying these algorithms to simulations that realistically mimic the evolution
of SARS-CoV-2, we find that simulated trees frequently only deviate from
maximally parsimonious trees locally, through simple structures consisting of
the same mutation appearing independently on sister branches.


  The emergence of collective cooperation in competitive environments is a
well-known phenomenon in biology, economics and social systems. While most
evolutionary game models focus on the evolution of strategies for a fixed game,
how strategic decisions co-evolve with the environment has so far mostly been
overlooked. Here, we consider a game selection model where not only the
strategies but also the game can change over time following evolutionary
principles. Our results show that co-evolutionary dynamics of games and
strategies can induce novel collective phenomena, fostering the emergence of
cooperative environments. When the model is taken on structured populations the
architecture of the interaction network can significantly amplify pro-social
behaviour, with a critical role played by network heterogeneity and the
presence of clustered groups, distinctive features observed in real-world
populations. By unveiling the link between the evolution of strategies and
games for different structured populations, our model sheds new light on the
origin of social dilemmas ubiquitously observed in real-world social systems.


  Molecular sciences address a wide range of problems involving molecules of
different types and sizes and their complexes. Recently, geometric deep
learning, especially Graph Neural Networks, has shown promising performance in
molecular science applications. However, most existing works often impose
targeted inductive biases to a specific molecular system, and are inefficient
when applied to macromolecules or large-scale tasks, thereby limiting their
applications to many real-world problems. To address these challenges, we
present PAMNet, a universal framework for accurately and efficiently learning
the representations of three-dimensional (3D) molecules of varying sizes and
types in any molecular system. Inspired by molecular mechanics, PAMNet induces
a physics-informed bias to explicitly model local and non-local interactions
and their combined effects. As a result, PAMNet can reduce expensive
operations, making it time and memory efficient. In extensive benchmark
studies, PAMNet outperforms state-of-the-art baselines regarding both accuracy
and efficiency in three diverse learning tasks: small molecule properties, RNA
3D structures, and protein-ligand binding affinities. Our results highlight the
potential for PAMNet in a broad range of molecular science applications.


  In this paper, we consider the infection-age-dependent Kermack--McKendrick
model in which host individuals are distributed in a continuous state space. To
provide a mathematical foundation for the heterogeneous model, we develop a
$L^1$-framework to formulate basic epidemiological concepts.
  First, we show the mathematical well-posedness of the basic model under
appropriate conditions allowing the unbounded parameters with non-compact
domain. Next we define the basic reproduction number and prove the pandemic
threshold theorem. We then present a systematic procedure to compute the
effective reproduction number and the herd immunity threshold. Finally we give
some illustrative examples and concrete results by using the separable mixing
assumption.


  Many infectious diseases are comprised of multiple strains with examples
including Influenza, tuberculosis, and Dengue virus. The time evolution of such
systems is linked to a complex landscape shaped by interactions between
competing strains. Possible long-term dynamics include the extinction of less
competitive strains, convergence to multi-strain steady-states, or
self-sustained oscillations. This work considers a two-strain epidemic model in
which the strains can interact indirectly via the immunity response generated
following infections, and in which this immune response wanes with time. In
particular, we focus on scenarios where the rate of waning immunity is
significantly faster than the rate of demographic turnover. The first key
result of this study is the explicit computation of the steady states of the
nonlinear system of seven equations. Following this result, we take advantage
of the separation of time scales in the problem and use perturbation methods to
analyze the stability of the fixed points. In particular, we establish the
conditions under which the system gives rise to the coexistence of the two
strains and whether coexistence is attained via convergence to an endemic
steady-state or via self-sustained oscillations. Our study unveils two
parameter regimes of distinct qualitative behavior of the system and
characterizes the separatrix between them. Within the first regime, the system
gives rise to oscillatory coexistence for all feasible conditions. In the
second regime, the system's behavior is governed by a solution to a quadratic
equation, potentially resulting in the convergence to a multi-strain endemic
equilibrium or the persistence of oscillatory coexistence.


  On August 9-10, 2023, a workshop was convened at the Pacific Northwest
National Laboratory (PNNL) in Richland, WA that brought together a group of
internationally recognized experts in metabolomics, natural products discovery,
chemical ecology, chemical and biological threat assessment, cheminformatics,
computational chemistry, cloud computing, artificial intelligence, and novel
technology development. These experts were invited to assess the value and
feasibility of a grand-scale project to create new technologies that would
allow the identification and quantification of all small molecules, or to
decode the molecular universe. The Decoding the Molecular Universe project
would extend and complement the success of the Human Genome Project by
developing new capabilities and technologies to measure small molecules
(defined as non-protein, non-polymer molecules less than 1500 Daltons) of any
origin and generated in biological systems or produced abiotically. Workshop
attendees 1) explored what new understanding of biological and environmental
systems could be revealed through the lens of small molecules; 2) characterized
the similarities in current needs and technical challenges between each science
or mission area for unambiguous and comprehensive determination of the
composition and quantities of small molecules of any sample; 3) determined the
extent to which technologies or methods currently exist for unambiguously and
comprehensively determining the small molecule composition of any sample and in
a reasonable time; and 4) identified the attributes of the ideal technology or
approach for universal small molecule measurement and identification. The
workshop concluded with a discussion of how a project of this scale could be
undertaken, possible thrusts for the project, early proof-of-principle
applications, and similar efforts upon which the project could be modeled.


  We study the problem of a target search by a Brownian particle subject to
stochastic resetting to a pair of sites. The mean search time is minimized by
an optimal resetting rate which does not vary smoothly, in contrast with the
well-known single site case, but exhibits a discontinuous transition as the
position of one resetting site is varied while keeping the initial position of
the particle fixed, or vice-versa. The discontinuity vanishes at a "liquid-gas"
critical point in position space. This critical point exists provided that the
relative weight $m$ of the further site is comprised in the interval
$[2.9028...,8.5603...]$. When the initial position follows the resetting point
distribution, a discontinuous transition also exists for the optimal rate as
the distance between the resetting points is varied, provided that $m$ exceeds
the critical value $m_c=6.6008...$ This setup can be mapped onto an
intermittent search problem with switching diffusion coefficients and
represents a minimal model for the study of distributed resetting.


  Bulk RNA sequencing (RNA-seq) has revolutionized gene expression analysis,
yet struggles with cellular heterogeneity. Traditional methods lack the ability
to examine diverse cell types simultaneously, while single-cell RNA sequencing
(scRNA-seq) is costly and complex, especially for tissues like the brain.
Recent methodologies have emerged to estimate cell type proportions from
RNA-seq samples, leveraging scRNA-seq matrices. Nevertheless, existing
deconvolution approaches face challenges, including being black-box methods
with unclear feature importance and not adequately addressing the
distributional shift between bulk and scRNA-seq data. This work presents
Sweetwater, an interpretable data-driven deconvolution model. Using an
autoencoder-based approach with real and simulated bulk samples, Sweetwater
creates a common low-dimensional embedding, minimizing platform-specific
variations. Moreover, interpretation analysis reveals Sweetwater's
effectiveness in identifying cell type marker genes, offering a transparent and
powerful tool for dissecting intricate cellular landscapes.


  Allostery, the phenomenon by which the perturbation of a molecule at one site
alters its behavior at a remote functional site, enables control over
biomolecular function. Allosteric modulation is a promising avenue for drug
discovery and is employed in the design of mechanical metamaterials. However, a
general principle of allostery, i.e. a set of quantitative and transferable
"ground rules", remained elusive. It is neither a set of structural motifs nor
intrinsic motions. Focusing on elastic network models, we here show that an
allosteric lever--a mode-coupling pattern induced by the perturbation--governs
the directional, source-to-target, allosteric communication: a structural
perturbation of an allosteric site couples the excitation of localized hard
elastic modes with concerted long range soft-mode relaxation. Perturbations of
non-allosteric sites instead couple hard and soft modes uniformly. The
allosteric response is shown to be generally non-linear and non-reciprocal, and
allows for minimal structural distortions to be efficiently transmitted to
specific changes at distant sites. Allosteric levers exist in proteins and
"pseudoproteins"--networks designed to display an allosteric response.
Interestingly, protein sequences that constitute allosteric transmission
channels are shown to be evolutionarily conserved. To illustrate how the
results may be applied in drug design, we use them to successfully predict
known allosteric sites in proteins.


  Given the increasing availability of RNA-seq data and its complex and
heterogeneous nature, there has been growing interest in applying AI/machine
learning methodologies to work with such data modalities. However, because
omics data is characterized by high dimensionality and low sample size (HDLSS),
current attempts at integrating AI in this domain require significant human
guidance and expertise to mitigate overfitting. In this work we look at how
transfer learning can be improved to learn from small RNA-seq sample sizes
without significant human interference. The strategy is to gain general prior
knowledge about a particular domain of data (e.g. RNA-seq data) by pre-training
on a general task with a large aggregate of data, then fine-tuning to various
specific, downstream target tasks in the same domain. Because previous attempts
have shown traditional transfer learning failing on HLDSS, we propose to
improve performance by using Guided Transfer Learning (GTL). Collaborating with
Robots Go Mental, the AI we deploy here not only learns good initial parameters
during pre-training, but also learns inductive biases that affect how the AI
learns downstream tasks. In this approach, we first pre-trained on recount3
data, a collection of over 400,000 mouse RNA-seq samples sourced from thousands
of individual studies. With such a large collection, patterns of expression
between the ~30,000 genes in mammalian systems were pre-determined. Such
patterns were sufficient for the pre-trained AI agent to efficiently learn new
downstream tasks involving RNA-seq datasets with very low sample sizes and
performed notably better on few-shot learning tasks compared to the same model
without pre-training.


  Cell growth and gene expression, essential elements of all living systems,
have long been the focus of biophysical interrogation. Advances in single-cell
methods have invigorated theoretical studies into these processes. However,
until recently, there was little dialog between the two areas of study. Most
theoretical models for gene regulation assumed gene activity to be oblivious to
the progression of the cell cycle between birth and division. But there are
numerous ways in which the periodic character of all cellular observables can
modulate gene expression. The molecular factors required for transcription and
translation increase in number during the cell cycle, but are also diluted due
to the continuous increase in cell volume. The replication of the genome
changes the dosage of those same cellular players but also provides competing
targets for regulatory binding. Finally, cell division reduces their number
again, and so forth. Stochasticity is inherent to all these biological
processes, manifested in fluctuations in the synthesis and degradation of new
cellular components as well as the random partitioning of molecules at each
cell division. The notion of gene expression as stationary is thus hard to
justify. In this review, we survey the emerging paradigm of cell-cycle
regulated gene expression, with an emphasis on the global expression patterns
rather than gene-specific regulation. We discuss recent experimental reports
where cell growth and gene expression were simultaneously measured in
individual cells, providing first glimpses into the coupling between the two.
While the experimental findings, not surprisingly, differ among genes and
organisms, several theoretical models have emerged that attempt to reconcile
these differences and form a unifying framework for understanding gene
expression in growing cells.


  Analytical studies of network epidemiology almost exclusively focus on the
extreme situations where the time scales of network dynamics are well separated
(longer or shorter) from that of epidemic propagation. In realistic scenarios,
however, these time scales could be similar, which has profound implications
for epidemic modeling (e.g., one can no longer reduce the dimensionality of
epidemic models). We build a theory for the critical behavior of
susceptible-infected-susceptible (SIS) epidemics in the vicinity of the
critical threshold on the activity-driven model of temporal networks. We find
that the persistence of links in the network leads to increasing recovery rates
reducing the threshold. Dynamic correlations (coming from being close to
infected nodes increases the likelihood of infection) drive the threshold in
the opposite direction. These two counteracting effects make epidemic
criticality in temporal networks a remarkably complex phenomenon.


  Metagenomics has led to significant advancements in many fields. Metagenomic
analysis commonly involves the key tasks of determining the species present in
a sample and their relative abundances. These tasks require searching large
metagenomic databases containing information on different species' genomes.
Metagenomic analysis suffers from significant data movement overhead due to
moving large amounts of low-reuse data from the storage system to the rest of
the system. In-storage processing can be a fundamental solution for reducing
data movement overhead. However, designing an in-storage processing system for
metagenomics is challenging because none of the existing approaches can be
directly implemented in storage effectively due to the hardware limitations of
modern SSDs. We propose MetaStore, the first in-storage processing system
designed to significantly reduce the data movement overhead of end-to-end
metagenomic analysis. MetaStore is enabled by our lightweight and cooperative
design that effectively leverages and orchestrates processing inside and
outside the storage system. Through our detailed analysis of the end-to-end
metagenomic analysis pipeline and careful hardware/software co-design, we
address in-storage processing challenges for metagenomics via specialized and
efficient 1) task partitioning, 2) data/computation flow coordination, 3)
storage technology-aware algorithmic optimizations, 4) light-weight in-storage
accelerators, and 5) data mapping. Our evaluation shows that MetaStore
outperforms the state-of-the-art performance- and accuracy-optimized software
metagenomic tools by 2.7-37.2$\times$ and 6.9-100.2$\times$, respectively,
while matching the accuracy of the accuracy-optimized tool. MetaStore achieves
1.5-5.1$\times$ speedup compared to the state-of-the-art metagenomic
hardware-accelerated tool, while achieving significantly higher accuracy.


  The genome sequence contains the blueprint for governing cellular processes.
While the availability of genomes has vastly increased over the last decades,
experimental annotation of the various functional, non-coding and regulatory
elements encoded in the DNA sequence remains both expensive and challenging.
This has sparked interest in unsupervised language modeling of genomic DNA, a
paradigm that has seen great success for protein sequence data. Although
various DNA language models have been proposed, evaluation tasks often differ
between individual works, and might not fully recapitulate the fundamental
challenges of genome annotation, including the length, scale and sparsity of
the data. In this study, we introduce BEND, a Benchmark for DNA language
models, featuring a collection of realistic and biologically meaningful
downstream tasks defined on the human genome. We find that embeddings from
current DNA LMs can approach performance of expert methods on some tasks, but
only capture limited information about long-range features. BEND is available
at https://github.com/frederikkemarin/BEND.


  Multimodal single-cell technologies enable the simultaneous collection of
diverse data types from individual cells, enhancing our understanding of
cellular states. However, the integration of these datatypes and modeling the
interrelationships between modalities presents substantial computational and
analytical challenges in disease biomarker detection and drug discovery.
Established practices rely on isolated methodologies to investigate individual
molecular aspects separately, often resulting in inaccurate analyses. To
address these obstacles, distinct Machine Learning Techniques are leveraged,
each of its own kind to model the co-variation of DNA to RNA, and finally to
surface proteins in single cells during hematopoietic stem cell development,
which simplifies understanding of underlying cellular mechanisms and immune
responses. Experiments conducted on a curated subset of a 300,000-cell time
course dataset, highlights the exceptional performance of Echo State Networks,
boasting a remarkable state-of-the-art correlation score of 0.94 and 0.895 on
Multi-omic and CiteSeq datasets. Beyond the confines of this study, these
findings hold promise for advancing comprehension of cellular differentiation
and function, leveraging the potential of Machine Learning.


  Least squares estimation of phylogenies is an established family of methods
with good statistical properties. State-of-the-art least squares phylogenetic
estimation proceeds by first estimating a distance matrix, which is then used
to determine the phylogeny by minimizing a squared-error loss function. Here,
we develop a method for least squares phylogenetic inference that does not rely
on a pre-estimated distance matrix. Our approach allows us to circumvent the
typical need to first estimate a distance matrix by forming a new loss function
inspired by the phylogenetic likelihood score function; in this manner,
inference is not based on a summary statistic of the sequence data, but
directly on the sequence data itself. We use a Jukes-Cantor substitution model
to show that our method leads to improvements over ordinary least squares
phylogenetic inference, and is even observed to rival maximum likelihood
estimation in terms of topology estimation efficiency. Using a Kimura
2-parameter model, we show that our method also allows for estimation of the
global transition/transversion ratio simultaneously with the phylogeny and its
branch lengths. This is impossible to accomplish with any other distance-based
method as far as we know. Our developments pave the way for more optimal
phylogenetic inference under the least squares framework, particularly in
settings under which likelihood-based inference is infeasible, including when
one desires to build a phylogeny based on information provided by only a subset
of all possible nucleotide substitutions such as synonymous or non-synonymous
substitutions.


  We propose HydraScreen, a deep-learning approach that aims to provide a
framework for more robust machine-learning-accelerated drug discovery.
HydraScreen utilizes a state-of-the-art 3D convolutional neural network,
designed for the effective representation of molecular structures and
interactions in protein-ligand binding. We design an end-to-end pipeline for
high-throughput screening and lead optimization, targeting applications in
structure-based drug design. We assess our approach using established public
benchmarks based on the CASF 2016 core set, achieving top-tier results in
affinity and pose prediction (Pearson's r = 0.86, RMSE = 1.15, Top-1 = 0.95).
Furthermore, we utilize a novel interaction profiling approach to identify
potential biases in the model and dataset to boost interpretability and support
the unbiased nature of our method. Finally, we showcase HydraScreen's capacity
to generalize across unseen proteins and ligands, offering directions for
future development of robust machine learning scoring functions. HydraScreen
(accessible at https://hydrascreen.ro5.ai) provides a user-friendly GUI and a
public API, facilitating easy assessment of individual protein-ligand
complexes.


  In the past decade, topological data analysis (TDA) has emerged as a powerful
approach in data science. The main technique in TDA is persistent homology,
which tracks topological invariants over the filtration of point cloud data
using algebraic topology. Although knot theory and related subjects are a focus
of study in mathematics, their success in practical applications is quite
limited due to the lack of localization and quantization. We address these
challenges by introducing knot data analysis (KDA), a new paradigm that
incorporating curve segmentation and multiscale analysis into the Gauss link
integral. The resulting multiscale Gauss link integral (mGLI) recovers the
global topological properties of knots and links at an appropriate scale but
offers multiscale feature vectors to capture the local structures and
connectivities of each curve segment at various scales. The proposed mGLI
significantly outperforms other state-of-the-art methods in benchmark protein
flexibility analysis, including earlier persistent homology-based methods. Our
approach enables the integration of artificial intelligence (AI) and KDA for
general curve-like objects and data.


  Splice sites play a crucial role in gene expression, and accurate prediction
of these sites in DNA sequences is essential for diagnosing and treating
genetic disorders. We address the challenge of splice site prediction by
introducing DeepDeCode, an attention-based deep learning sequence model to
capture the long-term dependencies in the nucleotides in DNA sequences. We
further propose using visualization techniques for accurate identification of
sequence motifs, which enhance the interpretability and trustworthiness of
DeepDeCode. We compare DeepDeCode to other state-of-the-art methods for splice
site prediction and demonstrate its accuracy, explainability and efficiency.
Given the results of our methodology, we expect that it can used for healthcare
applications to reason about genomic processes and be extended to discover new
splice sites and genomic regulatory elements.


  Phylogenies depicting the evolutionary history of genetically heterogeneous
subpopulations of cells from the same cancer i.e., cancer phylogenies, provide
useful insights about cancer development and inform treatment. Cancer
phylogenies can be reconstructed using data obtained from bulk DNA sequencing
of multiple tissue samples from the same cancer. We introduce Orchard, a fast
algorithm that reconstructs cancer phylogenies using point mutations detected
in bulk DNA sequencing data. Orchard constructs cancer phylogenies
progressively, one point mutation at a time, ultimately sampling complete
phylogenies from a posterior distribution implied by the bulk DNA data. Orchard
reconstructs more plausible phylogenies than state-of-the-art cancer phylogeny
reconstruction methods on 90 simulated cancers and 14 B-progenitor acute
lymphoblastic leukemias (B-ALLs). These results demonstrate that Orchard
accurately reconstructs cancer phylogenies with up to 300 mutations. We then
introduce a simple graph based clustering algorithm that uses a reconstructed
phylogeny to infer unique groups of mutations i.e., mutation clusters, that
characterize the genetic differences between cancer cell populations, and show
that this approach is competitive with state-of-the-art mutation clustering
methods.


  With the rapid advancement of computational techniques, Molecular Dynamics
(MD) simulations have emerged as powerful tools in biomedical research,
enabling in-depth investigations of biological systems at the atomic level.
Among the diverse range of simulation software available, LAMMPS (Large-scale
Atomic/Molecular Massively Parallel Simulator) has gained significant
recognition for its versatility, scalability, and extensive range of
functionalities. This literature review aims to provide a comprehensive
overview of the utilization of LAMMPS in the field of biomedical applications.
This review begins by outlining the fundamental principles of MD simulations
and highlighting the unique features of LAMMPS that make it suitable for
biomedical research. Subsequently, a survey of the literature is conducted to
identify key studies that have employed LAMMPS in various biomedical contexts,
such as protein folding, drug design, biomaterials, and cellular processes. The
reviewed studies demonstrate the remarkable contributions of LAMMPS in
understanding the behavior of biological macromolecules, investigating
drug-protein interactions, elucidating the mechanical properties of
biomaterials, and studying cellular processes at the molecular level.
Additionally, this review explores the integration of LAMMPS with other
computational tools and experimental techniques, showcasing its potential for
synergistic investigations that bridge the gap between theory and experiment.
Moreover, this review discusses the challenges and limitations associated with
using LAMMPS in biomedical simulations, including the parameterization of force
fields, system size limitations, and computational efficiency. Strategies
employed by researchers to mitigate these challenges are presented, along with
potential future directions for enhancing LAMMPS capabilities in the biomedical
field.


  Filamentous viruses like influenza and torovirus often display systematic
bends and arcs of mysterious physical origin. We propose that such viruses
undergo an instability from a cylindrically symmetric to a toroidally curved
state. This "toro-elastic" state emerges via a spontaneous symmetry breaking
under prestress, induced via short range spike protein interactions and
magnified by the filament's surface topography. Once surface stresses become
sufficiently large, the filament buckles and the toroidal, curved state
constitutes a soft mode that can propagate through the filament's material
frame around a "mexican-hat" potential. In the mucus of our airways, glycan
chains are omnipresent that influenza's spike proteins can bind to and cut. We
show that when coupled to such a non-equilibrium chemical reaction, the curved
toro-elastic state can attain a spontaneous rotation for sufficiently strong
enzymatic activity, leading to a whole body reshaping propulsion similar to --
but different from -- eukaryotic flagella and spirochetes.


  Microorganisms can create engineered materials with exquisite structures and
living functionalities. Although synthetic biology tools to genetically
manipulate microorganisms continue to expand, the bottom-up rational design of
engineered living materials still relies on prior knowledge of
genotype-phenotype links for the function of interest. Here, we utilize a
high-throughput directed evolution platform to enhance the fitness of whole
microorganisms under selection pressure and identify novel genetic pathways to
program the functionalities of engineered living materials. Using
Komagataeibacter sucrofermentans as a model cellulose-producing microorganism,
we show that our droplet-based microfluidic platform enables the directed
evolution of these bacteria towards a small number of cellulose overproducers
from an initial pool of 40'000 random mutants. Sequencing of the evolved
strains reveals an unexpected link between the cellulose-forming ability of the
bacteria and a gene encoding a protease complex responsible for protein
turnover in the cell. The ability to enhance the fitness of microorganisms
towards specific phenotypes and to discover new genotype-phenotype links makes
this high-throughput directed evolution platform a promising tool for the
development of the next generation of engineered living materials.


  Diffusion generative models have emerged as a powerful framework for
addressing problems in structural biology and structure-based drug design.
These models operate directly on 3D molecular structures. Due to the
unfavorable scaling of graph neural networks (GNNs) with graph size as well as
the relatively slow inference speeds inherent to diffusion models, many
existing molecular diffusion models rely on coarse-grained representations of
protein structure to make training and inference feasible. However, such
coarse-grained representations discard essential information for modeling
molecular interactions and impair the quality of generated structures. In this
work, we present a novel GNN-based architecture for learning latent
representations of molecular structure. When trained end-to-end with a
diffusion model for de novo ligand design, our model achieves comparable
performance to one with an all-atom protein representation while exhibiting a
3-fold reduction in inference time.


  Chemical reactions can be understood as transformations of multigraphs
(molecules) that preserve vertex labels (atoms) and degrees (sums of bonding
and non-bonding electrons), thereby implying the atom-atom map of a reaction.
The corresponding reaction mechanism is often described by an electron pushing
diagram that explains the transformation by consecutive local relocations of
invidudal edges (electron pairs). Here, we show that every degree-preserving
map between multigraphs, and thus every atom-atom map, can be generated by
cyclic electron pushing. Moreover, it is always possible to decompose such an
explanation into electron pushing diagrams involving only four electron pairs.
This in turn implies that every reaction can be decomposed into a sequence of
elementary reactions that involve at most two educt molecules and two product
molecules. Hence, the requirement of a mechanistic explantion in terms of
electron pushing and small imaginary transition states does not impose a
combinatorial constraint on the feasibility of hypothetical chemical reactions.


  A novel method combining the ensemble refinement by maximum entropy principle
and the force field fitting approach is presented. Its formulation allows to
continuously interpolate in between these two methods, which can thus be
interpreted as two limiting cases. A cross-validation procedure enables to
correctly assess the relative weight of both of them, distinguishing scenarios
where the combined approach is meaningful from those in which either ensemble
refinement or force field fitting separately prevails. The efficacy of their
combination is examined for a realistic case study of RNA oligomers. Within the
new scheme, molecular dynamics simulations are integrated with experimental
data provided by nuclear-magnetic-resonance measures. We show that force field
corrections are in general superior when applied to the appropriate force field
terms, but are automatically discarded by the method when applied to
inappropriate force field terms.


  The origin of life must have been preceded by Darwin-like evolutionary
dynamics that could propagate it. How did that adaptive dynamics arise? And
from what prebiotic molecules? Using evolutionary invasion analysis, we develop
a universal framework for describing any origin story for evolutionary
dynamics. We find that cooperative autocatalysts, i.e. autocatalysts whose
per-unit reproductive rate grows as their population increases, have the
special property of being able to cross a barrier that separates their initial
degradation-dominated state from a growth-dominated state with evolutionary
dynamics. For some model parameters, this leap to persistent propagation is
likely, not rare. We apply this analysis to the Foldcat Mechanism, wherein
peptides fold and help catalyze the elongation of each other. Foldcats are
found to have cooperative autocatalysis and be capable of emergent evolutionary
dynamics.


  The COVID-19 pandemic has highlighted the need to upgrade systems for
infectious disease surveillance and forecasting and modeling of the spread of
infection, both of which inform evidence-based public health guidance and
policies. Here, we discuss requirements for an effective surveillance system to
support decision making during a pandemic, drawing on the lessons of COVID-19
in the U.S., while looking to jurisdictions in the U.S. and beyond to learn
lessons about the value of specific data types. In this report, we define the
range of decisions for which surveillance data are required, the data elements
needed to inform these decisions and to calibrate inputs and outputs of
transmission-dynamic models, and the types of data needed to inform decisions
by state, territorial, local, and tribal health authorities. We define actions
needed to ensure that such data will be available and consider the contribution
of such efforts to improving health equity.


  With the emergence of advanced spatial transcriptomic technologies, there has
been a surge in research papers dedicated to analyzing spatial transcriptomics
data, resulting in significant contributions to our understanding of biology.
The initial stage of downstream analysis of spatial transcriptomic data has
centered on identifying spatially variable genes (SVGs) or genes expressed with
specific spatial patterns across the tissue. SVG detection is an important task
since many downstream analyses depend on these selected SVGs. Over the past few
years, a plethora of new methods have been proposed for the detection of SVGs,
accompanied by numerous innovative concepts and discussions. This article
provides a selective review of methods and their practical implementations,
offering valuable insights into the current literature in this field.


  What determines biodiversity in nature is a prominent issue in ecology,
especially in biotic resource systems that are typically devoid of
cross-feeding. Here, we show that by incorporating pairwise encounters among
consumer individuals within the same species, a multitude of consumer species
can self-organize to coexist in a well-mixed system with one or a few biotic
resource species. The coexistence modes can manifest as either stable steady
states or self-organized oscillations. Importantly, all coexistence states are
robust to stochasticity, whether employing the stochastic simulation algorithm
or individual-based modeling. Our model quantitatively illustrates species
distribution patterns across a wide range of ecological communities and can be
broadly used to explain biodiversity in many biotic resource systems.


  This paper seeks to study the evolution of the COVID-19 pandemic based on
daily published data from Worldometer website, using a time-dependent SIR
model. Our findings indicate that this model fits well such data, for different
chosen periods and different regions. This well-known model, consisting of
three disjoint compartments, susceptible , infected , and removed , depends in
our case on two time dependent parameters, the infection rate $\beta(t)$ and
the removal rate $\rho(t)$. After deriving the model, we prove the local
exponential behavior of the number of infected people, be it growth or decay.
Furthermore, we extract a time dependent replacement factor $\sigma_s(t)
={\beta(t)}s(t)/{\rho(t) }$, where $s(t)$ is the ratio of susceptible people at
time $t$. In addition, $i(t)$ and $r(t)$ are respectively the ratios of
infected and removed people, based on a population of size $N$, usually assumed
to be constant. Besides these theoretical results, the report provides
simulations on the daily data obtained for Germany, Italy, and the entire
World, as collected from Worldometer over the period stretching from April 2020
to June 2022. The computational model consists of the estimation of $\beta(t)$,
$\rho(t)$ and $s(t)$ based on the time-dependent SIR model. The validation of
our approach is demonstrated by comparing the profiles of the collected $i(t),
r(t)$ data and those obtained from the SIR model with the approximated
parameters. We also consider matching the data with a constant-coefficient SIR
model, which seems to be working only for short periods. Thus, such model helps
understanding and predicting the evolution of the pandemics for short periods
of time where no radical change occurs.


  Periodic travelling waves (PTW) are a common solution type of partial
differential equations. Such models exhibit multistability of PTWs, typically
visualised through the Busse balloon, and parameter changes typically lead to a
cascade of wavelength changes through the Busse balloon. In the past, the
stability boundaries of the Busse balloon have been used to predict such
wavelength changes. Here, motivated by anecdotal evidence from previous work,
we provide compelling evidence that the Busse balloon provides insufficient
information to predict wavelength changes due to a delayed loss of stability
phenomenon. Using two different reaction-advection-diffusion systems, we relate
the delay that occurs between the crossing of a stability boundary in the Busse
balloon and the occurrence of a wavelength change to features of the essential
spectrum of the destabilised PTW. This leads to a predictive framework that can
estimate the order of magnitude of such a time delay, which provides a novel
``early warning sign'' for pattern destabilization. We illustrate the
implementation of the predictive framework to predict under what conditions a
wavelength change of a PTW occurs.


  Density-dependent population growth is a feature of large carnivores like
wolves ($\textit{Canis lupus}$), with mechanisms typically attributed to
resource (e.g. prey) limitation. Such mechanisms are local phenomena and rely
on individuals having access to information, such as prey availability at their
location. Using over four decades of wolf population and range expansion data
from Wisconsin (USA) wolves, we found that the population not only exhibited
density dependence locally but also at landscape scale. Superficially, one may
consider space as yet another limiting resource to explain landscape-scale
density dependence. However, this view poses an information puzzle: most
individuals do not have access to global information such as range-wide habitat
availability as they would for local prey availability. How would the
population "know" when to slow their range expansion? To understand observed
large-scale spatial density dependence, we propose a reaction-diffusion model,
first introduced by Fisher and Kolmogorov, with a "travelling wave" solution,
wherein the population expands from a core range that quickly achieves local
carrying capacity. Early-stage acceleration and later-stage deceleration of
population growth can be explained by early elongation of an expanding frontier
and a later collision of the expanding frontier with a habitat boundary. Such a
process does not require individuals to have global density information. We
illustrate our proposal with simulations and spatial visualizations of wolf
recolonization in the western Great Lakes region over time relative to habitat
suitability. We further synthesize previous studies on wolf habitat selection
in the western Great Lakes region and argue that the habitat boundary appeared
to be driven by spatial variation in mortality, likely associated with human
use of the landscape.


  The epidemiology of pandemics is classically viewed using geographical and
political borders; however, these artificial divisions can result in a
misunderstanding of the current epidemiological state within a given region. To
improve upon current methods, we propose a clustering algorithm which is
capable of recasting regions into well-mixed clusters such that they have a
high level of interconnection while minimizing the external flow of the
population towards other clusters. Moreover, we analyze and identify so called
core clusters, clusters that retain their features over time (temporally
stable) and independent of the presence or absence of policy measures. In order
to demonstrate the capabilities of this algorithm, we use US county-level
cellular mobility data to divide the country into such clusters. Herein, we
show a more granular spread of SARS-CoV-2 throughout the first weeks of the
pandemic. Moreover, we are able to identify areas (groups of counties) that
were experiencing above average levels of transmission within a state, as well
as pan-state areas (clusters overlapping more than one state) with very similar
disease spread. Therefore, our method enables policymakers to make more
informed decisions on the use of public health interventions within their
jurisdiction, as well as guide collaboration with surrounding regions to
benefit the general population in controlling the spread of communicable
diseases.


  Mutualistic interactions, where species interact to obtain mutual benefits,
constitute an essential component of natural ecosystems. The use of ecological
networks to represent the species and their ecological interactions allows the
study of structural and dynamic patterns common to different ecosystems.
However, by neglecting the temporal dimension of mutualistic communities,
relevant insights into the organization and functioning of natural ecosystems
can be lost. Therefore, it is crucial to incorporate empirical phenology -- the
cycles of species' activity within a season -- to fully understand the effects
of temporal variability on network architecture. In this paper, by using two
empirical datasets together with a set of synthetic models, we propose a
framework to characterize phenology on ecological networks and assess the
effect of temporal variability. Analyses reveal that non-trivial information is
missed when portraying the network of interactions as static, which leads to
overestimating the value of fundamental structural features. We discuss the
implications of our findings for mutualistic relationships and intra-guild
competition for common resources. We show that recorded interactions and
species' activity duration are pivotal factors in accurately replicating
observed patterns within mutualistic communities. Furthermore, our exploration
of synthetic models underscores the system-specific character of the mechanisms
driving phenology, increasing our understanding of the complexities of natural
ecosystems.


  Advances in high-throughput sequencing technology have led to significant
progress in measuring gene expressions at the single-cell level. The amount of
publicly available single-cell RNA-seq (scRNA-seq) data is already surpassing
50M records for humans with each record measuring 20,000 genes. This highlights
the need for unsupervised representation learning to fully ingest these data,
yet classical transformer architectures are prohibitive to train on such data
in terms of both computation and memory. To address this challenge, we propose
a novel asymmetric encoder-decoder transformer for scRNA-seq data, called
xTrimoGene$^\alpha$ (or xTrimoGene for short), which leverages the sparse
characteristic of the data to scale up the pre-training. This scalable design
of xTrimoGene reduces FLOPs by one to two orders of magnitude compared to
classical transformers while maintaining high accuracy, enabling us to train
the largest transformer models over the largest scRNA-seq dataset today. Our
experiments also show that the performance of xTrimoGene improves as we scale
up the model sizes, and it also leads to SOTA performance over various
downstream tasks, such as cell type annotation, perturb-seq effect prediction,
and drug combination prediction. xTrimoGene model is now available for use as a
service via the following link: https://api.biomap.com/xTrimoGene/apply.


  Molecular docking, a key technique in structure-based drug design, plays
pivotal roles in protein-ligand interaction modeling, hit identification and
optimization, in which accurate prediction of protein-ligand binding mode is
essential. Conventional docking approaches perform well in redocking tasks with
known protein binding pocket conformation in the complex state. However, in
real-world docking scenario without knowing the protein binding conformation
for a new ligand, accurately modeling the binding complex structure remains
challenging as flexible docking is computationally expensive and inaccurate.
Typical deep learning-based docking methods do not explicitly consider protein
side chain conformations and fail to ensure the physical plausibility and
detailed atomic interactions. In this study, we present DiffBindFR, a full-atom
diffusion-based flexible docking model that operates over the product space of
ligand overall movements and flexibility and pocket side chain torsion changes.
We show that DiffBindFR has higher accuracy in producing native-like binding
structures with physically plausible and detailed interactions than available
docking methods. Furthermore, in the Apo and AlphaFold2 modeled structures,
DiffBindFR demonstrates superior advantages in accurate ligand binding pose and
protein binding conformation prediction, making it suitable for Apo and
AlphaFold2 structure-based drug design. DiffBindFR provides a powerful flexible
docking tool for modeling accurate protein-ligand binding structures.