  The prediction of protein-ligand binding affinity is of great significance
for discovering lead compounds in drug research. Facing this challenging task,
most existing prediction methods rely on the topological and/or spatial
structure of molecules and the local interactions while ignoring the
multi-level inter-molecular interactions between proteins and ligands, which
often lead to sub-optimal performance. To solve this issue, we propose a novel
global-local interaction (GLI) framework to predict protein-ligand binding
affinity. In particular, our GLI framework considers the inter-molecular
interactions between proteins and ligands, which involve not only the
high-energy short-range interactions between closed atoms but also the
low-energy long-range interactions between non-bonded atoms. For each pair of
protein and ligand, our GLI embeds the long-range interactions globally and
aggregates local short-range interactions, respectively. Such a joint
global-local interaction modeling strategy helps to improve prediction
accuracy, and the whole framework is compatible with various neural
network-based modules. Experiments demonstrate that our GLI framework
outperforms state-of-the-art methods with simple neural network architectures
and moderate computational costs.


  Single-cell RNA sequencing (scRNA-seq) provides a high throughput,
quantitative and unbiased framework for scientists in many research fields to
identify and characterize cell types within heterogeneous cell populations from
various tissues. However, scRNA-seq based identification of discrete cell-types
is still labor intensive and depends on prior molecular knowledge. Artificial
intelligence has provided faster, more accurate, and user-friendly approaches
for cell-type identification. In this review, we discuss recent advances in
cell-type identification methods using artificial intelligence techniques based
on single-cell and single-nucleus RNA sequencing data in vision science.


  Local and global inference methods have been developed to infer structural
contacts from multiple sequence alignments of homologous proteins. They rely on
correlations in amino-acid usage at contacting sites. Because homologous
proteins share a common ancestry, their sequences also feature phylogenetic
correlations, which can impair contact inference. We investigate this effect by
generating controlled synthetic data from a minimal model where the importance
of contacts and of phylogeny can be tuned. We demonstrate that global inference
methods, specifically Potts models, are more resilient to phylogenetic
correlations than local methods, based on covariance or mutual information.
This holds whether or not phylogenetic corrections are used, and may explain
the success of global methods. We analyse the roles of selection strength and
of phylogenetic relatedness. We show that sites that mutate early in the
phylogeny yield false positive contacts. We consider natural data and realistic
synthetic data, and our findings generalise to these cases. Our results
highlight the impact of phylogeny on contact prediction from protein sequences
and illustrate the interplay between the rich structure of biological data and
inference.


  A key goal in studies of ecology and evolution is understanding the causes of
phenotypic diversity in nature. Most traits of interest, such as those relating
to morphology, life-history, immunity and behaviour are quantitative, and
phenotypic variation is driven by the cumulative effects of genetic and
environmental variation. The field of quantitative genetics aims to quantify
the additive genetic component of this trait variance (i.e. the
"heritability"), often with the underlying assumption that trait variance is
driven by many loci of infinitesimal effects throughout the genome. This
approach allows us to understand the evolutionary potential of natural
populations and can be extended to examine the genetic covariation with fitness
to predict responses to selection. Therefore, quantitative genetic studies are
fundamental to understanding evolution in the wild. Over the last two decades,
there has been a wealth of studies investigating trait heritabilities and
genetic correlations, but these were initially limited to long-term studies of
pedigreed populations or common-garden experiments. However, genomic
technologies have since allowed quantitative genetic studies in a more diverse
range of wild systems and has increased the opportunities for addressing
outstanding questions in ecology and evolution. In particular, genomic studies
can uncover the genetic basis of fitness-related quantitative traits, allowing
a better understanding of their evolutionary dynamics. We organised this
special issue to highlight new work and review recent advances at the cutting
edge of "Wild Quantitative Genomics". In this Editorial, we will present some
history of wild quantitative genetic and genomic studies, before discussing the
main themes in the papers published in this special issue and highlighting the
future outlook of this dynamic field.


  We demonstrate the applicability of model-agnostic algorithms for
meta-learning, specifically Reptile, to GNN models in molecular regression
tasks. Using meta-learning we are able to learn new chemical prediction tasks
with only a few model updates, as compared to using randomly initialized GNNs
which require learning each regression task from scratch. We experimentally
show that GNN layer expressivity is correlated to improved meta-learning.
Additionally, we also experiment with GNN emsembles which yield best
performance and rapid convergence for k-shot learning.


  As antibiotic-resistant bacterial strains are rapidly spreading worldwide,
infections caused by these strains are emerging as a global crisis causing the
death of millions of people every year. Antimicrobial Peptides (AMPs) are one
of the candidates to tackle this problem because of their potential diversity,
and ability to favorably modulate the host immune response. However,
large-scale screening of new AMP candidates is expensive, time-consuming, and
now affordable in developing countries, which need the treatments the most. In
this work, we propose a novel active machine learning-based framework that
statistically minimizes the number of wet-lab experiments needed to design new
AMPs, while ensuring a high diversity and novelty of generated AMPs sequences,
in multi-rounds of wet-lab AMP screening settings. Combining recurrent neural
network models and a graph-based filter (GraphCC), our proposed approach
delivers novel and diverse candidates and demonstrates better performances
according to our defined metrics.


  Virus infection is one of the major diseases that seriously threaten human
health. To meet the growing demand for mining and sharing data resources
related to antiviral drugs and to accelerate the design and discovery of new
antiviral drugs, we presented an open-access antiviral drug resource and
machine learning platform (VDDB), which, to the best of our knowledge, is the
first comprehensive dedicated resource for experimentally verified potential
drugs/molecules based on manually curated data. Currently, VDDB highlights 848
clinical vaccines, 199 clinical antibodies, as well as over 710,000 small
molecules targeting 39 medically important viruses including SARS-CoV-2.
Furthermore, VDDB stores approximately 3 million records of pharmacological
data for these collected potential antiviral drugs/molecules, involving 314
cell infection-based phenotypic and 234 target-based genotypic assays. Based on
these annotated pharmacological data, VDDB allows users to browse, search and
download reliable information about these collects for various viruses of
interest. In particular, VDDB also integrates 57 cell infection- and 117
target-based associated high-accuracy machine learning models to support
various antivirals identification-related tasks, such as compound activity
prediction, virtual screening, drug repositioning and target fishing. VDDB is
freely accessible at http://vddb.idruglab.cn.


  A novel framework has recently been proposed for designing the molecular
structure of chemical compounds with a desired chemical property, where design
of novel drugs is an important topic in bioinformatics and chemo-informatics.
The framework infers a desired chemical graph by solving a mixed integer linear
program (MILP) that simulates the computation process of a feature function
defined by a two-layered model on chemical graphs and a prediction function
constructed by a machine learning method. A set of graph theoretical
descriptors in the feature function plays a key role to derive a compact
formulation of such an MILP. To improve the learning performance of prediction
functions in the framework maintaining the compactness of the MILP, this paper
utilizes the product of two of those descriptors as a new descriptor and then
designs a method of reducing the number of descriptors. The results of our
computational experiments suggest that the proposed method improved the
learning performance for many chemical properties and can infer a chemical
structure with up to 50 non-hydrogen atoms.


  While Neutral Theory famously describes the number of discrete genetic
differences in populations, we consider the number of genetic backgrounds under
which such differences are observed - setting limits to the generalizability of
their effects. This allow us to determine which population structures and
diversity rates have maximal effect generalization across (1) environmental and
(2) genetic variation, and to demonstrate that they correspond asymptotically
to those of populations under (1) natural selection and (2) drift. At the same
time, these results suggest distinct limits to the predictability of fitness
and evolution across evolutionary regimes. We employ both broad time,
large-scale genome sequencing datasets (including whole-genome autocorrelation
calculations), and fine time-scale barcoding experiments.


  Surges that have been observed at different periods in the number of COVID-19
cases are associated with the emergence of multiple SARS-CoV-2 (Severe Acute
Respiratory Virus) variants. The design of methods to support laboratory
detection are crucial in the monitoring of these variants. Hence, in this
paper, we develop a semi-automated method to design both forward and reverse
primer sets to detect SARS-CoV-2 variants. To proceed, we train deep
Convolution Neural Networks (CNNs) to classify labelled SARS-CoV-2 variants and
identify partial genomic features needed for the forward and reverse Polymerase
Chain Reaction (PCR) primer design. Our proposed approach supplements existing
ones while promoting the emerging concept of neural network assisted primer
design for PCR. Our CNN model was trained using a database of SARS-CoV-2
full-length genomes from GISAID and tested on a separate dataset from NCBI,
with 98\% accuracy for the classification of variants. This result is based on
the development of three different methods of feature extraction, and the
selected primer sequences for each SARS-CoV-2 variant detection (except
Omicron) were present in more than 95 \% of sequences in an independent set of
5000 same variant sequences, and below 5 \% in other independent datasets with
5000 sequences of each variant. In total, we obtain 22 forward and reverse
primer pairs with flexible length sizes (18-25 base pairs) with an expected
amplicon length ranging between 42 and 3322 nucleotides. Besides the feature
appearance, in-silico primer checks confirmed that the identified primer pairs
are suitable for accurate SARS-CoV-2 variant detection by means of PCR tests.


  While the distribution patterns of cold-water corals, such as Paragorgia
arborea, have received increasing attention in recent studies, little is known
about their in situ activity patterns. In this paper, we examine polyp activity
in P. arborea using machine learning techniques to analyze high-resolution time
series data and photographs obtained from an autonomous lander cluster deployed
in the Stjernsund, Norway. An interactive illustration of the models derived in
this paper is provided online as supplementary material. We find that the best
predictor of the degree of extension of the coral polyps is current direction
with a lag of three hours. Other variables that are not directly associated
with water currents, such as temperature and salinity, offer much less
information concerning polyp activity. Interestingly, the degree of polyp
extension can be predicted more reliably by sampling the laminar flows in the
water column above the measurement site than by sampling the more turbulent
flows in the direct vicinity of the corals. Our results show that the activity
patterns of the P. arborea polyps are governed by the strong tidal current
regime of the Stjernsund. It appears that P. arborea does not react to shorter
changes in the ambient current regime but instead adjusts its behavior in
accordance with the large-scale pattern of the tidal cycle itself in order to
optimize nutrient uptake.


  In this article, I put forward the idea that the neoplastic process (NP) has
deep evolutionary roots and make specific predictions about the connection
between cancer and the formation of the first embryo, which allowed for the
evolutionary radiation of metazoans. My main hypothesis is that the NP is at
the heart of cellular mechanisms responsible for animal morphogenesis and,
given its embryological basis, also at the center of animal evolution. It is
thus understood that NP-associated mechanisms are deeply rooted in evolutionary
history and tied to the formation of the first animal embryo. In my
consideration of these arguments, I expound on how cancer biology is perfectly
intertwined with evolutionary biology. I describe essential cellular components
of unicellular holozoans that served as a basis for the formation of the
neoplastic functional module (NFM) and its subsequent exaptation, which brought
forth two great biophysical revolutions within the first embryo. Finally, I
examine the role of Physics in the modeling of the NFM and its contribution to
morphogenesis to reveal the totipotency of the zygote.


  Feed-forward dynamics, which is well-known to have several important
implications in nonlinear dynamical systems, frequently occurs in gene
expression motifs, and has been well explored experimentally and
mathematically. However, dependency of the components of a genetic circuit upon
its host, due to the requirement for resources like ribosome, ATP,
transcription factors, tRNA, etc., and related effects are of utmost
importance, which is commonly ignored in mathematical models. In a
resource-limited environment, two apparently unconnected genes can compete for
resources for their respective expression and may exhibit indirect regulatory
connection; an emergent response thus arises in the system completely because
of resource competition. In this work, we have shown how the responses of the
feed-forward loop (FFL), a well-studied regulatory genetic motif, can be
recreated considering the resource competition in a three-gene pathway.
Exploring the genetic system with temporal as well as spatiotemporal stability
analysis, interesting transient and steady-state responses have been observed.
The genetic motifs explored in this paper show many of the characteristic
features of the conventional FFL structure, like response delay and pulse
generation. Most interestingly, in a two-dimensional cellular arrangement,
characteristic pattern formation under a concentration gradient of input
signals have also been observed. This study pinpoints a larger area of research
and exploration in synthetic and cellular systems, which will reveal novel
controlling ideas and unique behavioral changes in the system for its context
dependencies.


  While noise is generally associated with uncertainties and often has a
negative connotation in engineering, living organisms have evolved to adapt to
(and even exploit) such uncertainty to ensure the survival of a species or
implement certain functions that would have been difficult or even impossible
otherwise. In this article, we review the role and impact of noise in systems
and synthetic biology, with a particular emphasis on its role in the genetic
control of biological systems, an area we refer to as Cybergenetics. The main
modeling paradigm is that of stochastic reaction networks, whose applicability
goes beyond biology, as these networks can represent any population dynamics
system, including ecological, epidemiological, and opinion dynamics networks.
We review different ways to mathematically represent these systems, and we
notably argue that the concept of ergodicity presents a particularly suitable
way to characterize their stability. We then discuss noise-induced properties
and show that noise can be both an asset and a nuisance in this setting.
Finally, we discuss recent results on (stochastic) Cybergenetics and explore
their relationships to noise. Along the way, we detail the different technical
and biological constraints that need to be respected when designing synthetic
biological circuits. Finally, we discuss the concepts, problems, and solutions
exposed in the article; raise criticisms and concerns about current ideas and
approaches; suggest current (open) problems with potential solutions; and
provide some ideas for future research directions.


  Switch-like responses arising from bistability have been linked to cell
signaling processes and memory. Revealing the shape and properties of the set
of parameters that lead to bistability is necessary to understand the
underlying biological mechanisms, but is a complex mathematical problem. We
present an efficient approach to determine a basic topological property of the
parameter region of multistationary, namely whether it is connected or not. The
connectivity of this region can be interpreted in terms of the biological
mechanisms underlying bistability and the switch-like patterns that the system
can create.
  We provide an algorithm to assert that the parameter region of
multistationarity is connected, targeting reaction networks with mass-action
kinetics. We show that this is the case for numerous relevant cell signaling
motifs, previously described to exhibit bistability. However, we show that for
a motif displaying a phosphorylation cycle with allosteric enzyme regulation,
the region of multistationarity has two distinct connected components,
corresponding to two different, but symmetric, biological mechanisms. The
method relies on linear programming and bypasses the expensive computational
cost of direct and generic approaches to study parametric polynomial systems.
This characteristic makes it suitable for mass-screening of reaction networks.


  The adaptive immune system constantly remodels its lymphocyte repertoire for
better protection against future pathogens. Its ability to improve antigen
recognition on the fly relies on somatic mutation and selective expansion of B
lymphocytes expressing high-affinity antigen receptors. However, this Darwinian
process inside an individual appears ineffective, hitting a modest ceiling of
antigen-binding affinity. Experiment began to reveal that evolving B cells
physically extract antigens from presenting cells and that the extraction level
dictates clonal expansion; this challenges the prevailing assumption that the
equilibrium constant of receptor-antigen binding determines selective
advantage. We present a theoretical framework to explore whether, and how, such
tug-of-war antigen extraction impacts the quality and diversity of an evolved B
cell repertoire. We find that the apparent ineffectiveness of clonal selection
can be a direct consequence of the non-equilibrium nature of antigen
recognition. Our theory predicts that the physical strength of antigen
tethering under tugging forces sets the affinity ceiling. Meanwhile, the model
showed that, intriguingly, cells can use force variability to diversify binding
phenotype without compromising fitness, thus remaining plastic under resource
constraint. These results suggest that active probing of receptor quality via a
molecular tug of war during antigen recognition limit the potency of response
to the current antigen, but confer adaptive benefit for protection against
future variants. Importantly, a saddle point in the fitness landscape of B cell
phenotype evolution emerges from the tug-of-war setting, which rationalizes
multiple key phenomenology and puts forward a role of active physical dynamics
in immune adaptation.


  The rate at which individual bacterial cells grow depends on the
concentrations of cellular components such as ribosomes and proteins. These
concentrations continuously fluctuate over time and are inherited from mother
to daughter cells leading to correlations between the growth rates of cells
across generations. Division sizes of cells are also stochastic and correlated
across generations due to a phenomenon known as cell size regulation.
Fluctuations and correlations from both growth and division processes affect
the population dynamics of an exponentially growing culture. Here, we provide
analytic solutions for the population dynamics of cells with continuously
fluctuating growth rates coupled with a generic model of cell-size regulation.
We show that in balanced growth, the effects of growth and division processes
decouple; the population growth rate only depends on the single-cell growth
rate process, and the population cell size distribution only depends on the
model of division and cell size regulation. The population growth rate is
always higher than the average single-cell growth rate, and the difference
increases with growth rate variability and its correlation time. This
difference also sets the timescale for the population to reach its steady
state. We provide analytical solutions for oscillations in population growth
rate and traveling waves in size distribution during this approach to the
steady state.


  To better describe the spread of a disease, we extend a discrete time
stochastic SIR-type epidemic model of Tuckwell and Williams. We assume the
dependence on time of the number of daily encounters and include a parameter to
represent a possible quarantine of the infectious individuals. We provide an
analytic description of this Markovian model and investigate its dynamics. Both
a diffusion approximation and the basic reproduction number are derived.
Through several simulations, we show how the evolution of a disease is affected
by the distribution of the number of daily encounters and its dependence on
time.


  Molecular pretraining, which learns molecular representations over massive
unlabeled data, has become a prominent paradigm to solve a variety of tasks in
computational chemistry and drug discovery. Recently, prosperous progress has
been made in molecular pretraining with different molecular featurizations,
including 1D SMILES strings, 2D graphs, and 3D geometries. However, the role of
molecular featurizations with their corresponding neural architectures in
molecular pretraining remains largely unexamined. In this paper, through two
case studies -- chirality classification and aromatic ring counting -- we first
demonstrate that different featurization techniques convey chemical information
differently. In light of this observation, we propose a simple and effective
MOlecular pretraining framework with COmplementary featurizations (MOCO). MOCO
comprehensively leverages multiple featurizations that complement each other
and outperforms existing state-of-the-art models that solely relies on one or
two featurizations on a wide range of molecular property prediction tasks.


  Cryogenic electron microscopy (cryo-EM) provides a unique opportunity to
study the structural heterogeneity of biomolecules. Being able to explain this
heterogeneity with atomic models would help our understanding of their
functional mechanisms but the size and ruggedness of the structural space (the
space of atomic 3D cartesian coordinates) presents an immense challenge. Here,
we describe a heterogeneous reconstruction method based on an atomistic
representation whose deformation is reduced to a handful of collective motions
through normal mode analysis. Our implementation uses an autoencoder. The
encoder jointly estimates the amplitude of motion along the normal modes and
the 2D shift between the center of the image and the center of the molecule .
The physics-based decoder aggregates a representation of the heterogeneity
readily interpretable at the atomic level. We illustrate our method on 3
synthetic datasets corresponding to different distributions along a simulated
trajectory of adenylate kinase transitioning from its open to its closed
structures. We show for each distribution that our approach is able to
recapitulate the intermediate atomic models with atomic-level accuracy.


  The binding complexes formed by proteins and small molecule ligands are
ubiquitous and critical to life. Despite recent advancements in protein
structure prediction, existing algorithms are so far unable to systematically
predict the binding ligand structures along with their regulatory effects on
protein folding. To address this discrepancy, we present NeuralPLexer, a
computational approach that can directly predict protein-ligand complex
structures solely using protein sequence and ligand molecular graph inputs.
NeuralPLexer adopts a deep generative model to sample the 3D structures of the
binding complex and their conformational changes at an atomistic resolution.
The model is based on a diffusion process that incorporates essential
biophysical constraints and a multi-scale geometric deep learning system to
iteratively sample residue-level contact maps and all heavy-atom coordinates in
a hierarchical manner. NeuralPLexer achieves state-of-the-art performance
compared to all existing methods on benchmarks for both protein-ligand blind
docking and flexible binding site structure recovery. Moreover, owing to its
specificity in sampling both ligand-free-state and ligand-bound-state
ensembles, NeuralPLexer consistently outperforms AlphaFold2 in terms of global
protein structure accuracy on both representative structure pairs with large
conformational changes (average TM-score=0.93) and recently determined
ligand-binding proteins (average TM-score=0.89). Case studies reveal that the
predicted conformational variations are consistent with structure determination
experiments for important targets, including human KRAS$^\textrm{G12C}$,
ketol-acid reductoisomerase, and purine GPCRs. Our study suggests that a
data-driven approach can capture the structural cooperativity between proteins
and small molecules, showing promise in accelerating the design of enzymes,
drug molecules, and beyond.


  The extraction of sequence patterns from a collection of functionally linked
unlabeled DNA sequences is known as DNA motif discovery, and it is a key task
in computational biology. Several deep learning-based techniques have recently
been introduced to address this issue. However, these algorithms can not be
used in real-world situations because of the need for labeled data. Here, we
presented RL-MD, a novel reinforcement learning based approach for DNA motif
discovery task. RL-MD takes unlabelled data as input, employs a relative
information-based method to evaluate each proposed motif, and utilizes these
continuous evaluation results as the reward. The experiments show that RL-MD
can identify high-quality motifs in real-world data.


  Retrosynthetic planning aims to devise a complete multi-step synthetic route
from starting materials to a target molecule. Current strategies use a
decoupled approach of single-step retrosynthesis models and search algorithms,
taking only the product as the input to predict the reactants for each planning
step and ignoring valuable context information along the synthetic route. In
this work, we propose a novel framework that utilizes context information for
improved retrosynthetic planning. We view synthetic routes as reaction graphs
and propose to incorporate context through three principled steps: encode
molecules into embeddings, aggregate information over routes, and readout to
predict reactants. Our approach is the first attempt to utilize in-context
learning for retrosynthesis prediction in retrosynthetic planning. The entire
framework can be efficiently optimized in an end-to-end fashion and produce
more practical and accurate predictions. Comprehensive experiments demonstrate
that by fusing in the context information over routes, our model significantly
improves the performance of retrosynthetic planning over baselines that are not
context-aware, especially for long synthetic routes. Code is available at
https://github.com/SongtaoLiu0823/FusionRetro.


  Marine plankton play a crucial role in carbon storage, oxygen production,
global climate, and ecosystem function. Planktonic ecosystems are embedded in a
Lagrangian patches of water that are continuously moving, stretching, and
diluting. These processes drive inhomegeneities on a range of scales, with
implications for the integrated ecosystem properties, but are hard to
characterize. We present a theoretical framework which accounts for all these
aspects; tracking the water patch hosting a drifting ecosystem along with its
physical, environmental, and biochemical features. The model resolves patch
dilution and internal physical mixing as a function of oceanic strain and
diffusion. Ecological dynamics are parameterized by an idealized nutrient and
phytoplankton population and we specifically capture the propagation of the
biochemical spatial variances to represent within-patch heterogeneity. We find
that, depending only on the physical processes to which the water patch is
subjected, the plankton biomass response to a resource perturbation can vary
several fold. This work indicates that we must account for these processes when
interpreting and modelling marine ecosystems and provides a framework with
which to do so.


  Inverse molecular design is critical in material science and drug discovery,
where the generated molecules should satisfy certain desirable properties. In
this paper, we propose equivariant energy-guided stochastic differential
equations (EEGSDE), a flexible framework for controllable 3D molecule
generation under the guidance of an energy function in diffusion models.
Formally, we show that EEGSDE naturally exploits the geometric symmetry in 3D
molecular conformation, as long as the energy function is invariant to
orthogonal transformations. Empirically, under the guidance of designed energy
functions, EEGSDE significantly improves the baseline on QM9, in inverse
molecular design targeted to quantum properties and molecular structures.
Furthermore, EEGSDE is able to generate molecules with multiple target
properties by combining the corresponding energy functions linearly.


  We formulate and analyse a robust mathematical model of the dynamics of
gonorrhoea incorporating passive immunity and control. Our results show that
the disease-free and endemic equilibria of the model are both locally and
globally asymptotically stable. A sensitivity analysis of the model shows that
the dynamics of the model is variable and dependent on waning rate, control
parameters and interaction of the latent and infected classes. In particular,
the lower the waning rate, the more the exponential decrease in the passive
immunity but the susceptible population increases to the equilibrium and wanes
asymptotically due to the presence of the control parameters and restricted
interaction of the latent and infected classes.


  Parkinson's disease (PD) is the second most common neurodegenerative disorder
worldwide, yet there is no disease-modifying therapy up to this date. The
biological complexity underlying PD hampers the investigation of the principal
contributors to its pathogenesis. In this context, mechanistic models grounded
in molecular-level knowledge provide virtual labs to uncover the primary events
triggering PD onset and progression and suggest promising therapeutic targets.
Multiple modeling efforts in PD research have focused on the pathological role
of $\alpha$-synuclein ($\alpha$syn), a presynaptic protein that emerges from
the intricate molecular network as a crucial driver of neurodegeneration. Here,
we collect the advances in mathematical modeling of $\alpha$syn homeostasis,
focusing on aggregation and degradation pathways, and discussing potential
modeling improvements and possible implications in PD therapeutic strategy
design.


  The ability to computationally generate novel yet physically foldable protein
structures could lead to new biological discoveries and new treatments
targeting yet incurable diseases. Despite recent advances in protein structure
prediction, directly generating diverse, novel protein structures from neural
networks remains difficult. In this work, we present a new diffusion-based
generative model that designs protein backbone structures via a procedure that
mirrors the native folding process. We describe protein backbone structure as a
series of consecutive angles capturing the relative orientation of the
constituent amino acid residues, and generate new structures by denoising from
a random, unfolded state towards a stable folded structure. Not only does this
mirror how proteins biologically twist into energetically favorable
conformations, the inherent shift and rotational invariance of this
representation crucially alleviates the need for complex equivariant networks.
We train a denoising diffusion probabilistic model with a simple transformer
backbone and demonstrate that our resulting model unconditionally generates
highly realistic protein structures with complexity and structural patterns
akin to those of naturally-occurring proteins. As a useful resource, we release
the first open-source codebase and trained models for protein structure
diffusion.


  Electron cryo-microscopy (cryo-EM) produces three-dimensional (3D) maps of
the electrostatic potential of biological macromolecules, including proteins.
Along with knowledge about the imaged molecules, cryo-EM maps allow de novo
atomic modelling, which is typically done through a laborious manual process.
Taking inspiration from recent advances in machine learning applications to
protein structure prediction, we propose a graph neural network (GNN) approach
for automated model building of proteins in cryo-EM maps. The GNN acts on a
graph with nodes assigned to individual amino acids and edges representing the
protein chain. Combining information from the voxel-based cryo-EM data, the
amino acid sequence data and prior knowledge about protein geometries, the GNN
refines the geometry of the protein chain and classifies the amino acids for
each of its nodes. Application to 28 test cases shows that our approach
outperforms the state-of-the-art and approximates manual building for cryo-EM
maps with resolutions better than 3.5 \r{A}.


  Background: In the metagenome assembly of a microbiome community, we may
think abundant species would be easier to assemble due to their deeper
coverage. However, this conjucture is rarely tested. We often do not know how
many abundant species we are missing and do not have an approach to recover
these species.
  Results: Here we proposed k-mer based and 16S RNA based methods to measure
the completeness of metagenome assembly. We showed that even with PacBio
High-Fidelity (HiFi) reads, abundant species are often not assembled as high
strain diversity may lead to fragmented contigs. We developed a novel algorithm
to recover abundant metagenome-assembled genomes (MAGs) by identifying circular
assembly subgraphs. Our algorithm is reference-free and complement to standard
metagenome binning. Evaluated on 14 real datasets, it rescued many abundant
species that would be missing with existing methods.
  Conclusions: Our work stresses the importance of metagenome completeness
which is often overlooked before. Our algorithm generates more circular MAGs
and moves a step closer to the complete representation of microbiome
communities.


  Predicting the responses of a cell under perturbations may bring important
benefits to drug discovery and personalized therapeutics. In this work, we
propose a novel graph variational Bayesian causal inference framework to
predict a cell's gene expressions under counterfactual perturbations
(perturbations that this cell did not factually receive), leveraging
information representing biological knowledge in the form of gene regulatory
networks (GRNs) to aid individualized cellular response predictions. Aiming at
a data-adaptive GRN, we also developed an adjacency matrix updating technique
for graph convolutional networks and used it to refine GRNs during
pre-training, which generated more insights on gene relations and enhanced
model performance. Additionally, we propose a robust estimator within our
framework for the asymptotically efficient estimation of marginal perturbation
effect, which is yet to be carried out in previous works. With extensive
experiments, we exhibited the advantage of our approach over state-of-the-art
deep learning models for individual response prediction.


  Chan, Durrett, and Lanchier introduced a multitype contact process with
temporal heterogeneity involving two species competing for space on the
d-dimensional integer lattice. Time is divided into two seasons. They proved
that there is an open set of the parameters for which both species can coexist
when their dispersal range is sufficiently large. Numerical simulations
suggested that three species can coexist in the presence of two seasons. The
main point of this paper is to prove that this conjecture is incorrect. To do
this we prove results for a more general ODE model and contrast its behavior
with other related systems that have been studied in order to understand the
competitive exclusion principle.


  We analyze 14,651 HIV1 reverse transcriptase (HIV RT) sequences from the
Stanford HIV Drug Resistance Database labeled with treatment regimen in order
to study the evolution this enzyme under drug selection in the clinic. Our goal
is to identify distinct sectors of HIV RT's sequence space that are undergoing
evolution as a way to identify individual selections and/or evolutionary
solutions. We utilize Uniform Manifold Approximation and Projection (UMAP), a
graph-based dimensionality reduction technique uniquely suited for the
detection of non-linear dependencies and visualize the results using an
unsupervised clustering algorithm based on density analysis. Our analysis
produced 21 distinct clusters of sequences. Supporting the biological
significance of these clusters, they tend to represent phylogenetically related
sequences with strong correspondence to distinct treatment regimens. Thus, this
method for visualization of areas of HIV RT undergoing evolution can help infer
information about selective pressures, although it is correlative. The mutation
signatures associated with each cluster may represent the higher-order
epistatic context facilitating these evolutionary pathways, information that is
generally not accessible by other types of mutational co-dependence analyses.


  As the sequencing costs are decreasing, there is great incentive to perform
large scale association studies to increase power of detecting new variants.
Federated association testing among different institutions is a viable solution
for increasing sample sizes by sharing the intermediate testing statistics that
are aggregated by a central server. There are, however, standing challenges to
performing federated association testing. Association tests are known to be
confounded by numerous factors such as population stratification, which can be
especially important in multiancestral studies and in admixed populations among
different sites. Furthermore, disease etiology should be considered via
flexible models to avoid biases in the significance of the genetic effect. A
rising challenge for performing large scale association studies is the privacy
of participants and related ethical concerns of stigmatization and
marginalization. Here, we present dMEGA, a flexible and efficient method for
performing federated generalized linear mixed model based association testing
among multiple sites while underlying genotype and phenotype data are not
explicitly shared. dMEGA first utilizes a reference projection to estimate
population-based covariates without sharing genotype dataset among sites. Next,
dMEGA uses Laplacian approximation for the parameter likelihoods and decomposes
parameter estimation into efficient local-gradient updates among sites. We use
simulated and real datasets to demonstrate the accuracy and efficiency of
dMEGA. Overall, dMEGA's formulation is flexible to integrate fixed and random
effects in a federated setting.


  Exact solution of the Susceptible-Infectious-Recovered (SIR) epidemic model
is derived, and various properties of solution are obtained directly from the
exact solution. It is shown that there exists an exact solution of an initial
value problem for SIR differential system and that the parametric form of the
exact solution satisfies some linear differential system.


  Even though existence of non-convergent evolution of the states of
populations in ecological and evolutionary contexts is an undeniable fact,
insightful game-theoretic interpretations of such outcomes are scarce in the
literature of evolutionary game theory. As a proof-of-concept, we tap into the
information-theoretic concept of relative entropy in order to construct a
game-theoretic interpretation for periodic orbits in a wide class of
deterministic discrete-time evolutionary game dynamics, primarily investigating
the two-player two-strategy case. Effectively, we present a consistent
generalization of the evolutionarily stable strategy -- the cornerstone of the
evolutionary game theory -- and aptly term the generalized concept: information
stable orbit. The information stable orbit captures the essence of the
evolutionarily stable strategy in that it compares the total payoff obtained
against an evolving mutant with the total payoff that the mutant gets while
playing against itself. Furthermore, we discuss the connection of the
information stable orbit with the dynamical stability of the corresponding
periodic orbit.


  The 1000 Genomes Project provides sequencing data on 3,202 samples from 26
populations spanning five continental regions with no access or use
restrictions. The kgp R package provides consistent and comprehensive metadata
about samples and populations in the 1000 Genomes Project and other population
sequencing data in the International Genome Sample Resource collection. The kgp
package is distributed via the Comprehensive R Archive Network (CRAN) at
https://cran.r-project.org/package=kgp. Source code is available on GitHub at
https://github.com/stephenturner/kgp. Further documentation is online at
https://stephenturner.github.io/kgp/.


  The aim of this paper is to design a deep learning-based model to predict
proximal femoral strength using multi-view information fusion. Method: We
developed new models using multi-view variational autoencoder (MVAE) for
feature representation learning and a product of expert (PoE) model for
multi-view information fusion. We applied the proposed models to an in-house
Louisiana Osteoporosis Study (LOS) cohort with 931 male subjects, including 345
African Americans and 586 Caucasians. With an analytical solution of the
product of Gaussian distribution, we adopted variational inference to train the
designed MVAE-PoE model to perform common latent feature extraction. We
performed genome-wide association studies (GWAS) to select 256 genetic variants
with the lowest p-values for each proximal femoral strength and integrated
whole genome sequence (WGS) features and DXA-derived imaging features to
predict proximal femoral strength. Results: The best prediction model for fall
fracture load was acquired by integrating WGS features and DXA-derived imaging
features. The designed models achieved the mean absolute percentage error of
18.04%, 6.84% and 7.95% for predicting proximal femoral fracture loads using
linear models of fall loading, nonlinear models of fall loading, and nonlinear
models of stance loading, respectively. Compared to existing multi-view
information fusion methods, the proposed MVAE-PoE achieved the best
performance. Conclusion: The proposed models are capable of predicting proximal
femoral strength using WGS features and DXA-derived imaging features. Though
this tool is not a substitute for FEA using QCT images, it would make improved
assessment of hip fracture risk more widely available while avoiding the
increased radiation dosage and clinical costs from QCT.


  The experimental measurement of correlation functions and critical exponents
in disordered systems is key to testing renormalization group (RG) predictions.
We mechanically unzip single DNA hairpins with optical tweezers, an
experimental realization of the diffusive motion of a particle in a
one-dimensional random force field, known as the Sinai model. We measure the
unzipping forces $F_w$ as a function of the trap position $w$ in equilibrium
and calculate the force-force correlator $\Delta_m(w)$, its amplitude, and
correlation length, finding agreement with theoretical predictions. We study
the universal scaling properties since the effective trap stiffness $m^2$
decreases upon unzipping. Fluctuations of the position of the base pair at the
unzipping junction $u$ scales as $u \sim m^{-\zeta}$, with a roughness exponent
$ \zeta=1.34\pm0.06$, in agreement with the analytical prediction $\zeta =
\frac{4}{3}$. Our study provides a single-molecule test of the functional RG
approach for disordered elastic systems in equilibrium.


  Background: Drug synergy occurs when the combined effect of two drugs is
greater than the sum of the individual drugs' effect. While cell line data
measuring the effect of single drugs are readily available, there is relatively
less comparable data on drug synergy given the vast amount of possible drug
combinations. Thus, there is interest to use computational approaches to
predict drug synergy for untested pairs of drugs.
  Methods: We introduce a Graph Neural Network (GNN) based model for drug
synergy prediction, which utilizes drug chemical structures and cell line gene
expression data. We use information from the largest drug combination database
available (DrugComb), combining drug synergy scores in order to construct high
confidence benchmark datasets.
  Results: Our proposed solution for drug synergy predictions offers a number
of benefits: 1) It utilizes a combination of 34 distinct drug synergy datasets
to learn on a wide variety of drugs and cell lines representations. 2) It is
trained on constructed high confidence benchmark datasets. 3) It learns
task-specific drug representations, instead of relying on generalized and
pre-computed chemical drug features. 4) It achieves similar or better
prediction performance (AUPR scores ranging from 0.777 to 0.964) compared to
state-of-the-art baseline models when tested on various benchmark datasets.
  Conclusions: We demonstrate that a GNN based model can provide
state-of-the-art drug synergy predictions by learning task-specific
representations of drugs.


  We have gained access to vast amounts of multi-omics data thanks to Next
Generation Sequencing. However, it is challenging to analyse this data due to
its high dimensionality and much of it not being annotated. Lack of annotated
data is a significant problem in machine learning, and Self-Supervised Learning
(SSL) methods are typically used to deal with limited labelled data. However,
there is a lack of studies that use SSL methods to exploit inter-omics
relationships on unlabelled multi-omics data. In this work, we develop a novel
and efficient pre-training paradigm that consists of various SSL components,
including but not limited to contrastive alignment, data recovery from
corrupted samples, and using one type of omics data to recover other omic
types. Our pre-training paradigm improves performance on downstream tasks with
limited labelled data. We show that our approach outperforms the
state-of-the-art method in cancer type classification on the TCGA pan-cancer
dataset in semi-supervised setting. Moreover, we show that the encoders that
are pre-trained using our approach can be used as powerful feature extractors
even without fine-tuning. Our ablation study shows that the method is not
overly dependent on any pretext task component. The network architectures in
our approach are designed to handle missing omic types and multiple datasets
for pre-training and downstream training. Our pre-training paradigm can be
extended to perform zero-shot classification of rare cancers.


  Motivation: The rapid growth of metagenomics sequencing data makes
metagenomics increasingly dependent on computational and statistical methods
for fast and efficient analysis. Consequently, novel analysis tools for
big-data metagenomics are constantly emerging. One of the biggest challenges
for researchers occurs in the analysis planning stage: selecting the most
suitable metagenomics software tool to gain valuable insights from sequencing
data. The building process of data analysis pipelines is often laborious and
time-consuming since it requires a deep and critical understanding of how to
apply a particular tool to complete a specified metagenomics task.
  Results: We have addressed this challenge by using machine learning methods
to develop a classification system of metagenomics software tools into 13
classes (11 semantic annotations of EDAM and two virus-specific classes) based
on the descriptions of the tools. We trained three classifiers (Naive Bayes,
Logistic Regression, and Random Forest) using 15 text feature extraction
techniques (TF-IDF, GloVe, BERT-based models, and others). The manually curated
dataset includes 224 software tools and contains text from the abstract and the
methods section of the tools' publications. The best classification
performance, with an Area Under the Precision-Recall Curve score of 0.85, is
achieved using Logistic regression, BioBERT for text embedding, and text from
abstracts only. The proposed system provides accurate and unified
identification of metagenomics data analysis tools and tasks, which is a
crucial step in the construction of metagenomics data analysis pipelines.


  To successfully manage marine fisheries using an ecosystem-based approach,
long-term predictions of fish stock development considering changing
environmental conditions are necessary. Such predictions can be provided by
end-to-end ecosystem models, which couple existing physical and biogeochemical
ocean models with newly developed spatially-explicit fish stock models.
Typically, Individual-Based Models (IBMs) and models based on
Advection-Diffusion-Reaction (ADR) equations are employed for the fish stock
models. In this paper, we present a novel fish stock model called SPRAT for
end-to\hyp{}end ecosystem modeling based on Population Balance Equations (PBEs)
that combines the advantages of IBMs and ADR models while avoiding their main
drawbacks. SPRAT accomplishes this by describing the modeled ecosystem
processes from the perspective of individuals while still being based on
partial differential equations. We apply the SPRAT model to explore a
well-documented regime shift observed on the eastern Scotian Shelf in the 1990s
from a cod-dominated to a herring-dominated ecosystem. Model simulations are
able to reconcile the observed multitrophic dynamics with documented changes in
both fishing pressure and water temperature, followed by a predator-prey
reversal that may have impeded recovery of depleted cod stocks. We conclude
that our model can be used to generate new hypotheses and test ideas about
spatially interacting fish populations, and their joint responses to both
environmental and fisheries forcing.


  The stochastic reaction network in which chemical species evolve through a
set of reactions is widely used to model stochastic processes in physics,
chemistry and biology. To characterize the evolving joint probability
distribution in the state space of species counts requires solving a system of
ordinary differential equations, the chemical master equation, where the size
of the counting state space increases exponentially with the type of species,
making it challenging to investigate the stochastic reaction network. Here, we
propose a machine-learning approach using the variational autoregressive
network to solve the chemical master equation. Training the autoregressive
network employs the policy gradient algorithm in the reinforcement learning
framework, which does not require any data simulated in prior by another
method. Different from simulating single trajectories, the approach tracks the
time evolution of the joint probability distribution, and supports direct
sampling of configurations and computing their normalized joint probabilities.
We apply the approach to representative examples in physics and biology, and
demonstrate that it accurately generates the probability distribution over
time. The variational autoregressive network exhibits a plasticity in
representing the multimodal distribution, cooperates with the conservation law,
enables time-dependent reaction rates, and is efficient for high-dimensional
reaction networks with allowing a flexible upper count limit. The results
suggest a general approach to investigate stochastic reaction networks based on
modern machine learning.


  Glaucoma is the leading cause of irreversible blindness in people over the
age of 60, accounting for 6.6 to 8% of all blindness in 2010, but there is
still much to be learned about the genetic origins of the eye disease. With the
modern development of Next-Generation Sequencing (NGS) technologies, scientists
are starting to learn more about the genetic origins of Glaucoma. This research
uses differential expression (DE) and gene ontology (GO) analyses to study the
genetic differences between mice with severe Glaucoma and multiple control
groups. Optical nerve head (ONH) and retina data samples of genome-wide RNA
expression from NCBI (NIH) are used for pairwise comparison experimentation. In
addition, principal component analysis (PCA) and dispersion visualization
methods are employed to perform quality control tests of the sequenced data.
Genes with skewed gene counts are also identified, as they may be marker genes
for a particular severity of Glaucoma. The gene ontologies found in this
experiment support existing knowledge of Glaucoma genesis, providing confidence
that the results were valid. Future researchers can thoroughly study the gene
lists generated by the DE and GO analyses to find potential activator or
protector genes for Glaucoma in mice to develop drug treatments or gene
therapies to slow or stop the progression of the disease. The overall goal is
that in the future, such treatments can be made for humans as well to improve
the quality of life for human patients with Glaucoma and reduce Glaucoma
blindness rates.


  In this paper we study the appearance of bifurcations of limit cycles in an
epidemic model with two types of aware individuals. All the transition rates
are constant except for the alerting decay rate of the most aware individuals
and the rate of creation of the less aware individuals, which depend on the
disease prevalence in a non-linear way. For the ODE model, the numerical
computation of the limit cycles and the study of their stability are made by
means of the Poincar\'e map. Moreover, sufficient conditions for the existence
of an endemic equilibrium are also obtained. These conditions involve a rather
natural relationship between the transmissibility of the disease and that of
awareness. Finally, stochastic simulations of the model under a very low rate
of imported cases are used to confirm the scenarios of bistability (endemic
equilibrium and limit cycle) observed in the solutions of the ODE model.


  We propose a model to describe the adaptation of a phenotypically structured
population in a $H$-patch environment connected by migration, with each patch
associated with a different phenotypic optimum, and we perform a rigorous
mathematical analysis of this model. We show that the large-time behaviour of
the solution (persistence or extinction) depends on the sign of a principal
eigenvalue, $\lambda_H$, and we study the dependency of $\lambda_H$ with
respect to $H$. This analysis sheds new light on the effect of increasing the
number of patches on the persistence of a population, which has implications in
agroecology and for understanding zoonoses; in such cases we consider a
pathogenic population and the patches correspond to different host species. The
occurrence of a springboard effect, where the addition of a patch contributes
to persistence, or on the contrary the emergence of a detrimental effect by
increasing the number of patches on the persistence, depends in a rather
complex way on the respective positions in the phenotypic space of the optimal
phenotypes associated with each patch. From a mathematical point of view, an
important part of the difficulty in dealing with $H\ge 3$, compared to $H=1$ or
$H=2$, comes from the lack of symmetry. Our results, which are based on a fixed
point theorem, comparison principles, integral estimates, variational
arguments, rearrangement techniques, and numerical simulations, provide a
better understanding of these dependencies. In particular, we propose a precise
characterisation of the situations where the addition of a third patch
increases or decreases the chances of persistence, compared to a situation with
only two patches.


  Unlike vision and language data which usually has a unique format, molecules
can naturally be characterized using different chemical formulations. One can
view a molecule as a 2D graph or define it as a collection of atoms located in
a 3D space. For molecular representation learning, most previous works designed
neural networks only for a particular data format, making the learned models
likely to fail for other data formats. We believe a general-purpose neural
network model for chemistry should be able to handle molecular tasks across
data modalities. To achieve this goal, in this work, we develop a novel
Transformer-based Molecular model called Transformer-M, which can take
molecular data of 2D or 3D formats as input and generate meaningful semantic
representations. Using the standard Transformer as the backbone architecture,
Transformer-M develops two separated channels to encode 2D and 3D structural
information and incorporate them with the atom features in the network modules.
When the input data is in a particular format, the corresponding channel will
be activated, and the other will be disabled. By training on 2D and 3D
molecular data with properly designed supervised signals, Transformer-M
automatically learns to leverage knowledge from different data modalities and
correctly capture the representations. We conducted extensive experiments for
Transformer-M. All empirical results show that Transformer-M can simultaneously
achieve strong performance on 2D and 3D tasks, suggesting its broad
applicability. The code and models will be made publicly available at
https://github.com/lsj2408/Transformer-M.


  Predicting the binding structure of a small molecule ligand to a protein -- a
task known as molecular docking -- is critical to drug design. Recent deep
learning methods that treat docking as a regression problem have decreased
runtime compared to traditional search-based methods but have yet to offer
substantial improvements in accuracy. We instead frame molecular docking as a
generative modeling problem and develop DiffDock, a diffusion generative model
over the non-Euclidean manifold of ligand poses. To do so, we map this manifold
to the product space of the degrees of freedom (translational, rotational, and
torsional) involved in docking and develop an efficient diffusion process on
this space. Empirically, DiffDock obtains a 38% top-1 success rate (RMSD<2A) on
PDBBind, significantly outperforming the previous state-of-the-art of
traditional docking (23%) and deep learning (20%) methods. Moreover, while
previous methods are not able to dock on computationally folded structures
(maximum accuracy 10.4%), DiffDock maintains significantly higher precision
(21.7%). Finally, DiffDock has fast inference times and provides confidence
estimates with high selective accuracy.


  A phylogenetic variety is an algebraic variety parameterized by a statistical
model of the evolution of biological sequences along a tree. Understanding this
variety is an important problem in the area of algebraic statistics with
applications in phylogeny reconstruction. In the broader area of algebra
statistics, there have been important theoretical advances in computing certain
invariants associated to algebraic varieties arising in applications. Beyond
the dimension and degree of a variety, one is interested in computing other
algebraic degrees, such as the maximum likelihood degree and the Euclidean
distance degree. Despite these efforts, the current literature lacks explicit
computations of these invariants for the particular case of phylogenetic
varieties. In our work, we fill this gap by computing these invariants for
phylogenetic varieties arising from the simplest group-based models of
nucleotide substitution Cavender-Farris-Neyman model, Jukes-Cantor model,
Kimura 2-parameter model, and the Kimura 3-parameter model on small
phylogenetic trees with at most 5 leaves.


  The group of large aquatic mammals has representatives being the largest
living beings on earth, surpassing the weight and size of dinosaurs. In this
paper, we present some empirical evidence and a mathematical model to argue
that fat accumulation in marine mammals triggers a series of metabolic events
that result in these animals' increased size. Our study starts by analysing 43
ontogenetic trajectories of species of different types and sizes. For instance,
the analyses include organisms with asymptotic mass from 27g (Taiwan field
mouse) to $2.10^{7}$g (grey whale). The available data allows us to determine
all available species' ontogenetic parameters (catabolism and anabolism
constant, scaling exponent and asymptotic mass). The analyses of those data
show a minimisation of catabolism and scaling exponent in marine mammals
compared to other species analysed. We present a possible explanation for this,
arguing that the large proportion of adipose tissue in these animals can cause
this minimisation. That is because adipocytes have different scaling properties
in comparison to non-adipose (typical) cells, expressed in reduced energetic
demand and lower metabolism. The conclusion is that when we have an animal with
a relatively large amount of adipose tissue, as is the case of aquatic mammals,
the cellular metabolic rate decreases compared to other animals with the same
mass but with proportionally smaller fat tissue. A final consequence of this
cause-effect process is the increase of the asymptotic mass of these mammals.


  Spot instances offer a cost-effective solution for applications running in
the cloud computing environment. However, it is challenging to run long-running
jobs on spot instances because they are subject to unpredictable evictions.
Here, we present Spot-on, a generic software framework that supports
fault-tolerant long-running workloads on spot instances through checkpoint and
restart. Spot-on leverages existing checkpointing packages and is compatible
with the major cloud vendors. Using a genomics application as a test case, we
demonstrated that Spot-on supports both application-specific and transparent
checkpointing methods. Compared to running applications using on-demand
instances, it allows the completion of these workloads for a significant
reduction in computing costs. Compared to running applications using
application-specific checkpoint mechanisms, transparent checkpoint-protected
applications reduce runtime by up to 40%, leading to further cost savings of up
to 86%.


  Leveraging artificial intelligence for automatic retrosynthesis speeds up
organic pathway planning in digital laboratories. However, existing deep
learning approaches are unexplainable, like "black box" with few insights,
notably limiting their applications in real retrosynthesis scenarios. Here, we
propose MechRetro, a chemical-mechanism-driven graph learning framework for
interpretable retrosynthetic prediction and pathway planning, which learns
several retrosynthetic actions to simulate a reverse reaction via elaborate
self-adaptive joint learning. By integrating chemical knowledge as prior
information, we design a novel Graph Transformer architecture to adaptively
learn discriminative and chemically meaningful molecule representations,
highlighting the strong capacity in molecule feature representation learning.
We demonstrate that MechRetro outperforms the state-of-the-art approaches for
retrosynthetic prediction with a large margin on large-scale benchmark
datasets. Extending MechRetro to the multi-step retrosynthesis analysis, we
identify efficient synthetic routes via an interpretable reasoning mechanism,
leading to a better understanding in the realm of knowledgeable synthetic
chemists. We also showcase that MechRetro discovers a novel pathway for
protokylol, along with energy scores for uncertainty assessment, broadening the
applicability for practical scenarios. Overall, we expect MechRetro to provide
meaningful insights for high-throughput automated organic synthesis in drug
discovery.


  Mixing maternal and paternal genomes is the base of plant sexual
reproduction, but some so-called 'haploid inducer lines' lead to the formation
of seeds bearing well-developed embryos with solely the maternal genome. A
recent study adds a new piece to the puzzle of this enigmatic in planta haploid
embryo induction process.


  Human activity is leading to changes in the mean and variability of climatic
parameters in most locations around the world. The changing mean has received
considerable attention from scientists and climate policy makers. However,
recent work indicates that the changing variability, that is, the amplitude and
the temporal autocorrelation of deviations from the mean, may have greater and
more imminent impact on ecosystems.
  In this paper, we demonstrate that changes in climate variability alone could
drive cyclic predator-prey ecosystems to extinction via so-called phase tipping
(P-tipping), a new type of instability that occurs only from certain phases of
the predator-prey cycle. We construct a mathematical model of a variable
climate and couple it to two self-oscillating paradigmatic predator-prey
models. Most importantly, we combine realistic parameter values for the Canada
lynx and snowshoe hare with actual climate data from the boreal forest. In this
way, we demonstrate that critically important species in the boreal forest have
increased likelihood of P-tipping to extinction under predicted changes in
climate variability, and are most vulnerable during stages of the cycle when
the predator population is near its maximum. Furthermore, our analysis reveals
that stochastic resonance is the underlying mechanism for the increased
likelihood of P-tipping to extinction.


  Sheer cliffs on mountains and in deep canyons are among the worlds most
iconic landmarks but our understanding of the enigmatic flora that lives in
these vertical rock landscapes remains fragmented. In this article, I review
and synthesize recent studies on the evolution of specialization onto bare rock
and its consequences for plant diversification. Putative adaptations commonly
associated with growth on bare rock include specialized root structures, stress
tolerant leaf traits, and reduced dispersibility. Fitness tradeoffs are a
principal explanation for edaphic specialization, but adaptation to bare
environments stands apart as a precursor environment to specialization in other
stressful habitats, such as chemically harsh soils. In species level
phylogenies, many rock specialist plants are evolutionarily isolated and form
the sister lineage to congeners found on decomposed substrates, suggesting that
specialization onto bare rock may be an evolutionary trap or dead end. In other
cases, rock specialists form diverse clades, suggesting archipelago speciation
or ecological release enabled by innovations in stress tolerance. Today, in
environments severely impacted by poor management and introduced organisms,
such as islands, cliffs serve as important refuges for stress tolerant
endangered plants. New technology for climbing safety and drone assisted
reconnaissance open new possibilities for the discovery and conservation of
cliff plants in these landscapes but come with the inherent risks of increasing
their accessibility. The future success of efforts to save plants from
extinction may depend on our understanding of the unique and resilient flora
endemic to cliffy places.


  Food webs are complex ecological networks whose structure is both
ecologically and statistically constrained, with many network properties being
correlated with each other. Despite the recognition of these invariable
relationships in food webs, the use of the principle of maximum entropy
(MaxEnt) in network ecology is still rare. This is surprising considering that
MaxEnt is a statistical tool precisely designed for understanding and
predicting many different types of constrained systems. Precisely, this
principle asserts that the least-biased probability distribution of a system's
property, constrained by prior knowledge about that system, is the one with
maximum information entropy. Here we show how MaxEnt can be used to derive many
food-web properties both analytically and heuristically. First, we show how the
joint degree distribution (the joint probability distribution of the numbers of
prey and predators for each species in the network) can be derived analytically
using the number of species and the number of interactions in food webs.
Second, we present a heuristic and flexible approach of finding a network's
adjacency matrix (the network's representation in matrix format) based on
simulated annealing and SVD entropy. We built two heuristic models using the
connectance and the joint degree sequence as statistical constraints,
respectively. We compared both models' predictions against corresponding null
and neutral models commonly used in network ecology using open access data of
terrestrial and aquatic food webs sampled globally. We found that the heuristic
model constrained by the joint degree sequence was a good predictor of many
measures of food-web structure, especially the nestedness and motifs
distribution. Specifically, our results suggest that the structure of
terrestrial and aquatic food webs is mainly driven by their joint degree
distribution.


  The alignment of biological sequences such as DNA, RNA, and proteins, is one
of the basic tools that allow to detect evolutionary patterns, as well as
functional/structural characterizations between homologous sequences in
different organisms. Typically, state-of-the-art bioinformatics tools are based
on profile models that assume the statistical independence of the different
sites of the sequences. Over the last years, it has become increasingly clear
that homologous sequences show complex patterns of long-range correlations over
the primary sequence as a consequence of the natural evolution process that
selects genetic variants under the constraint of preserving the
functional/structural determinants of the sequence. Here, we present a new
alignment algorithm based on message passing techniques that overcomes the
limitations of profile models. Our method is based on a new perturbative
small-coupling expansion of the free energy of the model that assumes a linear
chain approximation as the $0^\mathrm{th}$-order of the expansion. We test the
potentiality of the algorithm against standard competing strategies on several
biological sequences.


  Inverse protein folding, the process of designing sequences that fold into a
specific 3D structure, is crucial in bio-engineering and drug discovery.
Traditional methods rely on experimentally resolved structures, but these cover
only a small fraction of protein sequences. Forward folding models like
AlphaFold offer a potential solution by accurately predicting structures from
sequences. However, these models are too slow for integration into the
optimization loop of inverse folding models during training. To address this,
we propose using knowledge distillation on folding model confidence metrics,
such as pTM or pLDDT scores, to create a faster and end-to-end differentiable
distilled model. This model can then be used as a structure consistency
regularizer in training the inverse folding model. Our technique is versatile
and can be applied to other design tasks, such as sequence-based protein
infilling. Experimental results show that our method outperforms
non-regularized baselines, yielding up to 3% improvement in sequence recovery
and up to 45% improvement in protein diversity while maintaining structural
consistency in generated sequences. Code is available at
https://github.com/IBM/AFDistill


  Epistatic interactions between mutations add substantial complexity to
adaptive landscapes, and are often thought of as detrimental to our ability to
predict evolution. Yet, patterns of global epistasis, in which the fitness
effect of a mutation is well-predicted by the fitness of its genetic
background, may actually be of help in our efforts to reconstruct fitness
landscapes and infer adaptive trajectories. Microscopic interactions between
mutations, or inherent nonlinearities in the fitness landscape, may cause
global epistasis patterns to emerge. In this brief review, we provide a
succinct overview of recent work about global epistasis, with an emphasis on
building intuition about why it is often observed. To this end, we reconcile
simple geometric reasoning with recent mathematical analyses, using these to
explain why different mutations in an empirical landscape may exhibit different
global epistasis patterns - ranging from diminishing to increasing returns.
Finally, we highlight open questions and research directions.


  Lampreys are jawless fish that evolved about 550 million years ago at the
base of the vertebrate line. Modern lampreys contain a corticoid receptor (CR),
the common ancestor of the glucocorticoid receptor (GR) and mineralocorticoid
receptor (MR), which first appear in cartilaginous fish, such as sharks. Until
recently, 344 amino acids at the amino terminus of adult lamprey CR were not
present in the lamprey CR sequence in GenBank. A search of the recently
sequenced lamprey germline genome identified two CR sequences, CR1 and CR2,
containing the 344 previously un-identified amino acids at the amino terminus.
CR1 also contains a novel four amino acid insertion in the DNA-binding domain
(DBD). We studied corticosteroid activation of CR1 and CR2 and found their
strongest response was to 11-deoxycorticosterone and 11-deoxycortisol, the two
circulating corticosteroids in lamprey. Based on steroid specificity, both CRs
are close to elephant shark MR and distant from elephant shark GR. HEK293 cells
transfected with full-length CR1 or CR2 and the MMTV promoter have about 3-fold
higher steroid-mediated activation compared to HEK293 cells transfected with
these CRs and the TAT3 promoter. Deletion of the amino-terminal domain (NTD) of
lamprey CR1 and CR2 to form truncated CRs decreased transcriptional activation
by about 70% in HEK293 cells transfected with MMTV, but increased transcription
by about 6-fold in cells transfected with TAT3, indicating that the promoter
has an important effect on NTD regulation of CR transcription by
corticosteroids.


  Spatially resolved transcriptomics (SRT) has evolved rapidly through various
technologies, enabling scientists to investigate both morphological contexts
and gene expression profiling at single-cell resolution in parallel. SRT data
are complex and multi-modal, comprising gene expression matrices, spatial
information, and often high-resolution histology images. Because of this
complexity and multi-modality, sophisticated computational algorithms are
required to accurately analyze SRT data. Most efforts in this domain have been
made to utilize conventional machine learning and statistical approaches,
exhibiting sub-optimal results due to the complicated nature of SRT datasets.
To address these shortcomings, researchers have recently employed deep learning
algorithms including various state-of-the-art methods mainly in spatial
clustering, spatially variable gene identification, and alignment. While great
progress has been made in developing deep learning-based models for SRT data
analysis, further improvement is still needed to create more biologically aware
models that consider aspects such as phylogeny-aware clustering or the analysis
of small histology image patches. Additionally, strategies for batch effect
removal, normalization, and handling overdispersion and zero inflation patterns
of gene expression are still needed in the analysis of SRT data using deep
learning methods. In this paper, we provide a comprehensive overview of these
deep learning methods, including their strengths and limitations. We also
highlight new frontiers, current challenges, limitations, and open questions in
this field. Also, we provide a comprehensive list of all available SRT
databases that can be used as an extensive resource for future studies.


  Shape-based virtual screening is widely employed in ligand-based drug design
to search chemical libraries for molecules with similar 3D shapes yet novel 2D
chemical structures compared to known ligands. 3D deep generative models have
the potential to automate this exploration of shape-conditioned 3D chemical
space; however, no existing models can reliably generate valid drug-like
molecules in conformations that adopt a specific shape such as a known binding
pose. We introduce a new multimodal 3D generative model that enables
shape-conditioned 3D molecular design by equivariantly encoding molecular shape
and variationally encoding chemical identity. We ensure local geometric and
chemical validity of generated molecules by using autoregressive fragment-based
generation with heuristic bonding geometries, allowing the model to prioritize
the scoring of rotatable bonds to best align the growing conformational
structure to the target shape. We evaluate our 3D generative model in tasks
relevant to drug design including shape-conditioned generation of chemically
diverse molecular structures and shape-constrained molecular property
optimization, demonstrating its utility over virtual screening of enumerated
libraries.


  During the Noachian, Mars' crust may have provided a favorable environment
for microbial life. The porous brine-saturated regolith would have created a
physical space sheltered from UV and cosmic radiations and provided a solvent,
while the below-ground temperature and diffusion of a dense reduced atmosphere
may have supported simple microbial organisms that consume H2 and CO2 as energy
and carbon sources and produce methane as a waste. On Earth, hydrogenotrophic
methanogenesis was among the earliest metabolisms but its viability on early
Mars has never been quantitatively evaluated. Here we present a probabilistic
assessment of Mars' Noachian habitability to H2-based methanogens, and quantify
their biological feedback on Mars' atmosphere and climate. We find that
subsurface habitability was very likely, and limited mainly by the extent of
surface ice coverage. Biomass productivity could have been as high as in early
Earth's ocean. However, the predicted atmospheric composition shift caused by
methanogenesis would have triggered a global cooling event, ending potential
early warm conditions, compromising surface habitability and forcing the
biosphere deep into the Martian crust. Spatial projections of our predictions
point to lowland sites at low-to-medium latitudes as good candidates to uncover
traces of this early life at or near the surface.


  One of the main concerns of Anfinsen was to reveal the connection between the
amino acid sequence and their biologically active conformation. This search
gave rise to two crucial questions in structural biology, namely, why the
proteins fold and how a sequence encodes its folding. As to the why, he
proposes a plausible answer, namely, at a given milieu a protein folds into its
functional form -- native state -- because such structure represents the lowest
free-energy minimum among all feasible conformations -- the thermodynamic
hypothesis or Anfinsen dogma. As to the how, this remains as an unsolved
challenge and, hence, this inquiry is examined here from a new perspective of
protein folding, namely, as an analytic whole -- a notion proposed by Leibnitz
and Kant. This new perspective forces us to discuss in detail why the
theoretical force-field-based approaches have failed in both their ability to
predict the three-dimensional structure of a protein accurately and in their
capacity to answer one of the most critical questions in structural biology:
how a sequence encodes its folding.


  Polyploidization is an evolutionary process by which a species acquires
multiple copies of its complete set of chromosomes. The reticulate nature of
the signal left behind by it means that phylogenetic networks offer themselves
as a framework to reconstruct the evolutionary past of species affected by it.
The main strategy for doing this is to first construct a so called
multiple-labelled tree and to then somehow derive such a network from it. The
following question therefore arises: How much can be said about that past if
such a tree is not readily available? By viewing a polyploid dataset as a
certain vector which we call a ploidy (level) profile we show that, among other
results, there always exists a phylogenetic network in the form of a beaded
phylogenetic tree with additional arcs that realizes a given ploidy profile.
Intriguingly, the two end vertices of almost all of these additional arcs can
be interpreted as having co-existed in time thereby adding biological realism
to our network, a feature that is, in general, not enjoyed by phylogenetic
networks. In addition, we show that our network may be viewed as a generator of
ploidy profile space, a novel concept similar to phylogenetic tree space that
we introduce to be able to compare phylogenetic networks that realize one and
the same ploidy profile. We illustrate our findings in terms of a publicly
available Viola dataset.


  Fragment-based drug discovery has been an effective paradigm in early-stage
drug development. An open challenge in this area is designing linkers between
disconnected molecular fragments of interest to obtain chemically-relevant
candidate drug molecules. In this work, we propose DiffLinker, an
E(3)-equivariant 3D-conditional diffusion model for molecular linker design.
Given a set of disconnected fragments, our model places missing atoms in
between and designs a molecule incorporating all the initial fragments. Unlike
previous approaches that are only able to connect pairs of molecular fragments,
our method can link an arbitrary number of fragments. Additionally, the model
automatically determines the number of atoms in the linker and its attachment
points to the input fragments. We demonstrate that DiffLinker outperforms other
methods on the standard datasets generating more diverse and
synthetically-accessible molecules. Besides, we experimentally test our method
in real-world applications, showing that it can successfully generate valid
linkers conditioned on target protein pockets.


  The complexity of mathematical models in biology has rendered model reduction
an essential tool in the quantitative biologist's toolkit. For stochastic
reaction networks described using the Chemical Master Equation, commonly used
methods include time-scale separation, the Linear Mapping Approximation and
state-space lumping. Despite the success of these techniques, they appear to be
rather disparate and at present no general-purpose approach to model reduction
for stochastic reaction networks is known. In this paper we show that most
common model reduction approaches for the Chemical Master Equation can be seen
as minimising a well-known information-theoretic quantity between the full
model and its reduction, the Kullback-Leibler divergence defined on the space
of trajectories. This allows us to recast the task of model reduction as a
variational problem that can be tackled using standard numerical optimisation
approaches. In addition we derive general expressions for the propensities of a
reduced system that generalise those found using classical methods. We show
that the Kullback-Leibler divergence is a useful metric to assess model
discrepancy and to compare different model reduction techniques using three
examples from the literature: an autoregulatory feedback loop, the
Michaelis-Menten enzyme system and a genetic oscillator.


  After the detection of coronavirus disease 2019 (Covid-19), caused by the
severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Wuhan, Hubei
Province, China in late December, the cases of Covid-19 have spiralled out
around the globe. Due to the clinical similarity of Covid-19 with other flulike
syndromes, patients are assayed for other pathogens of influenza like illness.
There have been reported cases of co-infection amongst patients with Covid-19.
Bacteria for example Streptococcus pneumoniae, Staphylococcus aureus,
Klebsiella pneumoniae, Mycoplasma pneumoniae, Chlamydia pneumonia, Legionella
pneumophila etc and viruses such as influenza, coronavirus,
rhinovirus/enterovirus, parainfluenza, metapneumovirus, influenza B virus etc
are identified as co-pathogens. In our current effort, we develop and analysed
a compartmental based Ordinary Differential Equation (ODE) type mathematical
model to understand the co-infection dynamics of Covid-19 and other influenza
type illness. In this work we have incorporated the saturated treatment rate to
take account of the impact of limited treatment resources to control the
possible Covid-19 cases. As results, we formulate the basic reproduction number
of the model system. Finally, we have performed numerical simulations of the
co-infection model to examine the solutions in different zones of parameter
space.


  Single-cell RNA-sequencing (scRNA-seq) has become a routinely used technique
to quantify the gene expression profile of thousands of single cells
simultaneously. Analysis of scRNA-seq data plays an important role in the study
of cell states and phenotypes, and has helped elucidate biological processes,
such as those occurring during development of complex organisms and improved
our understanding of disease states, such as cancer, diabetes, and COVID, among
others. Deep learning, a recent advance of artificial intelligence that has
been used to address many problems involving large datasets, has also emerged
as a promising tool for scRNA-seq data analysis, as it has a capacity to
extract informative, compact features from noisy, heterogeneous, and
high-dimensional scRNA-seq data to improve downstream analysis. The present
review aims at surveying recently developed deep learning techniques in
scRNA-seq data analysis, identifying key steps within the scRNA-seq data
analysis pipeline that have been advanced by deep learning, and explaining the
benefits of deep learning over more conventional analysis tools. Finally, we
summarize the challenges in current deep learning approaches faced within
scRNA-seq data and discuss potential directions for improvements in deep
algorithms for scRNA-seq data analysis.


  In this article, we introduce a random (directed) graph model for the
simultaneous forwards and backwards description of a rather broad class of
Cannings models with a seed bank mechanism. This provides a simple tool to
establish a sampling duality in the finite population size, and obtain a
path-wise embedding of the forward frequency process and the backward ancestral
process. Further, it allows the derivation of limit theorems that generalize
celebrated results by M\"ohle to models with seed banks, and where it can be
seen how the effect of seed banks affects the genealogies. The explicit
graphical construction is a new tool to understand the subtle interplay of seed
banks, reproduction and genetic drift in population genetics.


  Generating new molecules is fundamental to advancing critical applications
such as drug discovery and material synthesis. Flows can generate molecules
effectively by inverting the encoding process, however, existing flow models
either require artifactual dequantization or specific node/edge orderings, lack
desiderata such as permutation invariance, or induce discrepancy between the
encoding and the decoding steps that necessitates post hoc validity correction.
We circumvent these issues with novel continuous normalizing E(3)-equivariant
flows, based on a system of node ODEs coupled as a graph PDE, that repeatedly
reconcile locally toward globally aligned densities. Our models can be cast as
message-passing temporal networks, and result in superlative performance on the
tasks of density estimation and molecular generation. In particular, our
generated samples achieve state-of-the-art on both the standard QM9 and
ZINC250K benchmarks.


  In silico prediction of the ligand binding pose to a given protein target is
a crucial but challenging task in drug discovery. This work focuses on blind
flexible selfdocking, where we aim to predict the positions, orientations and
conformations of docked molecules. Traditional physics-based methods usually
suffer from inaccurate scoring functions and high inference cost. Recently,
data-driven methods based on deep learning techniques are attracting growing
interest thanks to their efficiency during inference and promising performance.
These methods usually either adopt a two-stage approach by first predicting the
distances between proteins and ligands and then generating the final
coordinates based on the predicted distances, or directly predicting the global
roto-translation of ligands. In this paper, we take a different route. Inspired
by the resounding success of AlphaFold2 for protein structure prediction, we
propose E3Bind, an end-to-end equivariant network that iteratively updates the
ligand pose. E3Bind models the protein-ligand interaction through careful
consideration of the geometric constraints in docking and the local context of
the binding site. Experiments on standard benchmark datasets demonstrate the
superior performance of our end-to-end trainable model compared to traditional
and recently-proposed deep learning methods.


  Molecular evolution is often conceptualised as adaptive walks on rugged
fitness landscapes, driven by mutations and constrained by incremental fitness
selection. It is well known that epistasis shapes the ruggedness of the
landscape's surface, outlining their topography (with high-fitness peaks
separated by valleys of lower fitness genotypes). However, within the strong
selection weak mutation (SSWM) limit, once an adaptive walk reaches a local
peak, natural selection restricts passage through downstream paths and hampers
any possibility of reaching higher fitness values. Here, in addition to the
widely used point mutations, we introduce a minimal model of sequence
inversions to simulate adaptive walks. We use the well known NK model to
instantiate rugged landscapes. We show that adaptive walks can reach higher
fitness values through inversion mutations, which, compared to point mutations,
allows the evolutionary process to escape local fitness peaks. To elucidate the
effects of this chromosomal rearrangement, we use a graph-theoretical
representation of accessible mutants and show how new evolutionary paths are
uncovered. The present model suggests a simple mechanistic rationale to analyse
escapes from local fitness peaks in molecular evolution driven by (intragenic)
structural inversions and reveals some consequences of the limits of point
mutations for simulations of molecular evolution.


  Reducing the cost of sequencing genomes provided by next-generation
sequencing technologies has greatly increased the number of genomic projects.
As a result, there is a growing need for better assembly and assembly
validation methods. One promising idea is to use heterogeneous data in assembly
projects. Optical Mapping (OM) is beneficial in validating genomic assemblies,
correction and scaffolding. Single raw OM read describes a DNA molecule's long
fragment, up to 1Mbp. Raw OM data from the same genome could be assembled to
create consensus maps that span an entire chromosome.
  The assembly process is computationally hard because of the large number of
errors in input data.
  This work describes a new algorithm and computer program to assemble OM reads
without a reference genome. In our algorithm, we explored binary representation
for genome maps. We focused on the efficiency of data structures and algorithms
and scale on parallel platforms. The algorithm consists of several steps, of
which the most important are : (1) conversion of the restriction maps into
binary strings, (2) detection of overlaps between restriction maps, (3)
determining the layout of restriction maps set, (4) creation of consensus
genomic maps. Our algorithm deals with optical mapping data with low error
levels but fails with high-level error reads.
  We developed a software library, console application and module for Python
language. The approach presented in this paper proved to be faster than a
dynamic programming approach and performed well on error-free data. It could be
used as a step of \textit{de~novo} assembly pipelines or to detect
misassemblies.The software is freely available in a public repository under GNU
LGPL v3 license (https://sourceforge.net/p/binary-genome-maps/code).


  Searching for extrasolar biosignatures is important to understand life on
Earth and its origin. Astronomical observations of exoplanets may find such
signatures, but it is difficult and may be impossible to claim unambiguous
detection of life by remote sensing of exoplanet atmospheres. Here, another
approach is considered: collecting grains ejected by asteroid impacts from
exoplanets in the Milky Way and then traveling to the Solar System. The optimal
grain size for this purpose is around 1 $\mu$m, and though uncertainty is
large, about $10^5$ such grains are expected to be accreting on Earth every
year, which may contain biosignatures of life that existed on their home
planets. These grains may be collected by detectors placed in space, or
extracted from Antarctic ice or deep-sea sediments, depending on future
technological developments.


  Antibodies comprise the most versatile class of binding molecules, with
numerous applications in biomedicine. Computational design of antibodies
involves generating novel and diverse sequences, while maintaining structural
consistency. Unique to antibodies, designing the complementarity-determining
region (CDR), which determines the antigen binding affinity and specificity,
creates its own unique challenges. Recent deep learning models have shown
impressive results, however the limited number of known antibody
sequence/structure pairs frequently leads to degraded performance, particularly
lacking diversity in the generated sequences. In our work we address this
challenge by leveraging Model Reprogramming (MR), which repurposes pretrained
models on a source language to adapt to the tasks that are in a different
language and have scarce data - where it may be difficult to train a
high-performing model from scratch or effectively fine-tune an existing
pre-trained model on the specific task. Specifically, we introduce ReprogBert
in which a pretrained English language model is repurposed for protein sequence
infilling - thus considers cross-language adaptation using less data. Results
on antibody design benchmarks show that our model on low-resourced antibody
sequence dataset provides highly diverse CDR sequences, up to more than a
two-fold increase of diversity over the baselines, without losing structural
integrity and naturalness. The generated sequences also demonstrate enhanced
antigen binding specificity and virus neutralization ability. Code is available
at https://github.com/IBM/ReprogBERT


  In drug discovery, aqueous solubility is an important pharmacokinetic
property which affects absorption and assay availability of drug. Thus, in
silico prediction of solubility has been studied for its utility in virtual
screening and lead optimization. Recently, machine learning (ML) methods using
experimental data has been popular because physics-based methods like quantum
mechanics and molecular dynamics are not suitable for high-throughput tasks due
to its computational costs. However, ML method can exhibit over-fitting problem
in a data-deficient condition, and this is the case for most chemical property
datasets. In addition, ML methods are regarded as a black box function in that
it is difficult to interpret contribution of hidden features to outputs,
hindering analysis and modification of structure-activity relationship. To deal
with mentioned issues, we developed Bayesian graph neural networks (GNNs) with
the self-attention readout layer. Unlike most GNNs using self-attention in node
updates, self-attention applied at readout layer enabled a model to improve
prediction performance as well as to identify atom-wise importance, which can
help lead optimization as exemplified for three FDA-approved drugs. Also,
Bayesian inference enables us to separate more or less accurate results
according to uncertainty in solubility prediction task We expect that our
accurate, reliable and interpretable model can be used for more careful
decision-making and various applications in the development of drugs.